12.09级Ch8.4单步法进一步讨论.pptVIP

* §4 单步法的进一步讨论 4.1 收敛性与相容性 本节内容 理论分析：收敛性、相容性、稳定性、 提高精度（的途径）：外推法、变步长 对单步法： 问题 求解 当局部截断误差为 即 (a) (4.1)（单步法）是 p 阶方法. 时，且 关于变量u,v 满足Lipschitz条件，则数值方法 y(x)为精确解. 收敛性 结论 阶单步法是收敛的. (1.1)的准确解为y(x). ？ 事实上, 当 时， , 定义4 设单步法(4.1)生成 (1.1),(1.2)的数值解: ， ，则称数值（单步）方法(4.1)是收敛的. 若对任意固定的 都有 成立. (a) 结合（a）式，则结论 p阶方法的局部截断误差一般为： 相容性 将 在点 展开： 由 代入上式，得 从而 的必要条件是: 所以 得 定义5 结论 说明 方法是否相容，即方法与方程是否配套。当数值方法 两边对 取极限，得 数值方法(4.1)称为与初值问题(1.1),(1.2)是相容的，若 (4.6)式成立. 阶单步法全是相容的。前面讨论的方法全是相 容的方法。 (4.1)与(1.1),(1.2)相容时，对固定的k，差分方程（4.7） 首先，设数值方法(4.1)是相容的. 4.2 稳定性 其次，设 是(4.1)的准确数值解，y(x)是(1.1),(1.2) 持有界，总的误差能被控制，但当 时, 数值解 检验稳定性的方法 预备知识 的准确解. 计算机求解时，每步都有可能有舍入误差（迭代误差 与计算机引起的误差），因此设 是(4.1)的计算解. 定的. 定义 若在计算过程中某部引入的舍入误差，在以后的计算 （或传播）中被压缩、衰减，则称(4.1)是数值稳定的；若在传播中 误差被放大，则称(4.1)是数值不稳定的；若在传播中舍入误差保 可能不随k增大而趋于零，约定这种情况下方法(4.1)也是数值不稳 原因 (1) 方法简单. (2)方程(1.1)可局部线性化为 . 使用检验方程 ，其中 为常数. 令 则 事实上 , 则有 (2)方程(1.1)可局部线性化为 . 令 则 因此有试验方程： 若 ; 对该微分方程由 若 此时，称试验方程(4.8)是稳定的. 用数值方法(4.1)解试验方程，得 的有理函数或多项式函数. 其中 该微分方程的解 用数值方法(4.1)解试验方程，得 的有理函数或多项式函数. 其中 若方法(4.1)是p阶的，则 试验方程： 则 的一个逼近. 数值解的稳定性 数值解的稳定性 定义6 若计算 时有误差 ，而以后的计算全部准确，则 若在(4.9)中 ，则称方法(4.1)是绝对稳定的. 举例（几种单步法的稳定性） 显式Euler方法： 计算公式： 因此， 绝对稳定区域为 绝对稳定区间为 . o 由 即 得 将产生误差 讨论(4.1)的数值稳定性,只要讨论 是否小于1即可. 定区域；绝对稳定区域与实轴的交称为绝对稳定区间. 在复平面上，变量 满足 的区域称为(4.1)的绝对稳 单位开圆 隐式Euler方法 计算公式： 因此， 梯形方法 计算公式: 因此， 对任意的 有 . 说明隐式Euler方法对任意步长h是 稳定的. 对任意的h，有 ， 说明梯形方法对任意步长h是稳定的. 预估－校正Euler方法 经典Runge-Kutta方法 绝对稳定区间: . 绝对稳定区间为 . 说明 通过以上的例子分析可以看出，有些方法不是对任意h 稳定性的好坏依次为 都是绝对收敛的，必须取足够小的步长，使 落在绝对稳定