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

求解超越方程组

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
方程组:
# a# _8 Y& o1 P: S0=-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3))
" I! {. P+ n7 O1 Z3 ]3 J+ b0=-0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3))
6 E; o6 Y  M( \( p0=-cos(x(1))*sin(x(2))0 Y5 [3 o* H+ J: ?4 V! A

. L6 U1 L' m1 W) a0 f  U, c3 q% S代码如下:
% M  W) \( q4 [4 G' `: S* G9 c& W7 j5 t/ v

0 P4 A3 w* m: Z! z& p! m, P4 y5 nfunction F = myfun(x)
! C: Z- U9 L4 ]9 iF = [-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));# r% S2 @( l& u, t  V5 Y
    -0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));
/ M9 d7 ]) a8 H    -cos(x(1))*sin(x(2))];( H  W6 u! p% f* K
x0=[0;0;0];
( `5 Y& y) i5 W; j8 [3 goptions=optimset('Dispaly','iter');
5 R8 R1 z* I; I. B# m- |1 P[x,fval]=fsolve(@myfun,x0,options)
  i" p! M( T. E) D* j2 Oend/ l! ~. }) t' X9 k; `7 l
9 n3 ^8 r$ T) i* {9 N
; y; ?! i1 }3 ]7 j4 h8 ~/ P

# d$ {5 b6 @2 n6 _" F! O出问题语句在函数段。7 f' J0 g8 Z( U$ o8 E
错误信息:
, I. M8 u( y( x; U6 ^6 Z" ^- d>> myfun
7 i+ U" `* ~2 v输入参数的数目不足。7 U7 }# \. l* g( D; H( M! C
  T7 b/ t, l5 ]% x4 `5 R  j7 i
出错 myfun (line 2)) M# q9 X6 D# a% A8 ^. j% k
F =/ a& ], ~2 M% [) \
[-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));
' d1 a$ Y. Q  z# O8 Y+ g>>1 Q$ t# }& I8 h

, T/ C. ]2 E1 l6 M* b; D7 S( u这应该怎么修改?
; B) x3 H. u/ ]. e8 a1 F* N( U- y

该用户从未签到

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));
$ o% s  V1 X4 ^2 Y    -0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));
  z8 f9 w3 E5 p5 [2 P    -cos(x(1))*sin(x(2))];
, T- G% h: d$ U: t1 P' zx0=[0;0;0];
! `+ h( f' j9 a. Soptions=optimset('Display','iter');
  r# S  c/ H4 ~* X$ C[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-7-27 14:36 , Processed in 0.109375 second(s), 23 queries , Gzip On.

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

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

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