声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1359|回复: 7

[综合讨论] 求解一个超越方程组

[复制链接]
发表于 2008-9-28 18:46 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
求解一个超越方程组,要求每给一个rx(41.35<rx<61.35,间隔1取值),就得到对应的(t,R)。大侠们,救命啊!
clear all
syms a b c d e A B K H P Ri R rx t C1 C2 C3 C4 C5 C6 C7 C8 dR dt M1 M2 M3 M4
a=20*pi/180
c=pi/6
d=(4+269/60)*pi/180
A=40
K=5.5
H=0
B=10
P=5
Ri=20
C1=sin(c)*cos(d)*sin(e)-sin(d)*cos(e)
C2=cos(c)*cos(b)*sin(e)+sin(c)*sin(b)*sin(d)*sin(e)+sin(b)*cos(d)*cos(e)
C3=(A-K*e)*sin(e)+H*cos(e)
C4=sin(c)*cos(d)*cos(e)+sin(d)*sin(e)-cos(c)*cos(d)/tan(a)
C5=cos(c)*cos(b)*cos(e)+sin(c)*sin(b)*sin(d)*cos(e)-sin(b)*cos(d)*sin(e)+(sin(c)*cos(b)-cos(c)*sin(b)*sin(d))/tan(a)
C6=(A-K*e)*cos(e)-H*sin(e)+(B-P*e)/tan(a)-Ri
C7=((A-K*e)*sin(c)*sin(d)+H*cos(d)-(B-P*e)*cos(c)*sin(d)+rx*(cos(d)*sin(e)-sin(c)*sin(d)*cos(e)))/tan(a)+(Ri-rx)*cos(c)*sin(d)
C8=((A-K*e)*cos(c)+(B-P*e)*sin(c)-rx*cos(c)*cos(e))/tan(a)-(Ri-rx)*sin(c)
b=atan(C7/C8)+pi
R=(C3*C4-C1*C6)/(C1*C5-C2*C4)
t=(C2*C6-C3*C5)/(C1*C5-C2*C4)
dR=diff(R)
dt=diff(t)
(cos(c)*cos(d)*dR/dt+sin(c)*cos(b)-cos(c)*sin(d)*sin(b))*P+(sin(c)*cos(d)*dR/dt-cos(c)*cos(b)-sin(c)*sin(d)*sin(b))*(sin(d)*t-cos(d)*R*sin(b)-H-K)+(-sin(d)*dR/dt-cos(d)*sin(b))*(sin(c)*cos(d)*t+cos(c)*R*cos(b)+sin(c)*sin(d)*R*sin(b)+A-K*e)=0

求解超越方程组.txt

1.16 KB, 下载次数: 12

回复
分享到:

使用道具 举报

发表于 2008-9-29 19:36 | 显示全部楼层

回复 楼主 longdebaby 的帖子

你这个怎么还含有dR/dt,那他就是微分方程了啊?
超越方程的求解可以试一下fsolve等
 楼主| 发表于 2008-10-4 21:50 | 显示全部楼层
问题就是这样啊,我解了很久都没有解出来。所以就到论坛来求救了啊,大侠们救命啊!!!!!
 楼主| 发表于 2008-10-4 21:51 | 显示全部楼层
我试了fsolve,可是老是提示我出错,说超出范围。
 楼主| 发表于 2008-10-6 20:18 | 显示全部楼层
请教高手啊
 楼主| 发表于 2008-10-9 14:46 | 显示全部楼层
自己顶起!
发表于 2008-10-11 06:58 | 显示全部楼层
改用数值解吧
 楼主| 发表于 2008-10-13 14:52 | 显示全部楼层
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-26 05:44 , Processed in 0.059054 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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