期权二叉树定价程序求助啊

我想算期权收敛的价格  但是为什么s和a的值是错误的啊 显示NaN
最后展示期权收敛的那个价格是用disp(price)吗？

对于 factorial(n) 来说，如果 n 太大了，MATLAB 出来的结果就是 Inf，所以最好的办法就是把

a=factorial(N)/(factorial(N-j)*factorial(j))*p^(j)*(1-p)^(N-j)*f;

换成

a=binopdf(j,N,p)*f;

就可以避免出现 NaN 的问题。我已经测试过，出来的 price 没问题

刚刚贴的图好像看不见。
S0=31.39;
K=31;
sigma=0.10911818;
T=1/12;
n=500;
r=0.23;
price=zeros(1,n);
for N=1:n;
    dt=T/N;
    u=exp(sigma*sqrt(dt));
    d=1/u;
    p=(exp(r*dt)-d)/(u-d);
    s=0;
    for j=0:N;
        f=max(S0*u^(j)*d^(N-j)-K,0);
        a=factorial(N)/(factorial(N-j)*factorial(j))*p^(j)*(1-p)^(N-j)*f;
        s=s+a;
    end
    price(1,N)=s*exp(-r*dt);
end
plot(price)

自己顶一下

我跑出来的图形是这样的，属于正常的吗？
正常的图形难道不该是这样的吗？

cmwei333 发表于 2016-6-1 06:41
对于 factorial(n) 来说，如果 n 太大了，MATLAB 出来的结果就是 Inf，所以最好的办法就是把

a=factoria ...

我跑出来的图形是图一这样的，属于正常的吗？

正常的图形难道不该是图二这样的吗？

成哲宇 发表于 2016-6-1 10:13
我跑出来的图形是图一这样的，属于正常的吗？

正常的图形难道不该是图二这样的吗？

我觉得

price(1,N)=s*exp(-r*dt);

应该是

price(1,N)=s*exp(-r*T);

出来的图其实差不多

cmwei333 发表于 2016-6-1 11:13
我觉得

price(1,N)=s*exp(-r*dt);

对啊 我看书上的收敛图都是上下对称的，所以我不知道是不是我的数据的问题。不过还是很感谢你之前的回答啦，学到了一个新函数哈哈