计算方法——常微分方程初值问题的数值解法.pptVIP

* 计 算 方 法 华中科技大学数学与统计学院 第六章 常微分方程初值 问题的数值解法 计算方法课程组 §6.1 基本离散方法 §6.2 Runge-Kutta方法 §6.3 线性多步法 §6.4 收敛性与稳定性 §6 常微分方程数值解法 基本思想: 根据微分中值定理有： 这里 称为在区间 上函数 的平均斜率。 §6.2 龙格 - 库塔方法 建立高精度的单步递推格式。 问题：计算近似值 的关键是如何选择算法确定平均斜率 ？ 龙格 - 库塔方法的基本思想 考察Euler格式： 它简单的取点 的斜率值 作为平均斜率 , 精度较低。 考察改进的Euler格式： 它可以理解为用 与 两个点的斜率值 和 取算数平均斜率 。 龙格-库塔方法 这个处理过程启示我们，如果设法在 内多预报几个点的斜率值，然后将它们加权平均作为平均斜率 ，则有可能构造出更高精度的计算格式。 二阶方法 推广改进的Euler格式，对于区间 内任意给定一点 ： 用 和 两个点的斜率值 和 加权平均得到 ： 即令 （ 为待定参数） 取 问题是如何取 先用Euler格式预报 ： 将 和 值代入 ，得到 二阶龙格-库塔方法 基本形式 二阶龙格-库塔方法 接下来问题在于如何选取参数值，使其具有较高精度。 考查近似关系式： 不管 取值如何，它都具有一阶精度。可以适当的选取参数 ，使其具有二阶精度。 为简化处理，令 ，上式可以化简为： 令它对于 准确成立，可得： 二阶龙格-库塔方法 满足 这一条件的二阶基本形式 统称为二阶的龙格-库塔格式。 二阶龙格-库塔方法 特别地，当 时，有 改进欧拉格式 三阶库塔格式 用三个点 ， ， 的斜率的加权平均得到平均斜率 。其形式： 取 用Euler格式预报得： ，可得 （3） 和 已知，用线性插值可得 三阶库塔格式 问题在于如何选取 值。 同二阶方法，得近似关系式： 它对 成立，令它对于 准确成立，可得 三阶库塔格式 综上所述，这样设计出的差分格式具有形式： 这种三阶格式称作库塔格式。 四阶经典龙格-库塔格式 继续前面的推导过程，设法在区间 内多预报几个点的斜率值，然后将它们加权平均作为平均斜率，即可设计出更高精度的单步格式。这类格式统称为龙格-库塔格式。 实际计算中常用的为四阶经典格式： 按照数学学习一般方法，继续上面的做法，经过较复杂的数学推理和计算，可以导出常用的四阶龙格－库塔格式。 四阶精度的计算框图： 四阶经典龙格-库塔格式 算例： 分别用Euler公式，改进的Euler公式，经典4阶R-K 公式计算一阶常微分方程初值问题。 并与准确解 比较。 解：Euler公式, 改进的Euler公式取步长h=0.1, 经典4阶R-K公式取步长h=0.2。 4阶R-K公式： X Euler公式 改进Euler公式 4阶R-K公式 准确值 0.0 1.0000 1.0000 1.0000 1.0000 0.1 1.1000 1.0959 1.0954 0.2 1.1918 1.1841 1.1832 1.1832 0.3 1.2774 1.2662 1.2649 0.4 1.3582 1.3434 1.3417 1.3416 0.5 1.4351 1.4164 1.4142 0.6 1.5090 1.4860 1.4833 1.4832 0.7 1.5803 1.5525 1.5492 0.8 1.6498 1.6153 1.6125 1.6125 0.9 1.7178 1.6782 1.6733 1.0 1.7848 1.7379 1.7321 1.7321 计算结果见下表： 分析：单从每一步看，步长越小，截断误差越小；但随着步长的缩小，在一定求解范围内所要完成的步数就增加了，步数的增加不但引起计算量的增大，而且可能导致舍入误差的严重积累。因此，同积分的数值计算一样，微分方程的数值解法也需要选择步长。在选择步长时，需要考虑的问题： 如何衡量和检验计算结果的精度； 如何依据所判定的精度来处理步长。 龙格 - 库塔方法 小结