|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
方程组:7 s2 s) ]3 z6 l
0=-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3))
_: q( ~- |! v0=-0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3))! K) I% J% c, l
0=-cos(x(1))*sin(x(2))$ l! `% ]7 T P) v, q3 z) D
, T- r' U. L& P5 F3 P7 u% g9 f- A
代码如下:
9 b* W, D/ O4 Z" @' J
" U8 D, s% Q! R0 N/ m" H" h% u3 {
function F = myfun(x)5 ~7 G; \: ^5 p- Y7 w
F = [-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3)); s8 j" q2 h% l0 U# S7 L3 e- |
-0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));) u8 w6 G: Y. P0 u; ?% ^
-cos(x(1))*sin(x(2))];
. N ], b c, Y9 q; b' v% t' \+ P/ Ux0=[0;0;0];
. r3 x4 A0 z: g* X3 Doptions=optimset('Dispaly','iter');
; [. G- Z$ t/ _3 v' I' f[x,fval]=fsolve(@myfun,x0,options)
+ Y% |: W, G* a7 C; [! R9 Rend, g. E' u& A; ~- s" A$ v
! s3 X. J* r8 i2 Z1 v
# q: @4 S2 v( U6 E$ I# \
- E) m3 R1 u" T: h9 N
出问题语句在函数段。
) m) u" N/ N% f% x& U错误信息:' R3 i" p4 t5 d- y1 O1 C3 A
>> myfun9 ~$ |- A) y+ \9 ]4 R; z
输入参数的数目不足。
$ ? ?/ e* ]0 E
1 _9 J6 I7 n: y# \$ `1 J9 M, ]+ r' T出错 myfun (line 2)
1 o7 ]7 ?: c) b% s5 CF =5 E. Y. a$ |9 l
[-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));
+ b" I1 K! i$ h. v! v9 u, l>>
6 ~% @. E. b" A( P- v2 ~1 l0 ~% O5 ^9 }$ h
这应该怎么修改?
8 z/ v$ n P, Y; w |
|