北京化工大学2003-2004学年第一学期 《计算化学》期末试卷答案 9 1.计算原理（化学原理和计算方法）(20分) 反应 A+B C+D B+C->D+E 其动力学方程式为一个常微分方程组： dIA]=-[A][B]+(CID] At d[B]-[A][B]+k.(CID]-&,[BJC] dIC]=IA]IB]-&(CID]-&,[B]C] d d[D]-k[A]B]-&[CID]+[BC] dt d回=kBC] 利用Run e-Kutta法解此方程组，可得不同时刻[A],B,[C],D,E值 2.程序框图 30分) 开始 输入：[A,B,[C,D,E初值 时间间隔H,计算精度EPS 计算：积分步长H0 调用变步长Runge--Kuta积分法子程序解动力学方程组计算时刻t 的[A,B,[C],D,E值 (其中调用计算动力学方程组的子程序计算DY值) 列表输出：0-200min内每隔10min所对应的[A,B,[C],D],E]值 北京化工大学 2003-2004 学年第一学期 《计算化学》期末试卷答案 A 1．计算原理（化学原理和计算方法）（20 分） 反应 k1 k2 A+B C+D B C D+E 3k + ⎯⎯→ 其动力学方程式为一个常微分方程组： 1 2 1 23 1 23 1 23 3 d[A] [A][B] [C][D] d d[B] [A][B] [C][D] [B][C] d d[C] [A][B] [C][D] [B][C] d d[D] [A][B] [C][D] [B][C] d d[E] [B][C] d k k t kkk t kkk t kkk t k t =− + =− + − = −− = −+ = 利用 Runge-Kutta 法解此方程组，可得不同时刻[A]，[B]，[C]，[D]，[E]值。 2．程序框图（30 分） 开始 输入：[A]，[B]，[C]，[D]，[E]初值 时间间隔 H， 计算精度 EPS 计算：积分步长 H0 调用变步长 Runge-Kutta 积分法子程序解动力学方程组计算时刻 t 的[A]，[B]，[C]，[D] ，[E]值 （其中调用计算动力学方程组的子程序计算 DY 值） 列表输出：0-200min 内每隔 10min 所对应的[A]，[B]，[C]，[D] ，[E]值