|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
程序如下function P=Price(V,r,vol_i,t00,t_i,K_i,init)* R! u, Q# T. N% Y2 o; q: O
format long
* E" h0 |& G# l3 ft=t_i;
- r7 D! G. q* F! M! r* W. IK=K_i;
7 j. Z/ R4 h9 i% Yvol=vol_i;7 Z; w' d3 ` o: ` d# Y5 W
n=length(t);% computes the dimension of the problem; w7 e9 u8 X/ o, O! w6 {$ M
qq=init;
: |% ?( P i# ~5 ~C=zeros(n,1);! | e( ?1 `8 Z% h8 w- U( d
initial=K(n,1)*qq; C4 [9 N# i- m7 A* ? |
Vbar=zeros(n,1);
$ {# l E2 _$ A6 n3 e3 H. P9 _& ~N=zeros(n,n);
2 a+ _& Z# \; I) W; JN1=zeros(n,n);
8 B& O: h/ s, Ga=zeros(n,n);4 Q5 C* O/ e- P0 o |$ a# P
a0=-inf.* ones(n,n);2 Y# t/ I! }7 i. F2 G
b=zeros(n,n);5 |# R0 g1 P' h
b0=-inf.* ones(n,n);
1 P/ n# A/ z( c6 _h=zeros(n,n);
' c- c- O! f1 ?1 qhl=zeros(n,n);
$ E7 W/ z$ L( y1 ^Vbar(n,1)=K(n,1);6 \8 W# k: Z' C: j
for s=n:-1:2
; R. h' e m) d$ D current=s-1;* w8 |2 b8 n8 q- d8 E+ M [
t0=t(current,1);
* ?/ w; j4 @! u1 }( | options=optimset('Display','off');# f: M: q& O9 |' e" @. K$ h# R0 g
Vbar(s-1,1)=fzero('@Price_ss',initial,options,s,Vbar,r,vol,t0,t,K); %在这一步出错$ }. L; w/ _) T4 s
initial=Vbar(s-1,1)*qq;
$ R1 n& z( o- @" ]0 l[c1,h1,al,b1,N1,w2,wl]=Price_s(V,s,Vbar,r,vol,t0,t,K);
7 @; P3 T9 A( {) [C(s,1)=c1;: M6 `) U P4 E v" m- D
h(s, =hl;# E9 h3 S, v6 K/ T# J5 Z& w
a(s, =al;
$ {9 Q/ `, V* z# V9 S8 sb(s, =b1;4 q! v1 ?! Z% ~
N(s,:)=N1;- b, S, K6 X p4 E7 w; r
end l- k/ P9 V5 p: j+ h* ]
k=1;- z: b+ Y1 w( Q
[c2,h2,a2,b2,N2,w4,w3]=Price_s(V,k,Vbar,r,vol,t00,t,K);7 S7 c% T, J: P4 h6 K
C(k,1)=c2;
2 j2 @& V9 Z# y( _2 Lh(k,:)=h2;
' U6 B+ I4 g; n1 d1 [a(k,:)=a2;
$ H8 c7 F; s7 M6 ]9 k& Ab(k,:)=b2;+ r. ^3 t. g" |; n: _
N(k,:)=N2;
+ L* g `& w2 m4 l7 RP=C(k,1);
3 F9 i) c+ A4 J. L4 \% [
0 L+ k$ c& j6 P& m/ V; W! C
- H1 T! p! X+ S/ S) ~$ y$ J/ }报错2 O. y) I/ K6 l% y7 Q
未定义与 'function_handle' 类型的输入参数相对应的运算符 '=='。% }9 ?1 O" C0 Y, g" b
出错 fzero (line 314)3 j) ?9 D9 m8 ]0 }1 M) X! J' ?
if fx == 0
3 l, g% c" g8 E6 M4 c* U出错 Price (line 24)
0 H4 ^5 G* h& ~( N9 Q7 v Vbar(s-1,1)=fzero('@Price_ss',initial,options,s,Vbar,r,vol,t0,t,K);
6 [* Y+ D1 [- E& E6 l, a0 V6 S8 L% s( v9 G
请大神指教8 R/ \0 ~2 l: h
|
|