我正在使用randn
函数生成一个高斯随机变量x1
,其均值为0,方差为1.生成x1
后,我生成另一个高斯随机变量{{1} }使用均值 m 和方差σ 2 ,即:x
。
我想创建一个log-normal随机变量x = m + sqrt(σ)·x1
并绘制其PDF。我想我可以用直方图函数来做到这一点:
y = ex
这是我得到的图表:
我认为我的结果是错误的,但我不知道如何修复它。
非常感谢你! :)
答案 0 :(得分:1)
lognormal distribution可以有一个很长的尾巴(即使y
prob(y) > 0
bins = linspace(0,5,500);
n = hist( y, bins );
bar( bins(1:end-1), n(1:end-1)/N ); axis([0 bins(end-1) 0 1]); % discard last bin that has the "tail"
并且不会消失。
如果你想将你的pdf情节与例如wiki web page的情节进行比较,你需要忽略尾巴。
{{1}}