数值分析93龙格库塔方法.pptxVIP

第九章 常微分方程数值解 §2 龙格 - 库塔法 /* Runge-Kutta Method */;复习;一、再看Taylor方法; 注：应用泰勒展开式构造公式时，希望近似值的泰勒展开式与精确值的泰勒展式有尽可能多的项相同，即：两者关于h的展式中的系数尽可能多的项相同，这样局部截断误差与h的更高次同阶，从而获得高精度的近似公式。 ;建立高精度的单步递推格式。;首先希望能确定系数 ?1、?2、p，使得到的算法格式有2阶精度，即在 的前提假设下，使得 ;Step 3: 将 yi+1 与 y( xi+1 ) 在 xi 点的泰勒展开作比较;其中?i ( i = 1, …, m )，?i ( i = 2, …, m ) 和 ?ij ( i = 2, …, m; j = 1, …, i?1 ) 均为待定系数，确定这些系数的步骤与前面相似。 ;? 最常用为四级4阶经典龙格-库塔法 /* Classical Runge-Kutta Method */ ：;注： ? 龙格-库塔法的运算在于计算 Ki 的值，即计算 f 的值。Butcher 于1965年给出??计算量与可达到的最高精度阶数的关系：