第三部分 龙格-库塔方法.pptVIP

内江师范学院数学与信息科学学院 吴开腾 制作 第三部分 龙格－库塔方法 2、分析 一、设计思想 1、微分中值定理 如果 在 上连续，在 内可导，则在 内至少存在一点 ，有 于是有 其中 称作区间 上的平均斜率。 问题：计算近似值 的关键是如何选择算法确定平均斜率 ？ 差商的定义： 加权平均斜率。 取点 的斜率 作为区间 上的平均斜率，精度低。 3、欧拉格式 4、改进的欧拉格式 于是： 其中 取点 和 的斜率值 和 的算术平均值作为区间 上的平均斜率. 5、Runge-Kutta算法的基本思想 在区间 上取几个点的斜率，然后计算它们的平均斜率 ，以此可以构造出更高精度的格式，这就是Runge-Kutta 算法的基本思想。 二、二阶Runge-Kutta 有 令 和 处的斜率值分别为 和 。其中： 条件 问题 如何确定 ，使得所构造的格式具有二阶精度? 令 其中 为待定参数。 由改进的欧拉格式有 处的预报值 ： 又由原方程（1），通过计算可以得到 处的另一个斜率值： 其中 为待定参数。 于是，有二阶Runge-Kutta格式 定理： 当 时，Runge-Kuta格式（12）具有二阶精度。 特别地，当 时，有 改进欧拉格式 三、三阶Runge-Kutta 除 和 外，再增加一个点 ，于是有: 令 其中 为待定参数， 为同（12）式。 又由方程得到 处的另一个斜率值： 利用 的平均斜率来预报 处的值 特别地，当 时，有 定理： 适当选取参数 ，可以保证Runge-Kuta格式（14）具有三阶精度。这类格式统称为三阶龙格－库塔格式，下列库塔格式就是其中之一。 注释1：可以用Taylor展示证明格式（14）具有三阶精度，并且还可以用类似的方法得到四阶及其以上的更高阶精度的Runge-Kutta格式。 四、四阶Runge-Kutta 按照数学学习一般方法，继续上面的做法，经过较复杂的数学推理和计算，可以导出四阶龙格－库塔格式，这里有一个较有用的四阶格式。 注释3：四阶精度的计算框图。 例3 求解以下初值问题 解：分析 1、 确定步长h＝0.2，由四阶龙格－库塔格式，有 具体形式如下 注意比较例1的步长增加了1倍。