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

求解超越方程组

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
方程组:
9 }* b8 G: ~4 ]2 p' \3 T# K0=-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3))
8 e( A5 U, u/ S8 E0=-0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3))
3 K  a4 C2 {0 F4 y0=-cos(x(1))*sin(x(2))
: ~( A& M1 k, S* ^) j. S: U/ U" H4 ?, x$ L- ^& c$ @
代码如下:
# g" t( L" b8 F7 n# g  s
% J9 T+ W5 r5 P
; a- ^! P4 N4 r9 n! J$ Kfunction F = myfun(x)' z3 V$ m% `: _+ P
F = [-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));
+ ?7 Q6 r: G& K+ M3 S' y    -0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));
2 ^, ]  P! Q$ e1 |7 h* p    -cos(x(1))*sin(x(2))];+ _6 o+ t: }$ S6 [# G( e* k# k. R/ @
x0=[0;0;0];
* X$ b9 k7 k  o$ Hoptions=optimset('Dispaly','iter');$ G. s- A1 G- \
[x,fval]=fsolve(@myfun,x0,options)
1 d2 S" q' h1 v0 y0 K/ N' r+ Dend6 h* y* Z% Z' j# ]5 ^) Y* N

; f. {+ j1 u3 V# h
4 n- l4 s: d7 E* L5 e
! v8 d' S7 W" K出问题语句在函数段。; f) I1 i# a: I* x# h* ^. i
错误信息:
$ N$ O3 C, c$ H. B' s* C>> myfun
; a! D/ y# `- R. Z输入参数的数目不足。; L- }* M& w0 g4 \
& S6 G$ u) G) k. ~
出错 myfun (line 2)
# S; D0 Q! P7 j- zF =
& d! Y' ~$ c$ o: e. y: r- F; ~7 ?6 R[-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));
6 Q3 Z7 k% _+ ^5 j  V2 i>>
! b4 P5 P6 ]: h6 ]  w; |' @/ T! F: L3 c6 e
这应该怎么修改?
) Q& O0 {& f1 `( w9 s4 D8 }

该用户从未签到

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));0 E" g5 o  V  A  h2 _- R: q
    -0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));
" c" H4 P" [! o+ {. Y    -cos(x(1))*sin(x(2))];
: F' H8 V6 w8 d' X' Fx0=[0;0;0];
# S0 k, i3 Z" S, b2 g5 Zoptions=optimset('Display','iter');
0 c9 _: W% c8 w) Z0 G2 \[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 19:12 , Processed in 0.171875 second(s), 23 queries , Gzip On.

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

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

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