8.Runge-Kutta方法.PDFVIP

第八章常微分方程数值解法 8.2 Runge-Kutta 方法 8.2.1 Runge-Kutta方法的基本思想 8.2.2 几类显式Runge-Kutta方法 第八章常微分方程数值解法 8.2.1 Runge-Kutta方法的基本思想 显式Euler方法是最简单的单步法，它是一阶的，它可以看作Talylor展开后 ( ) y x 取前两项。因此，得到高阶方法的一个直接想法是用Talylor展开，如果能计算 的高阶导数，则可写出p阶方法的计算方法 ′ h2 ′′ hp (p ) y n +1 y n +hy n + y n +L+ y n ， 2 p ! j (j ) ( )( ) y y x ( ) 其中 是 的近似值， 若将f x ，y ，∂f ∂x ，∂f ∂y ，L， n n j 0，1，2，Lp 。 分别记成 则对于二阶和三阶导数可表示为 f ，f x ，f y ，L， ′′ y f x +f y f ， y ′′′ f xx +2f xy f +f x f y +f yy f 2 +f y2 f 。 第八章常微分方程数值解法 ( ) 这个方法并不实用，因为一般情况下，求f x ，y 的导数相当麻烦。从计算高 阶导数的公式知道，方法的截断误差提高一阶，需要增加的计算量很大。但是 f 由此启发我们用区间上若干个点的导数 ，而不是高阶导数，将它们作线性 组合得到平均斜率，将其与解的Taylor展开相比较，使前面若干项吻合，从而 得到具有一定阶的方法。这就是Runge-Kutta方法的基本思想，其一般形式为 L y y +h∑λK ， n+1 n i i i 1 (8.2.1) K f (x y ) ， ， 1 n n ⎛ i−1 ⎞ ⎜ ⎟ K f x +c h