研究生数值分析(21,22)数值积分.pptVIP

* * 数值分析 第6章 数值积分 主讲老师: 雷鸣 在高等数学中，计算定积分 根据微积分学基本定理，若被积函数f(x)在区间[a，b]上连续，只要找到f(x)的一个原函数F(x)便可利用牛顿—莱布尼兹公式 第6章 数值积分 §1 求积公式及其代数精度 求得积分值。 但是在工程技术和科学研究当中，往往遇到如下困难，而不能使用牛顿—莱布尼兹公式。 为 1、 找不到用初等函数表示的原函数 例如： 等等 原函数的计算复杂性大大超过被积函数。 2、 虽然找到了原函数，但因表达式过于复杂而不便于计算 例如： 3、f(x) 是由测量和计算得到的函数列表，即给出的是 f(x) 的一张数据表。 可惜的是ξ值不易找到，因而难以求出 f(ξ)的准确值。 由于这些困难，我们必须研究积分的数值计算问题。回顾积分中值定理 但若能对 f(ξ)提供一种近似算法，也可以得到一种数值积分公式。 若取 ，则得到 如取 ，则得到 如取 ，则得到 以上三个公式分别称为左矩形公式，中矩形公式和右矩形公式。 由定积分的定义 可以得到定积分的一个近似计算公式 为求积结点， 它们均与f(x)的具体形式无关。 进一步，我们设想更一般的求积公式为 ① 称 为求积系数， 这类数值积分的方法通常称为机械求积法，主要有插值型和外推型两种。它们均是应用被积函数 f(x)在一些节点上的函数值得线性组合得出积分的近似值。 于是求积分值的问题就转化为计算被积函数在节点处的函数值的问题。对于求积公式①，关键在于确定求积系数 。 设给定一组节点 且已知函数 f(x)在这些节点上的值为 ，则可作f(x)的n次插值多项式 （1）构造数值积分公式的基本方法 其中 为节点 的基本插值多项式。 用 近似代替被积函数 f(x) ，则得 若记 我们称由系数式确定 的数值积分公式① 立即可得数值积分公式①。 为插值型求积公式。 记插值型求积公式 的截断误差为R[ f ] ，则有 （2）插值型求积公式的截断误差与代数精度 ② 由余项公式②可以看出，如果f(x)是一个n次多项式，则 即 为了使求积公式能对更多的积分具有较好的实际计算意义，就应要求它尽可能多的被积函数f(x)都能准确地成立。为此我们引入代数精度的概念。 求积公式是准确成立的。 其中 定义 如果一个数值积分公式 对于次数≤m 的代数多项式均能准确地成立，但至少对一个 m+1 次多项式不能准确地成立，则称该数值积分公式具有 m 次代数精度。 关于代数精度有如下结论： 定理1 含有n+1个节点的插值型数值积分公式 其中 至少具有 n 次代数精度。 证 当 f(x)为任何次数不高于n的多项式时， 根据插值型求积公式的截断误差不断式 此时，R n = 0，因而等式 成立，其中 根据代数精度的定义，可知定理的结论成立。 证毕 定理2 数值积分公式 具有m次代数精度的充分必要条件为该公式对 精确成立，而对 不精确成立。 证 记 必要性 设 对任一 m 次多项式精确成立， 但对某一m+1次多项式 不准确成立，其中 是一个 m 次多项式。 由于 是特殊的次数 小于等于 m 的多项式，所以求积公式是准确成立的。另外由于 及 得 充分性 设 。 任一 m 次多项式，可表示为 于是 即求积公式对任一 m 次多项式是准确成立的，再注意到求积公式对 是不精确成立的，充分性证毕。