Как использовать Multiprecision Computing Toolbox с normcdf в Matlab - PullRequest
1 голос
/ 08 апреля 2019

У меня та же проблема, что и здесь https://comp.soft -sys.matlab.narkive.com / HV95Az6b / как вычислить нормальную кумулятивную вероятность с чрезвычайно высокой точностью

Я загружаю пакет 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;

Добавление mp () после 0.5 не работает, возможно, мне также следует отредактировать исходный кодэрф (или написать свой эрф)?

Заранее спасибо!

...