第七章常微分方程数值解 课件.ppt

7.4 龙格 - 库塔（ Runge-Kutta ）法 7.4.1 龙格 - 库塔 (Runge-Kutta) 法的基本思想 Euler 公式可改写成 ? ? ? ? ? ? ? ) , ( 1 1 1 i i i i y x f K hK y y 则 y i+1 的表达式 y ( x i+1 ) 与的 Taylor 展开式的前两项 完全相同 , 即局部截断误差为 。 改进的 Euler 公式又可改写成 ) ( 2 h O ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ) , ( ) , ( ) ( 2 1 1 2 1 2 1 1 hK y x f K y x f K K K h y y i i i i i i 上述两组公式在形式上有一个共同点 : 都是用 f(x,y) 在某些点上值的线性组合得出 y(x i+1 ) 的近似 值 y i+1 , 而且增加计算的次数 f ( x , y ) 的次数 , 可提高截 断误差的阶。如欧拉公式 : 每步计算一次 f ( x , y ) 的值 , 为一阶方法。改进欧拉公式需计算两次 f ( x , y ) 的值， 它是二阶方法。它的局部截断误差为 。 ) ( 3 h O 于是可考虑用函数 f ( x , y ) 在若干点