教案二Runge-Kutta方法.PDFVIP

教案二 Runge-Kutta方法 基本内容提要: 1. Runge-Kutta方法的基本思想, 2. 二级二阶Runge-Kutta方法, 3. 三级三阶Runge-Kutta方法, 4. 四级四阶经典Runge-Kutta方法. 教学目的和要求： 1. 掌握二级二阶Runge-Kutta方法, 2. 掌握经典Runge-Kutta方法. 教学重点： 1. Runge-Kutta方法的基本思想, 2. 二级二阶Runge-Kutta方法, 3. 四级四阶经典Runge-Kutta方法. 教学难点： 1. Runge-Kutta方法的基本思想, 2. 经典Runge-Kutta方法. 课程类型： 新知识理论课. 教学方法： 结合提问，以讲授为主. 教学过程： 问题引入 1 显式Euler方法是最简单的单步法, 它是一阶的, 可以看作Taylor展开后取前面的两项. 因 此, 得到高阶方法的一个直接思想是用Taylor展开, 如果能计算y(x)的高阶导数, 则可写出p阶 方法的公式 h2 hp y = y + hy0 + y00 + · · · + y(p), n+1 n n 2 n p! n (j) (j) ?f ?f 其中yn 可以看作y (xn)的近似值(j = 0, 1, · · · , p). 若将f(x, y), ?x , ?y , · · · , 分别记成f, fx, fy, · · · , 则对于二阶和三阶导数可表示为: 00 y = fx + fyf, 000 2 2 y = fxx + 2fxyf + fxfy + fyyf + fy f. 这个方法并不实用, 因为一般情况下, 求f(x, y)的导数相当麻烦. 从计算高阶导数的公式知道, 方法的截断误差提高一阶, 需要增加的计算量很大. 但是由此启发我们用区间上的若干个点 的f导数, 而不是高阶导数, 将它们作线性组合得到平均斜率, 将其与解的Taylor展开相比较, 使得前面若干项相吻合, 从而得到具有一定阶的方法. Runge-Kutta 方法 Runge-Kutta方法的基本思想 Runge-Kutta方法的一般形式为 XL 0 yn+1 = yn + h λiKi i=1 K1 = f(xn, yn), Xi?1 Ki = f(xn + cih, yn + cih aijKj), i = 2, 3, · · · , L, j=1 PL Pi?1 其中, ci ≤ 1, i=1 = 1, j=1 aij = 1. 它的局部截断误差是