[所有分类]第五章 常微分方程的数值解法.pptVIP

第五章 常微分方程的数值解法 主要内容： 1、引言 2、欧拉方法 3、龙格－库塔方法 4、单步法的收敛性和稳定性 5、线性多步法 6、一阶方程组与高阶方程 第一节 引言 ●在常微分方程课程里面讨论的是一些典型方程求解解析解的基本方法。 ●然而在生产实践和科学研究中遇到的微分方程往往比较复杂，在很多情况下，不能给出解的解析表达式；有时候即时能用解析表达式来表示，又因为计算量太大而不实用，有时候一些是已经有了求解的基本方法的典型方程，但实际使用时也是有困难的。 ●以上情况说明用求解解析解的基本方法来求微分方程的解往往是不适宜的，甚至很难办到。 ●实际问题中，对于求解微分方程，一般只要求得到解的若干个点上的近似值或者解的便于计算的近似表达式。 ●本章研究微分方程的数值解法，而且着重讨论微分方程中最简单的一类问题——一阶方程的初值问题。 第一节 引言 1、一阶方程的初值问题 假定上式在区间[a, b]上存在唯一且足够 光滑的解y(x)。 ●所谓数值解法就是寻求解y(x)在一系列离散点，也称为节点处的值： 要计算出解函数 y(x) 在一系列节点 a = x0 x1… xn= b 处的近似值 第一节 引言 ●节点间距，即步长为： 通常采用等距节点，即hi = h (常数) ●等间距节点 ●在这些节点上采用离散化方法（通常用数值积分、微分、泰勒展开等）将上述初值问题化成关于离散变量的相应问题。把这个相应问题的解yn作为y(xn)的近似值。这样求得的yn就是上述初值问题在节点xn上的数值解。一般说来，不同的离散化导致不同的方法。 第二节 欧拉方法 一、欧拉法Euler 1、向前差商近似导数 第二节 欧拉方法 2、举例 例1 用欧拉法求初值问题 当h = 0.02时在区间[0, 0.10]上的数值解。 解：根据欧拉公式可以得到： 此外，可以得到方程的真解： 第二节 欧拉方法 求解过程如下： 第二节 欧拉方法 3、欧拉方法的几何意义 第二节 欧拉方法 从上述几何意义上得知，由Euler法所得的折线明显偏离了积分曲线，可见此方法非常粗糙即误差太大。 4、欧拉法的局部截断误差 （1）截断误差定义 在假设 yi = y(xi)，即第 i 步计算是精确的前提下，考虑的截断误差 Ri+1 = y(xi+1) ? yi+1 ，称为局部截断误差 如图所示：APi+1即为欧拉方法在xi+1点的截断误差 （2）如果某种方法的局部截断误差是 则称该方法具有p阶精度 第二节 欧拉方法 （3）则截断误差的大小？ 写出y(xn+1)的泰勒展开式： 由欧拉方法可以得到： 则上面两个公式相减得到： 第二节 欧拉方法 二、改进的欧拉法 一阶方程的初值问题与如下积分方程是等价的： 当x = x1时 可以借助于数值积分，求y(x1)的值 1、用矩形公式 第二节 欧拉方法 可以推导出： 用矩形法计算右端的积分与用欧拉法计出的结果完全相同 2、用梯形公式 则可以推导出： 第二节 欧拉方法 梯形公式的截断误差： 第二节 欧拉方法 和欧拉公式相比较，梯形公式在计算yi+1时候也只用到前一步的值yi，但是若yi已知，将yi带入公式求解时候，一般不能直接得到yi+1，而需要通过其他方法（比如迭代法）求解，所以梯形公式被称为隐式公式。 3、改进的欧拉方法 梯形公式是隐式的，一般用迭代法求解，计算量较大。实际中常将欧拉公式和梯形公式联合使用，先用欧拉公式得出一个y(xi+1)的近似值 称为预估值，然后对预估值使用梯形公式对它进行精确化，得到较为精确的近似值yi+1，称之为校正值，计算公式为： 这样的预估校正系统称为改进的欧拉方法。 第二节 欧拉方法 为了便于编写程序，常将上面的公式改写为如下式： 第二节 欧拉方法 4、举例 在区间[0, 1.5]上，取h = 0.1，求解。 解：（1）用欧拉法计算公式如下： （2）用改进欧拉法计算公式如下： 本题的精确解为 ， 可用来检验近似解 的精确程度。 计算结果如表： 第二节 欧拉方法 P93 例5-2 第二节 欧拉方法 5、欧拉两步公式 中心差商 第三节 龙格－库塔方法 一、龙格库塔法的基本思想 1、平均斜率 考察差商： 根据微分中值定理：在闭区间[a, b]上连续，开区间(a, b)上可导，则至少存在(a, b)上的一点，使得下式成立： 根据上面公式可以得到： 第三节 龙格－库塔方法 因此只要对 K* 提供一种算法，就可以求得数值解，根据该观点对欧拉法及改进的欧拉进行分析。 2、基于平均斜率对欧拉法和改进的欧拉法进行分析 （1）在欧拉公式中，是取了一个