|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
程序如下function P=Price(V,r,vol_i,t00,t_i,K_i,init)
8 c! f9 h% w$ Fformat long
7 |" f* C# Z( h1 \2 jt=t_i;
6 M c) u6 `3 fK=K_i;
+ ~6 ^( m4 M% Q; a9 z% Ovol=vol_i;+ m. o2 Z o! K' f
n=length(t);% computes the dimension of the problem# i# I( b4 z/ q9 D% l" Z/ j
qq=init;
. ~4 m/ @# z2 p6 H, `9 D$ mC=zeros(n,1);
: }! `* \6 s4 L1 `& ~- F* dinitial=K(n,1)*qq;
6 e1 P2 L. N9 \' o% f9 X9 EVbar=zeros(n,1);
% g! C b v; I& B+ M7 \N=zeros(n,n);) R4 v7 c* q. D% T
N1=zeros(n,n);
1 m! u1 V( R$ U( p& I$ Aa=zeros(n,n);
( f2 c9 }8 s9 K( `- l6 |& F3 \* ha0=-inf.* ones(n,n);
6 W( z# ]7 [! S8 A% ?( o% n6 g3 {b=zeros(n,n);& Q% h* |8 [1 P& c% D( ~1 G* l% q
b0=-inf.* ones(n,n);7 E* W( ~- ]1 E X* `
h=zeros(n,n);+ c o3 @; M: C: h) f2 ?' Y
hl=zeros(n,n);1 I# s* V, M) f$ ?/ |& ^9 ^ F* n
Vbar(n,1)=K(n,1);
4 e' K9 {2 L( |" }for s=n:-1:2
1 {8 k* c" y% _( ]# D current=s-1;
% ~8 e( s- o5 D" v t0=t(current,1);5 B) A+ f0 }! @" ?1 B6 I
options=optimset('Display','off');
7 f- @% b8 I1 `: Y; y Vbar(s-1,1)=fzero('@Price_ss',initial,options,s,Vbar,r,vol,t0,t,K); %在这一步出错. Y& {( x$ l+ b
initial=Vbar(s-1,1)*qq;
8 e" j+ b1 p, p+ g) z[c1,h1,al,b1,N1,w2,wl]=Price_s(V,s,Vbar,r,vol,t0,t,K);( U' W1 |9 c9 O) ~( R0 f5 ]0 `( H
C(s,1)=c1;
( s& }9 K# e0 uh(s, =hl;4 _2 Y3 x6 }* B3 V
a(s, =al;
0 Z; j1 q8 Y2 B$ P4 M2 V$ N7 Lb(s, =b1;
" I1 v I2 _, `N(s,:)=N1;: I. @) _7 E" U9 `( C2 {
end
6 b& \% u. }+ @' e4 l- tk=1;9 P: u$ A1 Y5 s# z
[c2,h2,a2,b2,N2,w4,w3]=Price_s(V,k,Vbar,r,vol,t00,t,K);3 T* j( K7 P4 n# k! U }8 n& r+ x
C(k,1)=c2;
+ w X- B8 F# G# k+ l3 Rh(k,:)=h2;) u6 p+ l6 g2 e1 d. Z/ `
a(k,:)=a2;+ P* _; w3 S: c( }. y
b(k,:)=b2;) R: u1 z1 d, u5 N. y {8 p
N(k,:)=N2;
1 `2 R4 K( V+ D' q$ b6 P' sP=C(k,1);1 A% q+ b/ p- m; ]3 ]
; b7 C& e9 n ~" |/ t3 l
D; R8 T( S0 n4 j" s! j
报错9 F' X) J& w% \# o7 M: ~
未定义与 'function_handle' 类型的输入参数相对应的运算符 '=='。7 a% o& Q, @( r- _
出错 fzero (line 314)
* J' Q! }0 {. V( p if fx == 0. d; E5 y4 k4 w/ I
出错 Price (line 24)
" {6 O% G* n/ \4 T" d$ _# H Vbar(s-1,1)=fzero('@Price_ss',initial,options,s,Vbar,r,vol,t0,t,K);
! A+ S1 u% q' E; v) t
' Y ?1 v; z% p; t" e请大神指教* @& G! ~1 x0 }; S' t- o t9 ^* u% m2 Q
|
|