小鸭子.docVIP

第九章 常微分方程初值问题的数值解法 第一部分 内容提要 数值解的一般概念 常微分方程初值问题的数值解是指通过一定的近似方法得出准确解在一列离散点上的近似值。数值解的特征是步进式，即在点的近似值是由等若干点处的近似值的信息给出的递推公式。若依赖于前面步的值，则称为步法；称为单步法。 利用在的精确解借助某种算法计算出，则称为该方法的局部截断误差。如果一个算法的局部截断误差是，则称该方法是阶的；而利用数值解得到的与微分方程的精确解之差称为整体截断误差，即是该数值方法的误差。 对于固定的，取，用某种算法得到，如有=0，则称该方法是收敛的。注意，因是固定的，随着，数值解的步数。 在实际计算时由于舍入误差不可避免，实际得到数值解是，稳定性即研究是否随着计算步骤的增加而增加。通常所提的稳定性是通过模型方程来讨论的。若当某一步有舍入误差时，在以后的计算中误差不会逐步扩大，则称这种稳定性为绝对稳定性。 简单单步法及其收敛性、稳定性 Euler法的局部截断误差为，整体截断误差为，即一阶收敛。对于模型问题，当时，Euler法是数值稳定的。 隐式Euler法的误差与Euler法相同，但是无条件稳定：即对任意步长，隐式Euler法都是稳定的。 梯形法的误差比Euler法高一阶，也是无条件稳定的。 改进Euler法是一种预测－校正方法：——Euler法预