|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
方程组:
& Z& @ }. b3 I5 p0=-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3))% b [* d7 x; J& P" f" U$ f6 O
0=-0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3))0 A& F0 c& u( i( L" N- C! M& U
0=-cos(x(1))*sin(x(2))5 p' m$ Y% r% c. M: M
5 n* N6 k. h/ ?, Z+ {代码如下: h9 c1 q" \. [' K a6 O
* |: M0 z; z0 p/ [8 e/ l+ H
8 E" [6 f4 b( P# J2 X! \& {
function F = myfun(x)$ b- X( E8 S) `& \* o, p
F = [-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));
2 p/ X, Y8 y7 H/ q8 [- i7 o -0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));
( C8 ?. u Q/ M9 s -cos(x(1))*sin(x(2))];8 |) A) k$ O5 V, W
x0=[0;0;0];
( I( t W" [. W3 j2 T$ Roptions=optimset('Dispaly','iter');
/ E" W- l; n& |; C0 x[x,fval]=fsolve(@myfun,x0,options)" J$ Z+ `. J; i5 {, n( J/ ^# i
end" _6 T; x9 T! f, H' }) R
% o0 l7 _% s# A4 g! ~7 h1 ^
' Q2 ~. }' \2 c. y' k' a, j( c9 m3 n' Q4 e$ X1 H7 K
出问题语句在函数段。
( C+ u+ O" r4 E错误信息:$ F# w( D1 I% L& z: Y
>> myfun$ M2 B7 Y8 R7 e$ Q$ G8 u& G
输入参数的数目不足。4 x: W% M3 b) a3 V4 G" q
3 D) d3 T, D1 O% N7 R. M出错 myfun (line 2)4 m% l" ~) T6 L: M( F9 R0 Z
F =
7 U |$ M" X, j, @) g' V+ g$ p[-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));, U* X5 q. B8 @+ c( ^% q4 P
>>2 F+ V9 S# v/ z+ A Z$ ~; h
4 `6 U% J) H: a" x# N
这应该怎么修改?
7 y- k- V( [) _! | |
|