※ 引述《abc76522000 (Ray)》之铭言： : 各位大大您好 我是Matlab初学者 : 而近在使用Matlab上遇到了一些问题，不知各位大大能否给点建议，谢谢您!!!! : 问题如下: : 1. : 如何使b精准的算出确切值(因为fzero算出得解是近似值) : clear : U=3; L=-3; : for i=1:300 : for a=i : f=@(b)1/3*norminv(1/2*normcdf((U-((a/100)-0.01))/b) + ... : 1/2*normcdf((((a/100)-0.01)-L)/b))-2; : b=fzero(f,1); : end : end 你要让norminv(x) = 2/3 ? 等同於 求 quantile of normal cdf = 2/3 所以就是在解 x = normcdf(2/3) = 0.7453731 因此，你只要专心解 1/2*normcdf((U-((a/100)-0.01))/b)+1/2*normcdf((((a/100)-0.01)-L)/b) =0.7453731 再两边都乘上2 normcdf((U-((a/100)-0.01))/b) + normcdf((((a/100)-0.01)-L)/b) = 1.490746 解这个会比你解上面那个会简单一点 不过你要确切值是不可能，我丢到symbolic toolbox去算 syms b solve(normcdf((U-((a/100)-0.01))/b) + normcdf((((a/100)-0.01)-L)/b) == 1.490746, b) 跑出warning: Warning: Cannot solve symbolically. Returning a numeric approximation instead. > In solve (line 303) 所以你只能用fzero求而已 : 2. : 如何精准的算出s，使其不会产生inf : clc : U=3; L=-3; : mu=0.023638; : sigma=0.362759; : s=1/3*norminv(1/2*normcdf((U-mu)./sigma)+1/2*normcdf((mu-L)./sigma)); (U-mu)./sigma = 8.2048 (mu-L)./sigma = 8.3351 normcdf(4.0) = 0.999968328758 normcdf(6.5) = 0.999999999960 你可能需要够长的位数才有办法去解出来不是无限大的s.... 程式只是工具别想靠程式帮你直接解决全部问题... 去了解一下你要解的问题比较重要 --

※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 140.109.73.190 ※ 文章网址: https://webptt.com/cn.aspx?n=bbs/MATLAB/M.1446088885.A.859.html