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

求解超越方程组

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
方程组:
& Z& @  }. b3 I5 p0=-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3))% b  [* d7 x; J& P" f" U$ f6 O
0=-0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3))0 A& F0 c& u( i( L" N- C! M& U
0=-cos(x(1))*sin(x(2))5 p' m$ Y% r% c. M: M

5 n* N6 k. h/ ?, Z+ {代码如下:  h9 c1 q" \. [' K  a6 O
* |: M0 z; z0 p/ [8 e/ l+ H
8 E" [6 f4 b( P# J2 X! \& {
function F = myfun(x)$ b- X( E8 S) `& \* o, p
F = [-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));
2 p/ X, Y8 y7 H/ q8 [- i7 o    -0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));
( C8 ?. u  Q/ M9 s    -cos(x(1))*sin(x(2))];8 |) A) k$ O5 V, W
x0=[0;0;0];
( I( t  W" [. W3 j2 T$ Roptions=optimset('Dispaly','iter');
/ E" W- l; n& |; C0 x[x,fval]=fsolve(@myfun,x0,options)" J$ Z+ `. J; i5 {, n( J/ ^# i
end" _6 T; x9 T! f, H' }) R

% o0 l7 _% s# A4 g! ~7 h1 ^
' Q2 ~. }' \2 c. y' k' a, j( c9 m3 n' Q4 e$ X1 H7 K
出问题语句在函数段。
( C+ u+ O" r4 E错误信息:$ F# w( D1 I% L& z: Y
>> myfun$ M2 B7 Y8 R7 e$ Q$ G8 u& G
输入参数的数目不足。4 x: W% M3 b) a3 V4 G" q

3 D) d3 T, D1 O% N7 R. M出错 myfun (line 2)4 m% l" ~) T6 L: M( F9 R0 Z
F =
7 U  |$ M" X, j, @) g' V+ g$ p[-0.3+cos(x(1))*cos(x(2))*cos(x(3))-sin(x(1))*sin(x(3));, U* X5 q. B8 @+ c( ^% q4 P
>>2 F+ V9 S# v/ z+ A  Z$ ~; h
4 `6 U% J) H: a" x# N
这应该怎么修改?
7 y- k- V( [) _! |

该用户从未签到

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));$ m+ y0 @5 u" h* _3 V1 r3 p
    -0.4+cos(x(1))*cos(x(2))*sin(x(3))+sin(x(1))*cos(x(3));; L& \# m' P4 l, L- D
    -cos(x(1))*sin(x(2))];, V: J+ Z1 p# R' V& O6 _4 H
x0=[0;0;0];/ Q2 y. B+ o! c/ V# Z( S: t
options=optimset('Display','iter');' Z  J1 m! |6 g
[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 17:25 , Processed in 0.171875 second(s), 23 queries , Gzip On.

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

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

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