数值分析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