|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
方程组:0 ?2 K' `) n$ {7 ]
0=-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3))4 Z2 T; v1 Q I" u; J
0=-0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3))6 l2 h& n7 @. X8 k" c9 ]
0=-cos(x(1))*sin(x(2))1 u- n. Q! d7 L
5 I) z9 ~& E7 d9 M- m9 _) E
代码如下:# I7 r* N9 Z8 ~# Q: O+ ]! Q
2 ] D* |: }- G% i2 ?
; M( J3 u2 p$ E% ^' l6 Z8 n! ^3 S1 ~' R& M
function F = myfun(x). A9 _1 O$ H' u* M0 J
F = [-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));7 V2 {3 I% t$ T! C- _
-0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));
3 P& u |& D" `' N: p+ g, a -cos(x(1))*sin(x(2))];
0 U" n5 X" C: o2 a& ]x0=[0;0;0];
+ K8 c3 k* d7 o! G* Aoptions=optimset('Dispaly','iter');5 O/ W O: S! p- \; s& ^
[x,fval]=fsolve(@myfun,x0,options)( M! G; J* U/ \8 ^( }
end. X, K3 O# ^0 B0 a9 |% F
6 \2 b) D; y, n: G% r! y% @/ S
2 e$ i, ^' c% v% N1 \3 y% L
# T1 Q" X3 I+ l+ e出问题语句在函数段。; p, [! } Q' {: t- r& v, ~2 Z
错误信息:
/ k: b6 c+ ~& `$ e9 o `% W7 J- N" B>> myfun* B6 T( r3 Z2 i/ p7 L/ @5 C* }/ T- [* d
输入参数的数目不足。" `% b: `6 B6 w# M; I' L% K' _
- R/ t7 S! |* P# U9 H( b& J. X
出错 myfun (line 2)
/ | b1 Y) T$ |/ K; l' n, u8 hF = F- N; L3 Z0 l6 ?7 v- t
[-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));
D8 E- B0 q- n4 F8 o1 U( B7 G>>0 Y9 L7 w* t! t) H& d
" F) q5 l* s$ B( [1 E- u0 O这应该怎么修改?
6 r3 X9 B. M- ~ |
|