YI)=Y④+RK()*H/6 50 CONTINUE RETURN END 4.运行结果。(20分)（每列5分） N=4T0= 0T1= 20.0H0=10.0 EPS=.100E-05 Y0.=1,2,3.4.N t/min H/min [A]/molL-1 [B]molL-1 [ClmolL-1 [DYmolL-I 00 10.00 100000 000000 000000 3 49574E-01 293546E-02 210712E-02 412842E-0 912842E 316478E-( 74E.0 881574E.01 2927 4000 35380 26432E02 5000 119775E-01 1 328818E-01 828818E-01 238091F-02 1473 3520- 215166E-02 172274E01 80.00 1.25 194817E0 766942E.01 21529 90.00 A0815E01 749815E.01 233968E.01 234109E-01 7341 148510E-02 251040E-01 11000 219667E.01 719667E.01 266696E-01 125 2063s7E01 125538E.02 281090E-01 1300n 104061G01 504n61E.n1 115833E02 30435sD01 14000 1856g1B01 66Q101 107107E02 150 00 175126E01 672126E0 00233703 317950E-01 160 00 3 160233B01 663232E01 921084E.03 1700n 153198E-01 65310gB01 856407E.03 33Q23gF01 1g000 125 144605E01 644A605E.01 707530.03 347320C.01 190.00 125 136759E01 636750D.01 .743818E-03 355803E-01 200.00 1.25 129343E-01 629343E-01 694677E-03 363711E-01 C 1.计算原理（化学原理和计算方法）(20分) 液相反应：A+B→R 该反应的动力学方程为 =-der=kca"cn (1) dt (1)式中，广为化学反应速率，c为生成物浓度，1为反应时间，m和n为反应级数， 6
6 40 Y(I)=Y(I)+RK(I)*H/6 50 CONTINUE RETURN END 4.运行结果。（20 分）（每列 5 分） N= 4 T0= .0 T1= 200.0 H0=10.0 EPS= .100E-05 Y(I),I=1,2,3,4,...,N --------------------------------------------------------------------------------------------------------------------- t/min H/min [A]/molL-1 [B]/molL-1 [C]/molL-1 [D]/molL-1 --------------------------------------------------------------------------------------------------------------------- .00 10.00 .500000E-01 .100000 .000000 .000000 10.00 .31 .449574E-01 .949574E-01 .293546E-02 .210712E-02 20.00 .31 .412842E-01 .912842E-01 .316478E-02 .555106E-02 30.00 .63 .381574E-01 .881574E-01 .292780E-02 .891484E-02 40.00 1.25 .353802E-01 .853802E-01 .264232E-02 .119775E-01 50.00 1.25 .328818E-01 .828818E-01 .238091E-02 .147373E-01 60.00 1.25 .306210E-01 .806210E-01 .215166E-02 .172274E-01 70.00 1.25 .285669E-01 .785669E-01 .195149E-02 .194817E-01 80.00 1.25 .266942E-01 .766942E-01 .177611E-02 .215297E-01 90.00 1.25 .249815E-01 .749815E-01 .162170E-02 .233968E-01 100.00 1.25 .234109E-01 .734109E-01 .148510E-02 .251040E-01 110.00 1.25 .219667E-01 .719667E-01 .136371E-02 .266696E-01 120.00 1.25 .206357E-01 .706357E-01 .125538E-02 .281090E-01 130.00 1.25 .194061E-01 .694061E-01 .115833E-02 .294355E-01 140.00 1.25 .182681E-01 .682681E-01 .107107E-02 .306609E-01 150.00 1.25 .172126E-01 .672126E-01 .992337E-03 .317950E-01 160.00 1.25 .162322E-01 .662322E-01 .921084E-03 .328468E-01 170.00 1.25 .153198E-01 .653198E-01 .856407E-03 .338238E-01 180.00 1.25 .144695E-01 .644695E-01 .797539E-03 .347330E-01 190.00 1.25 .136759E-01 .636759E-01 .743818E-03 .355803E-01 200.00 1.25 .129343E-01 .629343E-01 .694677E-03 .363711E-01 ------------------------------------------------------------------------------------------------------------------------ C 1.计算原理（化学原理和计算方法）（20 分） 液相反应： A + B → R 该反应的动力学方程为 R R AB d d c m n r kc c t =− = （1） （1）式中， Rr 为化学反应速率，cR为生成物浓度，t 为反应时间，m 和 n 为反应级数
k为反应速率常数，dc山为生成物浓度随时间的变化率。 (5分) 因c0=C,故存在CA=CB,于是(1)式可改写为 =-der =ke (2) d 只要求得待定参数k、m和n,该动力学方程即可确定。欲确定k、m十，首先要建立 ~Ca相应的数据组。 (5分) (2)式两边取对数：n=ln(-d)=nk+(m+m)nca d 令y=ln,x=lnca,a=lnk,b=m+n,得y=a+bx (5分) 对于cR(,)对应的反应物的浓度cA()可由式(3)求出 CA()=CA0-CR(1) (3) 将计算所得不同时间1的反应物浓度ca数据通过插值和差分法求出反应速率。计算 nR和lnca后，用线性回归子程序计算反应级数n和反应速率常数k。 (5分) 2.程序框图(30分) 开始 输入：数据点数N, 反应温度TO, 时间间隔H 时间t与cR的实验数据T④，CR① (3分) CA(D)=CAO-CR(I)(I=1,N)( X=T①，YI=CR0I=1,N)(5分) 调用中心差分子程序计算反应速率R(W)(dcd)(5分) X①=nCA,Y①=In-R(=1,N)(5分)
7 k 为反应速率常数，dcR/dt 为生成物浓度随时间的变化率。 （5 分） 因 A0 B0 c c = ，故存在 A B c c = ，于是（1）式可改写为 R R A d d c m n r kc t + =− = （2） 只要求得待定参数 k、m 和 n，该动力学方程即可确定。欲确定 k、m＋n，首先要建立 Rr ~ Ac 相应的数据组。 （5 分） （2）式两边取对数： R R A d ln ln( ) ln ( )ln d c r k mn c t =− = ++ 令 R A y = = = = + =+ ln , ln , ln , , r x c a k b m n y a bx 得 （5 分） 对于 cR（ti）对应的反应物的浓度 cA（ti）可由式（3）求出 A A0 R () () i i ct c ct = − （3） 将计算所得不同时间 t 的反应物浓度 cA数据通过插值和差分法求出反应速率 Rr 。计算 ln Rr 和 lncA后，用线性回归子程序计算反应级数 n 和反应速率常数 k。 （5 分） 2.程序框图（30 分） 开始 输入：数据点数 N, 反应温度 TO, 时间间隔 H 时间 t 与 cR的实验数据 T(I), CR(I) （I（3 分） CA(I)=CA0-CR(I) （I=1,N） （4 分） 调用中心差分子程序计算反应速率 R(W)（dcR/dt） （5 分） X(I)=ln[CA(I)],Y(I)=ln[-R(I)] (I=1,N) （5 分） X(I)=T(I),Y(I)=CR(I) (I=1,N) （5 分）