江西财经大学数学与决策科学系制作华长生8.pptVIP

华长生制作 * 江西财经大学 数学与决策科学系 制作：华长生 第四章 数值积分与数值微分 4.3 Romberg算法 4.3 Romberg算法 综合前几节的内容,我们知道 梯形公式,Simpson公式,Cotes公式的代数精度分别为 1次,3次和5次 复合梯形、复合Simpson、复合Cotes公式的收敛阶分别为 2阶、4阶和6阶 无论从代数精度还是收敛速度,复合梯形公式都是较差的 有没有办法改善梯形公式呢? 一、复合梯形公式的递推化 各节点为 复合梯形(Trapz)公式为 --------(1) --------(2) --------(3) 则由(1)(2)(3)式,有 因此(1)(2)(3)式可化为如下递推公式 (4)------- 上式称为递推的梯形公式 递推梯形公式加上一个控制精度,即 可成为自动选取步长的复合梯形公式 具体的方法请同学们完成 思考 二、外推加速公式 由复合梯形公式的余项公式 可得 由(3)式 复合Simpson公式 --------(5) --------(6) 因此由复合Simpson公式的余项 可得 即 当然 令 自己 证明 --------(6) --------(7) --------(8) 即 当然 同样由复合Cotes公式的余项 得 令 --------(9) 外推 加速 公式 以上整个过程称为Romberg算法 将 上 述 结 果 综 合 后 * * * * *