24小时热门版块排行榜    

查看: 1288  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 面上合作单位盖章 +3 ssyjh 2026-08-27 3/150 2026-08-27 09:59 by 阿拉贡aragon
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +4 a089 2026-08-26 5/250 2026-08-27 09:24 by Howard28
[基金申请] 基金不中,共勉 +7 eulota 2026-08-26 7/350 2026-08-27 09:02 by hsding2026
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 17/850 2026-08-27 08:20 by 大萍1987
[基金申请] 国自然面上复盘~欢迎讨论 +7 晴天加油 2026-08-26 8/400 2026-08-27 07:53 by longfie172
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 科研孤儿太难了 +21 我4大白菜 2026-08-20 22/1100 2026-08-26 11:29 by nlgza
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见