作者yellow41535 (多多)
看板MATLAB
標題[運算] 詢問用matlab求解normcdf裡面的x
時間Mon Dec 26 21:13:52 2016
版上各位大大您們好,最近剛學matlab,遇到了一些問題,想詢問一下,謝謝您們
我的程式碼如下,一直出現錯誤訊息,好像是normcdf那邊出了問題,但看了許久無法解
決,所以想來求助版上各位大神,謝謝您們
clear
clc
gc=0.0184;
gd=0.0184;
c=0.0379;
d=0.12;
w=0.15;
r=1;
p=0.98;
h=2.25;
b0=2;
rf=0.3;
f=@(x) p*exp(gd-r*gc+((r^2)*(c^2)*(1-(w^2))/2))*((1+x)/x)*exp((d^2)-(2*r*w*d*c
)+((r^2)*(w^2)*(c^2)))+b0*p*(h*((1+x)/x)*exp(gd+((d^2)/2))-rf)+(1-h)*(((1+x)/x
)*exp(gd+((d^2)/2))*normcdf(((-log10(rf)+log10((1+x)/x))+gd+(d^2))/d)-(rf*(nor
mcdf(((-log10(rf)+log10((1+x)/x))+gd)/d))));
x=fzero(f,1);
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 110.28.12.27
※ 文章網址: https://www.ptt.cc/bbs/MATLAB/M.1482758034.A.EAF.html
推 ksr310135: 跟normcdf沒什麼關係,你這個函數在x>0的時候都是正的 01/01 03:40
→ ksr310135: 而x<=0的時候,log10((1+x)/x)就會出現虛數 01/01 03:42