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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 工科数学 » 求助化学反应产物生成解析解
汕头大学海洋科学接受调剂
51/1返回列表
查看: 1358  |  回复: 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的回帖

peterflyer

木虫之王 (文学泰斗)

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

1.png
...

这个问题刚刚在私信中我已经讲的很清楚了。应该可以求得解析解的。
赞一下(1人) 回复此楼
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的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085408光电信息工程专硕355一志愿长春光机所调剂 +6 王ymaa 2026-04-13 13/650 2026-04-14 11:33 by 王ymaa
[考研] 求调剂 +8 小聂爱学习 2026-04-11 11/550 2026-04-14 09:23 by 1587964
[考研] 305求调剂 +8 玛卡巴卡boom 2026-04-11 8/400 2026-04-14 09:04 by pengliang8036
[考研] 288求调剂，一志愿华南理工大学071005 +24 ioodiiij 2026-04-08 24/1200 2026-04-14 08:42 by nana_zhong
[考研] 291 求调剂 +35 化工2026届毕业� 2026-04-09 35/1750 2026-04-14 08:34 by 科宇老师
[考研] 求调剂 +3 我爱高数高数爱� 2026-04-12 3/150 2026-04-14 01:00 by 王珺璞
[考研] 一志愿中南大学 0855 机械 286 求调剂 +11 不会吃肉 2026-04-12 11/550 2026-04-13 21:59 by bljnqdcc
[考研] 材料考研调剂 +29 云木达达 2026-04-11 31/1550 2026-04-13 13:32 by lyh鲁老师
[考研] 314求调剂 +24 wakeluofu 2026-04-09 25/1250 2026-04-13 08:58 by lhj2009
[考研] 22408 352分求调剂 +5 努力的夏末 2026-04-09 5/250 2026-04-12 19:17 by wj165256
[考研] 305求调剂 +6 77Qi 2026-04-07 6/300 2026-04-12 02:30 by 秋豆菜芽
[考研] 296求调剂 +14 汪！？！ 2026-04-08 15/750 2026-04-11 20:28 by dongdian1
[考研] 275求调剂 +9 1624447980 2026-04-08 10/500 2026-04-11 10:20 by Delta2012
[考研] 0854调剂 +4 长弓傲 2026-04-09 4/200 2026-04-11 09:18 by 猪会飞
[考研] 材料与化工调剂 +12 否极泰来2026 2026-04-10 13/650 2026-04-11 00:28 by wangjihu
[考研] 一志愿京区985，085401电子信息，本科电子信息 +3 阳光开朗的男孩 2026-04-10 3/150 2026-04-10 16:29 by sophia_93
[考研] 301求调剂 +5 149. 2026-04-10 5/250 2026-04-10 15:45 by 柴小白
[考研] 278求调剂 +27 范婷娜 2026-04-07 31/1550 2026-04-09 20:49 by zhouxiaoyu
[考研] 283电子信息求调剂 +4 三石WL 2026-04-08 4/200 2026-04-09 10:21 by wp06
[考研] 求调剂 +13 柒luck 2026-04-07 13/650 2026-04-08 22:46 by 猪会飞
信息提示
关闭
请填处理意见
关闭