|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
方程组:
& C- |& n6 x4 {& Q& y# D0=-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3))+ `. N$ o2 L6 s+ S, a
0=-0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3))
* R& Q; t1 n9 G0 n0=-cos(x(1))*sin(x(2))3 d. |' f7 U/ U; w7 r
7 J5 _* p+ f. R D. n
代码如下:
" d+ J* C7 Z2 o0 x
% w0 S6 x4 G9 {/ w9 i6 {2 H9 w* J" }/ x2 W ?4 ^! [0 O& o& A3 w8 y
function F = myfun(x)
8 ^4 ~, ^' n# z% Q1 cF = [-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));
8 y# h g& |; c: @. G -0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));' D& f* Z/ `& e
-cos(x(1))*sin(x(2))];
- I7 ~* i# P) L/ b9 I& _9 e6 S9 gx0=[0;0;0];$ H9 I) b9 m3 y* _7 j5 p& W3 g, E% d
options=optimset('Dispaly','iter');; [" r- R) ?* J* e( h9 d$ M7 C
[x,fval]=fsolve(@myfun,x0,options)
. S3 ~5 A7 \1 k9 `9 c1 O7 z, xend# r' @4 A2 w8 Q v
8 \" R, v( O5 ]
5 P. p4 w! ^' w. J, A* V6 l
, E; ?/ H2 W; ^出问题语句在函数段。
, L- i0 ^+ k4 l# R7 t' y, c4 q* |错误信息:
5 |9 ~9 `6 a- k# y2 N5 |>> myfun
6 v" h/ ]' y2 X) x5 I C1 ?输入参数的数目不足。
: _' d7 H3 Q0 Z F
- g& V U3 c& M出错 myfun (line 2)9 f/ o: \0 W) `: d7 Z
F =
6 z8 J2 E: L3 N, h[-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));& T& p" q: u) D1 |9 U0 l
>>
& i# r& r" X1 q: i# i% A3 R0 p, e+ D
4 A2 p$ C- N6 { R: d这应该怎么修改?5 E& ]8 y. H) Q, Z/ f( Z& }& R
|
|