24小时热门版块排行榜    

查看: 1973  |  回复: 16
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

taxue1988

银虫 (小有名气)

[求助] 求助循环反应反应速率-线性微分方程组求解

X生成Y的速率是K1,Y生成Z的速率是K2,Z生成X的速率是K3。得到线性微分方程组:
dX/dt=-K1*X+K3*Z
dY/dt=-K1*Y+K1*X
dZ/dt=-K1*Z+K2*Y
t是反应时间
求X,Y,Z表达式

[ Last edited by taxue1988 on 2014-2-28 at 16:14 ]
回复此楼

» 猜你喜欢

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

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

taxue1988

银虫 (小有名气)

引用回帖:
6楼: Originally posted by feixiaolin at 2014-03-01 19:31:48
4#补正
于是,x=-K3*exp(-K1*z),  y=-K2*exp(-K1*z)
代入 y'=-K1*y+K1*x
K1*K2*exp(-K1*z)*z'=K1*K2*exp(-K1*z)-K1*K3*exp(-K1*z),  
K2*z'=K2-K3
z=(K2-K3)/K2*t
最后
x=-K3*exp(-K1*(K2-K3)/K2*t),  
y=- ...

感谢你耐心的回答,金币先送上。但对解析式(z=(K2-K3)/K2*t)并不满足原方程(dZ/dt=-K1*Z+K2*Y),求释疑
8楼2014-03-01 23:38:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 17 个回答

feixiaolin

荣誉版主 (文坛精英)

优秀版主

2楼2014-02-28 23:06:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

taxue1988

银虫 (小有名气)

送红花一朵
引用回帖:
2楼: Originally posted by feixiaolin at 2014-02-28 23:06:03
http://muchong.com/bbs/viewthread.php?tid=6650922&fpage=2&target=self&page=1

能帮忙写一下解析解吗?链接中是代码,我需要用解析解取拟合数据。数学比较差,一碰到复根就晕菜了
3楼2014-03-01 18:07:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与,应助指数 +1
taxue1988: 金币+100, ★★★很有帮助 2014-03-01 23:38:49
jiangyun2: 金币+2, 谢谢回帖交流,欢迎常来~~ 2014-03-04 09:25:01
原式dX/dt=-K1*X+K3*Z
dY/dt=-K1*Y+K1*X
dZ/dt=-K1*Z+K2*Y
做拉普拉斯变换,有
(s+K1)*X+K3*Z=0
(s+K1)*Y+K1*X=0
(s+K1)*Z+K2*Y=0
直接用辣死你变换求解有困难,间接法:
将z视为 x, y的函数,
X=-K3/(s+K1)*Z
Y=-K2/(s+K1)*Y
于是,x=-K3*exp(-K1*z),  y=-K2*exp(-K1*z)
代入 z'=-K1*z + K2*y,有
z'=-K1*z  -K2^2*exp(-K1*z)
z=z1+z2, 裂项处理
z=exp(-K1*t) + ln[-K1*(c1+K2^2*t)]/K1

可参考 http://www.wolframalpha.com/inpu ... 80%99%5Bx%5D%3D%3Da*exp%28-b*y%5Bx%5D%29%5D&dataset=
4楼2014-03-01 18:46:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +28 医学老男孩 2026-08-13 60/3000 2026-08-15 19:54 by 医学老男孩
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +8 majunge000 2026-08-11 10/500 2026-08-15 18:55 by Equinoxhua
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +7 Ldrop2023 2026-08-13 7/350 2026-08-15 15:57 by mzhh2000
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
请填处理意见