第三章数值积分与常微分方程的数值解法.pptVIP

第三章数值积分与常微分方程的数值解法 第一页，共四十页，2022年，8月28日 3-1-1 梯形法——方法原理 基本思想： 复化求积，即从近似计算为出发点，用有限项 的求和计算来代替从而求出定积分的近似值。 定步长： 求f(x)在[a,b]上的定积分 x y y=f(x) a b xk-1 xk h Ik h——步长 第二页，共四十页，2022年，8月28日 变步长： N个区间，h，T1 2N个区间，h/2，T2 |T2-T1|EPS （xk-1,xk） xk-1/2 其中： 3-1-1 梯形法——方法原理 第三页，共四十页，2022年，8月28日 例：Debye-Einstein公式推导得到计算固体热容的公式为 其中： ?D为Debye温度，R为气体常数8.314JK-1mol-1 已知固体的Debye温度如下： Pb Ag Cu Al Fe KCl NaCl C ?D/*K 88 215 315 398 420 227 281 1910 求在50，100，298.15，500，1500K时，各固体的热容。 3-2-1-1 Simpson法——问题的提出 第四页，共四十页，2022年，8月28日 求积分 Simpson法是把积分区间分割成有限个小区间，在每个小区间上采用二次抛物线来近似被积函数f(x)的图形，近似求出小区间的面积，然后再将有限个小区间相加得到被积函数的近似值。 x y y=f(x) xi-1 xi+1 xi y=g(x) h h Si 定步长: 3-2-1-2 Simpson法——方法原理 第五页，共四十页，2022年，8月28日 变步长： 其中： 3-2-1-2 Simpson法——方法原理 第六页，共四十页，2022年，8月28日 判据： 3-2-1-2 Simpson法——方法原理 第七页，共四十页，2022年，8月28日 3-2-1-3 Simpson法——程序框图 Simp(A,B,EPS,S2,F) N=1, H=B-A, S1=0, T1=H*(F(A)+F(B))/2 DO K=1，N S=0 S=S+F(A+(K-1/2)*H) T2=(T1+H*S)/2, S2=T2+(T2-T1)/3, D=|S2-S1| |S2|1 D=|(S2-S1)/S2| DEPS No RETURN Yes No N=N*2 H=H/2 T1=T2 S1=S2 Yes 第八页，共四十页，2022年，8月28日 3-2-1-4 Simpson法——应用示例 开始 输入： Debye温度T（5），精度EPS，温度THETA 输出：固体的热容Cv 结束 调用Simpson积分法子程序计算式右方积分值S2 计算：XM=THETA/T(I) (I=1,N) 输入：积分上下限A=10-4，B=XM B=0 Yes No 固体的热容Cv=9R/XM**3*S2 显示程序 显示输出 第九页，共四十页，2022年，8月28日 3-1-3 –1 离散点数据的求积——方法原理 实验时，得不到变量间的关系式，只测量到（xi,yi） 的离散点数据。 x y a b 方法 ： 1.用插值程序求任意 点的函数值。 一元三点Lagrange插值： 2．用Simpson求积程 序计算[a,b]区间中离 散点下的面积。 第十页，共四十页，2022年，8月28日 Simp(M,A,B,X,Y,EPS,S2) N=1, H=B-A, （1）； S1=0, T1=H*(F(A)+F(B))/2 DO K=1，N S=0 （2）； S=S+F(A+(K-1/2)*H) T2=(T1+H*S)/2, S2=T2+(T2-T1)/3, D=|S2-S1| |S2|1 D=|(S2-S1)/S2| DEPS No RETURN Yes No N=N*2 H=H/2 T1=T2 S1=S2 Yes (1)调用Lagrange一元三点插值F(A),F(B) (2)调用Lagrange一元三点插值F(A+(K-1/2)*H) 3-1-3 –2 离散点数据的求积——程序框图 第十一页，共四十页，2022年，8月28日 纯气体的逸度由定义： （1） 代入（1） 积分，并取极低压力下气体视为理想气体，得 逸度： φ为逸度系数 （2） 例1： 实际气体逸度的计算 已知p~Vm数据 3-1-3 –3 离散点数据的求积——应用示例 第十二页，共四十页，2022年，8月28日 开始 输入：数据点数N，精度EPS，温度T 压力p和摩尔体积Vm的实验数据X(I),Y(I) （I=1，N） 输出：B, FI, FF 结束 调用离散点求积子程序计算（2）式右方积分值S