找回密码
 注册
关于网站域名变更的通知
查看: 417|回复: 4
打印 上一主题 下一主题

使用fzero出现的问题

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2020-8-14 15:17 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

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

该用户从未签到

2#
发表于 2020-8-14 15:43 | 只看该作者
帮你顶一下

该用户从未签到

3#
发表于 2020-8-17 10:45 | 只看该作者
楼主发一下完整代码
( }/ f5 S0 j: C; r
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-11-24 15:10 , Processed in 0.140625 second(s), 23 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表