第九章 常微分方程数值解.ppt

通过计算结果的比较可以看出，改进的Euler方法 的计算精度比Euler方法要高。 §3 龙格 - 库塔法 /* Runge-Kutta Method */ 建立高精度的单步递推格式。 单步递推法的基本思想是从 ( xi , yi ) 点出发， 欧拉法及其各种变形所能达到的最高精度为2阶。 以某一斜率沿直线达到 点。 ? 考察改进的欧拉法，可以将其改写为： 斜率 一定取K1 K2 的平均值吗？ 步长一定是一个h 吗？ 将改进欧拉法推广为： ) , ( ) , ( ] [ 1 2 1 2 2 1 1 1 phK y ph x f K y x f K K K h y y i i i i i i + + = = + + = + l l 首先希望能确定系数 ?1、?2、p，使得到的算法格式有2阶精度，即在 的前提假设下，使得 Step 1: 将 K2 在 ( xi , yi ) 点作 Taylor 展开 Step 2: 将 K2 代入第1式，得到 Step 3: 将 yi+1 与 y( xi+1 ) 在 xi 点的泰勒展开作比较 要求 ，则必须有： 这里有 个未知数， 个方程。 3 2 存在无穷多个解。所有满足上式的格式统称为2阶龙格 - 库塔格式。 注意到， 就是改进的欧拉法。 为获得更高的精度，应该如何进一步推广？ ) ... , ( ... ... ) , ( ) , ( ) , ( ] ... [ 1 1 2 2 1 1 2 32 1 31 3 3 1 21 2 2 1 2 2 1 1 1 - - + + + + + + = + + + = + + = = + + + + = m m m m m m i m i i i i i i m m i i hK hK hK y h x f K hK hK y h x f K hK y h x f K y x f K K K K h y y b b b a b b a b a l l l 其中?i ( i = 1, …, m )，?i ( i = 2, …, m ) 和 ?ij ( i = 2, …, m; j = 1, …, i?1 ) 均为待定系数，确定这些系数的步骤与前面相似。 考虑一阶常微分方程初值问题 将区域[a,b]进行分划： 若 则 n级显式Runge-Kutta方法 二阶Runge-Kutta方法 取n=2 记 由此得 另一方面 为使局部截断误差为 ,应取 改进的Euler方法 取 中点方法 取 二阶Heun方法 取 二级Runge-Kutta方法不超过二阶 记 则 因此局部截断误差只能达到 三级Runge-Kutta方法 取n=3 * * （Numerical Methods for Ordinary Differential Equations ） 问题驱动：蝴蝶效应 洛伦兹吸引子(Lorenz attractor)是由MIT大学的气象学家E dward Lorenz在1963年给出的，他给出第一个混沌现象——蝴 蝶效应。 图10.1.1蝴蝶效应示意图 洛伦兹方程是大气流体动力学模型的一个简化的常微分方程组： 该方程组来源于模拟大气对流，该模型除了在天气预报中有显 著的应用之外，还可以用于研究空气污染和全球侯变化。洛伦 兹借助于这个模型，将大气流体运动的强度x与水平和垂直方 向的温度变化y和z联系了起来。参数 称为普兰特数， 是规范 化的瑞利数， 和几何形状相关。洛伦兹方程是非线性方程组， 无法求出解析解，必须使用数值方法求解上述微分方程组。洛 伦兹用数值解绘制结果图10.1.1，并发现了混沌现象。 §1 引 言 微分方程数值解一般可分为：常微分方程数值解和偏微分 方程数值解。自然界与工程技术中的许多现象，其数学表达式 可归结为常微分方程（组）的定解问题。一些偏微分方程问题 也可以转化为常微分方程问题来（近似）求解。Newton最早采 用数学方法研究二体问题，其中需要求解的运动方程就是常微 分方程。许多著名的数学家，如 Bernoulli（家族），Euler、 Gauss、Lagrange和Laplace等，都遵循历史传统，研究重要 的力学问题的数学模型，在这些问题中，许多是常微分方程的 求解。作为科学史上的一段佳话