常微分方程§4.3.ppt

* 1. 显式龙格-库塔法的一般形式 2. 二阶显式R-K方法 3. 三阶与四阶显式R-K方法 4. 变步长的龙格-库塔方法 §4.3 龙格-库塔(Runge-Kutta)方法 1.显式龙格-库塔法的一般形式 此时增量函数为 可知，若要使公式阶数提高，就必须使右端积分的数值求积 　由上节知欧拉法的局部截断误差为Tn?1 ?O(h2)，即方法为p=1阶，若用改进欧拉法，它可表为 它比欧拉法的?(xn , yn , h) ? f(xn , yn)，增加了计算一个右函数f的值，可望 p=2．若要使得到的公式阶数 p更大， ?就必须包含更多的f 值．实际上从方程(4.1)等价的积分形式 公式精度提高，它必然要增加求积节点，为此可将(4.16)的右端用求积公式表示为 一般来说，点数 r 越多，精度越高，上式右端相当于增量函数?(x , y , h)，为得到便于计算的显式方法，可类似于改进欧拉法(4.14), (4.15)将公式表示为 其中 这里ci , ?i , ?ij均为常数． (4.17)和(4.18)称为r 级显式龙格-库塔(Runge-Kutta)法，简称R-K方法． 　当r =1， ?(xn , yn , h) ? f(xn , yn)时，就是欧拉法，此时方法的阶为p =1．当r =2时，改进欧拉法(4.14), (4.15)就是其中的一种，下面将证明阶p=2．要