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

求解超越方程组

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
方程组:0 ?2 K' `) n$ {7 ]
0=-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3))4 Z2 T; v1 Q  I" u; J
0=-0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3))6 l2 h& n7 @. X8 k" c9 ]
0=-cos(x(1))*sin(x(2))1 u- n. Q! d7 L
5 I) z9 ~& E7 d9 M- m9 _) E
代码如下:# I7 r* N9 Z8 ~# Q: O+ ]! Q
2 ]  D* |: }- G% i2 ?
; M( J3 u2 p$ E% ^' l6 Z8 n! ^3 S1 ~' R& M
function F = myfun(x). A9 _1 O$ H' u* M0 J
F = [-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));7 V2 {3 I% t$ T! C- _
    -0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));
3 P& u  |& D" `' N: p+ g, a    -cos(x(1))*sin(x(2))];
0 U" n5 X" C: o2 a& ]x0=[0;0;0];
+ K8 c3 k* d7 o! G* Aoptions=optimset('Dispaly','iter');5 O/ W  O: S! p- \; s& ^
[x,fval]=fsolve(@myfun,x0,options)( M! G; J* U/ \8 ^( }
end. X, K3 O# ^0 B0 a9 |% F
6 \2 b) D; y, n: G% r! y% @/ S

2 e$ i, ^' c% v% N1 \3 y% L
# T1 Q" X3 I+ l+ e出问题语句在函数段。; p, [! }  Q' {: t- r& v, ~2 Z
错误信息:
/ k: b6 c+ ~& `$ e9 o  `% W7 J- N" B>> myfun* B6 T( r3 Z2 i/ p7 L/ @5 C* }/ T- [* d
输入参数的数目不足。" `% b: `6 B6 w# M; I' L% K' _
- R/ t7 S! |* P# U9 H( b& J. X
出错 myfun (line 2)
/ |  b1 Y) T$ |/ K; l' n, u8 hF =  F- N; L3 Z0 l6 ?7 v- t
[-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));
  D8 E- B0 q- n4 F8 o1 U( B7 G>>0 Y9 L7 w* t! t) H& d

" F) q5 l* s$ B( [1 E- u0 O这应该怎么修改?
6 r3 X9 B. M- ~

该用户从未签到

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));4 a0 q8 k9 Z3 _" g" k$ j
    -0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));
6 m9 t8 N& u7 z! {    -cos(x(1))*sin(x(2))];
$ c3 o/ W+ e6 ^+ e: rx0=[0;0;0];2 `% c& i; x* b
options=optimset('Display','iter');
5 \* A( r7 F7 q; x. O3 j, P[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-6-25 04:29 , Processed in 0.078125 second(s), 23 queries , Gzip On.

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

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

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