|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
方程组:
# a# _8 Y& o1 P: S0=-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3))
" I! {. P+ n7 O1 Z3 ]3 J+ b0=-0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3))
6 E; o6 Y M( \( p0=-cos(x(1))*sin(x(2))0 Y5 [3 o* H+ J: ?4 V! A
. L6 U1 L' m1 W) a0 f U, c3 q% S代码如下:
% M W) \( q4 [4 G' `: S* G9 c& W7 j5 t/ v
0 P4 A3 w* m: Z! z& p! m, P4 y5 nfunction F = myfun(x)
! C: Z- U9 L4 ]9 iF = [-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));# r% S2 @( l& u, t V5 Y
-0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));
/ M9 d7 ]) a8 H -cos(x(1))*sin(x(2))];( H W6 u! p% f* K
x0=[0;0;0];
( `5 Y& y) i5 W; j8 [3 goptions=optimset('Dispaly','iter');
5 R8 R1 z* I; I. B# m- |1 P[x,fval]=fsolve(@myfun,x0,options)
i" p! M( T. E) D* j2 Oend/ l! ~. }) t' X9 k; `7 l
9 n3 ^8 r$ T) i* {9 N
; y; ?! i1 }3 ]7 j4 h8 ~/ P
# d$ {5 b6 @2 n6 _" F! O出问题语句在函数段。7 f' J0 g8 Z( U$ o8 E
错误信息:
, I. M8 u( y( x; U6 ^6 Z" ^- d>> myfun
7 i+ U" `* ~2 v输入参数的数目不足。7 U7 }# \. l* g( D; H( M! C
T7 b/ t, l5 ]% x4 `5 R j7 i
出错 myfun (line 2)) M# q9 X6 D# a% A8 ^. j% k
F =/ a& ], ~2 M% [) \
[-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));
' d1 a$ Y. Q z# O8 Y+ g>>1 Q$ t# }& I8 h
, T/ C. ]2 E1 l6 M* b; D7 S( u这应该怎么修改?
; B) x3 H. u/ ]. e8 a1 F* N( U- y |
|