我正在使用绘制概率密度函数的图表:
y = zeros(1,10000);
for j=1:10000
r = rand(100,1);
for i=1:100
y(j) = y(j) + r(i) - 0.5;
end
y(j) = y(j)/sqrt(100);
end
[n,x] = hist(y,100);
plot(x,n/10000/diff(x(1:2)));
hold on;
但我也想打印理论部分。目前为止,最好的方法似乎是以下方式:
plot(x,normpdf(x,0,1),'r');
但这与实际情况完全不符。我在这里漏掉了什么?目前我的图表如下所示。蓝色是实际情况,红色是理论情况。
![enter image description here](https://istack.dev59.com/g1jdS.webp)
plot(x,normpdf(x,0,1/pi),'r');
。 - Maksim Gorkiy