1.5线性多步法汇编.pptVIP

（二） 任一显式多步法和隐式多步法都可搭配成预——校算法 (四) 多步法计算的问题 1. 需附加初值 为保持精度，计算附加初值时要节点加密或采用同样精度的Runge-Kutta法 2. 如何选择阶p(或步数k): 要考虑光滑性和稳定性及总的工作量；高阶多步法的绝对稳定域小， 3. 如何选择步长h: 可根据精度要求由粗到细逐渐调整，计算过程中可以改变步长，但计算过程复杂 对于多步法而言，表头的值式非常重要的，如果表头的值精度低，则用精度再高的多步格式也无济于事，要得到高精度的表头值，除了用高阶单步法外，也可从较低阶单步法如欧拉法，配合缩短步长以及外推技巧( 外推法)等计算得高精度的开始值，例如设我们利用欧拉法计算得开始点 处的值为 ，由于欧拉法整体截断误差为 ，故可写成 (五) Richardson 外推法 利用 为计算步长，则得 则 因此，取 作为 的近似值，误差将提高一阶为 ，如果用 作为的近似值，整体截断误差为 ，则有 利用 为计算步长，有 消去 项，则得 以它作为 的近似值的截断误差 ，精度提高一阶，这就是利用外推技巧( 外推法)提高格式精度的方法。 谢 谢 作业： 作业要求：写出程序，列表或用图形显示结果，并给出图或表所说明的结果 P20-22，例2.1,2.2和习题 3 * 机动 目录 上页 下页 返回 结束 1.5 线性多步法 1.5 线性多步法 一、线性差分方程 称序列 满足的方程 为 阶线性差分方程，序列 是差分方程的解。 为确定差分解，需给定 个初值 引入记号 (1.40) 依次递推， 可用 的一阶、二阶直到 阶差分表示，所以差分方程(1.40)的最高阶为k.即(1.40) 称为k阶差分方程。它是k阶线性常微分方程的离散化。二者有许多平行的性质。 二、线性多步法 欧拉法和梯形法的回顾 求未知函数在 的近似值 ，基本思想：对于积分表达式 将被积函数 用水平直线 或连接 ， 两点的直线 代替。 单步法：计算节点 的近似值 仅用到前一节点的值 ，故从初值 可算出以后各节点的值。 用它近似 中的 ，则得到积分近似值 像欧拉格式一样，我们可得近似求解格式如下： 曲线，如我们利用 Lagrange 插值可得经过 然而为了近似 上的曲线 也可用多点插值 上面的格式与欧拉格式 不同之处在于增加了包括 的两项,即由 计算 ，与欧拉格式仅由前面一点 计算 的这种单步法不同，格式(1.51)称为多步法。 为提高精度，构造线性多步法，一般形式为 其中 (1.41) 用(1.41)计算 需用到前 k 个节点的值 故称为多步法或 k 步法。又因(1.41)关于 是线性的，故称为线性多步法。 注意：需附加初值 的解 在 处的近似值 ，或者说给出表头 现在，设已给出常微分方程初值问题 (一) 数值积分法 适当取k+1个节点，用 的k次Lagrange插值多项式 近似 将方程 写成积分形式，比如在 上积分得 (1.42) 不同的插值节点导出不同的多步法 (1)Adams外插法 （显式多步法） 取 为节点，构造 其中 是插值余项，代入(1.42)得 舍去余项 (1.43) (1.44) 下面给出(1.44)的具体形式。由插值节点等距，且被插值 点 靠近最后一个节点 ，故用牛顿向后插值 公式 引进记号 则 因此 (1.45) 式中 式(1.45)就是著名的 外插公式, 即为Euler法 定义生成函数 两端展开成幂级数有 例如 则Adams 外插公式可改写为 其中 的值见表2.2