常微分方程的差分方法.pptVIP

为例对欧拉方法做稳定性分析。 先分析显式欧拉格式，对于（5 .1）有 （5.1） 这里简单地以模型方程（称为试验方程） 将上式反复递推后可得 * 第94页，共147页。 要使 有界，其充要条件为： 或 式中 * 第95页，共147页。 即 由于 λ 0，故有 可见如欲保证算法的稳定，显式欧拉格式的步长 h 的选取要受到（5.2）的限制。 的绝对值越大，则限制的h值就越小。 （5.2） * 第96页，共147页。 例如初值问题 的精确解为： 如取步长h = 0.05，作显式欧拉格式计算时，有 * 第97页，共147页。 按上式计算所得的结果列于 表7-8中。从表7-8中可见，这时显式欧拉方法是不稳定的。究其原因，就是由于h不符合（5.2）的规定： 再分析隐式欧拉格式，对于（5.1）有 由于 ，则恒有 * 第98页，共147页。 表7-8 * 第99页，共147页。 （5.3） 因此，隐式欧拉格式是绝对稳定（无条件稳定）的（对任何h 0）。 对于初值问题 故恒有 * 第100页，共147页。 同样取步长h = 0.05，作隐式欧拉格式计算时，有 按上式计算所得结果也列于表7－8中。从表中的数值可见，同样的步长 h = 0.05 ，隐式欧拉方法确实稳定的。因为这时它对任何h 0都是绝对稳定的。 稳定区域：能保证算法稳定的 的取值范围，称为稳 定区域。 对于（5.1）有 * 第101页，共147页。 显然，稳定区域越大，意味着该算法的稳定性越好。 由（5.2）可知，显式欧拉方法的稳定区域为： 由（5.3）可知，隐式欧拉方法的稳定区域为： * 第102页，共147页。 可见隐式欧拉方法的稳定性比显式欧拉方法的好。 可用与上述类似的方法来分析 k 步阿达姆斯显式和隐式方法的稳定性。对于典型方程，它们的稳定区域如分别记为 和 ，则其数值如表7-9所示，由表7－9可见： （1） 阿达姆斯隐式方法的稳定区域都比同阶 的显式方法的大，这是隐式方法最大的 优点； （2） k 越大（步数、阶数越大）时稳定区域就 越小，只有隐式欧拉方法是绝对稳定的。 * 第103页，共147页。 将其代入（4.2），且利用 ，可得： 顿向后插值多项式为 其余项为 * 第62页，共147页。 （4.3） 当 有界，且 时，可以估计出上式中 此即为截断误差项。于是可得 （4.4） * 第63页，共147页。 即为显式欧拉格式，其局部断截误差为 当k=0时： 当k=1时： * 第64页，共147页。 ……..（4.5） * 第65页，共147页。 式中 将（4.5）代入（4.4）后，可得： （4.6） 式中 * 第66页，共147页。 代入（4.6）后得 其中 * 第67页，共147页。 对应于不同的k值下的 和 值可算出，分别列于表7-3中。 从表7-3可知： 当k=2时 局部截断误差为 * 第68页，共147页。 表7-3 * 第69页，共147页。 当k=3时 （4.7） 其中k=3时的（4.7）称为阿达姆斯四步显示方法，因为它要用到前面四个节点上的 值，是一种常用的多步算法，其精度为四阶。 （4.2）中的 也可用牛顿内插多项式来逼近。设由 * 第70页，共147页。 共k+1个数据来构造一个牛顿内插多项式 ，则与上面类似推导可得： 局部截断误差为 对于不同k值下的 和 值可算出，分别列于表7-4中。 * 第71页，共147页。 从表7-4可知： 当k=0时 即为隐式欧拉格式，其局部截断误差为： 当k=1时 * 第72页，共147页。 即为梯形格式，其局部截断误差为 当k=2时 * 第73页，共147页。 ……（4.8） 当k=3时 其中k=3时的（4.8）称为三步四阶阿达姆斯隐式算法，它只用到前面三个节点上的 值，但其精度却为四阶。 * 第74页，共147页。 表7-4 *