计算机数值方法第五章 常微分方程数值解法.ppt

Step 2: 将 K2 代入第1式，得到 Step 3: 将 yi+1 与 y( xi+1 ) 在 xi 点的泰勒展开作比较 要求 ，则必须有： 这里有 个未知数， 个方程。 3 2 所有满足上式的格式统称为2阶龙格 - 库塔格式。 注意到， 就是改进的欧拉法。 ? 最常用为四阶经典龙格-库塔法 ： 注： ? 龙格-库塔法的主要运算在于计算 Ki 的值，即计算 f 的值。Butcher 于1965年给出了计算量与可达到的最高精度阶数的关系： 7 5 3 可达到的最高精度 6 4 2 每步须算Ki 的个数 ? 由于龙格-库塔法的导出基于泰勒展开，故精度主要受解函数的光滑性影响。对于光滑性不太好的解，最好采用低阶算法而将步长h 取小。 例 3 用四阶龙格－库塔法求解初值问题 解 因为x0=0，y0=1，h=0.2 ， 用标准四阶龙格－库塔公式计算，并与解析解 的结果进行比较，如表5－3： 表 5―3 §4 线性多步法 (2) 基于数值积分的构造法 将 在 上积分，得到 只要近似地算出右边的积分