北京化工大学2009一一2010学年第一学期 将表中所列实验数据变形后，代入多元线性回归子程序（其中调用高斯消去法子程序解正 规方程组)，通过计算确定方程中的各个系数、加、值和复相关系数.(5分) 通过计算的b、b、b:值得到安托尼方程中的参数A,B和C,A=b1,C=-b,=brAC. 《计算化学》期末考试试卷标准答案 (5 第一都分笔试都分（共30分）】 2程序框图(20分) 一、判断题，正确的在括号中划“√”，错误的划“x”(每题1.5分，共15分) 1.×2.×3.4.5.￥6.￥7.89.×10.￥ 开始 二、填空题（每空0.5分，共15分） 输入：数据点数N, 1.(1)函数在区间1ab1内连续，(2)函数在a,b两点之值Ra)与b)洞号. 2.利用矩阵的初等变换将系数矩阵A约化为一个上三角形矩阵，然后再进行求解.(1)选 苯酚的饱和蒸汽压P(I)和温度T(I)(4分) 主元，(2)消元，(3)回代。 3.函数关系，相关关系。插值法，回归分析。 ngp=y,T=x1, 1gp=x2(4分) 4.差商：中心差商. 5.x)函数形式已知，但其积分不能表示成初等函数的闭合形式：x)函数形式未知，但其 离散数据表己给出。 调用多元线性回归子程序计算回归方程的系数b、b、b:值 6.(1)用插值程序求任意点的函数值，2.用Simpson求积程序计算[a,b]区间中离散点下 及复相关系数R(其中调用高斯消去法子程序解正规方程组) (4分) 的面积。 7.实对称矩阵的本征值和本征向量：量子化学中的分子轨道相关的性质。 8。初值。边值。扩散、反应、传质、传热。 A=b1,C=b2,B=-bo+A*C (4 9.ChemOffice:Origin:Gaussian:Materials Studio. 10.线性方程组的系数矩阵为主对角线占优势。 输出：参数A,B和C(4分) 11.三角形：四面体：不用求导，甚至在没有目标函数表达式的情况下也可使用：色谱过 程的顺序优化。 结束 第二部分上机题（共70分） A 3.源程序(20分) 1.计算原理（化学原理和计算方法）(15分) program main (10分) 苯酚的饱和蕊汽压(p)和温度(T)关系遵从安托尼方程 DOUBLE PRECISION XV.S,A,X,BO,QT,U,QE,SE,R,F,X1,AA,BB,CC dimension XV(7)S(7),XI(6,7) B COMMON /AX/A(6,7),X(6,7) Ig(p'/kPa)=A- (T/K)+C OPEN(5,FILE-2010TAI.TXT,STATUS-OLD) OPEN(6,FILE-2010TA2.TXT,STATUS-UNKNOWN) READ(5,*)M,N 将两边进行变形T1gp=(AC-B)+AT-C1gp(省略单位) M1=M+1 READ(5,)((X(I,J),J=1,MI),I=1,N) 令ng,T=x,lg即=，b0=4AC-B,b=A,b=-C DO151=1N X1(LMI=X(L,I(LOGIO(X(1,2))) 得 y=bo+bx+b2xz (5分) X1L,1=XL,1) 15 X1(L,2=L0G10X1,2)
1 北京化工大学 2009——2010 学年第一学期 《计算化学》期末考试试卷标准答案 第一部分 笔试部分（共 30 分） 一、判断题，正确的在括号中划“3”，错误的划“2”（每题 1.5 分，共 15 分） 1. 2 2. 2 3. 3 4. 3 5. 2 6. 2 7. 3 8. 3 9. 2 10. 2 二、填空题（每空 0.5 分，共 15 分） 1．（1）函数在区间[a,b]内连续，（2）函数在 a、b 两点之值 f(a)与 f(b)同号。 2．利用矩阵的初等变换将系数矩阵 A 约化为一个上三角形矩阵，然后再进行求解。（1）选 主元，（2）消元，（3）回代。 3．函数关系，相关关系。插值法，回归分析。 4．差商；中心差商。 5．f(x)函数形式已知，但其积分不能表示成初等函数的闭合形式； f(x)函数形式未知，但其 离散数据表已给出。 6．（1）用插值程序求任意点的函数值，2．用 Simpson 求积程序计算[a,b]区间中离散点下 的面积。 7．实对称矩阵的本征值和本征向量；量子化学中的分子轨道相关的性质。 8．初值，边值。扩散、反应、传质、传热。 9．ChemOffice；Origin；Gaussian；Materials Studio。 10．线性方程组的系数矩阵为主对角线占优势。 11．三角形；四面体；不用求导，甚至在没有目标函数表达式的情况下也可使用；色谱过 程的顺序优化。 第二部分 上机题（共 70 分） A 1.计算原理（化学原理和计算方法）（15 分） 苯酚的饱和蒸汽压（p*）和温度（T）关系遵从安托尼方程 T C B p A + = − ( / K) lg( / kPa) * 将两边进行变形 * * T lg p = (AC − B) + AT − C lg p （省略单位） 令 Tlgp=y, T=x1, lgp=x2, b0=AC-B, b1=A, b2= -C 得 0 1 1 2 2 y = b + b x + b x （5 分） 将表中所列实验数据变形后，代入多元线性回归子程序（其中调用高斯消去法子程序解正 规方程组），通过计算确定方程中的各个系数 b0、b1、b2值和复相关系数。（5 分） 通过计算的 b0、b1、b2值得到安托尼方程中的参数 A，B 和 C，A＝b1，C＝-b2，B= b0-A*C。 （5 分） 2.程序框图（20 分） 3.源程序（20 分） program main （10 分） DOUBLE PRECISION XV,S,A,X,B0,QT,U,QE,SE,R,F,X1,AA,BB,CC dimension XV(7),S(7),X1(6,7) COMMON /AX/A(6,7),X(6,7) OPEN(5,FILE='2010TA1.TXT',STATUS='OLD') OPEN(6,FILE='2010TA2.TXT',STATUS='UNKNOWN') READ(5,*) M,N M1=M+1 READ(5,*) ((X(I,J),J=1,M1),I=1,N) DO 15 I=1,N X1(I,M1)=X(I,1)*(LOG10(X(I,2))) X1(I,1)=X(I,1) 15 X1(I,2)=LOG10(X(I,2)) 开始 输入：数据点数 N, 苯酚的饱和蒸汽压 P（I）和温度 T（I）（4 分） 调用多元线性回归子程序计算回归方程的系数 b0、b1、b2 值 及复相关系数 R（其中调用高斯消去法子程序解正规方程组） （4 分） 输出：参数 A，B 和 C （4 分） 结束 Tlgp=y, T=x1, lgp=x2（4 分） A＝b1，C＝b2，B= - b0+A*C（4 分）
D0181=1,N DO 10 K=I,N D018J=1,M1 10 XV(I)=XV(I+X(K,I) X(IJEXI() 20 XVOEXVUVDBLEONY CALLMLR(M,N,XV,S,BO,QT,U,QE,SE,R,F) D040=1,M AA=A1,MI)） S1=0 CC=-1*A2,M1) D030K=1,N BB=AA*CC-B0 30 SI=s+XK,IXV》*2 WRITE(6,20)M,N S(I)=DSQRT(S(1)) FORMAT(/4X.'INPUT DATA:,'Xnl'1X.'Xn2.....1X.'Xnm,Yn(m=1.'& D050=1,M 2m-1:2:T/7X.XU/T.IX.X2/gp.IX.Tgp-(AC+B ALIO(I)=0 +A*T-C1gp/1X,661H-V6X,XI',8X,'X2',8X,& D050K=1,N Y,7X,Ycale',7X,'DIFF/1X,66(IH-)) 50 ALIO(D=ALIO(D)+(X(K,I}XV(D))(X(K,MI)-XV(MI)) D040=1,N D060=1,MI YCALC=BO DO60 K=LN D030J=1,M 60 X(K,I(X(K,I)-XV(I)VS(I) YCALC=YCALC+A(J,MIYX(LJ) D0701=1,M DIFF=YCALC-X(I,M1) A(LI)=1.0 WRITEI6 50(X(DEL MILYCALC DIFF D070Jj=+1,M1 FORMAT(3X,I0G10.4) A(L,J)-0 WRITE(6,60)(XVI),1=1,M1) D070K=1,N 60 FORMAT(/IX,'XV,10G10.4) 70 A(IJ=A(IJ)+X(K,I)*X(K.J) WRTE(6,70)(SD,=1,M1) D080I=L,M-1 之 FORMAT(/IX,'SV,10G10.4) D080J=+1,M WRITE(6.S0) 80 AUIFA(LI 80 FORMAT(/5X,Regression equation is:'/1x,60(1H-)) WRITE(6,11)M WRITE(6,90)BO,(A(I,M1),I=1,M) 11 FORMAT(4X,'Normal cquation A(m,m+1).2x,'m=,i2/1x,60(IH-)) 90 FORMAT(IX,TIgp-,G10.4,+,G10.4,"T,+,& DO90ELM G10.4,*gp) 90 WRITE(6,33)(A(1J)JJ=I,MI) WRITE(6,100)QT,U,QE,SE,U/QT,R,F,AA,BB,CC CALL GS(A,M,MI,1.D-6) 100 FORMAT(5X'Q=,G18.6W5X,U=,G18.6/5X,Qe-,G18.6/5X& WRITEI6 22) Sc=,G18.65X,N=,G18.6,5X,R=,G18.6/5X,T=,& WRITE(6,33)(A(I,MI).I=1,M) G18.6,/5X,A=,G18.6/5X,'B=,G18.6,5X,C=,G18.6) B0=0.D0 END J=0.D0 D01001=1,M SUBROUTINE MLR(M,N,XV,S,BO,QT,U,QE,SE,R,F) (5分) A(LMIA(LMIYS(MIVS(D) DOUBLE PRECISION XV,S,A,X,BO,QT,U,QE,SE,R,F BO=B0+A(LMI)*XV(I) DIMENSION XV(M+1),S(M+1),ALIO(50) 100 【U=U+A(L,MI)ALI0I) COMMON /AX/A(6,7),X(6,7) BO=XV(MI)-BO MI=M+I QT=S(MI)S(MI) D0201=1,MI QE-QT-U XV(I)=0.DO SE-DSORT(QE/(N-M-1))
2 DO 18 I=1,N DO 18 J=1,M1 18 X(I,J)=X1(I,J) CALL MLR(M,N,XV,S,B0,QT,U,QE,SE,R,F) AA=A(1,M1) CC=-1*A(2,M1) BB=AA*CC -B0 WRITE(6,20) M,N 20 FORMAT(/4X,'INPUT DATA:','Xn1,'1X,'Xn2,....',1X,'Xnm,Yn(m=1,' & ,i2,',n=1,',i2,')'/7X,'X1/T',1X,'X2/lgp',1X,'Tlgp=(AC+B)& +A*T-Clgp'/1X,66(1H-)/6X,'X1',8X,'X2',8X, & 'Y',7X,'Ycalc',7X,'DIFF'/1X,66(1H-)) DO 40 I=1,N YCALC=B0 DO 30 J=1,M 30 YCALC=YCALC+A(J,M1)*X(I,J) DIFF=YCALC-X(I,M1) 40 WRITE(6,50) (X(I,J),J=1,M1),YCALC,DIFF 50 FORMAT(3X,10G10.4) WRITE(6,60) (XV(I),I=1,M1) 60 FORMAT(/1X,'Xv',10G10.4) WRITE(6,70) (S(I),I=1,M1) 70 FORMAT(/1X,'Sv',10G10.4) WRITE(6,80) 80 FORMAT(/5X,'Regression equation is:',/1x,60(1H-)) WRITE(6,90) B0,(A(I,M1),I=1,M) 90 FORMAT(1X,'Tlgp=',G10.4,'+',G10.4,'*T','+', & G10.4,'*lgp') WRITE(6,100) QT,U,QE,SE,U/QT,R,F,AA,BB,CC 100 FORMAT(/5X,'Qt=',G18.6,/5X,'U=',G18.6,/5X,'Qe=',G18.6,/5X, & 'Se=',G18.6,/5X,'N=',G18.6,/5X,'R=',G18.6,/5X,'F=', & G18.6,/5X,'A=',G18.6,/5X,'B=',G18.6,/5X,'C=',G18.6) END SUBROUTINE MLR(M,N,XV,S,B0,QT,U,QE,SE,R,F) （5 分） DOUBLE PRECISION XV,S,A,X,B0,QT,U,QE,SE,R,F DIMENSION XV(M+1),S(M+1),ALI0(50) COMMON /AX/A(6,7),X(6,7) M1=M+1 DO 20 I=1,M1 XV(I)=0.D0 DO 10 K=1,N 10 XV(I)=XV(I)+X(K,I) 20 XV(I)=XV(I)/DBLE(N) DO 40 I=1,M1 S(I)=0 DO 30 K=1,N 30 S(I)=S(I)+(X(K,I)-XV(I))**2 40 S(I)=DSQRT(S(I)) DO 50 I=1,M ALI0(I)=0 DO 50 K=1,N 50 ALI0(I)=ALI0(I)+(X(K,I)-XV(I))*(X(K,M1)-XV(M1)) DO 60 I=1,M1 DO 60 K=1,N 60 X(K,I)=(X(K,I)-XV(I))/S(I) DO 70 I=1,M A(I,I)=1.0 DO 70 J=I+1,M1 A(I,J)=0 DO 70 K=1,N 70 A(I,J)=A(I,J)+X(K,I)*X(K,J) DO 80 I=1,M-1 DO 80 J=I+1,M 80 A(J,I)=A(I,J) WRITE(6,11) M 11 FORMAT(4X,'Normal equation A(m,m+1)',2x,'m=',i2/1x,60(1H-)) DO 90 I=1,M 90 WRITE(6,33) (A(I,J),J=1,M1) CALL GS(A,M,M1,1.D-6) WRITE(6,22) WRITE(6,33) (A(I,M1),I=1,M) B0=0.D0 U=0.D0 DO 100 I=1,M A(I,M1)=A(I,M1)*S(M1)/S(I) B0=B0+A(I,M1)*XV(I) 100 U=U+A(I,M1)*ALI0(I) B0=XV(M1)-B0 QT=S(M1)*S(M1) QE=QT-U SE=DSQRT(QE/(N-M-1))
R=DSQRT(1.DO-QE/QT) Regression equation is (5分) F=N-M-】)LU/OEM FORMAT(/4X,'Roots of mormal equation ,5X,'A(i,m+1),i=1,m'& T1gp=2077. +6.081 *T+110.2 /1x,601H- FORMAT(1X,10G12.4) R= .999999 RETURN Fs 686329. END 6.08138 B= 1407.22 SUBROUTINE GS(A.N.M,EPS) (5分) = -110.160 (10分) DOUBLE PRECISION EPS,A DIMENSION A(6,7) 1.计算原理（化学原理和计算方法）(15分) DO 50K=LN 放射性衰变系列 P=0.0 D0101=K,N 2U4→20Th点→2Ra—→… IFOPGE ABS(AOLKINGOTO I0 P=ABS(A(LK)) 衰变方式服从连串反应动力学， L=I 其动力学方程式为一个常微分方程组： CONTINUE IF(P.LT.EPS)STOPNO SOLUTION IF(L.EQ.K)GOTO 30 d凹=-k,U叮 DO20 EK M (10分) TMP=A(L,J) A(LJA(KJ) d=k心]-kT A(K JETMP CONTINUE P=L0A(KK) 利用RngK法解此方程组。可符2过U和20Th这两种核素在10a内每释2a DO40 K+LM A(KJ)=A(K.J)*P 的质量变化 (5分) DO40 BK+IN A(1,JFA(IJ)-A(L,KA(K,J) 2.程序框图(20分》 CONTINUE 50 CONTINUE 开始 DOMK=2N I=M-K DO 60 J=I+1,N 输入：[U,Th初值234,230 A(L,M)=A(L,M)-A(LJ)*A(J,M) 时间间隔H, 计算精度PS (5分) 60 CONTINUE RETURN 计算：积分步长H0(5分) END 4.运行结果。(15分) 调用变步长Runge--Kuta法子程序解动力学方程组计算时刻t的U, [Th值 (其中调用计算动力学方程组的子程序计算DY值)(5分)
3 R=DSQRT(1.D0-QE/QT) F=(N-M-1)*U/QE/M 22 FORMAT(/4X,'Roots of mormal equation',5X,'A(i,m+1),i=1,m' & /1x,60(1H-)) 33 FORMAT(1X,10G12.4) RETURN END SUBROUTINE GS(A,N,M,EPS) （5 分） DOUBLE PRECISION EPS,A DIMENSION A(6,7) DO 50 K=1,N P=0.0 DO 10 I=K,N IF(P.GE.ABS(A(I,K)))GOTO 10 P=ABS(A(I,K)) L=I 10 CONTINUE IF(P.LT.EPS)STOP'NO SOLUTION' IF(L.EQ.K)GOTO 30 DO 20 J=K,M TMP=A(L,J) A(L,J)=A(K,J) A(K,J)=TMP 20 CONTINUE 30 P=1.0/A(K,K) DO 40 J=K+1,M A(K,J)=A(K,J)*P DO 40 I=K+1,N A(I,J)=A(I,J)-A(I,K)*A(K,J) 40 CONTINUE 50 CONTINUE DO 60 K=2,N I=M-K DO 60 J=I+1,N A(I,M)=A(I,M)-A(I,J)*A(J,M) 60 CONTINUE RETURN END 4.运行结果。（15 分） Regression equation is: （5 分） ------------------------------------------------------------ Tlgp=-2077. + 6.081 *T+ 110.2 *lgp R= .999999 F= 686329. A= 6.08138 B= 1407.22 C= -110.160 （10 分） B 1.计算原理（化学原理和计算方法）（15 分） 放射性衰变系列 234 92 U⎯⎯→k1 23090Th⎯⎯→k2 22684 Ra ⎯⎯→…… 衰变方式服从连串反应动力学。 其动力学方程式为一个常微分方程组： [U] [Th] d d[Th] [U] dd[U] 1 2 1 k k t k t = − = − （10 分） 利用 Runge-Kutta 法解此方程组，可得 U 23492 和 Th 23090 这两种核素在 10 a 内每隔 2 a 的质量变化。 （5 分） 2.程序框图（20 分） 开始 输入：[U]，[Th]初值 234，230 时间间隔 H， 计算精度 EPS （5 分） 计算：积分步长 H0 （5 分） 调用变步长Runge-Kutta法子程序解动力学方程组计算时刻t的[U]， [Th]值 （其中调用计算动力学方程组的子程序计算 DY 值） （5 分）
D060=,N 60 YO(J)Y2(J) WRITE(6,22YT0,H,(Y2(I),I=1,N) 80 CONTINUE 列表输出：两种核素在10a内每隔2a的质量变化U,[Th) WRITE(6,44) (5分) 11 FORMAT(/IX,N=,I2,2X,T0=,F6.1,2X,T1=,F8.1,2X,& H0=,G8.3,1XEPs=,G10.3 结束 1X,661H-/7X,/a,5X,H/a,10X,", 12X.1h/1X.661H-》 3.源程序(20分) 22 F0RMAT1X2F8.2,2X,2G15.6) PROGRAMMAIN (10分) 33 FORMAT(/IX,'FAILED TO FIND STABLE SOLUTION IN MAX IT/) IMP鬥LICIT REAL◆8A-H.O-Z 44 FORMAT(IX66H-1 DIMENSION YO(10),YI(10),Y2(10) END OPEN(6,FILE-2010TESTB.TXT,STATUS-UNKNOWN) DATA Y0/234.0230.0.80.0/ SUBROUTINE FCYDY (5分) N=2 IMPLICIT REAL*8(A-H,O-Z) 10=0.0 REAL◆8KI,K2 T1=10.0 DIMENSION YU0)DY(10) H0=2.0 DATA K12.8E-6/,K2/8.7E-6d EPs=1.0E-6 DY(IE-KIY(1) H=Ho DY2)=k1*Y1yK2*Y(2) WRITE(6,11)N,TO,TI,HO,EPS RETURN WRITE(6,22)TO,H,(YO(I),I=1,N) END M=IDINT(CTI-TOVHO) D08010=1,M SUBROUTINE RK4(XO,YO,Y,H,N,M) (5分 K=l IMPLICIT REAL8(A-HO-Z) CALL RK4(TO,YO,Y1,HO,N,K) DIMENSION YO(N),Y(N) D030KK=1,10 DIMENSION Y2(20)DY(20)RK(20) K=K+K X=X0 H=HO/K DOSELN CALLRK4(TO,Y0,Y2,H,N,K) 5 Y(I=YO(I) ES=0.0 D050L=1,M D010J=1N CALL F(YDY) 10 ES=ES+ABS((Y2(J)-Y1(J)/Y2(J)) D0101=1,N IF(ES.LT.EPS)GOTO 50 RK(IFDY(D) D020J=1,N 10 Y2(I=Y(IH+0.5*H*DY(I) Y10=Y2U) CALL F(Y2,DY) 30 CONTINUE D0201=1,N WRITE(6,33) RK(IFRK(I)+2*DY(I) STOP 20 Y2(I=Y(IH+0.5H+DY(D) T0=T0+H0 CALL F(Y2,DY)
4 3.源程序（20 分） PROGRAM MAIN （10 分） IMPLICIT REAL *8(A-H,O-Z) DIMENSION Y0(10),Y1(10),Y2(10) OPEN(6,FILE='2010TESTB.TXT',STATUS='UNKNOWN') DATA Y0/234.0,230.0,8*0.0/ N=2 T0=0.0 T1=10.0 H0=2.0 EPS=1.0E-6 H=H0 WRITE(6,11)N,T0,T1,H0,EPS WRITE(6,22)T0,H,(Y0(I),I=1,N) M=IDINT((T1-T0)/H0) DO 80 I0=1,M K=1 CALL RK4(T0,Y0,Y1,H0,N,K) DO 30 KK=1,10 K=K+K H=H0/K CALL RK4(T0,Y0,Y2,H,N,K) ES=0.0 DO 10 J=1,N 10 ES=ES+ABS((Y2(J)-Y1(J))/Y2(J)) IF(ES.LT.EPS)GOTO 50 DO 20 J=1,N 20 Y1(J)=Y2(J) 30 CONTINUE WRITE(6,33) STOP 50 T0=T0+H0 DO 60 J=1,N 60 Y0(J)=Y2(J) WRITE(6,22)T0,H,(Y2(I),I=1,N) 80 CONTINUE WRITE(6,44) 11 FORMAT(/1X,'N=',I2,2X,'T0=',F6.1,2X,'T1=',F8.1,2X, & 'H0=',G8.3,1X,'EPS=',G10.3//& 1X,66(1H-)/7X,'t/a',5X,'H/a',10X,'U', & 12X,'Th' /1X,66(1H-)) 22 FORMAT(1X,2F8.2,2X,2G15.6) 33 FORMAT(/1X,'FAILED TO FIND STABLE SOLUTION IN MAX IT'/) 44 FORMAT(1X,66(1H-)) END SUBROUTINE F(Y,DY) （5 分） IMPLICIT REAL *8(A-H,O-Z) REAL *8 K1,K2 DIMENSION Y(10),DY(10) DATA K1/2.8E-6/,K2/8.7E-6/ DY(1)=-K1*Y(1) DY(2)=k1*Y(1)- K2*Y(2) RETURN END SUBROUTINE RK4(X0,Y0,Y,H,N,M) （5 分） IMPLICIT REAL *8(A-H,O-Z) DIMENSION Y0(N),Y(N) DIMENSION Y2(20),DY(20),RK(20) X=X0 DO 5 I=1,N 5 Y(I)=Y0(I) DO 50 L=1,M CALL F(Y,DY) DO 10 I=1,N RK(I)=DY(I) 10 Y2(I)=Y(I)+0.5*H*DY(I) CALL F(Y2,DY) DO 20 I=1,N RK(I)=RK(I)+2*DY(I) 20 Y2(I)=Y(I)+0.5*H*DY(I) CALL F(Y2,DY) 列表输出：两种核素在 10 a 内每隔 2 a 的质量变化[U]，[Th] （5 分） 结束
D030=1,N RK(IERK(IH2.DY(I) 2. 程序框图(20分) 30 Y2(I)=Y(I)+HDY(I) CALL F(Y2.DY) 开始 D0401=1,N RK(I=RK(D+DY(D) Y(I=Y(I)+RK(I*H/6 输入：x1,x11,r1,r2 (4分)N) CONTINUE RETuRN X2=1-x1;x22=1-x11 (4分) END 4运行结果(15分，每列7.5分) 调用GAUSS消去法计算A,B(4分) N=2T0= 0T1= 10.0H0=2.00 EPS=.100E-05 常数A,B及x1=0200,x2=0.800代入原方程组计算出相应的常数 t H/a U r11,22(4分) ,00 2.00 234.000 230.000 2.00 1.00 233999 229.997 4.00 229.995 3. 结束 233.99 源程序(20分 6.00 1.00 233.996 229.992 PROGRAM MAIN (10分) 8.00 100 233.995 229.989 DIMENSION A(2 3 10.00 1.00 233.993 229.987 DOUBLE PRECISION A,EPS,X11,X21,R11,R21,X12,X22,R12,R22 INTEGER I,N,M OPEN(6,FILE-2010TESTC.Txt.STATUS-UNKNOWN) 计算原理（化学原理和计算方法）(15分) N=2 按照玛格物斯方程，二组分液态混合物中两个组分的活度系数，上与浓度，2间存在 M=N+1 如下关系： EPS=1.0D-06 X11=0.486 In=[A+2(B-A)xx (1) X21=1-X11 R11=0.813 R21=0.782 n5=[B+2(A-B)x] X12=0.200 22=1-X12 将上述方程组(1)变形为： A(1,M)=LOG(R11) A(2,M)=LOG(R21) lnr=[x(1-2x】A+2xx2B (2) A41,1=X21*2)(1-2*X11) Inr =2x,XA+[x(1-2x2)]B AM2,1=2*X21*X112 AM1,2=2*X11*X21**2 将=0.486， 1=0.813,=0.782,=1-1 代入(2)。买用高临消去法星方程阁 A(2,2=(X11*2)*(1-2X21) 可得到4，B的值.然后将x1=0.200,=1-1代入(1)，可求得相应的丙酮和氯仿的活度系数 WRITE(6,YORIGINAL LINEAR EQUATIONS: n,乃的值. WRITE(6.11(A(JEI.N+11EI.N)
5 DO 30 I=1,N RK(I)=RK(I)+2*DY(I) 30 Y2(I)=Y(I)+H*DY(I) CALL F(Y2,DY) DO 40 I=1,N RK(I)=RK(I)+DY(I) 40 Y(I)=Y(I)+RK(I)*H/6 50 CONTINUE RETURN END 4.运行结果（15 分，每列 7.5 分） N= 2 T0= .0 T1= 10.0 H0=2.00 EPS= .100E-05 ------------------------------------------------------------------ t/a H/a U Th ------------------------------------------------------------------ .00 2.00 234.000 230.000 2.00 1.00 233.999 229.997 4.00 1.00 233.997 229.995 6.00 1.00 233.996 229.992 8.00 1.00 233.995 229.989 10.00 1.00 233.993 229.987 ------------------------------------------------------------------ C 1． 计算原理（化学原理和计算方法）（15 分） 按照玛格勒斯方程，二组分液态混合物中两个组分的活度系数 r1，r2与浓度 x1，x2间存在 如下关系： ㏑ 2 1 12 r A B Ax x =+ − [ 2( ) ] （1） ㏑ 2 2 21 r B A Bx x =+ − [ 2( ) ] 将上述方程组（1）变形为： r x x A x x B r x x A x x B ln 2 [ (1 2 )] ln [ (1 2 )] 2 2 2 1 2 2 2 1 2 1 1 2 2 1 2 = + − = − + （2） 将 x1＝0.486， r1＝0.813， r2＝0.782，x2＝1-x1 代入（2），采用高斯消去法解方程组， 可得到 A，B 的值。然后将 x1＝0.200，x2＝1-x1代入（1）,可求得相应的丙酮和氯仿的活度系数 r1，r2的值。 2． 程序框图（20 分） 3． 源程序（20 分） PROGRAM MAIN （10 分） DIMENSION A(2,3) DOUBLE PRECISION A,EPS,X11,X21,R11,R21,X12,X22,R12,R22 INTEGER I,N,M OPEN(6,FILE='2010TESTC.Txt',STATUS='UNKNOWN') N=2 M=N+1 EPS=1.0D-06 X11=0.486 X21=1-X11 R11=0.813 R21=0.782 X12=0.200 X22=1-X12 A(1,M)=LOG(R11) A(2,M)=LOG(R21) A(1,1)=(X21**2)*(1-2*X11) A(2,1)=2*X21*X11**2 A(1,2)=2*X11*X21**2 A(2,2)=(X11**2)*(1-2*X21) WRITE(6,*)'ORIGINAL LINEAR EQUATIONS:' WRITE(6,11)((A(I,J),J=1,N+1),I=1,N) 开始 输入：x1,x11, r1,r2 （4 分）N） X2=1-x1; x22=1-x11 （4 分） 调用 GAUSS 消去法计算 A，B （4 分） ） 输出：r11，r22 （4 分） 结束 常数 A，B 及 x1=0.200， x2=0.800 代入原方程组计算出相应的常数 r11，r22 （4 分）