Shaded error area of Std Dev.
4 次查看(过去 30 天)
显示 更早的评论
Hello Community,
I am trying to shade the area between 2 traces of mean+Std and mean-Std with the mean passes in between. I am using:
x = CAD;
mean_y = Pav;
std_y = std(P_all,0,2);
y1 = mean_y;
y2 = mean_y - std_y;
y3 = mean_y + std_y;
figure(1)
plot(x, y1)
hold on
plot(x, y2)
hold on
plot(x, y3)
patch([x fliplr(x)], [y1 fliplr(y2)], 'g')
patch([x fliplr(x)], [y2 fliplr(y3)], 'g')
hold off
However, I get the whole area shaded under the trace due to a horizontal line that appears when excuting
patch([x fliplr(x)], [y1 fliplr(y2)], 'g')
patch([x fliplr(x)], [y2 fliplr(y3)], 'g')
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1317365/image.jpeg)
Please advise if there is a better approach
0 个评论
采纳的回答
Star Strider
2023-3-8
编辑:Star Strider
2023-3-8
It would of course help to have your data.
However this call:
std_y = std(P_all,0,2);
implies to me that your data are column vectors, not row vectors.
If that is correct, then this would be more appropriate:
patch([x; flip(x)], [y1; flip(y2)], 'g')
patch([x; flip(x)], [y2; flip(y3)], 'g')
It is the same essential idea, simply recognising the orientation of the data.
The,flip works regardless of the orientation, while fliplr simply filps the column vectors in this instance (although is appropriate if that is what you want to do).
EDIT — (8 Mar 2023 at 2:28)
With the supplied data (they are column data) —
LD = load(websave('Mathworks','https://www.mathworks.com/matlabcentral/answers/uploaded_files/1317410/Mathworks.mat'))
CAD = LD.CAD;
P_all = LD.P_all;
Pav = LD.Pav;
x = (-360:0.5:359.5)';
mean_y = Pav;
std_y = std(P_all,0,2);
y1 = mean_y;
y2 = mean_y - std_y;
y3 = mean_y + std_y;
figure
plot(x,y1)
hold on
patch([x; flip(x)], [y1; flip(y2)], 'g')
patch([x; flip(x)], [y2; flip(y3)], 'g')
hold off
xlabel('CAD')
figure
plot(x,y1)
hold on
patch([x; flip(x)], [y1; flip(y2)], 'g')
patch([x; flip(x)], [y2; flip(y3)], 'g')
hold off
xlabel('CAD')
ylim([50 80])
title('Zoom to show detail')
The standad deviation is quite small, so the lower (optional) plot zooms in to show the details.
.
4 个评论
更多回答(0 个)
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!