EDA365电子论坛网
标题:
使用fzero出现的问题
[打印本页]
作者:
lupkpu9
时间:
2020-8-14 15:17
标题:
使用fzero出现的问题
程序如下function P=Price(V,r,vol_i,t00,t_i,K_i,init)
; }, [* A- f- ^3 q! p
format long
. e( x% i& W# R Q" D% A
t=t_i;
X% T/ ?' x% f0 t
K=K_i;
: ]# b$ g$ B1 Y0 M1 L
vol=vol_i;
( [: i. d% j/ q
n=length(t);% computes the dimension of the problem
i# ]% t; @7 ~6 f% Z& X3 K
qq=init;
* e4 l6 E, ]+ u6 z, j C
C=zeros(n,1);
' _. W. ^/ D2 Q2 ^7 {! {( G
initial=K(n,1)*qq;
% ]# I. w1 {3 ?) z3 n
Vbar=zeros(n,1);
2 u% t3 Y+ f5 }
N=zeros(n,n);
) _2 O4 C4 \; \0 l
N1=zeros(n,n);
1 N- S% ^6 d$ G, }
a=zeros(n,n);
4 F% @8 a2 s9 b4 P
a0=-inf.* ones(n,n);
) S6 Z! ~# J* y0 j5 W% m: l1 V
b=zeros(n,n);
" \4 z% F" Y4 H
b0=-inf.* ones(n,n);
+ F; K r& l+ m( b
h=zeros(n,n);
6 }* N0 {' W% V, `9 e- f
hl=zeros(n,n);
' X$ V) T V2 ?& X& F7 i6 q- ?
Vbar(n,1)=K(n,1);
5 Q9 w, N( E1 ]1 s8 f4 B
for s=n:-1:2
% u% _1 M! ?+ k+ f. C+ _5 f
current=s-1;
# Q U5 }& M2 @3 ^
t0=t(current,1);
* a" I- s0 c9 W& ? k0 [! G
options=optimset('Display','off');
4 R8 E0 t% I$ n/ R* s: K( n4 e( ~: E
Vbar(s-1,1)=fzero('@Price_ss',initial,options,s,Vbar,r,vol,t0,t,K);
%在这一步出错
9 A, ]+ h# t, G# e; [4 L6 }% Y* _
initial=Vbar(s-1,1)*qq;
" B: H l' A/ G, R7 ]: O
[c1,h1,al,b1,N1,w2,wl]=Price_s(V,s,Vbar,r,vol,t0,t,K);
# N! L8 v7 T. _# D/ B7 Z2 J' k# ^& O- d
C(s,1)=c1;
0 [- r# @' @4 \4 I
h(s,
=hl;
4 F" Y7 @! W( I2 ?* X5 w
a(s,
=al;
# Z! b) W* Y$ e/ Y X+ j% S$ l2 s/ W9 I
b(s,
=b1;
+ y7 f# I; p9 T) }3 N% c7 Q
N(s,:)=N1;
& V2 ^5 X/ P4 c6 n, [ Q
end
+ o8 W! B' c' m" _
k=1;
& o6 z" e2 j( p" j! E
[c2,h2,a2,b2,N2,w4,w3]=Price_s(V,k,Vbar,r,vol,t00,t,K);
' M: M6 h6 U2 T
C(k,1)=c2;
! f& v3 |. h O$ M% r0 a7 D7 S
h(k,:)=h2;
' o8 r' |7 s4 g. {: W( F' Y) [
a(k,:)=a2;
' h1 `# T7 z9 | }/ o8 u2 z' |
b(k,:)=b2;
7 H0 {8 e, L" }+ s% Q9 l) z2 |3 ?
N(k,:)=N2;
; X. I' q2 e( y, Y, H1 c! i( ]
P=C(k,1);
. x0 X1 e# l/ I! Y8 C2 U
) z: b5 J/ `# j$ B# j# z) e9 v, W
; @* t. A$ F% r4 t& Y# Y( W+ B/ B6 ^
报错
( v0 X$ G+ T$ L' N0 {4 F t
未定义与 'function_handle' 类型的输入参数相对应的运算符 '=='。
& i& e' J3 I8 }/ h6 O
出错 fzero (line 314)
) {( b3 w. h' S
if fx == 0
3 K i) O6 j9 B! S( R# C' S6 J
出错 Price (line 24)
; _* n7 a+ }) I- W- S
Vbar(s-1,1)=fzero('@Price_ss',initial,options,s,Vbar,r,vol,t0,t,K);
# w% Y9 T, `; S9 Y e
8 U9 `( P+ m& Q, d6 `- [
请大神指教
% s0 A" I. B' x; l* w1 ^' s4 V
作者:
shuddkk
时间:
2020-8-14 15:43
帮你顶一下
作者:
kekek
时间:
2020-8-17 10:45
楼主发一下完整代码
3 P, \2 l5 ~. [! j+ }# N
作者:
zzz.dan
时间:
2020-8-18 11:02
作者:
nkkopd
时间:
2020-8-18 13:15
？？？？？？？？？？？？？？？？？？
# L7 [, _8 G6 @& n, b
欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)
Powered by Discuz! X3.2
=hl;4 F" Y7 @! W( I2 ?* X5 w=al;# Z! b) W* Y$ e/ Y X+ j% S$ l2 s/ W9 I=b1;+ y7 f# I; p9 T) }3 N% c7 Q