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

3-2-2-1 Runge-Kutta法——问题的提出 例： 基元反应： k2 已知：[A]0=[B]0=1mol·L-1, [C]0=[D]0=[E]0=0 1mol·L-1, k1=1.0min-1， k2=0.25min-1， k3=0.5min-1， 求各组分浓度随时间的变化（0-20min）. Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. ● Euler法及其改进 ● Runge-Kutta法 ● 梯形法 ● Simpson法 ●离散点数据的求积 第三章 数值积分与常微分方程的数值解法 数值积分 常微分方程的数值解法 1.f(x)函数形式已知，但其积分不能表示成初等函数的闭合形式 2. f(x)函数形式未知，但其离散数据表已给出 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 3-1-1 梯形法——方法原理 基本思想： 复化求积，即从近似计算为出发点，用有限项 的求和计算来代替从而求出定积分的近似值。 定步长： 求f(x)在[a,b]上的定积分 x y y=f(x) a b xk-1 xk h Ik h——步长 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 变步长： N个区间，h，T1 2N个区间，h/2，T2 |T2-T1|EPS （xk-1,xk） xk-1/2 其中： 3-1-1 梯形法——方法原理 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 例：Debye-Einstein公式推导得到计算固体热容的公式为 其中： ?D为Debye温度，R为气体常数8.314JK-1mol-1 已知固体的Debye温度如下： ?D/*K Pb 88 Ag 215 Cu 315 Al 398 Fe 420 KCl 227 NaCl 281 C 1910 求在50，100，298.15，500，1500K时，各固体的热容。 3-2-1-1 Simpson法——问题的提出 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 求积分 Simpson法是把积分区间分割成有限个小区间，在每个小区间上采用二次抛物线来近似被积函数f(x)的图形，近似求出小区间的面积，然后再将有限个小区间相加得到被积函数的近似值。 x y y=f(x) xi-1 xi+1 xi y=g(x) h h Si 定步长: 3-2-1-2 Simpson法——方法原理 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 变步长： 其中： 3-2-1-2 Simpson法——方法原理 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 判据： 3-2-1-2 Simpson法——方法原理 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 3-2-1-3 Simpson法——程序框图 Simp(A,B,EPS,S2,F) N=1, H=B-A, S1=0, T1=H*(F(A)+F(B))/