How can I plot distribution with shades?

Question

Giorgi 2015-9-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/243822-how-can-i-plot-distribution-with-shades

评论： Kelly Kearney 2015-9-18

Hello Matlabers :)

I have timeseries that I forecast for 12 period for 1000 times so I have 12X1000 matrix now and I want to make chart like this

Can anyone help me?

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

请先登录，再进行评论。

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

Answer 1

Joseph Cheng 2015-9-18

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/243822-how-can-i-plot-distribution-with-shades#answer_192919

在 MATLAB Online 中打开

you'd do something like this where you would use fill to generate the ranges of the distribution

t = 0:.1:10;
x = 2*t.^2;
ulimit= 10*rand(size(x));
llimit= 10*rand(size(x));
hold on
h(1)=fill([t fliplr(t)],[x+3*ulimit fliplr(x-3*llimit)],'g','edgecolor','none')
h(2)=fill([t fliplr(t)],[x+2*ulimit fliplr(x-2*llimit)],'b','edgecolor','none')
h(3)=fill([t fliplr(t)],[x+ulimit fliplr(x-llimit)],'r','edgecolor','none')
set(h,'facealpha',.5)
 plot(t,x,'b');

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Kelly Kearney 2015-9-18

在 MATLAB Online 中打开

My boundedline.m can help automate the procedure Joseph outlined:

% Data
t = (1:1000)';
y = 2*t.^2;
err = bsxfun(@times, rand(1000,12), t*1000);
y = bsxfun(@plus, y, err);
% Calculate bounds and plot
ymed = median(y, 2);
prc = [0 10 25 75 90 100];
p = prctile(y, prc, 2);
ylo = bsxfun(@minus, ymed, p(:,1:length(prc)/2));
yhi = bsxfun(@minus, p(:,end:-1:length(prc)/2+1), ymed);
bnd = permute(cat(3,ylo,yhi), [1 3 2]);
[hl, hp] = boundedline(t, [ymed ymed ymed], bnd);