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

求解超越方程组

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
方程组:
/ a$ a3 v5 t8 c3 J* a0=-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3))
5 w; p9 P% h9 z4 }0=-0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3))! W$ v1 q4 M$ a( x, E. s
0=-cos(x(1))*sin(x(2))
. w" ~$ l" y) m7 d! \3 O: l( m" H% H' h
代码如下:& `" w$ m) B5 V2 G: P
' f2 A( m* T9 v, Q4 K4 t
1 d$ B1 u6 ~/ f$ Y9 h6 `
function F = myfun(x)
* ?* U7 O4 a0 o# _3 b' n+ z( k  w/ [F = [-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));+ q  ]. B$ R" \$ [: q* H8 @
    -0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));! o: }, r. b0 A1 a4 U" m* z5 p
    -cos(x(1))*sin(x(2))];
8 {) M( M: ]7 W2 o8 t1 s) v2 b" g7 r; Rx0=[0;0;0];0 U3 b, F; ^3 Y/ r
options=optimset('Dispaly','iter');! E- t9 z  c5 s7 r  S
[x,fval]=fsolve(@myfun,x0,options)
0 c( m8 N- x$ I3 oend( L/ Z* }) B3 T6 U

: c; M6 l- }3 X. J3 D3 D0 M
; I; O  P6 C2 F* Q. P# }' S& `) r& h4 ^/ G6 k: G3 n3 _4 q
出问题语句在函数段。
* R1 C( K$ o4 f, q, q错误信息:
3 h  U+ v" r( F1 @( Q* V>> myfun
- H5 }; c* _$ d8 J& S- A输入参数的数目不足。
9 ?7 l. Z/ `& E7 q+ R7 w) N! \) H2 B+ M; `/ l
出错 myfun (line 2)- [' A8 |9 p4 D; K% M3 ]8 h6 r
F =5 N' y( L8 i. k! L9 X& o% H- R
[-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));4 K8 `& o; l/ X# f
>>' ]; w1 C) |# U0 e1 q

& m+ Q6 R8 g) ?1 j这应该怎么修改?8 s. R; c/ |8 W* }, r; e

该用户从未签到

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));
5 F2 o6 L' x* x. b9 W    -0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));) T. g! J4 y  T0 ]. R7 u. W: S0 F
    -cos(x(1))*sin(x(2))];
6 d, m. R# K; f/ Z  P1 Ox0=[0;0;0];
/ o- B: W8 R) I3 y4 roptions=optimset('Display','iter');
% J) d5 ~& U  P) ]! S  F[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 13:51 , Processed in 0.156250 second(s), 23 queries , Gzip On.

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

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

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