帮你编写了一个求正态分布的分布函数的反函数的程序,请保存成inversNormCDF.m。基本思路是二分法求方程的解。
function x=inversNormCDF(y)
xup=-100;
xlow=100;
for i=1:100
xmiddle=(xup+xlow)/2;
if normcdf(xmiddle)>y
xlow=xmiddle;
elseif normcdf(xmiddle)
else
x=xmiddle;
return
end
end
x=(xup+xlow)/2;
保存好了以后,在matlab命令行里输入如下命令,可以得到你想要的结果。
x=inversNormCDF(0.6)
>> norminv(0.6,900,100)
ans =
925.3347