|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
方程组:
9 }* b8 G: ~4 ]2 p' \3 T# K0=-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3))
8 e( A5 U, u/ S8 E0=-0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3))
3 K a4 C2 {0 F4 y0=-cos(x(1))*sin(x(2))
: ~( A& M1 k, S* ^) j. S: U/ U" H4 ?, x$ L- ^& c$ @
代码如下:
# g" t( L" b8 F7 n# g s
% J9 T+ W5 r5 P
; a- ^! P4 N4 r9 n! J$ Kfunction F = myfun(x)' z3 V$ m% `: _+ P
F = [-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));
+ ?7 Q6 r: G& K+ M3 S' y -0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));
2 ^, ] P! Q$ e1 |7 h* p -cos(x(1))*sin(x(2))];+ _6 o+ t: }$ S6 [# G( e* k# k. R/ @
x0=[0;0;0];
* X$ b9 k7 k o$ Hoptions=optimset('Dispaly','iter');$ G. s- A1 G- \
[x,fval]=fsolve(@myfun,x0,options)
1 d2 S" q' h1 v0 y0 K/ N' r+ Dend6 h* y* Z% Z' j# ]5 ^) Y* N
; f. {+ j1 u3 V# h
4 n- l4 s: d7 E* L5 e
! v8 d' S7 W" K出问题语句在函数段。; f) I1 i# a: I* x# h* ^. i
错误信息:
$ N$ O3 C, c$ H. B' s* C>> myfun
; a! D/ y# `- R. Z输入参数的数目不足。; L- }* M& w0 g4 \
& S6 G$ u) G) k. ~
出错 myfun (line 2)
# S; D0 Q! P7 j- zF =
& d! Y' ~$ c$ o: e. y: r- F; ~7 ?6 R[-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));
6 Q3 Z7 k% _+ ^5 j V2 i>>
! b4 P5 P6 ]: h6 ] w; |' @/ T! F: L3 c6 e
这应该怎么修改?
) Q& O0 {& f1 `( w9 s4 D8 } |
|