偶阶求积公式的代数精度.ppt

第4章 数值积分与数值微分 4.1.2 代数精度的概念 4.1.3 插值型的求积公式 4.1.4 求积公式的收敛性与稳定性 4.2 牛顿-柯特斯公式 4.2.2 偶阶求积公式的代数精度 4.2.3 几种低阶求积公式的余项 4.3 复化求积公式 4.3.2 复化辛普森求积公式 4.4 龙贝格求积公式 4.4.2 龙贝格算法 4.4.3 理查森外推加速法 为研究辛普森公式(2.3)的余项 构造次数 不超过3的多项式 满足 （2.6） 其中 辛普森公式具有三次代数精度，对于这样构造出的三次式 应是准确的，即 对于多项式 ，其插值余项由第2章(5.11)得 由插值条件(2.6)，上式右端实际上等于按辛普森公式(2.3) 求得的积分值 ， 因此积分余项 故有 这时积分的核函数 在 上保号 类似的，对于柯特斯公式(2.4)，结果如下： （2.8） (非正)，再用积分中值定理有 （2.7） 复化求积的基本思想是把积分区间分成若干子区间(通 常是等分)，再在每个子区间上用低阶求积公式，目的是提 高精度. 4.3.1 复化梯形公式 将区间 划分为 等分， 在每个子区间 上 分点 采用梯形公式(1.1)，则得 （3.1） 记 称为复化梯形公式. （3.2） 由(2.5) ，其余项 由于 , 且 所以 使 于是复化梯形公式余项为 （3.3） 误差是 阶， 且当 时有 即复化梯形公式是收敛的. 将 改写为 此外， 的求积系数为正，由定理2知复化梯形公式是 稳定的. 只要 则当 时，上式右端括号内的两个 和式均收敛到积分 所以复化梯形公式(3.2)收敛. 将区间 分为 等分， 在每个子区间 上 若记 记 （3.5） 采用辛普森公式(2.3)， 则得 （3.4） 称为复化辛普森求积公式. 由(2.7)，其余项 于是当 时， （3.6） 误差阶为 ，显然是收敛的. 与复化梯形公式相似有 实际上，只要 则可得到收敛性， 即 此外，由于 中求积系数均为正数，故知复化辛普森 公式计算稳定. 例1 对于函数 ， 给出 的函数表 并估计误差. 解 将积分区间 划分为8等分， (见表4-2)， 计算积分 应用复化梯形法求得 试用复化梯形公式(3.2)及复化辛普森公式(3.5) 而如果将 分为4等分，应用复化辛普森法有 以上得到的两个结果 与 ，都需要提供9个点上的 同积分的准确值 比较， 接下来看误差估计 ，由于 所以有 函数值， 计算量基本相同，然而精度却差别很大. 只有两位有效数字. 复化梯形法的结果 于是 由(3.3)得复化梯形公式误差 对复化辛普森公式，由(3.6)得 4.4.1 梯形法的递推化 复化求积方法可提高求积精度，实际计算时若精度不 够可将步长逐次分半. 设将区间 分为 等分，共有 个分点， 如果将求积区间再二分一次，则分点增至 个， 我们将二分前后两个积分值联系起来加以考察. 用复化梯形公式求得该子区间上的积分值为 每个子区间 经过二分只增加了一个分点 这里 代表二分前的步长. 将每个子区间上的积分值相加得 例2 解 先对整个区间 使用梯形公式. 从而利用式(3.2)可导出下列递推公式 （4.1） 计算积分值 对于函数 定义它在 的值 而 由梯形公式 将区间二等分，求出中点的函数值 利用递推公式(4.1)，有 进一步二分求积区间，并计算新分点上的函数值 再利用式(4.1)，有 这样不断二分下去，计算结果见下表. 它表明用复化梯形公式计算积分 要达到7位有效数 字的精度需要二分区间10次，即要有分点1025个，计算量 很大. 梯形法计算简单但收敛慢，本节讨论如何提高收敛速 度以节省计算量. 根据复化