5
金融计算与编程
上海财经大学金融学院 曹志广
直接调用 c=latticeeucall(52,50,0.1,6/12,0.2,500)得到：c= 5.5644 直接调用 p=latticeeuput(52,50,0.1,6/12,0.2,100)得到：p=1.1308； 直接调用 p=latticeeuput(52,50,0.1,6/12,0.2,200)得到：p=1.1240； 直接调用 p=latticeeuput(52,50,0.1,6/12,0.2,500)得到：p=1.1259 由以上计算可以看出：随着二叉树阶段数的增加，即时间间隔 ∆t 的减少，二叉树模型 的计算结果与期权价格的解析解也逐步接近。
T T T − r∆t
u = eσ
, d = e −σ
∆t
p=
e r∆t − d = 0.5076 u−d
∆t
股票价格的运动如图所示，期权的二叉树图如图所示。
u = eσ
= 1.1224, d = e −σ
∆t
= 0.8909
2
金融计算与编程
上海财经大学金融学院 曹志广
股票价格的运动
股票价格的运动
金融计算与编程
上海财经大学金融学院 曹志广
lattice(i+1,j+1)=exp(-r*deltaT)*(p*lattice(i+2,j+2)+(1-p)*lattice(i+2,j+1)); end end price=lattice(1,1);
function [price,lattice]=latticeamcall(S0,K,r,T,sigma,N) deltaT=T/N;u=exp(sigma*sqrt(deltaT)); d=1/u;p=(exp(r*deltaT)-d)/(u-d); lattice=zeros(N+1,N+1); for j=0:N lattice(N+1,j+1)=max(0,(S0*(u^j)*(d^(N-j))-K)); end for i=N-1:-1:0 for j=0:i lattice(i+1,j+1)=max(exp(-r*deltaT)*(p*lattice(i+2,j+2)+(1-p)*lattice(i+2,j+1)),-K+S0*u^j*d^ (i-j)); end end price=lattice(1,1); function [price,lattice]=latticeamput(S0,K,r,T,sigma,N) deltaT=T/N;u=exp(sigma*sqrt(deltaT)); d=1/u;p=(exp(r*deltaT)-d)/(u-d); lattice=zeros(N+1,N+1); for j=0:N lattice(N+1,j+1)=max(0,-(S0*(u^j)*(d^(N-j))-K)); end for i=N-1:-1:0 for j=0:i lattice(i+1,j+1)=max(exp(-r*deltaT)*(p*lattice(i+2,j+2)+(1-p)*lattice(i+2,j+1)),K-S0*u^j*d^( i-j)); end end price=lattice(1,1);
