数值分析与9-3 龙格-库塔方法 .pptVIP

* 第九章 常微分方程数值解 §2 龙格 - 库塔法 /* Runge-Kutta Method */ 一、再看Taylor方法 二、Runge-Kutta 方法 复习 ? 一元函数的Taylor展开 ? 二元函数的Taylor展开 一、再看Taylor方法 一般可取公式为如下形式 注：应用泰勒展开式构造公式时，希望近似值的泰勒展开式与精确值的泰勒展式有尽可能多的项相同，即：两者关于h的展式中的系数尽可能多的项相同，这样局部截断误差与h的更高次同阶，从而获得高精度的近似公式。 建立高精度的单步递推格式。 单步递推法的基本思想是从 ( xi , yi ) 点出发，以某一斜率沿直线达到 ( xi+1 , yi+1 ) 点。欧拉法及其各种变形所能达到的最高精度为2阶。 ? 考察改进的欧拉法，可以将其改写为： 斜率 一定取K1 K2 的平均值吗？ 步长一定是一个h 吗？ 二、Runge-Kutta 方法 首先希望能确定系数 ?1、?2、p，使得到的算法格式有2阶精度，即在 的前提假设下，使得 Step 1: 将 K2 在 ( xi , yi ) 点作 Taylor 展开(用二元函数泰勒展开) 将改进欧拉法推广为： ) , ( ) , ( ] [ 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 Step 2: 将 K2 代入第1式，得到 Step 3: 将 yi+1 与 y( xi+1 ) 在 xi 点的泰勒展开作比较 要求 ，则必须有： 这里有 个未知数， 个方程。 3 2 存在无穷多个解。所有满足上式的格式统称为2阶龙格 - 库塔格式。 注意到， 就是改进的欧拉法。 问题: 为获得更高的精度，应该如何进一步推广？