How to add different colours to positive and negative regions of area plot?

Question

Kabitri Chattopadhyay 2015-9-22

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/244481-how-to-add-different-colours-to-positive-and-negative-regions-of-area-plot

评论： Kabitri Chattopadhyay 2015-9-23

fig1.jpg

Hi, I am trying to make an area plot with separate colours for positive and negative regions. I tried separating the positive and the negative parts and then plotted both one upon another:

for i=1:length(mismatch)
  if (mismatch(i) < 0)
      neg_mismatch(i)=mismatch(i);
      pos_mismatch(i)=0.0;
  else
      neg_mismatch(i)=0.0;
      pos_mismatch(i)=mismatch(i);
  end
end

The problem is there are overlap regions of positive and negative values around 0 values (see attachment) which makes the plot look ugly. I tried to interpolate values in between, but still no improvement!

     x_neg=(1:length(neg_mismatch))';
     y_neg=neg_mismatch;
     xi_neg=(1:0.01:length(neg_mismatch))';
     yi_neg = interp1(x_neg,y_neg,xi_neg,'cubic');
     x_pos=(1:length(pos_mismatch))';
     y_pos=pos_mismatch;
     xi_pos=(1:0.01:length(pos_mismatch))';
     yi_pos = interp1(x_pos,y_pos,xi_pos,'cubic');
     figure;
     hold on
     h1=area(yi_neg);
     h2=area(yi_pos);
     hold off

Can anyone please help!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

WAT 2015-9-22

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/244481-how-to-add-different-colours-to-positive-and-negative-regions-of-area-plot#answer_193343

在 MATLAB Online 中打开

The problem is that the area command is drawing lines from the actual data points to the zeroed out data, rather than connecting to the actual (but opposite signed) data point.

As a simple example, think of what would happen if you tried to do a positive/negative area plot for the simple two-point data set [-1,1].

There may be a more elegant solution, but one way is to insert a fake 0 data point wherever the data switch from positive to negative (or vice-versa):

 %%MAKE UP SOME DATA
 x = -10:.1:10;
 y = sinc(x)+0.2*random('normal',0,1,1,length(x));
 plot(y)
 %%TRY SPLITTING IT THE OLD WAY
 y1=y;y2=y;
 y1(y1<0)=0; % zero out negative entries
 y2(y2>0)=0; % zero out positive entries
 figure
 hold on
 h1 = area(x,y1,'FaceColor','b')
 h2 = area(x,y2,'FaceColor','r');
 %%NOW INSERT ZEROS WHERE LINE CROSSES THE Y AXIS
 % Insert 0s wherever data transitions from +/- or -/+
 count = 0; % counter for number of fake 0's inserted
 for ii=1:length(y)-1
   i = ii+count; % i is the index in our new, possibly extended vectors
   % find points with different signs
   if(y(i+1)*y(i) < 0) 
       % find y intercept
       x_fake = x(i) + abs(y(i))/(abs(y(i))+abs(y(i+1))) * (x(i+1)-x(i)); 
       % put fake data into x and y arrays
       x = [ x(1:i), x_fake , x(i+1:length(x))];
       y = [ y(1:i), 0 , y(i+1:length(y))];
       % increment counter
       count = count+1;
   end  
end
 % create y1 and y2, consisting of nonnegative and nonpositive entries
 % of y respectively
 y1=y;  y2=y;
 y1(y1<0) = 0;
 y2(y2>0) = 0;
 figure
 hold on
 area(x,y1,'FaceColor','b');
 area(x,y2,'FaceColor','r');
 hold off