数值分析--第08章常微分方程数值解法.pptVIP

数 值 分 析 Numerical Analysis 机械与汽车工程学院 主讲人：孔胜利 kongsl@spu.edu.cn 2012-09-01 第8章 常微分方程的数值解法 Euler方法 Runge-Kutta方法 线性多步法 边值问题的数值解法 常微分方程的数值解法 工程实际问题的解决往往需要求解常微分方程的定解问题。 对于初值问题 若函数 f (x, y) 满足 1）在矩形区域D内连续； 2）在区域D内对y满足Lipschitz条件，即存在常数L。对于D中的任意 x和y的两个值y1和y2，有 则此初值问题存在唯一的解 y = y (x) ，且满足初值条件。 解析方法只能求解一些特殊类型的方程，而应用问题中碰到的方程 主要靠数值解法。 常微分方程的数值解法 所谓数值解法，就是寻求y (x)在一系列离散节点 上的近似值 y1， y2, … yn , yn+1 ,…。 相邻两个节点的间距称为步长。 初值问题的数值解法有个基本特点，就是都采取“步进式”，即求解 过程顺着节点排列的次序一步一步地向前推进。描述这类算法，只 要给出用已知信息yn， yn-1, yn-2 ,…,计算 yn+1 的递推公式。 8.1 Euler方法 1、Euler公式的构造法 用差商代替微商（导数） 用向前差商代替微商得 或 我们用 y (xn) 的近似值即解在 xn 点的数值解 yn 代入上式，算出来的 y (xn+1) 的近似解yn +1，即得 称为Euler公式。 由于y (x0) = y0 已知，代入上式依次计算，可算成各离散点上的数值解。 Euler方法 Euler公式的构造法 后退的Euler公式 用向后差商代替微商得 称为后退的Euler公式 。 梯形公式 用中心差商代替微商得 称为梯形公式。 8.1 Euler方法 2、数值积分法 将y (xn)微分方程改写成积分的形式，并在 [xn, xn+1] 上积分，得 用左矩形公式计算右端的积分，得 从而 若用右矩形公式计算积分，可得后退的Euler公式； 若用梯形公式计算积分，可得梯形公式。 8.1 Euler方法 3、 Taylor展开法 将y (xn+1)在xn点展开，得 略去O ( h2 ) ，得 从而 为Euler公式。 8.1 Euler方法 若将y (xn)在xn+1点展开，得 略去O ( h2 ) ，得 从而 为后退的Euler公式。 结合Euler公式和后退的Euler公式即可得梯形公式。 例题 用Euler法解初值问题 取步长 ，计算到x=0.3（保留到小数点后4位）。 解： 因 ，故由Euler公式得 得到各节点处的数值解列表如下 几个基本概念 1、单步法和多步法 2、显式与隐式 3、自开始的与不是自开始的 单步法是自开始的，多步法不是自开始的 4、截断误差 局部截断误差和整体截断误差 8.2 Runge-Kutta方法 问题的提出 在构造数值积分时，由 ，围绕着如何用函数 值表示出 的近似值得出各种求积公式。 的近似值可以用若干 点处 之值的加权平均值来表示，用的点数越多代数精度越高。 二阶公式 三阶公式 四阶公式 经典R-K公式 例题 用经典R-K公式求解 解 由于 ，故求解公式为 具体计算过程为 即 同理，由 可得 以此类推，计算结果如下表 8.3 线性多步法 线性多步法的基本概念 求 时只用到前一步 的值，即为单步法； 事实上，在求 时，前面各点的数值解 均已经求 出来。如果充分利用这些现有的结果，便可获得较高精度的公式， 即我们构造形为 的数值公式，式中， 为与微分方程 无关的待定参数，即线性多步法。 线性多步法的主