常微分方程数值解(第六章)讲述.ppt

§6.4 单步方法的收敛性和稳定性 初值问题的数值解法是经过某种离散化过程导出的,因此需要对数值解法进行定性分析.本节主要讨论单步方法的收敛性与稳定性. 6.4.1 单步方法的收敛性 定义6.1 定理6.1 [证明] Lipschiz条件： |f(x,y2)-f(x,y1)|≤L|y2-y1| 等比数列公比为 在收敛性的讨论中,我们已假定差分方程是精确求解的,但实际情况并非如此.例如,初始数据可能存在误差,计算过程中也不可避免地产生计算舍入误差,这些误差的传播和积累都会影响到数值解.那么实际计算得出的数值解能否作为精确解的近似呢?这取决于计算误差是否可控制,这就是数值方法稳定性的问题. 定义6.2 6.4.2 单步方法的稳定性 （-2.78，0） 4 四阶R-K方法 （-2，0） 2 改进Euler方法 （-2.51，0） 3 三阶R-K方法 （-∞，0） 2 梯形方法 （-2，0） 2 二阶R-K方法 （-2，0） 1 Euler方法 确定区间 方法的阶数 方法 确定区间 方法的阶数 方法 综上所述，收敛性是反映差分公式本身的截断误差对数值解的影响；稳定性是反映计算过程中舍入误差对数值解的影响。单步显式方法的稳定性与步长密切相关，在一种步长下是稳定的差分公式，取大一点步长就可能是不稳定的，只有既收敛又稳定的差分公式才有实用价值。