MATLAB有条件期望似乎对对数正态不起作用

Question

我已经尝试了其他发行版（指数，正常等）的类似代码，在这种情况下，代码可行。但是，对于对数正态分布情况，Matlab计算数值积分的方式似乎存在问题。

代码：

function y=condint(sigma,mi,phihat)
y=quad(@phipr,phihat,110,0.000000000000001)
    function y=phipr(x)
      y=1./(1-logncdf(phihat)).*1/(2*pi*sigma^2).^(1/2).*exp(-((log(x)-mi).^2)./(2*sigma^2));
    end
end

没有明显的属性，即更高截止（Phihat）的条件期望值应更高。这是否是衰减非常快的分布的数值过程的通用问题？这是一个已知问题吗？谢谢您的帮助

P.S。在上面，我用作上限110，因为当使用参数时，MI = 0，Sigma = 0.25使用较高的值也会为CDF带来WIERD的结果。同样，这与梯形积分近似中的问题有关？

Answer 1

代码中的错误：logncdf（phihat）应为logncdf（phihat，mi，sigma）。

解决它。