数值积分li演示文件修改版.pptVIP

第4章 数值积分与数值微分 4.1.3 插值型的求积公式 4.1.4 求积公式的收敛性与稳定性 4.2牛顿-柯特斯公式(Newton-cotes) 4.2.1柯特斯系数 前一节讨论的插值型的求积公式没有对节点作任 何限制.事实上,对求积节点作适当限制或选择,可 以简化问题的复杂性和提高公式的性能,如余项、 代数精度 4.2.2 偶阶求积公式的代数精度 4.2.3 公式的余项 作 业(第九周) 这时有 ，即辛普森公式对次数不超过三次的多项式 均能准确成立，而它对 通常是不准确的，因此， 辛普森公式实际上具有三次代数精度. 另一方面，直接求积得 定理3 当阶 为偶数时，牛顿-柯特斯公式(2.1)至少有 次代数精度. （2.1） 证明 只要验证，当 为偶数时，牛顿-柯特斯公式对  的余项为零,即代数精度至少为n+1次. 由于这里 引进变换 并注意到 有 所以按余项公式有 若 为偶数，则 为整数， 再令 进一步有 因为被积函数 为奇函数，所以 4.1 数值积分概论 4.2 牛顿-柯特斯公式 4.3 复合求积公式 4.4 龙贝格求积公式 4.5 自适应积分方法 4.6 高斯求积公式 4.7 多重积分 4.8 数值微分 4.1 数值积分概论 4.1.1引言 (数值积分的必要性) 依据微积分基本定理，对于积分 只要找到被积函数 的原函数 ，便有下列牛顿-莱 布尼兹(Newton-Leibniz)公式： 但对于下列情形： （1）被积函数，诸如 等，找不到用 初等函数表示的原函数，或者即使能求得原函数但原函数的 表达式非常复杂，计算困难； （3）当 是由测量或数值计算给出的一张数据表 时，牛顿-莱布尼兹公式也不能直接运用. 因此有必要研究积分的数值计算问题. （2）被积函数能求得原函数,但原函数的 表达式非常复杂,计算困难,比如 ； 就是说，底为 而高为 的矩形面积恰等于所求 曲边梯形的面积 (图4-1). 图4-1 4.1.2 求积公式与代数精度 由积分中值定理知，在积分区间 内存在一点ξ， 成立 问题在于点ξ的具体位置一般是不知道的，因而难以 准确算出 的值. 将 称为区间 上的平均高度. 这样，只要对平均高度 提供一种算法，相应地便 获得一种数值求积方法. 用两端点“高度“ 与 的算术平均作为平均高度 的近似值，这样导出的求积公式 （1.1） 是梯形公式(几何意义参看图4-2). 图4-2 用区间中点 的“高度” 近似地取代平均 高度 ，则又可导出所谓中矩形公式(简称矩形公式) （1.2） 一般地，可以在区间 上适当选取某些节点 ， 然后用 加权平均得到平均高度 的近似值， （1.3） 式中 称为求积节点； 称为求积系数，亦称伴随节点 的权. 权 仅仅与节点 的选取有关，而不依赖于被积函数 的具体形式. 这样构造出的求积公式具有下列形式： 这类数值积分方法通常称为机械求积，其特点是将积 分求值问题归结为函数值的计算,基本思想就是用被积函数f(x) 在一些离散点处的函数值f(x0),f(x1),…,f(xn)的线性组合来计算 积分的近似值.这就避开了牛顿-莱布尼兹公式需要寻求原函数 的困难. 构造或确定一个求积公式,要讨论得问题有: (1)确定求积系数Ak和求积节点xk; (2)求积公式的误差估计和收敛性. 定义1 如果某个求积公式对于次数不超过 的多项式 均能准确地成立，但对于 次多项式就不准确成立， 则称该求积公式具有 次代数精度. 梯形公式(1.1)和矩形公式(1.2)均具有一次代数精度. 数值求积是近似方法，为保证精度，自然希望求积公 式对尽可能多的函数准确成立. P52魏尔斯特拉斯(weierstrass)定理:多项式函数可以任意 逼近连续的函数. 欲使求积公式(1.3)具有 次代数精度，则只要令它 对 都准确成立，就得到 （1.4） （1.3） 例1 给定形如