24小时热门版块排行榜    

查看: 1755  |  回复: 5

约定msn

新虫 (小有名气)

[求助] 求高手帮我解一阶常微分方程组,非常感谢!!

dy1/dt=-0.4y1-0.1y2
dy2/dt=0.2y1-0.3y2-0.2y3
dy3/dt=0.1y2-0.4y3-0.1y4
dy4/dt=0.2y3-0.3y4-0.2y5
dy5/dt=0.1y4-0.4y5-0.1y6
dy6/dt=0.2y5-0.3y6-0.2y7
dy7/dt=0.1y6-0.4y7-0.1y8
dy8/dt=0.2y7-0.3y8-0.2y9
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

约定msn

新虫 (小有名气)

2楼2013-11-30 11:32:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

LT[dy1/dt]=s*Y1(s)-y1(0);
仿此,做拉普拉斯变换,
求解出Yi(s),i=1-8,再逆变换回来即可。
3楼2013-11-30 13:29:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer


【答案】应助回帖

感谢参与,应助指数 +1
楼主给个初始条件吧,只有这样才能获得通解和特解。
  假如t=0, yi=Ci ,  i=1~8, 再记yi的拉氏变换为Ui,则获得如下的八元一次方程组:
    s*U1-C1=-0.4*U1-0.1*U2                                      (1)
      s*U2-C2=0.2*U1-0.3*U2-0.2*U3                             (2)
      s*U3-C3=0.1*U2-0.4*U3-0.1*U4                             (3)
      s*U4-C4=0.2*U3-0.3*U4-0.2*U5                              (4)
      s*U5-C5=0.1*U4-0.4*U5-0.1*U6                             (5)
      s*U6-C6=0.2*U5-0.3*U6-0.2*U7                              (6)
      s*U7-C7=0.1*U6-0.4*U7-0.1*U8                              (7)
      s*U8-C8=0.2*U7-0.3*U8-0.2*U9                              (8)
     求解上面的方程组,即可获得Ui (i=1~8)的表达式。由上可知,Ui 的表达式很复杂,但它们均是s的有理分式,当进行拉是反变换时,它们无不是由常数、exp(ai*t)、exp(-bi*t)、Cos(hi*t)、Sin(ki*t)等初等函数的组合。由于计算太过复杂,就不具体解出具体表达式了。但应该已不存在数学上的困难了。

完毕
4楼2013-11-30 22:20:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

约定msn

新虫 (小有名气)

引用回帖:
4楼: Originally posted by peterflyer at 2013-11-30 22:20:50
楼主给个初始条件吧,只有这样才能获得通解和特解。
  假如t=0, yi=Ci ,  i=1~8, 再记yi的拉氏变换为Ui,则获得如下的八元一次方程组:
    s*U1-C1=-0.4*U1-0.1*U2                                      (1)
   ...

恩 好的 就是过于复杂
5楼2013-12-01 09:55:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer


引用回帖:
5楼: Originally posted by 约定msn at 2013-12-01 09:55:27
恩 好的 就是过于复杂...

复杂是肯定的,谁让你是八元常微分方程组呢。幸亏全部系数都是常系数,还解得出,要是变系数的话,一般情况下就只有数值法一条路了。
6楼2013-12-01 11:30:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 约定msn 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 极限拉扯,如何抉择 +10 Isak0320 2026-10-12 14/700 2026-10-14 11:39 by thegarfield
[基金申请] 申请科研项目纯靠实力能成功吗? 200+3 David1982 2026-10-13 6/300 2026-10-14 11:14 by 北京莱茵编辑
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 12/600 2026-10-13 20:40 by Merryya
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见