数值分析72龙格—库塔方法剖析课件.pptVIP

7.3龙格—库塔方法对许多实际问题来说，欧拉公式与改进欧拉公式精度还不能满足要求，为此从另一个角度来分析这两个公式的特点，从而探索一条构造高精度方法的途径.

7.3.1显式龙格—库塔法的一般形式上节给出了显式单步法的表达式(2.11),其局部截断误差为O(hp+1)，对欧拉法T=O(h)，即方法为2n+1p=1阶，若用改进欧拉法(2.9)：此时增量函数为

它比欧拉法的?(x,y,h)=f(x,y),增加了计算一个函nnnn数f的值，期望p=2.若要使得到的公式阶数p更大,?就必须包含更多的f值.实际上从方程(1.1)等价的积分形式(2.4)，即若要使公式阶数提高，就必须使右端积分的数值求积公式精度提高，它必然要增加求积节点，为此可将(3.3)的右端用求积公式表示为

一般说来，点数r越多，精度越高，上式右端相当于增量函数?(x,y,h)，为得到便于计算的显式方法，nn可类似于改进欧拉法(3.1),(3.2)，将公式表示为其中这里c,λ,μ均为常数.(3.4)和(3.5)称为r级显式龙格-iiij库塔(Runge-Kutta)法,简称R-K方法.

注：当r=1,?(x,y,h)=f(x,y)时，就是欧拉法，此nnnn时方法的阶为p=1.当r=2时，改进欧拉法(3.1)是其中一种，下面将证明其阶p=2.要使公式(3.4),(3.5)具有更高的阶p，就要增加点数r.下面我们只就r=2推导R-K方法.并给出r=3,4时的常用公式，其推导方法与r=2时类似，只是计算较复杂.

补充：二元函数的泰勒公式一元函数的泰勒公式:推广多元函数泰勒公式

记号：(设下面涉及的偏导数连续):??表示?一般地,表示

定理：的某一邻域内有直为此邻域内任到n+1阶连续偏导数,一点,则有①②其中①称为f在点(x,y)的n阶泰勒公式,②称为其拉格00朗日型余项.

7.3.2二阶显式R-K方法对r=2的R-K方法，由(3.4),(3.5)式可得如下计算公式这里c,c,λ,μ均为待定常数，我们希望适当选取12221这些系数，使公式阶数p尽量高.根据局部截断误差定义，推导出(3.6)的局部截断误差为

这里y=y(x),y=y(x).为得到T的阶p，要将上nnn+1n+1n+1式各项在(x,y)处做泰勒展开，用二元泰勒展开，各nn项展开式为其中

将以上结果代入(3.7)，则有

要使公式(3.6)具有p=2阶，必须使即(3.9)的解是不唯一的.可令c=a≠0，则得2这样得到的公式称为二阶R-K方法.

如取a=1/2，则c=c=1/2,λ=μ=1.这就是改进12221的欧拉公式(3.1).

如取a=1，则c=0,c=1,λ=μ=1/2.得计算公式12221称为中点公式(变形的欧拉公式)，相当于数值积分的中矩形公式.也可以表示为

对r=2的R-K公式(3.6)能否使局部误差提高到O(h4)?为此需把K多展开一项，从(3.8)的看到展2开式中的项是不能通过选择参数消掉的.故r＝2的显式R-K方法的阶只能是p=2，而不能得到三阶公式.

7.3.3三阶与四阶显式R-K方法要得到三阶显式R-K方法，必须r=3.此时计算(3.4),(3.5)的公式表示为其中c,c,c及λ,μ,λ,μ,μ均为待定常数，公式12322133132(3.11)的局部截断误差为

只要K,K将按二元泰勒展开，使T＝O(h)，可得412n+1待定参数满足方程

这是8个未知数6个方程的方程组，解不是唯一的.可以得到很多公式.满足条件(3.12)的公式(3.11)统称为三阶R-K公式.下面只给出其中一个常见的公式.此公式称为三阶库塔方法.

继续上述过程，经过较复杂的数学演算，可以导出各种四阶R-K公式，下列经典公式是其中常用的一个：四阶R-K方法的每一步需要计算四次函数值f，可以证明其局部截断误差为O(h).5

然而值得指出的是，龙格-库塔方法的推导基于泰勒展开方法，因而它要求所求的解具有较好的光滑性质.反之，如果解的光滑性差，那么，使用龙格-库塔方法求得的数值解，其精度可能反而不如改进的欧拉方法.实际计算时，我们应当针对问题的具体特点选择合适的算法.

*7.3.4变步长的龙格-库塔方法单从每一步看，步长越小，截断误差就越小，但随着步长的缩小，在一定求解范围内所要完成的步数就增加了.步数的增加不但引起计算量的增大，而且可能导致舍入误差的严重积累.因此同积分的数值计算一样，微分方程的数