现代科学与工程计算-孟大志-第三篇、第10章.pptVIP

代入（10.4.3）式得 与二阶Taylor展开式 比较系数可知，要使计算格式（10.4.3）的截断误差为 , 只要 即可。满足这一条件的一簇格式（10.4.3）统称 为二阶R-K格式。若以 为自由参数，分别取 ，可依次得到如下三个著名公式： (1) 中点公式： (2) Heun公式： (3) 改进Euler公式： 10.4.3 四阶R-K公式 构造二级二阶R-K方法的思路和方法，可以推广到 。 构造出一簇三级三阶，四级四阶等R-K方法。最常用的是四级 四阶的经典R-K方法： 方法的阶越高，方法就越精确。Euler法和向后Euler法各 用了 两端的一个导数，它们都是一阶方法。而梯形公式 和改进Euler法都用了两端的两个 导数，它们的精度的阶提高 了，都是二阶方法。由于R-K方法的构造是基于Taylor展开式， 因此它要求所求的解具有较好的光滑性。如果光滑性较差，则 有可能使用R-K方法求得的数值解的精度反而不如改进Euler法。 这就要求在实际计算时，针对具体问题的特点选择适当的算法。 10.4.4 R-K法的稳定性 对 使用R-K法，可得 因此 误差方程是 于是R-K法的绝对稳定域为 它在复平面上的轨迹如图10.4.1所示。 例2 用四阶R-K方法求初值问题 , 取步长 ，求 。 解：取 ，计算四个组成部分 将其代入 得 即使用四阶R-K方法在 处的估计值为 。 10.5* 线性多步法（阅读材料） 通常可以采用线性多步法，它的构造思想是利用前面 多步信息进行计算，以获得较高精度的数值公式。 多步法中最常用的是线性多步法。 10.5.1一般迭代格式 设线性多步法的一般迭代格式为： 其中 均为常数， 。若 为显式， 为隐式。 构造线性多步法的主要途径有基于数值积分法和基于Tayl or展开方法。这里介绍基于Taylor展开式的构造方法，即将线 性多步公式在 处Taylor展开，与 在 处Taylor展开式相 比较，要求它们前面有尽可能多的项重合，由此确定参数 。 下面介绍几种常见的线性多步法。 10.5.2 Adams 格式 取 ，并令 ，由方程组 * * * 第十章 常微分方程初值问题的数值方法 10.1 问题驱动：蝴蝶效应 洛伦兹吸引子(Lorenz attractor)是由MIT大学的气象学家E dward Lorenz在1963年给出的，他给出第一个混沌现象——蝴 蝶效应。 图10.1.1蝴蝶效应示意图 洛伦兹方程是大气流体动力学模型的一个简 化的常微分方程组： 该方程组来源于模拟大气对流，该模型除了在天气预报中有显 著的应用之外，还可以用于研究空气污染和全球侯变化。洛伦 兹借助于这个模型，将大气流体运动的强度x与水平和垂直方 向的温度变化y和z联系了起来。参数 称为普兰特数， 是规范 化的瑞利数， 和几何形状相关。洛伦兹方程是非线性方程组， 无法求出解析解，必须使用数值方法求解上述微分方程组。洛 伦兹用数值解绘制结果图10.1.1，并发现了混沌现象。 10.2 常微分方程的数值解 微分方程数值解一般可分为：常微分方程数值解和偏微分 方程数值解。自然界与工程技术中的许多现象，其数学表达式 可归结为常微分方程（组）的定解问题。一些偏微分方程问题 也可以转化为常微分方程问题来（近似）求解。Newton最早采 用数学方法研究二体问题，其中需要求解的运动方程就是常微 分方程。许多著名的数学家，如Bernoulli（家族），Euler、G auss、Lagrange和Laplace等，都遵循历史传统，研究重要的 力学问题的数学模型，在这些问题中，许多是常微分方程的求 解。作为科学史上的一段佳话，海王星的发现就是通过对常微 分方程的近似计算得到的。本章主要介绍常微分 方程数值解的若干方法。 10.2.1 初边值问题 由于积分求出解析解的局限性，绝大多数的常微分方程的 初值或边值问题都很难给出解析解。因此，研究常微分方程的 数值解法具有重要意义。 常微分方程的数值解法主要分为两大部分： (1) 初值问题的数值解法； (2) 边值问题的数值解法。 这里只介绍初值问题，并且首先研究一阶常微分方程初值 问题的数值解法。其模型可以如下表示： （10.2.1） 1. 迭代格式的构造 (1) 构造思想：将连续的微分方程及初值条件离散为线性 方程组加以求解。由于离散化的出发点不同，产生出各种不同 的数值方法。基本方法有：有限差分法（数值微分）、有限体 积法（数值积分