Function: 計財系碩一 108071603 黃偉德 1. Latticeeurput.M

Download as docx, pdf, or txt
Download as docx, pdf, or txt
You are on page 1of 4

金融商品設計與評價

HW6
計財系碩一 108071603 黃偉德
1. LatticeEurPut.m
function [price, lattice] =
LatticeEurPut(S0,X,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 , X - S0*(u^j)*(d^(N-j)));
end
for i=N-1:-1:0
for j=0:i
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);

2. CompLatticeBLSPut.m
% CompLatticeBLS.m
S0 = 50;
X = 50;
r = 0.1;
sigma = 0.4;
T = 5/12;
N=50;
BlsP = blsprice_Put(S0,X,r,T,sigma);
LatticeP = zeros(1,N);
for i=(1:N)
LatticeP(i) = LatticeEurPut(S0,X,r,T,sigma,i);
end
plot(1:N, ones(1,N)*BlsP);
hold on;
plot(1:N, LatticeP);

3. LatticeAmCall vs.BLS
function [price, lattice] =
LatticeAmCall(S0,X,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))-X);
end
for i=N-1:-1:0
for j=0:i
lattice(i+1,j+1) = max( S0*u^j*d^(i-j)-X , ...
exp(-r*deltaT) *(p * lattice(i+2,j+2) + (1-
p) * lattice(i+2,j+1)));
end
end
price = lattice(1,1);

4. LatticeAmCall vs. LatticeEurCall


function [price, lattice] =
LatticeAmCall(S0,X,r,T,sigma,N,q)
deltaT = T/N;
u=exp(sigma * sqrt(deltaT));
d=1/u;
p=(exp((r-q)*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))-X);
end
for i=N-1:-1:0
for j=0:i
lattice(i+1,j+1) = max( S0*u^j*d^(i-j)-X , ...
exp(-r*deltaT) *(p * lattice(i+2,j+2) + (1-
p) * lattice(i+2,j+1)));
end
end
price = lattice(1,1);

 用 LatticeEurCall(50,50,0.1,5/12,0.4,5)算出來的 call price=6.3595

 若要讓 LatticeAmCall(50,50,0.1,5/12,0.4,5,q)-6.3595=0
q=3.6269e-06
%LatticeAmCall v.s LatticeEurCall
fzero(@(q) LatticeAmCall(50,50,0.1,5/12,0.4,5,q)-
LatticeEurCall(50,50,0.1,5/12,0.4,5)+0.00001, 1)

 當 q>3.6269e-06,AmCall 會大於 EurCall,代表會提早履約。

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy