数值计算方法第九章 常微分方程数值解法.pptVIP

其绝对稳定域为三、隐式Runge-Kutta 方法m级隐式R–K 方法的一般形式其中系数的确定方法同显式R–K 方法完全类似(1) 一级二阶的隐式中点方法：(2) 二级四阶的隐式R-K 方法：N级隐式R-K 法可以达到2N 阶缺陷：需要求解非线性方程（组）一、k步线性多步法§4 线性多步法与预估-校正格式/*Linear Mutistep Method and Predictor-Corrector Format*/ 所谓的线性多步法，指的是某一步解的公式不仅与前一步的值有关，而且与前面若干步解的值有关的方法。对初值问题两边积分得将换为节点取节点，构造的q+1 个点的Lagrange 插值多项式：多步显式公式其中记若函数值已知，则得r+1 步显式方法如时，可得二步显式阿达姆斯（Adams ）格式其中Adams 显式公式的局部截断误差：由Lagrange 插值余项知其中（第二积分中值定理）q+1 阶方法取节点，构造的q+1 个点的Lagrange 插值多项式：多步隐式公式其中记则得到r+1 步q+1 阶的隐式方法如时，可得二步隐式阿达姆斯（Adams ）格式梯形公式* 第九章常微分方程数值解法/*Numerical Method for Ordinary Differential Equations*/ 本章主要介绍一阶常微分方程初值问题的数值解法。初值问题及其数值解的概念§1 引言常用的一些解析解法：常数变易法、Lapalace 变换等分离变量法、变量代换、一阶常微分方程初值问题：对于初值问题，如果在下列区域内连续：（解的存在唯一性）且关于满足Lipschitz 条件，即存在常数，使则初值问题存在唯一解，且解是连续可微的。所谓数值解是指：在解的存在区间上取一系列点逐个求出的近似值等距节点：步长初值问题的解析解及其数值解的几何意义：初值问题的解表示过点的一条曲线初值问题的数值解表示一组离散点列可用拟合方法求该组数据的近似曲线积分曲线§2 Euler 方法Euler 方法的导出将在点处进行Taylor 展开略去项：然后用代替，即得称上述公式为向前Euler 公式。若将在点处进行Taylor 展开略去项：然后用代替，即得称上述公式为向后Euler 公式。向后Euler 公式为隐式格式，需要利用迭代法求解解：向前Euler 公式：例1：分别利用向前和向后Euler 方法求解初值问题的数值解（取步长为）向后Euler 公式：具体计算结果见教材表 。常微分方程数值解法的稳定性设一个数值方法以定步长求解实验方程得到线性差分方程的解。当时，若，则称该方法对步长为绝对稳定的；否则称为不稳定的。将数值方法应用于实验方程，若对一切都是绝对稳定的，则称区域为该方法的绝对稳定域。上述定义表明，若数值方法可使任何一步产生的误差在后面的计算中都能逐步削弱，则该方法为绝对稳定。例如，对于向前Euler 法：将其应用于实验方程当时，误差将逐步减弱，故此时方法稳定。向前Euler 法绝对稳定域：当因有误差变为时，则有单步方法的局部误差和阶单步法的一般形式隐式单步法通常称为增量函数显式单步法设是准确的，用某种方法计算时产生的截称为某方法在点的整体截断误差断误差，称为该方法的局部截断误差，即其中为自然数，则称该方法是阶的或具有阶精度。如果给定方法的局部截断误差为如果一个阶单步方法的局部截断误差为则称为该方法的局部截断误差的主项。如向前Euler 方法的局部截断误差一阶方法Euler 方法的误差分析对初值问题中的微分方程两端在区间上积分如果用左矩形公式计算右端积分，并令其中上述等式中如果用代替，即得向前Euler 格式。其局部截断误差为设关于和均满足Lipschitz 条件，即和其中而整体截断误差为注意到对于初值问题，如果关于满足（向前Euler 方法的整体截断误差）Lipschitz 条件，为对应的Lipschitz 常数，当时，向前Euler 方法的数值解一致收敛于初值问题的精确解，且整体截断误差满足估计式如果，Euler 方法的整体截断误差为一、Runge-Kutta 方法的基本思想§3 龙格-库塔（Runge-Kutta ）方法显式单步法的一般形式：R-K 方法是利用一些点的线性组合构造增量函数，使得相应方法的局部截断误差的阶数尽可能高。二阶Runge-Kutta 方法确定参数，使得与在点的Taylor 展开式有尽可能多的相同项。比较两式的相同项得方程组有无穷多解若取其一组解则得到改进的Euler 公式（二阶方法）若取其另一组解则得到二阶的Heun （休恩）公式（见教材）。二、显式Runge-Kutta 方法及其稳定性和设是一个正整数，代表使用函数值的个数，是一些特定的权因子（均为实数），则称下列方法（公式）为初值问题的m级显式Runge–Kutta 公式