当我在matlab中使用log-normal PDF时,结果与我的预期相符

时间:2016-03-03 07:07:47

标签: matlab

我正在学习论文

该论文提出了一个数字 enter image description here

该图显示了建筑物高度的CDF 并且该论文还提供了关于这个数字的详细信息

  

建筑物高度统计:本模型使用统计数据   在典型的建成区域中建立高度作为输入数据。一个合适的   通过与该市的地理数据进行比较来寻求形式   英国吉尔福德。概率密度函数是   选择适合数据的是未知的对数正态分布   参数:平均值p和标准差t。可以注意到   从图3中可以看出,它非常适合地理数据   参数p = 7.3m,t = 0.26的值。

它告诉平均值是7.3,标准差是0.26对吗?

然而,当我在matlab中通过添加代码

来尝试它们时
x=0:0.01:20;
meanValue = 7.3;
standardDeviation = 0.26;
y1 = logncdf(x,meanValue,standardDeviation);
plot(x,y1);

结果显示与图3不同 enter image description here

我试图重新阅读论文,以确保参数正确。 并在matlab上查看有关如何使用此方法的文档。 除模拟结果外,一切似乎都可以。

请帮我解决一下!感谢

1 个答案:

答案 0 :(得分:0)

如评论中所述,参数mu和sigma是相关正态分布的均值和标准推导,而不是对数正态分布。细节,特别是两者之间的联系是explained in the Wikipedia article

要根据均值和方差计算mu和sigma,公式在维基百科文章or here in the matlab syntax中给出:

m=7.3
t=0.26
v=t.^2;
%A lognormal distribution with mean m and variance v has parameters
mu = log((m^2)/sqrt(v+m^2));
sigma = sqrt(log(v/(m^2)+1));
%finally your code:
x=0:0.01:20;
y1 = logncdf(x,mu,sigma);
plot(x,y1);

enter image description here

这与您问题中的图表更接近,但您问题中的图表似乎是CDF的更高标准派生。从视觉上猜测参数形成你的情节,我会说它大致是t=5