常微分方程初值问题数值解法1.ppt

第9章 常微分方程初值问题数值解法 9.2 简单的数值方法与基本概念 9.2.2 梯形方法 9.2.3 单步法的局部截断误差与阶 9.3 龙格—库塔方法 9.3.2 二阶显式R-K方法 9.3.3 三阶与四阶显式R-K方法 对r=2的R-K方法，由(3.4),(3.5)式可得如下计算公式 这里 c1, c2, λ2, μ21 均为待定常数，我们希望适当选取这些系数，使公式阶数 p 尽量高. 根据局部截断误差定义，推导出(3.6)的局部截断误差为 其中 这里yn=y(xn), yn+1=y(xn+1). 为得到Tn+1的阶p，要将上式各项在(xn, yn)处做泰勒展开，由于f(x, y )是二元函数，故要用二元泰勒展开，各项展开式为 将以上结果代入(3.7)，则有 要使公式(3.6)具有p=2阶，必须使 即 (3.9)的解是不唯一的. 可令c2=a≠0，则得 这样得到的公式称为二阶R-K方法. 则由此可以看出在改进的欧拉公式中相当于取(xn,yn), (xn+1,yn+1)两点处斜率的平均值，近似代替平均斜率，其精度比欧拉公式提高了. 如取a=1/2，则c1= c2=1/2, λ2=μ21=1. 这就是改进的欧拉公式(3.1). 称为中点公式(变形的欧拉公式)，相当于数值积分的中矩形公式.也可以表示为 如取a=1，则c