数值分析Runge-Kutta方法.pptVIP

第三节 Runge-Kutta方法 7.2 Runge?Kutta方法 7.2.1 构造高阶单步法的直接方法 由Taylor公式： 当h充分小时，略去Taylor公式余项，并以yi、yi+1分别代替y(xi)、y(xi+1)，得到差分方程： (7.2-1) 其局部截断误差为： 当h充分小时，略去Taylor公式余项，并以yi、yi+1分别代替y(xi)、y(xi+1)，得到差分方程： (7.2-1) 其局部截断误差为： 即(7.2-1)为p阶方式，上述方式称为Taylor方式。 注：利用Tayler公式构造，不实用，高阶导数f (i)不易计算。 7.2.2 Runge?Kutta方法 1. 基本思想 因为 = y(xi) + hf (?，y(?)) = y(xi) + hK? 其中K? = f (?，y(?))称为y(x)在[xi，xi+1]上的平均斜率。 若取 K1 = f (xi，y(xi)) ——Euler公式 取 K2 = f (xi+1，y(xi+1)) —— 向后Euler公式 一阶精度 取 —— 梯形公式