由于给出的积分方程为标准的正态分布函数,所以可以用matlab自带的normcdf()来求解积分值。
function y =norm_fun(x) %norm_fun.m
y=normcdf(x,0,1);
end
>>x=0.5;y =norm_fun(x)
y=0.69146
也可以利用 quadgk()自适应数值积分函数(Gauss Kronrod quadrature)来求解积分值。
function y =GK_fun(x) %GK_fun.m
func=@(u)exp(-u.^2/2)/sqrt(2*pi);
y=quadgk(func,-inf,x);
end
>>x=0.5;y =GK_fun(x)
y=0.69146