工程数理建模第九章.ppt

* 因为 * 类似前面的处理方法，可以得到四级方法：r =4 局部截断误差 最常用的一种四阶方法：经典显式Runge-Kutta公式 四阶Runge-Kutta方法 解 例4 用经典的四阶Runge-Kutta 方法求解下列初值问题 。 经典的四阶Runge-Kutta公式： 四阶Runge-Kutta方法 四阶Runge-Kutta方法 四阶Runge-Kutta方法 解 向前Euler法步长取 例5：分别利用三种方法 求解初值问题的数值解 改进Euler法步长取 经典4阶R—K法步长取 计算结果为： 常微分方程数值解法 Euler方法 改进的Euler法 经典的R—K方法 准确值 0 0 0 0 0 0.1 0.096312 0.095123 000.2 0.183348 0.181193 000.3 0.262001 0.259085 000.4 0.333079 0.329563 000.5 0.397312 0.393337 00常微分方程数值解法 注： 1. 三种方法的步长取法有何特点？ 2. 比较三种方法在节点处的误差？ 3. 如在节点0.5处，三种方法的误差分别为 常微分方程数值解法 k步线性多步法 线性多步法 所谓的线性多步法，指的是某一步解的公式不仅与前一步的值有关，而且与前面若干步解的值有关的方法。 利用前面多步的信息，则可以期望得到较高的精度。 构造多步法的主要途径是基于数值积分的方法和基于Taylor展开法.前者是将常微分方程两端积分后利用插值求积公式得到，后者是利用局部截断误差定义和Taylor展开得到。 Linear Mutistep Method k步线性多步法的一般形式 其中 为 的近似， ， 为常数， 不全为零.由于上式给出了 之间的线性关系，故称为线性k步法. ,称为显式k步法， 称为隐式k步法。 可根据局部截断误差以及阶确定。 线性多步法 将 在 处进行Taylor展开 在 处局部截断误差 线性多步法 代入局部截断误差的表达式 其中 线性多步法 则有 选择 满足 从而该多步法是p阶方法 线性多步法 基于数值积分方法 积分可得 将原初值问题中的方程在区间 上积分 被积函数用在三个节点 上的Lagrange插值代替 线性多步法 称为米尔尼（Milne）方法：显式四步四阶方法 被积函数用在三个节点 上的Lagrange插值代替 用 表示 的近似值，记 ，从而得到 线性多步法 基于数值积分方法 称为Simpson方法：隐式两步四阶方法 将原初值问题中的方程在区间 上积分 被积函数用在三个节点 上的Lagrange插值代替，计算后可得 线性多步法 基于Taylor展开方法 考虑形如 的K步法，称为阿当姆斯（Adams）方法 为显式方法， 为隐式方法 线性多步法 K=3时，Adams显式三步三阶方法 K=3时，Adams隐式三步四阶方法 K=4时，Adams显式四步四阶方法 线性多步法 利用前面所述的Taylor展开方法可得： 基于Taylor展开方法构造线性多步法比较灵活，可以构造任意多步法公式。 例5 构造初值问题 试确定参数 使方法的局部截断误差的阶尽可能高，并求局部截断误差。 的显式二步公式 线性多步法 解 由局部截断误差定义并利用Taylor公式得到 线性多步法 令 得到 局部截断误差为 线性多步法 从而所求二步法为 显式多步法计算简单，但是其精度及计算的稳定性没有隐式方法好.隐式多步