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

使用fzero出现的问题

[复制链接]

该用户从未签到

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

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

该用户从未签到

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

该用户从未签到

3#
发表于 2020-8-17 10:45 | 只看该作者
楼主发一下完整代码
$ q& W) }; _) B6 `1 A
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

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

EDA365公众号

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

GMT+8, 2025-7-27 19:52 , Processed in 0.125000 second(s), 23 queries , Gzip On.

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

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

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