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

求解超越方程组

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
方程组:
& C- |& n6 x4 {& Q& y# D0=-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3))+ `. N$ o2 L6 s+ S, a
0=-0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3))
* R& Q; t1 n9 G0 n0=-cos(x(1))*sin(x(2))3 d. |' f7 U/ U; w7 r
7 J5 _* p+ f. R  D. n
代码如下:
" d+ J* C7 Z2 o0 x
% w0 S6 x4 G9 {/ w9 i6 {2 H9 w* J" }/ x2 W  ?4 ^! [0 O& o& A3 w8 y
function F = myfun(x)
8 ^4 ~, ^' n# z% Q1 cF = [-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));
8 y# h  g& |; c: @. G    -0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));' D& f* Z/ `& e
    -cos(x(1))*sin(x(2))];
- I7 ~* i# P) L/ b9 I& _9 e6 S9 gx0=[0;0;0];$ H9 I) b9 m3 y* _7 j5 p& W3 g, E% d
options=optimset('Dispaly','iter');; [" r- R) ?* J* e( h9 d$ M7 C
[x,fval]=fsolve(@myfun,x0,options)
. S3 ~5 A7 \1 k9 `9 c1 O7 z, xend# r' @4 A2 w8 Q  v
8 \" R, v( O5 ]
5 P. p4 w! ^' w. J, A* V6 l

, E; ?/ H2 W; ^出问题语句在函数段。
, L- i0 ^+ k4 l# R7 t' y, c4 q* |错误信息:
5 |9 ~9 `6 a- k# y2 N5 |>> myfun
6 v" h/ ]' y2 X) x5 I  C1 ?输入参数的数目不足。
: _' d7 H3 Q0 Z  F
- g& V  U3 c& M出错 myfun (line 2)9 f/ o: \0 W) `: d7 Z
F =
6 z8 J2 E: L3 N, h[-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));& T& p" q: u) D1 |9 U0 l
>>
& i# r& r" X1 q: i# i% A3 R0 p, e+ D
4 A2 p$ C- N6 {  R: d这应该怎么修改?5 E& ]8 y. H) Q, Z/ f( Z& }& R

该用户从未签到

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));
6 ~+ G6 e  o; F    -0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));' b: t  y. @- e3 B
    -cos(x(1))*sin(x(2))];
/ q: q5 ?: t" ~x0=[0;0;0];$ b# a# H, K- z; k: ~* @# ^
options=optimset('Display','iter');% z' j8 Y1 b! W( }0 j5 |
[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-8-16 06:39 , Processed in 0.125000 second(s), 23 queries , Gzip On.

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

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

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