称为绝对稳定域.PPT

9.4 单步法的收敛性与稳定性 9.4.2 绝对稳定性与绝对稳定域 9.5 线性多步法 9.5.1 线性多步法的一般公式 9.5.2 阿当姆斯显式与隐式公式 9.4.1 收敛性与相容性 数值解法的基本思想是，通过某种离散化手段将微分 方程（1.1）转化为差分方程，如单步法（2.10），即 （4.1） 它在 处的解为 ，而初值问题（1.1），（1.2）在 处的精确解为 ，记 称为整体截断 误差. 收敛性就是讨论当 固定且 时 的问题. 定义3 若一种数值方法（如单步法（4.1））对于固定 的 , 当 时有 ，其中 是（1.1），（1.2）的准确解，则称该方法是收敛的. 显然数值方法收敛是指 ，对单步 法（4.1）有下述收敛性定理： 定理1 假设单步法（4.1）具有 阶精度，且增量函数 关于 满足利普希茨条件 （4.2） 又设初值 是准确的，即 ，则其整体截断误差 （4.3） 证明 设以 表示取 用公式（4.1）求 得的结果，即 （4.4） 则 为局部截断误差，由于所给方法具有 阶 精度，按定义2，存在定数 ，使 又由式（4.4）与（4.1），得 利用假设条件（4.2），有 从而有 即对整体截断误差 成立下列递推关系式 （4.5） 反复递推，可得 （4.6） 再注意到当 时 最终得下列估计式 （4.7） 由此可以断定，如果初值是准确的，即 ，则（4.3） 式成立. 依据这一定理，判断单步法（4.1）的收敛性，归结为 验证增量函数 能否满足利普希茨条件（4.2）. 对于欧拉方法，由于其增量函数 就是 ，故当 关于 满足利普希茨条件时它是收敛的. 再考察改进的欧拉方法，其增量函数由（3.2）式 给出，这时有 假设 关于 满足利普希茨条件，记利普希茨常数为 ，则由上式推得 设限定 为定数），上式表明 关于 的利普希茨 常数 因此改进的欧拉方法也是收敛的. 类似地，不难验证其他龙格-库塔方法的收敛性. 定理1表明 时单步法收敛，并且当 是初值 问题（1.1），（1.2）的解，（4.1）具有 阶精度时，则 有展开式 所以 的充要条件是 ，而 ，于是可给出如下定义： 定义4 若单步法（4.1）的增量函数 满足 则称单步法（4.1）与初值问题（1.1），（1.2）相容. 以上讨论表明 阶方法（4.1）当 时与(1.1)，(1.2) 相容，反之相容方法至少是1阶的. 于是由定理1可知方法（4.1）收敛的充分必要条件是此 方法是相容的. 定义5 若一种数值方法在节点值 上大小为 的扰动， 于以后各节点值 上产生的偏差均不超过 ，则称 该方法是稳定的. 以欧拉法为例考察计算稳定性. 例4 考察初值问题 其准确解 是一个按指数曲线衰减得很快的函数， 如图9-3所示. 用欧拉法解方程 得 若取 ，则欧拉公式的 具体形式为 计算结果列于表9-4的第2列. 可以看到，欧拉方法的解 （图9-3中用×号标出）在准确值 的上下波动，计算过程明 显地不稳定. 但若取 则计算过程稳定. 图9-3 再考察后退的欧拉方法，取 时计算公式为 计算结果列于表9-4的第3列（图9-3中标以·号），这时计算 过程是稳定的. 例题表明稳定性不但与方法有关，也与步长 的大小 有关，当然也与方程中的 有关. 为了只考察数值方法本身. 通常只检验将数值方法用 于解模型方程的稳定性，模型方程为 （4.8） 其中 为复数，对一般方程可以通过局部线性