Romberg积分及初值问题.ppt

于是 由于 实际上 简单的数值方法与基本概念 1. 简单欧拉法（Euler） 2．后退的欧拉法 3．梯形法 4．改进Euler法 §2、初值问题的数值解法―单步法 * * 5.3 Romberg算法 综合前几节的内容,我们知道 梯形公式,Simpson公式,Cotes公式的代数精度分别为 1次,3次和5次 复合梯形、复合Simpson、复合Cotes公式的收敛阶分别为 2阶、4阶和6阶 无论从代数精度还是收敛速度,复合梯形公式都是较差的 有没有办法改善梯形公式呢? * 一、复合梯形公式的递推化 各节点为 复合梯形(Trapz)公式为 --------(1) --------(2) * --------(3) * 则由(1)(2)(3)式,有 * 因此(1)(2)(3)式可化为如下递推公式 (4)------- 上式称为递推的梯形公式 * 三种公式之间的关系 Romberg算法求解步骤 （3，-3）（-1,2）（0,1）（1，-1）（3，-4） * 二、外推加速公式 由复合梯形公式的余项公式 可得 由(3)式 * 复合Simpson公式 --------(5) --------(6) * 因此由复合Simpson公式的余项 可得 即 当然 令 自己 证明 --------(6) --------(7) * --------(8) 即 当然 同样由复合C