24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1578  |  回复: 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的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

引用回帖:
12楼: Originally posted by taxue1988 at 2014-03-03 17:33:05
您误会了,我是说拉普拉斯逆变换这一步,逆变换表可以查到K3/(s+K1)的变换形式,但K3/(s+K1)*Z就不会了...

X=-K3/(s+K1)*Z
X/Z=-K3/(s+K1)   InvLT[X/Z]= -K3*exp(-K1*t)
x= InvLT[X/Z]卷积z = -K3*exp(-K1*t)卷积z
我这里投机,取成了x= -K3*exp(-K1*z)
其实,写成   x= -K3*exp(-K1*t)卷积(at^2+bt+c)     都比x= -K3*exp(-K1*z) 更合理。
13楼2014-03-03 20:01:45
已阅   回复此楼   关注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的回帖
信息提示
请填处理意见