数值分析7.4.pptVIP

梯形法的递推化 龙贝格算法 我们可以在步长逐步分半过程中将粗糙的积分值逐步加工为精度较高的积分值 ： 或者说将收敛缓慢的梯形值序列 加工成收敛迅速的积分值序列 ，这种加速方法称为龙贝格算法。 理查森外推加速法 * * 第七章 微积分的数值计算方法 7.4 Romberg算法 7.4 Romberg算法 综合前几节的内容,我们知道 梯形公式,Simpson公式,Cotes公式的代数精度分别为 1次,3次和5次 复化梯形、复化Simpson、复化Cotes公式的收敛阶分别为 2阶、4阶和6阶 无论从代数精度还是收敛速度,复化梯形公式都是较差的 有没有办法改善复化梯形公式呢? 实际计算中，由于要事先给出一个合适的步长往往很困难，所以我们往往采用变步长的计算方案,即在步长逐步分半的过程中，反复利用复化求积公式进行计算，直到所求得的积分值满足精度要求为止. 复化求积方法可提高求积精度，实际计算时若精度不够可将步长逐次分半. 设 表示复化梯形求得的积分值. 如果将求积区间再二分一次，则分点增至 2n+1 个. 每个子区间 经过二分只增加了一个分点 例2 解 先对整个区间 使用梯形公式. 从而可导出下列递推公式 （4.1） 计算积分值 对于函数 定义它在 的值 而 由梯形公式 （3.2） 例2 将区间二等分，求出中点的函数值 利用递推公式，有 进一步二分求积区间，并计算新分点上的函数值 再利用式，有 这样不断二分下去，计算结果见下表. 它表明用复化梯形公式计算积分 要达到7位有效数 字的精度需要二分区间10次，即要有分点1025个，计算量 很大. 梯形法计算简单但收敛慢，本节讨论如何提高收敛速 度以节省计算量. 根据复化梯形公式的余项表达式 假定 则有 （3.3） 由此可见，只有二分前后的两个积分值 与 相 当接近，就可以保证计算结果 的误差很小. 积分近似值 的误差大致等于 作为积分近似值，可能是更好的结果？！ 事实上， 是复化Simpson公式 启示：可以用低阶的公式组合一个高阶公式 由梯形公式出发，将区间 逐次二分可提高求积 公式的精度，上述加速过程还可继续下去，设 若记 当区间 分为 等分时， 梯形公式余项可展成级数形式，即有下面定理. 并且有 则有 定理4 设 则有 （*） 其中系数 与 无关. 定理表明 是 阶， 代替 ，有 （7。8） 若用4乘(7.8)式，减去(*)式再除3记之为 在(*)中，若用 则得 这里 以及后面将出现的 均为与 无关的系数， 这样构造的 与积分值 近似的阶为 . （7.9） 比较(7.9)与(4.4)可知，这样构造的序列 就是辛普森公式序列 又根据(7.9)，有 （