5第五章常微分方程数值解法详解.ppt

1、 收敛性简介（以显式Euler格式为例 ） 如果对任意固定的 ，数值解 当 （ ）时趋于准确解 ，则称该算法是收敛的。 为理解收敛性的含义，仅考察下列模型方程 对于Euler格式，上述问题的Euler格式具有形式 从而数值解 因此当时 , 即Euler法收敛！ 返回 这样做的根据是： 1）对试验模型不收敛的方法，不可用； 2）一般的初始问题在其解的存在区域内，可局部线性化转化为试验方程。 单步法收敛的充要条件为收敛阶p≧1。 我们称算法是稳定的，如果在节点值 上有大小为 的扰动，而以后各节点值 上产生的偏差值均不超过 。 2、 稳定性简介 考察下列试验模型方程 进一步，若以后各节点产生的误差趋于零，则称算法是绝对稳定的。 1）对试验模型数值不稳定的方法，不可用； 2）一般的初始问题在其解的存在区域内，可局部线性化转化为试验方程。 对于Euler格式，上述模型方程的Euler格式为 当选择h使得 时，Euler格式绝对稳定！ 其绝对稳定区域为 （1）显式欧拉法 Numerical Analysis J