|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
方程组:% h5 }( K) f) p0 M/ G$ x* w
0=-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3)), r# D% b; z `' s
0=-0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3))
0 v. T. m' l( h" X' o' M0=-cos(x(1))*sin(x(2)); S) l( l B0 v4 f) v) u
2 w0 D# _3 u; C8 M U代码如下:
% Y0 j4 I' \1 L3 ?! n3 r4 V
7 U- N; B1 O$ O0 h- a5 C- I* `# D$ W8 E( {* p
function F = myfun(x)9 n, Z; V: ^. x# U; O" F
F = [-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));
! [& x5 r- ?$ v! P9 { -0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));: q. S6 \/ _# y6 x( r; Y
-cos(x(1))*sin(x(2))];0 A5 m# r: ^! [& i; e \$ N
x0=[0;0;0];. e2 b% s6 w* L4 s
options=optimset('Dispaly','iter');* J- V) w5 Y& }2 k
[x,fval]=fsolve(@myfun,x0,options)
# `4 N: N7 M9 F2 Tend
) z5 U2 k. y! u; C+ F5 Z; @$ N4 i$ u+ b) x
" C, y) G# n1 Z( O- r4 Q
/ Q/ i/ d, X5 V1 M3 s/ b# `
出问题语句在函数段。
0 D3 T) l) k. r7 [# ~: ^" Z7 g错误信息:2 V4 E1 T" E0 _% i' ` Y
>> myfun2 A5 m4 U, P4 H* \5 V
输入参数的数目不足。
2 w( l7 \" R4 J8 d5 g
; y7 p) E6 N% d2 v* k' \3 x0 O; P+ K出错 myfun (line 2)# o/ @" T( T7 {4 O% q% h: L
F =
& t% ^( B, Y9 n[-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));$ i' q/ V$ T+ F3 u
>>4 j* s9 f# {# x' \5 H
% i" ~1 W% U8 L5 s" y这应该怎么修改?# l3 u; m6 w5 Y, Z5 A( H
|
|