北京化工大学2002-2003学年第一学期 《计算化学》期末试卷答案（上机部分） 注：上机考试总分60分。 1.计算原理（化学原理和计算方法）(12分) 根据热力学知识，有 K 由上式得：KR=Kexp(S) 此题中，已知热烙与温度的关系：△H~T 可推导出识T关系，可根据离散点求积的方法即：用S1即sOn求积方法 RT2 (其中调用Lagrange插值法子程序计算任意点的函数值)计算S,然后通过S 求出773K下的平衡常数。 2.程序框图(18分) 开始 输入：热焓与温度△H~T的实验数据 X①.Y① XⅪ=X0,Y=YXI*XⅫ）=l,N 输入：积分上、下限T2,T1及KP1 两次调用Simpson积分法子程序计算面积S (其中调用Lagrange插值法子程序计算任意点的函数值) 计算KP2=KPI*EXP(S/R) 输出：KP2北京化工大学 2002-2003 学年第一学期 《计算化学》期末试卷答案（上机部分） 注：上机考试总分 60 分。 A 1．计算原理（化学原理和计算方法）（12 分） 根据热力学知识，有 ∫ = Δ = 2 1 1 2 2 ln T T p p dT S RT H K K 由上式得： exp( ) 2 1 K p = K p ⋅ S 此题中，已知热焓与温度的关系：ΔH ~ T 可推导出 T RT H ~ 2 Δ 关系，可根据离散点求积的方法即：用 Simpson 求积方法 （其中调用 Lagrange 插值法子程序计算任意点的函数值）计算 S，然后通过 S 求出 773K 下的平衡常数。 2．程序框图（18 分） 开始 输入：热焓与温度 ΔH ~ T 的实验数据 X(I),Y(I) （I=1，N） 输入：积分上、下限 T2,T1 及 KP1 两次调用 Simpson 积分法子程序计算面积 S （其中调用 Lagrange 插值法子程序计算任意点的函数值） 计算 KP2=KP1*EXP(S/R) 输出：KP2 结束 XI=X(I),Y(I)=Y(I)/(XI*XI) (I=1,N)