找回密码
 注册
关于网站域名变更的通知
查看: 510|回复: 4
打印 上一主题 下一主题

求解超越方程组

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2020-8-21 09:51 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
方程组:7 s2 s) ]3 z6 l
0=-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3))
  _: q( ~- |! v0=-0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3))! K) I% J% c, l
0=-cos(x(1))*sin(x(2))$ l! `% ]7 T  P) v, q3 z) D
, T- r' U. L& P5 F3 P7 u% g9 f- A
代码如下:
9 b* W, D/ O4 Z" @' J
" U8 D, s% Q! R0 N/ m" H" h% u3 {
function F = myfun(x)5 ~7 G; \: ^5 p- Y7 w
F = [-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));  s8 j" q2 h% l0 U# S7 L3 e- |
    -0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));) u8 w6 G: Y. P0 u; ?% ^
    -cos(x(1))*sin(x(2))];
. N  ], b  c, Y9 q; b' v% t' \+ P/ Ux0=[0;0;0];
. r3 x4 A0 z: g* X3 Doptions=optimset('Dispaly','iter');
; [. G- Z$ t/ _3 v' I' f[x,fval]=fsolve(@myfun,x0,options)
+ Y% |: W, G* a7 C; [! R9 Rend, g. E' u& A; ~- s" A$ v
! s3 X. J* r8 i2 Z1 v
# q: @4 S2 v( U6 E$ I# \
- E) m3 R1 u" T: h9 N
出问题语句在函数段。
) m) u" N/ N% f% x& U错误信息:' R3 i" p4 t5 d- y1 O1 C3 A
>> myfun9 ~$ |- A) y+ \9 ]4 R; z
输入参数的数目不足。
$ ?  ?/ e* ]0 E
1 _9 J6 I7 n: y# \$ `1 J9 M, ]+ r' T出错 myfun (line 2)
1 o7 ]7 ?: c) b% s5 CF =5 E. Y. a$ |9 l
[-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));
+ b" I1 K! i$ h. v! v9 u, l>>
6 ~% @. E. b" A( P- v2 ~1 l0 ~% O5 ^9 }$ h
这应该怎么修改?
8 z/ v$ n  P, Y; w

该用户从未签到

2#
发表于 2020-8-21 10:36 | 只看该作者
帮你顶一下

该用户从未签到

3#
发表于 2020-8-24 13:46 | 只看该作者
F = @(x) [-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));/ p2 V2 [9 O7 m* [+ \+ T# t
    -0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));3 q/ n" \$ n9 t# Z/ _
    -cos(x(1))*sin(x(2))];& G6 J# ?; ?+ p
x0=[0;0;0];; F& W0 G3 [% i! {5 ]7 t5 L0 G
options=optimset('Display','iter');3 U, u9 A. g7 ~' _4 I8 |7 ]) ?# u
[x,fval]=fsolve(F,x0,options)

该用户从未签到

4#
发表于 2020-8-25 10:51 | 只看该作者
大佬,问题解决了吗,我的方程tan(x*y)-x*y=x^3(tany-y),   y的取值范围为15度到90度。请问一下这个图怎么画啊?或者求x的解。请教一下大佬,万分感谢

该用户从未签到

5#
发表于 2020-8-25 11:21 | 只看该作者
来学习一下
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-11-24 22:08 , Processed in 0.171875 second(s), 23 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表