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

使用fzero出现的问题

[复制链接]

该用户从未签到

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

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

该用户从未签到

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

该用户从未签到

3#
发表于 2020-8-17 10:45 | 只看该作者
楼主发一下完整代码
. x6 ~) a% P8 s$ i- m2 c
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

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

EDA365公众号

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

GMT+8, 2025-6-24 03:22 , Processed in 0.078125 second(s), 23 queries , Gzip On.

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

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

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