第六讲常微分方程初值问题的数值解法2课题.pptVIP

* §6.3 Runge-Kutta法 考虑改进Euler法 通常写成如下的预报—校正公式 预报—校正公式是由梯形公式和Euler公式复合而成， 它具有二阶代数精度。 一、Runge-Kutta方法的基本思想 对于常微分方程的初值问题 的解 即 引入记号 则 改进的欧拉公式计算时进一步修改为预报—校正公式 推广上面的格式 其中 这就是龙格—库塔方法的思想。 如何确定其中的待定系数？ 通过Taylor展开式求上述式中的待定系数。 二、Runge-Kutta公式的推导 此时，有 利用二元函数的Taylor展开形式 再利用 得到 比较上面的两个式子，得 满足上式的解有无穷多，每一公式都是二阶精度。统称为二阶龙格—库塔方法。 若取 则正好是预报—校正公式。 若取 则得到 同样的推导方法可以构造出三阶、四阶的Runge-Kutta 公式。 一般常用的三阶Runge-Kutta公式为 局部截断误差为 ，三阶精度 实际中更常用的是下列四阶Runge-Kutta公式，又称为经典的龙格—库塔公式。 例1 用经典的龙格-库塔法计算上节例题，取步长 解： 从而 利用上面公式可计算出 欧拉法 预报-校正法 准确解 1 1.183216 1.341641 1.483240 1.612453 1.732051 1 1.184096 1.34360 1.485956 1.616976 1.737869 1 1.191818 1.358213 1.508966 1.649783 1.784770 0 0.2 0.4 0.6 0.8 1.0 1 1.18323 1.341667 1.483281 1.612513 1.732140 0 0.2 0.4 0.6 0.8 1.0 0 1 2 3 4 5 依次可计算出 ，其计算结果列于下表中。通过与上节例题结果相比可以发现，尽管步长增加一倍，但精度比预报-校正方法高。 例2 用Euler法（h=0.025），改进Euler法（h=0.05）和经典R-K方法（ h=0.1）计算初值问题 。 从计算结果看，在工作量大致相同的情况下，还是经典方法比其它两种方法的结果好得多．在 时，三种方法的误差分别是 。 计算结果如下表 解 设从 出发以 为步长，经过一步计算得 的近似值为 ，其截断误差为 ，即 三、变步长方法 问题：对于经典的龙格一库塔法，在计算过程中如何选择合适的步长才能满足精度？ 当 不大时， 可近似地看作常数，然后将步长折半，取 为步长，从 出发经两步计算求得 的近似值为 ，每一步的截断误差为 ，于是有 比较上述两式可得： 整理得： 这样，我们可以通过检查步长折半前后两次计算结果的偏差 来判定步长是否合适。 （1）对给定的精度 ，若 ，则将步长折半进行计算，直到 为止，这时取最终得到的 作为结果。 （2）如果 ，则反复将步长加倍，直到 为止，这时将步长折半一次，就得到所要的结果。 这种通过加倍或折半的手续处理步长的方法称为变步长方法，这种方法对于给定的精度能自动选取合适的步长，具有自适应能力，在实际运用中经常使用。 前面所讲的各方法在求 时，只用到前一步 的信息，即由 就可以计算出 ，这种方法称为单步法。如果在计算 时不仅要用到 ，而且还要用到 等等，这样的方法就称为多步法。由于多步法利用了更多的信息，因而能获得较高的精度。 一、Adams显示公式 对本章开始的微分方程，在 上两边积分得 如果被积函数 用插值点为 的线性插值函数代替，就得到改进的欧拉公式 §6.4 阿当姆斯（Adams）方法 问题： 通常插值多项式的次数越高，计算的值越准确，因此我们可以用更高次的插值多项式 逼近 ， 以得到较高