数值分析第九章常微分方程数值解法.ppt

* 以 yn, zn表示节点 xn 上 y(x), z(x) 的近似解，则其改进的Euler法具有形式 例 对于方程组 令 原方程组变为 * 预报 即 * 校正 即 * 关于高阶微分方程（或方程组）的初值问题，原则上总可以归结为一阶微分方程组来求解. 引进新变量 z=y?，即可化为一阶方程组的初值问题 从而可以运用前面介绍的算法来求解. ? 高阶微分方程的数值解法 例 对于下列二阶方程的初值问题 * ? 三阶Adams显式公式: F(x)用过节点 xn , xn?1 , xn ?2的2次插值函数 L2(x)代替 * 再离散化，得三阶Adams显式公式 * ? 三阶Adams显式公式截断误差 * ? 三阶Adams隐式公式: F(x)用过节点 xn+1 , xn , xn ?1的2次插值函数 L2(x)代替 * 再离散化，得三阶Adams隐式公式 * ? 三阶Adams隐式公式截断误差 * 上机作业 第342页 数值实验第1题 * * ? Milne公式 F(x)用过节点 xn , xn?1 , xn ?2的2次插值函数 L2(x)代替 再离散化，得Milne公式 * Milne公式是四阶四步显式方法, 截断误差为 * ? 基于Taylor展开方法能导出所有的常用线性多步法, 下面以线性二步法为例来说明这种方法 ? 线性二步法的一般格式 ? 适当选取参