9常微分方程数值解法.pptVIP

Numerical Analysis 第九章 常微分方程初值问题数值解法 9.2 简单的数值方法与基本概念 9.2.2 梯形方法 9.2.3 单步法的局部截断误差与阶 9.3 龙格—库塔方法 9.3.2 二阶显式R-K方法 9.3.3 三阶与四阶显式R-K方法 当r=1, ?(xn, yn, h)=f(xn, yn)时，就是欧拉法，此时方法的阶为p=1. 当r=2时，改进欧拉法(3.1),(3.1)是其中一种，下面将证明其阶p=2. 要使公式(3.4),(3.5)具有更高的阶p，就要增加点数r. 下面我们只就r=2推导R-K方法. 并给出 r=3,4 时的常用公式，其推导方法与r=2时类似，只是计算较复杂. 对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)的解是不唯一的