数值分析教教案15.docVIP

3.5龙格-库塔法 龙格-库塔方法的基本原理是：用方程 中函数在前一节点上取值的线性组合构造一个表示的近似值公式，从而避免了求时用的高阶导数。该方法有多种推导方法，这里用数值积分法推导。为此，先将微分方程略加变形得出：，对该式两边在相邻两节点和之间求积分，移项得出： (3-15) 采用不同的方法计算式(3-15)中定积分，便可得出数值积分的不同近似结果。如果用矩形求积公式计算，代入式(3-15)就可得出，这个结果和由Taylor公式得出的Euler公式结果完全一样。 3.5.1 二阶龙格-库塔公式 若用梯形求积公式计算式(3-15)中的积分，则有： 令，上式中里的用欧拉公式代换，则可得出，把和代入上式则有： 令：，则得出式(3-15)的一个近似结果： (3-16) 这就是二阶龙-格库塔公式，它的局部截断差为。式(3-16)右边的函数是和的线性组合，而和是把和的值代入函数得出的。这样计算时不再像用Taylor公式那样求的导数。 3.5.2 三阶龙格-库塔公式 若用抛物线(Simpson)求积公式计算式(3-15)中的积分，则有： 式中，而和都是未知的，可以用欧拉格式估算和。类似二阶龙格-库塔公式的推导，并令： ，， 把它们代入式(3-15)则得出三阶龙格-库塔公式，它的局部截断误差为。 （3-17） 3.5.3三阶龙格-库塔公式的MATLAB实