应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 工科数学 » 求助化学反应产物生成解析解
51/1返回列表
查看: 1270  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

树懒

银虫 (小有名气)

[求助] 求助化学反应产物生成解析解 已有1人参与

小弟想求教一个化学反应产物生成的浓度解析解
反应如下
                A   +    B     --------     C     +    B   ---------     D
                                      k1                               k2
t=0          a          b                     0                                0
t=t0        a-x       b-2y                x                                y
化合物A可以和B反应生成产物C反应速率常数为k1，C也可以和B反应生成D反应常数为k2。两个反应都不考虑反向的反应。
在反应刚开始时A的浓度为a，B的浓度为b。C和D都未生成
当反应进行到t0时产物C和D的浓度分别为x和y。
A和B的浓度分别a-x和b-2y
是否能够求得x和y随t变化的解析解？

[ Last edited by 树懒 on 2013-12-20 at 23:01 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : 未命名.bmp
    • 2013-12-20 22:55:54, 201.85 K

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-12-20 22:56:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Mr__Right

专家顾问 (著名写手)
引用回帖:
6楼: Originally posted by 277790320 at 2016-04-27 14:57:04
你好，请教一个问题：
我拟合一个动力学方程，但是这个方程得不到其解析解，我看一些资料上用差分的方法进行操作。由于本人数学功底较差，还请大神帮忙下，式中k K是未知常数，Ks已经求得，Ca=50-Cn。

1.png
...

最好另开一个新帖@ 想要找的答题者
赞一下 回复此楼
文章乃身外之物,要多考虑编辑、审稿人和读者的感受。
8楼2016-04-28 09:35:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

peterflyer

木虫之王 (文学泰斗)

peterflyer

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
feixiaolin: 金币+1 2013-12-21 09:48:49
树懒: 金币+5, ★★★很有帮助, 大神非常抱歉，题目中的条件我写得有一些错误的地方，能不能看一下帖子，帮我重新解一下 2013-12-21 11:33:12
依据题意有：
     dx/dt=k1*(a-x)                   （1）
        dy/dt=k2*x                        （2）
由（1），x=a-C1*e^(-k1*t)
由t=0时x=0，C1=a ,故 x=a*[1-e^(-k1*t)]
  将x代入（2），y=k2*a*t+k2*a/k1*e^(-k1*t)+C2
由t=0时y=0，C2=-k2/k1*a ,
  故 y=k2*a*t+k2*a/k1*[e^(-k1*t)-1]
赞一下 回复此楼
2楼2013-12-21 01:24:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

树懒

银虫 (小有名气)
楼上的大神非常抱歉，题目中的条件我写得有一些错误的地方

在t=t0的时候化合物A的浓度应该是等于

a-x-y而不是a-x

所以能不能麻烦您重新帮我解一次
赞一下 回复此楼
3楼2013-12-21 11:32:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
树懒: 金币+15, ★★★很有帮助, 大神，到这个程度我已经看不明白了，我让我师兄看过，他说没问题他再好好看一下。谢谢你了 2013-12-21 17:55:47
引用回帖:
3楼: Originally posted by 树懒 at 2013-12-21 11:32:24
楼上的大神非常抱歉，题目中的条件我写得有一些错误的地方

在t=t0的时候化合物A的浓度应该是等于

a-x-y而不是a-x

所以能不能麻烦您重新帮我解一次

楼主，你这一改，题目一下子复杂多了。现将解题过程写在下边，请你检查是否有误，有什么问题或不明白的，及时沟通和讨论解决。
解题过程：
dx/dt=k1*(a-x-y)                                                    （1）
dy/dt=k2*x                                                          （2）
x(0)=a,  y(0)=b ，t>0                                               (3)
使用拉普拉斯变换法解常微分方程组。设x、y对t的拉氏变换分别为U、V。
由（1）和（3）：s*U-a=k1/s-k1*U-k1*V
           整理后：（s+k1）*U+k1*V=k1/s                      (4)
由（2）和（3）：s*V-b=k2*U
            整理后：（s+k1）*U+k1*V=k1/s                             (5)
由（4）和（5）联立解得：
   U=k1*(1-b)/{(s+k1/2)^2+[4*k1*k2-k1^2]/4}                   (6)
    V=k1*k2*(1-b)/{s*[(s+k1/2)^2+(4*k1*k2-k1^2)/4]}+b/s   (7)

下面求U、V的拉氏反变换，求出x和y。
（1）若4*k1*k2-k1^2]/4<0，令 -m^2=4*k1*k2-k1^2]/4,  m>0
则由代数中的真有理分式的分解定理，（6）、（7）可改写为：
    U=k1*(1-b)/(2*m)*1/[s-(m-k1/2)]-k1*(1-b)/(2*m)*1/[s-(-k1/2-m)]
    V=α/s+β/[s-(m-k1/2)]+γ/[s-(-m-k1/2)]+b/s
其中α、β、γ可由比较系数法确定，经计算可得：
   α= k1*k2*(1-b)/{2*m^2*[(k1/2)^2-m^2]}
    β= -(k1+2*m)*α/(4*m)
    γ= -(α+β)
故：x=L^(-1)[U]= k1*(1-b)/(2*m)*e^{[(m-k1/2)]*t}- k1*(1-b)/(2*m)* e^{[(-m-k1/2)]*t}
Y= L^(-1)[V]=b+α + β*e^[-( k1/2-m)*t]+γ*e^[-(m+k1/2)*t]

（2）若4*k1*k2-k1^2]/4=0，
U= k1*(1-b)/(s+k1/2)^2= k1*(1-b)/[s-(-k1/2)]^2
V= k1*k2*(1-b)/{s*[(s+k1/2)^2]} +b/s
=b/s+ k1*k2*(1-b)*{4/k1^2*1/s-2/k1*1/[s+k1/2]^2-4/k1^2*1/(s+k1/2)}
x=L^(-1)[U]= k1*(1-b) *t*e^{[(-k1/2)*t}
Y= L^(-1)[V]=b+4*k2*(1-b)/k1-2*k2*(1-b)*t*e^[-(k1/2)*t]
                       -4*k2*(1-b)/k1* e^[-(k1/2)*t]

（3）若4*k1*k2-k1^2]/4>0  ,令m^2=4*k1*k2-k1^2]/4, m>0
U=k1*(1-b)/{(s+k1/2)^2+m^2}
  = k1*(1-b)/m*m/{[s-(-k1/2)]^2+m^2}                        
V=k1*k2*(1-b)/{s*[(s+k1/2)^2+m^2]} +b/s
=b/s+ k1*k2*(1-b)*{1/m^2*1/s+{[-1/m^2*(s+k1/2)]
                             +k1/(2*m^3) *m}/{[s-(-k1/2)]^2+m^2}}
=[b+ k1*k2*(1-b)/m^2]/s
     - k1*k2*(1-b)/m^2*(s+k1/2)/ {[s-(-k1/2)] ^2+m^2}
     + k1^2*k2*(1-b)/(2*m^3)* m/{[s-(-k1/2)]^2+m^2}
故有：
  x=L^(-1)[U]= k1*(1-b)/m*e^[(-k1/2)*t]*Sin(m*t)
  Y= L^(-1)[V]= b+ k1*k2*(1-b)/m^2
                       - k1*k2*(1-b)/m^2* e^[(-k1/2)*t]*Cos(m*t)
                      + k1^2*k2*(1-b)/(2*m^3)* e^[(-k1/2)*t]*Sin(m*t)

注：L^(-1)[U]表示U的拉氏反变换，其他依次类推。


解题完毕。
赞一下 回复此楼
4楼2013-12-21 16:27:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭