1

我和这里有同样的问题https://comp.soft-sys.matlab.narkive.com/HV95Az6b/how-to-compute-the-normal-cumulative-probability-at-extremely-precision

我下载了包https://www.advanpix.com/documentation/users-manual/ 并且知道如何在一个简单的表达式中使用它

mp('pi/4')

ans = 

    0.78539816339744830961566.....

但是我不知道如何使用这个mp。在 normcdf 函数上，我尝试编辑以下代码中的最后一行表达式，但仍然可以解决。

function uu = normcdf(xx, mu, sigma)
%NORMCDF Gaussian CDF without using the stats toolbox
%     uu = normcdf(xx, mu, sigma)
% Iain Murray, July 2012
if ~exist('mu', 'var')
    mu = 0;
end
if ~exist('sigma', 'var')
    sigma = 1;
end
uu = 0.5*erf((xx-mu)./(sigma*sqrt(2))) + 0.5;

在 0.5 之后添加 mp() 不起作用，也许我还应该编辑 erf 的源代码（或者，编写我自己的 erf）？

提前致谢！

4

0 回答 0