|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
程序如下function P=Price(V,r,vol_i,t00,t_i,K_i,init), v5 R T8 J! T- M4 X; q& j
format long! { ?/ N# \% d( h& o2 R2 k
t=t_i;2 `! \ R) f; ?
K=K_i;
& T6 P0 ]; S. G% x/ o4 Tvol=vol_i;* q4 e$ `) d1 p* g
n=length(t);% computes the dimension of the problem
) l0 c7 m$ O9 u& J. H9 D8 ~qq=init;: `7 F, z# x# c3 }6 c9 s) u) A9 T5 @
C=zeros(n,1);3 f+ w& O) h6 x- q" y8 M, z; I. P
initial=K(n,1)*qq;
' e Z6 k- `8 @! V; [Vbar=zeros(n,1);
! V9 r3 d# F6 j0 X# x1 L; j6 T; |+ pN=zeros(n,n);0 J! x" D0 @+ r K; B. \
N1=zeros(n,n);
5 |7 W" L' f& |" \: q G2 ?a=zeros(n,n);
" N4 X) z: y( W8 p7 wa0=-inf.* ones(n,n);
6 Z6 I- u) b: y* B1 S+ nb=zeros(n,n);
: I$ X: i" H8 W; E; e5 Z# w5 Cb0=-inf.* ones(n,n);
8 M$ h9 q; a( |h=zeros(n,n);
H2 c* ^) O' I9 ~5 yhl=zeros(n,n);# h( k! E* j+ o3 J, K; S
Vbar(n,1)=K(n,1);
! c/ m1 x* ], X! p, L' g- zfor s=n:-1:2
- s0 L! w8 d9 g* a current=s-1;
% d0 X6 B1 k6 t. s t0=t(current,1);# i/ q" _- ]# r9 R# _* R: ~
options=optimset('Display','off');
' ~3 S4 Z6 s2 w, C' }$ a Vbar(s-1,1)=fzero('@Price_ss',initial,options,s,Vbar,r,vol,t0,t,K); %在这一步出错
/ d( _/ w' K5 H5 ?7 v0 f" t initial=Vbar(s-1,1)*qq;
0 ^2 E0 _2 M5 L- R# b5 E[c1,h1,al,b1,N1,w2,wl]=Price_s(V,s,Vbar,r,vol,t0,t,K);
/ f8 R! T0 L- n6 }' u' ]C(s,1)=c1;
! F6 ]; w, b/ B, ^2 hh(s, =hl;
' v. R' o& r9 Y6 ra(s, =al;
4 H5 T6 ~* }. w) q# `) ^! Qb(s, =b1;( }; ^2 E% n$ x) f. L
N(s,:)=N1;
7 v4 X( f+ k0 v0 s0 L' E& N0 T9 oend
, i4 R/ n# ~! I) x6 T/ z2 mk=1;( Z0 h Z4 q# ]2 n
[c2,h2,a2,b2,N2,w4,w3]=Price_s(V,k,Vbar,r,vol,t00,t,K);
# S. m7 d) C9 W8 BC(k,1)=c2;
w0 b4 ?* z& c+ F- }h(k,:)=h2;. u/ d @1 h& G [( e4 L
a(k,:)=a2;
/ Y& b; R2 m! {! A! mb(k,:)=b2;! N; b3 K# }# c# ?: f9 i4 C# u
N(k,:)=N2;
' w# ]6 F% L# x) gP=C(k,1);: j- O" A1 i" L: s9 G5 P; L, g/ F
& s( q4 Q" K: {5 O
( N5 Z: R, X) S" i% n
报错
! O* r, Y& `; k' r( h- \- D9 j7 R未定义与 'function_handle' 类型的输入参数相对应的运算符 '=='。
6 P& [3 y# g: ? A2 @4 B出错 fzero (line 314)
3 k" T; n, G$ O% L3 |/ ~ if fx == 0
" }' R4 M! W1 C6 G* @1 c' ]5 V出错 Price (line 24)
- _- r, Q$ A+ k$ ]9 a8 o4 Q/ T Vbar(s-1,1)=fzero('@Price_ss',initial,options,s,Vbar,r,vol,t0,t,K);- X9 U6 K$ i) M
% h: m* h9 _; v请大神指教2 f8 B) {) Z6 v
|
|