24小时热门版块排行榜    

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

树懒

银虫 (小有名气)

[求助] 求助化学反应产物生成解析解 已有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

» 猜你喜欢

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

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

peterflyer

木虫之王 (文学泰斗)

peterflyer


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

1.png
...

这个问题刚刚在私信中我已经讲的很清楚了。应该可以求得解析解的。
7楼2016-04-28 09:11:07
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 328求调剂,英语六级551,有科研经历 +4 生物工程调剂 2026-03-17 8/400 2026-03-21 02:12 by JourneyLucky
[考研] 一志愿华南师大 070300(化学)304分求调剂 +3 0703武芊慧雪304 2026-03-18 3/150 2026-03-21 00:48 by JourneyLucky
[考研] 299求调剂 +4 shxchem 2026-03-20 5/250 2026-03-21 00:24 by 星空星月
[考研] 274求调剂 +10 S.H1 2026-03-18 10/500 2026-03-20 23:51 by JourneyLucky
[考研] 材料专硕英一数二306 +7 z1z2z3879 2026-03-18 7/350 2026-03-20 23:48 by JourneyLucky
[考研] 南昌大学材料专硕311分求调剂 +5 77chaselx 2026-03-20 5/250 2026-03-20 23:42 by lovewei0727
[考研] 295求调剂 +4 一志愿京区211 2026-03-18 6/300 2026-03-20 23:41 by JourneyLucky
[考研] 一志愿中海洋材料工程专硕330分求调剂 +8 小材化本科 2026-03-18 8/400 2026-03-20 23:16 by JourneyLucky
[考研] 考研调剂求学校推荐 +3 伯乐29 2026-03-18 5/250 2026-03-20 22:59 by JourneyLucky
[考研] 323求调剂 +3 洼小桶 2026-03-18 3/150 2026-03-20 22:54 by JourneyLucky
[考研] 265求调剂 +12 梁梁校校 2026-03-19 13/650 2026-03-20 21:01 by 无际的草原
[考研] 295复试调剂 +8 简木ChuFront 2026-03-19 8/400 2026-03-20 20:44 by zhukairuo
[考研] 一志愿华中农业071010,总分320求调剂 +3 困困困困坤坤 2026-03-20 3/150 2026-03-20 20:38 by 学员8dgXkO
[考研] 260求调剂 +3 朱芷琳 2026-03-20 3/150 2026-03-20 20:35 by 学员8dgXkO
[考研] 289求调剂 +6 怀瑾握瑜l 2026-03-20 6/300 2026-03-20 20:30 by 学员8dgXkO
[考研] 广西大学家禽遗传育种课题组2026年硕士招生(接收计算机专业调剂) +3 123阿标 2026-03-17 3/150 2026-03-20 15:58 by 飞行琦
[考研] 梁成伟老师课题组欢迎你的加入 +9 一鸭鸭哟 2026-03-14 11/550 2026-03-19 17:22 by !本暗一次!
[考研] 266求调剂 +5 阳阳哇塞 2026-03-14 10/500 2026-03-19 15:08 by 阳阳哇塞
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
信息提示
请填处理意见