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

求解超越方程组

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
方程组:# Q) H# Q+ M6 L
0=-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3))
7 P$ s' ?/ W3 \$ V: o( B) ^0=-0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3))) _9 B6 q% x% }, r% J
0=-cos(x(1))*sin(x(2))
- A: ^- N, J" m, f" W/ Y
' U2 V+ {$ R# L# E( {4 D代码如下:
/ U- q3 x( _7 Y: n! z0 d- M2 ^. _
. R2 `" v1 S, s6 {6 f+ B' `9 U2 T  s+ _  X' B
function F = myfun(x)
* v  x; [- x5 O4 ]F = [-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));+ _: w- x) Y$ s  U* }( V
    -0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));) J1 N6 g' u8 q$ U1 B, y" f
    -cos(x(1))*sin(x(2))];0 D2 x9 F6 A# C/ B
x0=[0;0;0];/ n% [- b5 S+ `5 E5 t; h. e; z
options=optimset('Dispaly','iter');; `0 B/ J3 f- v6 Y$ C
[x,fval]=fsolve(@myfun,x0,options)
; l8 n1 N+ o& Qend4 g) o4 Q& m5 }: m- |# D
3 ], E. A& F! a0 d
0 b1 P1 Y5 i" e4 b
! b5 K2 H2 d" \2 k* d
出问题语句在函数段。
7 j9 a, F6 h& @- k' _* _) u错误信息:1 j5 y+ ^# ~4 |! w# e
>> myfun- P4 @6 i- N, {% e2 E* I0 P
输入参数的数目不足。
) e- n# B* R7 O! R- z$ W' k- U( z7 W! A+ P* v! ^
出错 myfun (line 2)
% c3 K$ ~5 b0 s6 b0 V+ v; DF =( y$ O( J: m+ U0 G- h: v2 W( \
[-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));
7 q8 J# w9 b4 [+ z6 Q5 b>>
% `, P+ d. Q0 g4 K$ e
# V8 `+ ~, U# d这应该怎么修改?
9 w2 g" }' j# A1 F9 a+ W6 {7 k

该用户从未签到

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));
; n  C" R# Q' C/ i/ R    -0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));
, f6 y% b& [: _+ ]8 ^0 B, c    -cos(x(1))*sin(x(2))];8 J3 Z  |) i8 m# \! o% v  C+ b
x0=[0;0;0];
( v1 a- \% S$ ]* poptions=optimset('Display','iter');+ F3 c8 I6 I9 G% O; }  y: E
[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-5 13:07 , Processed in 0.140625 second(s), 23 queries , Gzip On.

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

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

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