2.程序框图(20分) D0100I=1,N X(ILOG(Y(D)) 开始 YIDLOG(ABS(V(DD) 100 CONTINUE 输入：数据点数N, 反应温度TO, 时间间隔H CALL PK(N,X.YA,B,R) =B 时间t与XA的实验数据T),XA(① (2分) SS=EXP(A) WRITE(6,11)A,B,R,S,SS 1 FORMAT/5XA=,D18.10M5X'B=,D18.10,5X,R=D18.10M& CA(D)=CAO (1-xA(D))(I=1.N)( 5X,=,D18.10,5Xk=,D18.10W END X(I=T(1).Y(1)=CA(I)(I=1.N)( SUBROUTINE CF(N,H,X,Y,Z) (4分) 调用中心差商子程序计算反应速率R(wW)(dcd山r)(4分) DIMENSION X(N),Y(N),Z(N) DOUBLE PRECISION X,Y.Z,H,T,YA,YB T=X(1)+H X(I)=In[CA(I)].Y(1)=In[-R(I)](I=I,N)( CALL LGRGZ(X,YN,T,YA) T=X(1H7H 调用一元线性回归子程序计算A,B CALLLGRG2(X,YN,T,YB) Z1=4*YA-3Y(1-YB2*H (反应级数N=B,速率常数KS=EXPA) (4分)N) D0101=2N-1 T-X(D)-H 输出：A,B,RN,KS (2分) CALL LGRG2(X,Y,N,T,YA) TX(IH+H 结束 CALLLGRG2(X,YN,T,YB) Z(I(YB-YAM2+H) 10 CONTINUE 3.源程序(20分) T=X(N)-H PROGRAM MAIN (10分) CALL LGRG2(X,Y,N,T,YA) DIMENSION X(5).Y(5)V(5) T=XN-2◆H DOUBLE PRECISION X,Y,V,A,B,R,S,SS,CA,H CALLLGRG2(X,Y.N,T,YB) OPEN(6,FILE-2010TESTE.TXT,STATUS-UNKNOWN) Z(N3Y(N)-4*YA+YBM(2*H) N=5 REIΠRN 1.0E-6 END CA=2.00 DATA X00306090120/ SUBROUTINE PK(N,X,Y,A,B,R) (3分) DATA Y0,0.46,0.63,0.72,0.77/ DOUBLE PRECISION X,Y,A,B,R DO 50 I=1,N DIMENSION X(N),Y(N) Y(IFCA*(1-Y(D)) sX=0.0 0 CONTINUE SY三0.0 CALL CF(N,H,X,Y,V) SXX=0.08 2.程序框图（20 分） 3.源程序（20 分） PROGRAM MAIN （10 分） DIMENSION X(5),Y(5),V(5) DOUBLE PRECISION X,Y,V,A,B,R,S,SS,CA,H OPEN(6,FILE='2010TESTE.TXT',STATUS='UNKNOWN') N=5 H=1.0E-6 CA=2.00 DATA X/0.0,3.0,6.0,9.0,12.0/ DATA Y/0.,0.46,0.63,0.72,0.77/ DO 50 I=1,N Y(I)=CA*(1-Y(I)) 50 CONTINUE CALL CF(N,H,X,Y,V) DO 100 I=1,N X(I)=LOG(Y(I)) Y(I)=LOG(ABS(V(I))) 100 CONTINUE CALL PK(N,X,Y,A,B,R) S=B SS=EXP(A) WRITE(6,11)A,B,R,S,SS 11 FORMAT(/5X,'A=',D18.10//5X,'B=',D18.10, 5X,'R=',D18.10// & 5X,'n=',D18.10,5X,'k=',D18.10/) END SUBROUTINE CF(N,H,X,Y,Z) （4 分） DIMENSION X(N),Y(N),Z(N) DOUBLE PRECISION X,Y,Z,H,T,YA,YB T=X(1)+H CALL LGRG2(X,Y,N,T,YA) T=X(1)+2*H CALL LGRG2(X,Y,N,T,YB) Z(1)=(4*YA-3*Y(1)-YB)/(2*H) DO 10 I=2,N-1 T=X(I)-H CALL LGRG2(X,Y,N,T,YA) T=X(I)+H CALL LGRG2(X,Y,N,T,YB) Z(I)=(YB-YA)/(2*H) 10 CONTINUE T=X(N)-H CALL LGRG2(X,Y,N,T,YA) T=X(N)-2*H CALL LGRG2(X,Y,N,T,YB) Z(N)=(3*Y(N)-4*YA+YB)/(2*H) RETURN END SUBROUTINE PK(N,X,Y,A,B,R) （3 分） DOUBLE PRECISION X,Y,A,B,R DIMENSION X(N),Y(N) SX=0.0 SY=0.0 SXX=0.0 开始 输入：数据点数 N, 反应温度 TO, 时间间隔 H 时间 t 与 xA 的实验数据 T(I), XA(I) （I（2 分） CA(I)=CA0（1－xA(I)）（I=1,N） （2 分） 调用中心差商子程序计算反应速率 R(W)（dcA/dt） （4 分） X(I)=ln[CA(I)], Y(I)=ln[-R(I)] (I=1,N) （4 分） 调用一元线性回归子程序计算 A，B （反应级数 N＝B, 速率常数 KS=EXP(A)） （4 分）N） 输出： A，B，R, N, KS （2 分） 结束 X(I)=T(I),Y(I)=CA(I) (I=1,N) （2 分）