EDA365电子论坛网
标题:
求解超越方程组
[打印本页]
作者:
zzz.dan
时间:
2020-8-21 09:51
标题:
求解超越方程组
方程组：
" @ o+ @( K& E% o0 ?4 |
0=-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3))
+ \2 `. K0 ^# K1 a: \) _5 D
0=-0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3))
: U0 l. V. E6 S1 O4 q" M
0=-cos(x(1))*sin(x(2))
! {1 I% C' Q. I' U
/ Y8 w, R, x8 w- b" ~! Y
代码如下：
$ P0 h: P @# j1 O3 G1 j
- ]$ n( w+ q$ a
$ q/ R& k1 y5 ?) a9 @- i
function F = myfun(x)
) s; b! U! t8 n @2 L) B* p
F = [-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));
- H7 j/ Y9 b; ]) L8 y- ^
-0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));
( `/ N# H- { I# M
-cos(x(1))*sin(x(2))];
" o% {$ m6 ?# d8 Y' n2 s
x0=[0;0;0];
$ G5 b8 S; O/ r- f3 J0 Q/ G! c. N
options=optimset('Dispaly','iter');
- G; T/ o' ^4 o: u/ M
[x,fval]=fsolve(@myfun,x0,options)
( Z4 n! H: r- `% @5 j' d
end
8 {/ P! \) R; I8 [/ c' [$ @* l% K
4 Y' N8 A( B3 f% m; L! Q0 F
" W9 ^# o/ g* Q$ x) Y% W" o0 j: V
^2 b: @( j) G5 j
出问题语句在函数段。
: {7 b k- M* v3 a9 \" N- x- d
错误信息：
2 `: Q( U; G5 Y! V7 i; A
>> myfun
4 D+ M3 M. x/ _$ a# f \3 ?/ ~
输入参数的数目不足。
$ w9 S8 s6 z0 y) w, r5 b* F
- I4 A( j1 {) P3 B* [8 }" o0 s
出错 myfun (line 2)
/ h' O ]. Y8 |: ^
F =
! Z; s& ~! T6 Y* A' }, E! G1 k, L
[-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));
: q" ]2 J5 Y9 D: b b0 Z! q' J# U
>>
- O: D9 Y. ?# y
% m: [: \, ?6 c7 s: N- z& x
这应该怎么修改？
5 ]- d, D4 c) K0 N9 w& C- M3 b
作者:
Uifhjvv
时间:
2020-8-21 10:36
帮你顶一下
作者:
llbnmo
时间:
2020-8-24 13:46
F = @(x) [-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));
/ |/ y8 \: |. ]: |1 U
-0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));
, w! c8 ]& L: Z! T
-cos(x(1))*sin(x(2))];
+ j6 H5 `0 C. E. n8 x
x0=[0;0;0];
. ^+ _* M& @2 Z) ^* W N3 w
options=optimset('Display','iter');
# y$ ]2 c# C- w+ X8 N
[x,fval]=fsolve(F,x0,options)
作者:
zaiyiaaaa
时间:
2020-8-25 10:51
大佬，问题解决了吗，我的方程tan(x*y)-x*y=x^3(tany-y), y的取值范围为15度到90度。请问一下这个图怎么画啊？或者求x的解。请教一下大佬，万分感谢
作者:
kekek
时间:
2020-8-25 11:21
来学习一下
欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)
Powered by Discuz! X3.2