|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
方程组:
/ a$ a3 v5 t8 c3 J* a0=-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3))
5 w; p9 P% h9 z4 }0=-0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3))! W$ v1 q4 M$ a( x, E. s
0=-cos(x(1))*sin(x(2))
. w" ~$ l" y) m7 d! \3 O: l( m" H% H' h
代码如下:& `" w$ m) B5 V2 G: P
' f2 A( m* T9 v, Q4 K4 t
1 d$ B1 u6 ~/ f$ Y9 h6 `
function F = myfun(x)
* ?* U7 O4 a0 o# _3 b' n+ z( k w/ [F = [-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));+ q ]. B$ R" \$ [: q* H8 @
-0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));! o: }, r. b0 A1 a4 U" m* z5 p
-cos(x(1))*sin(x(2))];
8 {) M( M: ]7 W2 o8 t1 s) v2 b" g7 r; Rx0=[0;0;0];0 U3 b, F; ^3 Y/ r
options=optimset('Dispaly','iter');! E- t9 z c5 s7 r S
[x,fval]=fsolve(@myfun,x0,options)
0 c( m8 N- x$ I3 oend( L/ Z* }) B3 T6 U
: c; M6 l- }3 X. J3 D3 D0 M
; I; O P6 C2 F* Q. P# }' S& `) r& h4 ^/ G6 k: G3 n3 _4 q
出问题语句在函数段。
* R1 C( K$ o4 f, q, q错误信息:
3 h U+ v" r( F1 @( Q* V>> myfun
- H5 }; c* _$ d8 J& S- A输入参数的数目不足。
9 ?7 l. Z/ `& E7 q+ R7 w) N! \) H2 B+ M; `/ l
出错 myfun (line 2)- [' A8 |9 p4 D; K% M3 ]8 h6 r
F =5 N' y( L8 i. k! L9 X& o% H- R
[-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));4 K8 `& o; l/ X# f
>>' ]; w1 C) |# U0 e1 q
& m+ Q6 R8 g) ?1 j这应该怎么修改?8 s. R; c/ |8 W* }, r; e
|
|