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

求解超越方程组

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
方程组:% h5 }( K) f) p0 M/ G$ x* w
0=-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3)), r# D% b; z  `' s
0=-0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3))
0 v. T. m' l( h" X' o' M0=-cos(x(1))*sin(x(2)); S) l( l  B0 v4 f) v) u

2 w0 D# _3 u; C8 M  U代码如下:
% Y0 j4 I' \1 L3 ?! n3 r4 V
7 U- N; B1 O$ O0 h- a5 C- I* `# D$ W8 E( {* p
function F = myfun(x)9 n, Z; V: ^. x# U; O" F
F = [-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));
! [& x5 r- ?$ v! P9 {    -0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));: q. S6 \/ _# y6 x( r; Y
    -cos(x(1))*sin(x(2))];0 A5 m# r: ^! [& i; e  \$ N
x0=[0;0;0];. e2 b% s6 w* L4 s
options=optimset('Dispaly','iter');* J- V) w5 Y& }2 k
[x,fval]=fsolve(@myfun,x0,options)
# `4 N: N7 M9 F2 Tend
) z5 U2 k. y! u; C+ F5 Z; @$ N4 i$ u+ b) x
" C, y) G# n1 Z( O- r4 Q
/ Q/ i/ d, X5 V1 M3 s/ b# `
出问题语句在函数段。
0 D3 T) l) k. r7 [# ~: ^" Z7 g错误信息:2 V4 E1 T" E0 _% i' `  Y
>> myfun2 A5 m4 U, P4 H* \5 V
输入参数的数目不足。
2 w( l7 \" R4 J8 d5 g
; y7 p) E6 N% d2 v* k' \3 x0 O; P+ K出错 myfun (line 2)# o/ @" T( T7 {4 O% q% h: L
F =
& t% ^( B, Y9 n[-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));$ i' q/ V$ T+ F3 u
>>4 j* s9 f# {# x' \5 H

% i" ~1 W% U8 L5 s" y这应该怎么修改?# l3 u; m6 w5 Y, Z5 A( H

该用户从未签到

5#
发表于 2020-8-25 11:21 | 只看该作者
来学习一下

该用户从未签到

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

该用户从未签到

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));1 N+ r: {3 K  d
    -0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));9 @  w4 O" q; M  i& C8 W4 v7 D9 B
    -cos(x(1))*sin(x(2))];
! F& D, i$ J5 J! fx0=[0;0;0];
) \/ f  G3 L. w7 |4 W4 z3 Z- ]! @options=optimset('Display','iter');
* p* x7 T# W. z! e! H/ c[x,fval]=fsolve(F,x0,options)

该用户从未签到

2#
发表于 2020-8-21 10:36 | 只看该作者
帮你顶一下
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

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

EDA365公众号

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

GMT+8, 2025-11-24 17:43 , Processed in 0.218750 second(s), 24 queries , Gzip On.

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

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

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