数值分析第9章-2.pptVIP

* * 结束 §9.3 龙格-库塔(Runge-Kutta)法 欧拉方法是显式的一步法,使用方便,但精度较低.本节将构造出高精度的显式一步法：龙格-库塔法,简称R-K法. 9.3.1 二阶R-K法 欧拉法的公式为： yi+1=yi+h f (xi ,yi) i=0,1,2,…,n-1 决定其精度的是函数f (xi ,yi).如能改进这个函数,就可能提高公式的精度.为此把公式改写成： yi+1=yi+h ? (xi ,yi,h) i=0,1,2,…,n-1 (9.10) 选择函数? (xi ,yi,h),一种方法是用若干个点的函数值的线性组合代替? (xi ,yi,h),如： * 结束 其中cj,aj,bjl是待定参数, aj和bjl满足 以上方法称为p级R-K法,选择cj,aj和bjl,可能使以上方法为p阶方法.显然欧拉法就是一阶R-K法. * 结束 二级R-K法的形式是: 此时 由二元函数的泰勒展开： 其中所有的偏导数都是它们在点(xi,y(xi))的值,下同 * 结束 又由于： 所以 代入(9.11) 代入(9.10) 而Taylor展开式 * 结束 二式相减,得局部截断误差 只要c1,c2,a2满足以上方程,就得到一个二阶的R-K法. 这是一个不定方程,有无穷多解.比如： * 结束 (1)取c1=c2=1/2,a2=1得 这实际上是(9.9)公式,即梯形公式的预估-校正公式只迭代一次的形式,通常称为改进的欧拉法. (2)取c1=0,c2=1,a2=1/2 得: 这公式又称中点公式.我们还可以构造其他的二阶R-K法. * 结束 9.3.2 四阶 R-K 法 用类似的方法可以确定三级和四级R-K法的参数,构造出三阶和四阶的R-K法.但最常用的是四阶R-K法,四阶R-K法也不只一个,下面给出的是最常用的四阶经典的R-K公式： * 结束 例3 用经典的四阶R-K法计算例2题目,取步长为0.2,且与准确值比较. 计算结果列入表9-3： 可见即使用h=0.2计算,也比一阶和二阶方法精度好得多 * 结束 §9.4 线性多步法 单步法只利用前一步的结果,只要给出初值,就能开始计算,但也因为它只利用前一步的值,为了提高精度就要计算一些非结点处的函数值,增加了计算量.R-K法就是通过这一途径提高精度的.下面介绍的线性多步法,在求yi+1时,不仅用到yi的值, 还用到前若干步的yi-1,…,yi-k的值,这些值都是已知的,因此可在计算量增加不多的情况下提高精度. 9.4.1 用待定系数法构造线性多步法 线性多步法的一般形式是： 或写为： * 结束 其中αj,βj(j=0,1,…,k)都是实常数,且αk≠0, ｜α0｜+｜β0｜≠0, fi+j= f (xi+j,yi+j), j=0,1,…,k, 由(9.15)可看出要计算yi+k,要利用它前面的k个值yi,yi+1,…,yi+k-1,又因为(9.15)关于yi+j和 f i+j都是线性组合,所以这一类方法都称为线性k步法.欧拉法,隐式欧拉法和梯形法都是线性一步法,欧拉中点公式是线性二步法. 把y(x+jh)和y’(x+jh)作Taylor展开： * 结束 其中 * 结束 若选择αj,βj,使C0=C1=…=Cp=0,Cp+1≠0,则 将x=xi代入上式,设αk=1,并注意到y?(xi+jh)=f(xi+jh,y(xi+jh)),可推出 即 * 结束 设yi=y(xi),yi+1=y(xi+1),…,yi+k-1=y(xi+k-1),记左端为 yi+k,并舍去最后两项,(9.17)变为： 就是一种p阶的线性k步方法,Cp+1hp+1y(p+1)(xi)称为局部截断误差的主项.当βk=0时,是显式方法,当βk≠0时是隐式方法. 下面构造几个实用的线性多步法公式 例4 形如： yi+4= -α0yi+h(β1fi+1+β2fi+2+β3fi+3)的线性4步法公式,试确定α0,β1,β2,β3并求其局部截断误差主项. * 结束 解: 由(9.15)知α1=α2=α3=0,α4=1,β0=β4=0,因为有4个待定系数,由(9.16)写出前4个方程： 解之,得α0=-1,β1=β3=8/3,β2=-4/3 故所求的公式为 * 结束 将α0,β1,β2,β3代入C4 再求C5 公式(9.19)称为米尔尼(Milne)公式,它的局部截断误差为： 它也可以写成 * 结束 例5 试确定下列公式的系数和局部截断误差 yi+1=a yi+b yi-2+h (c fi+1+d fi+e fi-1) 解 按(9.15)和(9.16)可知k=3,α0= - b,α1=0,α2= - a,α3=1, β0=0,β1=e,β2=d,β3=c,有五个未知参数,写出前五个方程： * 结束 解之得：a