常微分方程出值问题得数值解法.PPTVIP

* * 其插值余项为 基于数值积分的构造方法（续） 在（*）式中取k=0，并选择等距节点 作为插值节点，作函数F（x）的三次插值多项式 将 代入（*）式，得 * * 基于数值积分的构造方法（续） 略去上式右端第三项，得 对上式积分部分作变量代换 ，并注意到 则 * * 线性四步Adams显式公式 基于数值积分的构造方法（续） 其局部截断误差为 * * 基于数值积分的构造方法（续） 利用第二积分中值定理，得 由于插值多项式 是在 上作出的，而积分区间是 ，故上述线性四步Adams显式公式称为Adams外插公式 。 * * 基于数值积分的构造方法（续） 若在（*）式中取k=0，但选择等距节点 作为插值节点，作函数F（x）的三次插值多项式，仿照上面的做法，可以得到线性三步的Adams隐式公式及其局部截断误差 上式称为Adams内插公式。 由于Adams内插公式是隐式公式，故用它计算时也需用迭代法。通常把Adams外插公式与内插公式联合起来交替使用，先由前者提供初始值，再由后者进行修正，即 * * 基于数值积分的构造方法（续） 可以看出，无论四步Adams显式公式还是三步Adams隐式公式，均为四阶公式。类似于以上讨论，还可以采用数值积分方法得到其它p阶的线性多步法公式。 可以证明，当 时，迭代式收敛．若上式中的第二式只迭代一次，则可以得到Adams预测校正公式 * * 三、基于Taylor 展开的构造方法 Taylor展开法更具一般性。下面构造线性两步方法 利用线性多步法局部截断误差的定义，有 将上式在xn处作Taylor展开，并按h的升幂整理排列，得到 * * 基于Taylor 展开的构造方法（续） * * 基于Taylor 展开的构造方法（续） 令 求解上述方程组，得出 ，所得到的公式的局部截断误差至少为 。 * * 基于Taylor 展开的构造方法（续） 可以只要求前面的几个方程成立，例如要求前面的四个方程成立时，所得公式的局部截断误差至少为 。由于此时方程个数少于未知量个数，故此种情形下方程组有无穷多组解。此时方程的解可以写成 若取 ，则 ，从而得到一个线性两步隐式公式 其局部截断误差为 上式称为Simpson公式。该公式为四阶方法，它也可以由数值积分方法得到。 * * 基于Taylor 展开的构造方法（续） 若取 ，则 ，从而得到一个线性两步显式公式 其局部截断误差为 需要指出的是，上式不能用数值积分方法得到 。 可以按照上述方法构造其它的计算格式，也可以用类似的方法构造其它的线性多步方法，如前述的Adams公式等。 * * 基于Taylor 展开的构造方法（续） Remark1：用Taylor展开法比用数值积分法推导线性多步法更加灵活，推导出的方法也比积分法的更多。用数值积分法可以得到的方法，用Taylor展开法都可以得到，但用Taylor展开法得到的方法，用数值积分法却不一定能够得到。 Remark2：应用线性多步法求解初值问题时，开头几点处的数值要用别的方法先计算出来。一般选用与多步法同阶或更高阶的单步法，如Runge-Kutta法、Taylor展开法等。 Remark3：对线性隐式多步方法，除开头几点的函数值需要单独计算外，还需要迭代求解或采用预测-校正方法求解。 * * 数值算例 求解常微分方程初值问题 在[0,1]上的解，取步长h＝0.1。 计算结果如下图所示： 该问题的理论解为： * * * * * * * * 下面以二级Runge—Kutta公式为例进行具体推导。 对 要求适当选取系数 ，使当 时，上式的局部截断误差为 ，即成为二阶方法。 Runge—Kutta方法（续） * * Runge—Kutta方法（续） 按照局部截断误差的定义，有 * * Runge—Kutta方法（续） * * Runge—Kutta方法（续） 要使上式等于 ，只需满足 在上式中有四个未知量，三个方程，故可以得到无穷多组解，也就是可以得到无穷多个二级二阶Runge-Kutta公式．由于是二级方法，从而有 ，此时上式可以改写为 在上述方程中，选取不同的c2，即可获得