Matlab 欧式期权单步二叉树定价问题（风险中性）

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

      想请教大家一下关于写欧式看涨期权的二叉树定价问题，网上看到的二叉树定价代码是首先通过二叉树一步一步得到终端的股票价格，然后得到终端期权的payoff，再通过终端期权payoff与风险中性概率，一步一步地倒腾回去。
      但是，因为对于欧式期权而已，可以把期权的终端支付的期望直接通过无风险利率贴现回初始时刻，这样就一步到位了，我基于这个思想写了如下的Matlab代码：
function [CallP] = MYBSM(S0,K, r,T, sigma,M)
    dt = T/M;
    u = exp(sqrt(dt)*sigma);
    d = 1/u;
    S = zeros(M+1,1);
    NP = (exp(r*dt)-d)/(u-d);
    VectorP = zeros(M+1,1);
    for j = M:-1:0
            S(j+1,1) = S0*u^(j)*d^(M-j);
            VectorP(j+1,1) = NP^(j) * (1-NP)^(M-j)* nchoosek(M,j);
    end
    VectorC = max([S - K, zeros(M+1,1)],[],2);
    CallP = mean(VectorC .* VectorP) * exp(-r* T);
end


      这时候会出现问题，nchoosek(M,j)，当M很大时候，Matlab无法计算，但是二叉树定价在200期左右才收敛。请问大家有没有办法解决这个问题，如何通过风险中性一步计算出期权价格？
      谢谢大家！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：MATLAB matla atlab 欧式期权 Mat

代码倒数第二行有误，应该修改为CallP = sum(VectorC .* VectorP) * exp(-r* T);但是问题还是不变