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