如何在MATLAB中实现此功能？

[poster]

我正在使用具有我要使用的功能的书。但是我不认为我从函数中得到了正确的值。

这是书中的说明：

https://i.stack.imgur.com/yCns1.png

这是我在MATLAB中创建的函数：

function [ shortProbability ] = pShort( zkt, zktStar, short) if zkt > zktStar shortProbability = 0; else normalizer = 1/(1-exp(-short*zktStar)); shortProbability = normalizer * (short * exp(-short*zkt)); end end 我插入的值是：

zkt = 0:1:100 zktStar = 50; short = 0.01; 但是，我的图的行为不像我应该以的形式结束，这是这样的：

https://i.stack.imgur.com/G7z7M.png

我从图表中得到了这个，看起来很正确，但是我认为它没有被正确规范化：

https://i.stack.imgur.com/vxIeO.png

谁能帮我纠正此功能？



回答：

这就是指数分布 (http://en.wikipedia.org/wiki/Exponential_distribution) 。您可以从统计信息工具箱使用EXPPDF (http://www.mathworks.com/help/toolbox/stats/exppdf.html)和EXPCDF (http://www.mathworks.com/help/toolbox/stats/expcdf.html) ：

normalizer = 1 ./ ( expcdf(zktStar,1/short) - expcdf(0,1/short) ); shortProbability = exppdf(zkt, 1/short) * normalizer; 它应该等同于您所拥有的...

这是我将其与您先前问题中 (https://stackoverflow.com/questions/3644037/help-understanding-a-definitive-integral/3644076#3644076)的图组合后得到的结果：

https://i.stack.imgur.com/EP01w.png

为了确认，我们计算每条曲线下的面积（足够接近1 ）：

>> trapz(zkt,hitProbabilty) ans = 1 >> trapz(zkt,shortProbability) ans = 1.0077

更多&回答... (https://stackoverflow.com/questions/3646231)