第4章1-5节数值积分.pptVIP

* 3 . 几种低阶求积公式的余项 按余项公式，梯形公式的余项 这里积分的核函数 在区间 上保号(非正)， 应用积分中值定理，在 内存在一点 使 ， * * §3 复化求积公式 复化求积的基本思想是把积分区间分成若干子区间(通 常是等分)，再在每个子区间上用低阶求积公式，目的是提 高精度. 1. 复化梯形公式 分点 在每个子区间[xk,xk+1] (k=0,…n) 上 采用梯形公式，则得 将区间 划分为 等分， * 记 称为复化梯形公式. * 其余项 由于 , 且 所以 使 于是复化梯形公式余项为 * 误差是 阶， 且当 时有 即复化梯形公式是收敛的. * 此外， 的求积系数为正，由定理2知复化梯形公式是 稳定的. 只要 则当 时，上式均收敛到积分 所以复化梯形公式收敛. 将Tn 改写为 *  对复化梯形公式，还有 如果f(x)在[a,b]上有2r+2阶连续导数，余项 *  定义 设 * 2. 复化辛普森求积公式 记 将区间 分为 n 等分， n=2m, xk=a+kh,k=0,…,2m 在每个子区间 [x2k-2,x2k] 上用Simpson公式 * 称为复化辛普森求积公式. 于是当 时， 与复化梯形公式相似有 误差阶为 4 ，显然是收敛的. * 实际上，只要 则可得到收敛性， 即 此外，由于 Sn 中求积系数均为正数，故知复化辛普森 公式计算稳定 * 例2 对于函数 ， 给出 的函数表 并估计误差. 解 (见表4-2)， 计算积分 应用复化梯形法求得T8=0.9456909 试用复化梯形公式(及复化辛普森公式 将积分区间[0,1]划分为8等分， * 而如果将[0,1] 分为4等分，应用复化辛普森法有 S4=0.9460832 同积分的准确值 I=0.9460831比较， 接下来看误差估计 ，由于 所以有 只有两位有效数字,复化辛普森法有6位有效数字. 复化梯形法的结果 以上得到的两个结果 T8 与S4 ，都需要提供9个点上的 函值， 计算量基本相同，然而精度却差别很大. * 第4章 数值积分 * §1 引言 1 . 数值求积的基本思想 依据微积分基本定理，对于积分 只要找到被积函数 的原函数 ，便有下列牛顿-莱 布尼茨(Newton-Leibniz)公式： 但对于下列情形： * （1）被积函数，诸如 等等，找不到用 初等函数表示的原函数； （2）当 是由测量或数值计算给出的一张数据表. 这时，牛顿-莱布尼茨公式也不能直接运用. 因此有必要研究积分的数值计算问题. 由积分中值定理知，在积分区间 内存在一点ξ， 成立 * 就是说，底为 而高为 的矩形面积恰等于所求 曲边梯形的面积 (图4-1). 图4-1 * 问题在于点ξ的具体位置一般是不知道的，因而难以 准确算出 的值. 将 称为区间 上的平均高度. 这样，只要对平均高度 提供一种算法，相应地便 获得一种数值求积方法. 用两端点“高度“ 与 的算术平均作为平均高度 的近似值，这样导出的求积公式 是梯形公式(几何意义参看图4-2). * 图4-2 用区间中点 的“高度” 近似地取代平均 高度 ，则又可导出所谓中矩形公式(简称矩形公式) * 一般地，可以在区间 上适当选取某些节点 ， 然后用 加权平均得到平均高度 的近似值，这样 式中 称为求积节点； 称为求积系数,亦称伴随节点 的权. 权 仅仅与节点 的选取有关， 构造出的求积公式具有下列形式： 的具体形式. 而不依赖于被积函数 k A * 这类数值积分方法通常称为机械求积，其特点是将积 分求值问题归结为函数值的计算，这就避开了牛顿-莱布尼 茨公式需要寻求原函数的困难. * 2. 代数精度的概念 定义1 如果某个求积公式对于