24小时热门版块排行榜    

查看: 1333  |  回复: 4

cooking

铜虫 (小有名气)

[求助] 关于微分方程组的求解

请教虫友关于关于微分方程组的求解(见图1):
1.用什么软件运算才可以得到图2的解(图2中的a去掉或认为a=1),也就是说怎样才可以得到图2的解(不包括图2中的a,因为解中是没有a的,a是t=0时,A1=a)。
2.如果用matlab具体怎样操作?谢谢

图1.jpg



图2.jpg



[ Last edited by cooking on 2012-10-13 at 11:12 ]
回复此楼

» 猜你喜欢

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

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

guardian_qq

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
cooking: 金币+10, ★★★很有帮助, 谢谢! 2012-10-13 17:09:25
从图1得不到你要的图2结果(最多只能得到图2的第一式),因为你实际只给了一个边界条件A1(0)=a,你只有再补充3个边界条件(A2(0),A3(0),A4(0))才有可能,否则解出的A2A3A4必含有3个待定常数c1c2c3。
可以用Mathematics的DSolve命令求解。

» 本帖已获得的红花(最新10朵)

2楼2012-10-13 13:35:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooking

铜虫 (小有名气)

送鲜花一朵
引用回帖:
2楼: Originally posted by guardian_qq at 2012-10-13 13:35:21
从图1得不到你要的图2结果(最多只能得到图2的第一式),因为你实际只给了一个边界条件A1(0)=a,你只有再补充3个边界条件(A2(0),A3(0),A4(0))才有可能,否则解出的A2A3A4必含有3个待定常数c1c2c3。
可以用Mathemat ...

另外3个边界条件是A2(0)=0,A3(0)=0,A4(0)=0。请问具体怎样求解?谢谢!
3楼2012-10-13 17:10:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guardian_qq

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
cooking: 金币+40, ★★★★★最佳答案, 解决了!非常感谢! 2012-10-13 23:44:01
引用回帖:
3楼: Originally posted by cooking at 2012-10-13 17:10:45
另外3个边界条件是A2(0)=0,A3(0)=0,A4(0)=0。请问具体怎样求解?谢谢!...

我用的是Mathematics(我觉得符号运算不是Matlab的强项),直接输入:
Simplify[
   DSolve[
   {
      y1'[t] + k1*y1[t] == 0,
      y2'[t] - k1*y1[t] + k3*y2[t] == 0,
      y3'[t] - k3*y2[t] + k5*y3[t] - k6*y4[t] == 0,
      y4'[t] - k5*y3[t] + k6*y4[t] == 0,
      y1[0] == a, y2[0] == 0,
      y3[0] == 0, y4[0] == 0
    },
    {y1[t], y2[t], y3[t], y4[t]}, t
    ]
]
然后按shift+enter就能得到你要的结果了。
(上面为了让你看清才写成那样的格式,其实格式无所谓的)
4楼2012-10-13 18:02:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

醉乡常客

木虫 (知名作家)

混之

都是一级反应,难度不是想象中那么大。既然来数学版块求助,问得明确一些更容易得到有效的帮助。
混混,混混!(求助请注意礼貌!)
5楼2012-10-13 18:32:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cooking 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 上海工程技术大学激光智能制造课题组|2027级博士研究生招生公告 +7 水士口 2026-09-16 8/400 2026-09-16 13:47 by kkklens
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-16 06:27 by LH5NkK5Ud4bh
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 05:35 by LH5NkK5Ud4bh
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 05:15 by LH5NkK5Ud4bh
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-16 04:51 by LH5NkK5Ud4bh
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 04:47 by LH5NkK5Ud4bh
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 5/250 2026-09-16 03:59 by DgNGHc3h5tPl
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-13 7/350 2026-09-16 03:26 by DgNGHc3h5tPl
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-16 03:14 by DgNGHc3h5tPl
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-16 02:50 by DgNGHc3h5tPl
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +7 3n8v2C8RimXI 2026-09-13 7/350 2026-09-16 02:35 by DgNGHc3h5tPl
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 20:26 by zNcEhTcH7ihq
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 19:59 by CgyNCDVNhGVg
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 17:15 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:35 by CgyNCDVNhGVg
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 QUjhNVAcOSff 2026-09-13 9/450 2026-09-15 14:10 by VMDqgnlKimtZ
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:50 by 5BDX0d0WFp7t
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:25 by vZfe6xYu34yj
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见