红色填色的部分如何画出来，谢谢

re-us

下面的代码就是画一个高斯分布。如何把填色的部分画出来，谢谢

for i=1:201,
x(i)=(i-101)*0.05;
c=sqrt(2*pi);
f(i)=exp(-0.5*x(i)*x(i))/c;
end
plot(x,f)
axis([-5 5 0 0.5])
grid
title('Normal (Gaussian) Distribution')

[ 本帖最后由 ChaChing 于 2009-9-20 13:41 编辑 ]

VibrationMaster

for i=1:201,
x(i)=(i-101)*0.05;
c=sqrt(2*pi);
f(i)=exp(-0.5*x(i)*x(i))/c;
end
plot(x,f)
axis([-5 5 0 0.5])
grid
title('Normal (Gaussian) Distribution')
hold on
fill([2,2,5,2],[0,exp(-0.5*2^2)/c,exp(-0.5*5^2)/c,0],'r')

re-us

非常感谢 太子少傅
另外，最近我上穿jpg的文件总说不支持此类格式？奇怪ing

[ 本帖最后由 ChaChing 于 2009-9-21 09:53 编辑 ]

ChaChing

刚刚练了下, 去了for loop, 楼主参考下!
clc; clear
x=[-100:100]*0.05; f=exp(-0.5*x.*x)/sqrt(2*pi);
plot(x,f); axis([-5 5 0 0.5]); grid
title('Normal (Gaussian) Distribution')
hold on
ii=find(x==-2); fill([x(1:ii),-2,-5],[f(1:ii),0,0],'r')
ii=find(x==2);  fill([x(ii:end),5,2],[f(ii:end),0,0],'r')
hold off