王仁宏老师的《数值逼近》4.doc.docVIP

数值积分 教学目的及要求： 掌握Newton-Cotes 公式、Romberg 方法、Euler-Maclaurin公式、Gauss型求积公式等数值积分公式及方法。 §1．数值积分的一般概念 本章讨论定积分的近似计算问题。从微积分学中我们知道能够利用Newton-Leibni 公式 去计算的定积分是很少的。事实上，在实际问题中，我们常常无法利用初等函数去表出原函数例如，对于概率积分与椭圆积分 和 来说，我们便遇到了上述的困难。因此不能不考虑定积分的近似计算问题。 以下，我们所讨论的求积公式绝大多数具有如下形式： （1.1） 其中为求积公式的结点，为求积系数。通常，称右端的和为求积和； 又称 为求积误差。有时，也将求积公式写成 在（1.1）式中，[a,b]是实直线上的有限或无限的空间；函数是已知的固定的函数且常常是，以后我们将称它为权函数。此外，我们还假定积分 总是存在的，并且函数在点处是有定义的。 一般说来，求积公式（1.1）中的结点和系数可以按所希望的方式随意选取（除非是被积函数仅在一离散点集上是已知的，那时只好限制从离散点集中去选取了）。自然，我们总是希望通过和的选取使得在某种意义下求积误差尽可能地小。 概括来说，数值积分问题可分解为下述的三个主要问题： 求积公式的具体构造问题； 精确性程度的衡量标准问题； 余项估计问题（亦即，误差估计问题）。 为了解决第一个问题，我们必须考虑结点和求积系数的决定（或选择）问题。为了合理的解决第二个问题，我们将引进代数精度的概念。至于第三个问题，则主要是借助于内插多项式的余项估计公式来解决。 由第一章的Weierstrass多项式逼近定理可知，对于闭区间上的连续函数，都可以用多项式去一致地逼近它。换句话说，任一连续函数都可以用多项式作为它的最简单的近似函数，一般说来，多项式的次数取得越高，用它们来近似连续函数的程度也就越高。这自然使我们想到利用多项式的次数去规定求积公式的精确性程度（所谓代数精度）。 代数精度的概念是这样：就形如（1.1）式的求积公式来说，假如对，公式衡精确地成立（亦即），而当时公式不精确成立，则称公式（1.1）的代数精度为m。容易看出，m越大，则就一般的连续函数而言，公式（1.1）的右端数值与左端积分值的接近程度也就越高。事实上，当m越大时，用次数不高于m的多项式（例如p(x)）去近似f(x)亦就越好，即便越小，因而公式（1.1）的误差亦就越小。理由是， 由此可见，引进代数精度的概念作为衡量求积公式的精确性是十分自然的。 下面的定理说明了具有代数精度的求积公式的存在性。 定理1 对于任意给定的n个不同的结点，有常数使 得当是次数的多项式时求积公式（1.1）精确成立，亦即 证明 设已给定点，并且是函数f(x)在这些点上的Lagrange差值多项式，亦即 此处 于是 （1.2） 定义 则（1.2）式变为 但是，若是次数的多项式，则。这意味着故 证毕。 注意，上述定理并没有要求一定要属于区间[a,b].另外，定理只是断言了，当由（1.3）式决定时，公式（1.1）对于一切次数的多项式是精确的。换言之，定理1 说明了公式（1.1）的代数精度 以后，我们称求积系数由（1.3）式决定的求积公式为插值型求积公式。由于对次数不超过n-1次的任意多项式说来，所以n个结点的插值型求积公式的代数精度反之，容易证明，代数精度的n个结点的求积公式一定是插值型求积公式。事实上，因Lagrange插值基本多项式，而因求积公式的代数精度因此该求积公式对 必精确成立，即特别，当属于[a,b]时，我们称公式 （1.5） 为内插型求积公式，其中求积系数由下式确定： （1.6） §2．Newton-Cotes公式 设[a,b]是一有限区间，令依定理1，有常数使得求积公式 （2.1） 对于一切次数的多项式是精确的。