82第二节 龙格—库塔法.pptVIP

数学学院 信息与计算科学系 第二节 龙格-库塔法 　　前面所讨论的欧拉法与改进的欧拉法都是一步法，即计算yn+1时，只用到前一步的值，但精度却不一样，欧拉法是1阶的，而改进的欧拉法却是2阶的. 龙格-库塔(Runge-Kutta)法(简称R-K方法)是一类高精度的一步法，这类方法与泰勒展开级数法有着密切的关系. 一、泰勒级数法 　　设有初值问题 由泰勒展开式 (1) 若令 (2) 则 即公式(2)为k 阶方法. 二、龙格-库塔方法(R-K方法) 　　分析欧拉法 每步计算f 的值一次, 其截断误差为O(h2). R-K方法不是通过求导数的方法构造近似公式, 而是通过计算不同点上的函数值, 并对这些函数值作线性组合, 构造近似公式, 再把近似公式与解的泰勒展开式进行比较, 使前面的若干项相同 , 从而使近似公式达到一定的阶数. 每步计算f 的值两次, 其截断误差为O(h3). 　　分析预估-校正法 其中R1, R2, a, b 为待定常数. (3) 　　下面对预估-校正法进行改进, 将该公式写成更一般的形式 选择这些常数的原则是在 yn=y(xn) 的前提下，使 y(xn+1)-yn+1的截断误差的阶尽量高. 为此, 对k2作二元泰勒展开式( 这里已有k1=hf ) 其中f,fx,fy,…都是在(xn, yn)处的函数值. 将k1, k2代入yn+1得 只要4个参数满足 与泰勒展开式(2)进行比较,要使得 (4) 若取 ,即得前面预估-校正公式. 　　满足(4)式的R1,R2, a, b 可以有各种不同的取法,但不管如何取法，都要计算两次f 的值(即计算f 在两个不同点的函数值)，截断误差都是O(h3). 满足条件(4)的一族公式(3)统称为二阶龙格-库塔公式.可以证明不增加计算函数值的次数，公式的截断误差的阶数不会再提高了. 　　如果每步计算3次f 的值，可将公式写成下列形式 (5) 　　与二阶龙格-库塔公式的讨论方法类似，要使 只需8个参数满足 方程组包含6个方程，8个未知量，其解不唯一. (6) 满足条件(6)的一族公式(5)统称为三阶龙格—库塔公式. 　　一个比较简单的三阶龙格—库塔公式是 (7)式中13个待定常数需满足下列11个方程的方程组 (7) 　　截断误差O(h5)的四阶龙格—库塔公式是常用的公式，每步都要计算4 次f 的值. 它的一般形式是 　　最常用的四阶龙格-库塔公式是标准四阶龙格-库塔公式 和 吉尔公式 　　标准四阶龙格-库塔公式手算时采用表8-2所示的表格计算 表8-2 例1 用标准四阶龙格-库塔方法求解初值问题 解 计算过程和结果如表8-3所示 0.1584376 0.0849171 0.0794465 0.0787495 0.0744037 1.183229 1.267746 1.262676 1.340728 0.2 0.3 0.3 0.4 0.1832292 0.1000000 0.0918182 0.0908637 0.0843239 1.000000 1.100000 1.091818 1.181727 0 0.1 0.1 0.2 表 8-3 1.483281 0.6 0.1416245 0.0745394 0.0710094 0.0708400 0.0673253 1.341667 1.416026 1.412676 1.482627 0.4 0.5 0.5 0.6 接上图 因此有 对该例, 用几种不同的一步法计算的结果如下表8-4 表 8-4 1 1.095445 1.183216 1.264911 1.341641 1.414214 1.483240 1.549193 1.612452 1 1.183229 1.341667 1.483281 1.612514 1 1.095909 1.184097 1.266201 1.343360 1.416402 1.485956 1.552514 1.616475 1 1.1 1.191818 1.277438 1.358213 1.435133 1.508966 1.580388 1.649738 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 准确值 四阶R-K法 (h=0.2) 预估-校正法 (h=0.1) 欧拉法 (h=0.1) 　　由表8-4可见，虽然四阶龙格-库塔方法每步要计算四次 f 的值，但以h=0.2为步长的计算结果就有5 位有效数字，而欧拉法与预估计-校正方法以 h=0.1为步长的计算