|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
方程组:
( g- w$ U, W% B8 h, d) G0=-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3))8 M0 F: |( g6 o% V! e& L7 q+ X
0=-0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3))
7 D/ t9 h) L: s# \3 v8 m0=-cos(x(1))*sin(x(2))
; e7 g% ^+ i& r* e. g# c" ~# w
# o( |& t5 v; H8 e4 |8 N& j代码如下:' Q: ]6 Z6 ~, K4 n9 n* v0 C
3 I2 l' s* ^1 d* ?& `
/ }6 _+ Q1 l* F% D2 @
function F = myfun(x)
! U$ D& D0 Z0 mF = [-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));
# J0 \1 E& c$ s4 O$ A! H -0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));6 @) p& H) F9 U" T3 k1 ]
-cos(x(1))*sin(x(2))];3 w4 H3 C1 `: U, j, X# C$ u
x0=[0;0;0];( L. h& C1 e: M# W
options=optimset('Dispaly','iter');
0 m, `) n8 f9 W& x# K8 o q[x,fval]=fsolve(@myfun,x0,options)
# X! A e; D" b+ ^8 rend
, \/ P# w) L9 I* @+ {# k
+ X# a, L6 O4 D! u
5 o+ k3 H3 t* N* o H+ d3 T) o. c5 b( `" G2 G/ k
出问题语句在函数段。6 ?' ~' a) y/ p8 n; g) d' w* N. ^+ x
错误信息:
( d" N& p5 z/ E( l" J1 v>> myfun8 b" s6 c+ \9 r, N% m
输入参数的数目不足。# D' Q+ }- C/ _( D
5 _4 _4 x" P& [& ?出错 myfun (line 2)$ d4 d% a& L) V: u) B h
F =. d# x, q( }$ v, w! x0 m$ M2 T
[-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));1 {" j/ l4 C, a$ O8 x
>>
7 p& ?! U3 d' \+ X6 k0 \ S3 Z7 @/ Q3 }) T8 T
这应该怎么修改?: ]- b( b O. `. }+ T' A9 {6 s/ _
|
|