数值分析--第4章数值积分与数值微分[1].doc

第4章数值积分与数值微分
1 数值积分的基本概念
实际问题当中常常需要计算定积分。在微积分中,我们熟知,牛顿—莱布尼兹公式是计算定积分的一种有效工具,在理论和实际计算上有很大作用。对定积分,若在区间上连续,且的原函数为,则可计算定积分
似乎问题已经解决,其实不然。如
1)是由测量或数值计算给出数据表时,Newton-Leibnitz公式无法应用。
2)许多形式上很简单的函数,例如
等等,它们的原函数不能用初等函数的有限形式表示。
3)即使有些被积函数的原函数能通过初等函数的有限形式表示,但应用牛顿—莱布尼兹公式计算,仍涉及大量的数值计算,还不如应用数值积分的方法来得方便,既节省工作量,又满足精度的要求。例如下列积分
对于上述这些情况,都要求建立定积分的近似计算方法——数值积分法。
数值求积分的基本思想
根据以上所述,数值求积公式应该避免用原函数表示,而由被积函数的值决定。由积分中值定理:对,存在,有
表明,定积分所表示的曲边梯形的面积等于底为而高为的矩形面积(图4-1)。问题在于点的具体位置一般是不知道的,因而难以准确算出。我们将称为区间上的平均高度。这样,只要对平均高度提供一种算法,相应地便获得一种数值求积分方法。
如果我们用两端的算术平均作为平均高度的近似值,这样导出的求积公式
(4-1)
便是我们所熟悉的梯形公式(图4-2)。而如果改用区间中点的“高度”近似地取代平均高度,则可导出所谓中矩形公式(简称矩形公式)
(4-2)
更一般地,我们可以在区间上适当选取某些节点,然后用加权平均得到平均高度的近似值,这样构造出的求积公式具有下列形式:

图4-1 图4-2
(4-3)
式中称为求积节点;成为求积系数,亦称伴随节点的权。权仅仅与节点的选取有关,而不依赖于被积函数的具体形式。
这类由积分区间上的某些点上处的函数值的线性组合作为定积分的近似值的求积公式通常称为机械求积公式,它避免了Newton-Leibnitz公式寻求原函数的困难。对于求积公式(4-3),关键在于确定节点和相应的系数。
代数精度的概念
由Weierstrass定理可知,对闭区间上任意的连续函数,都可用多项式一致逼近。一般说来,多项式的次数越高,逼近程度越好。这样,如果求积公式对阶多项式精确成立,那么求积公式的误差仅来源于阶多项式对连续函数的逼近误差。因此自然有如下的定义
如果某个求积公式对于次数不超过的多项式均准确地成立,但对于次多项式就不准确成立,则称该求积公式具有次代数精度。
例1 判断求积公式
的代数精度。
解记
因为
所以求积公式具有5次代数精度。

最直接自然的一种想法是用在上的插值多项式代替,由于代数多项式的原函数是容易求出的,我们以在上的积分值作为所求积分的近似值,即
这样得到的求积分公式称为插值型求积公式。通常采用Lagrange插值。
设上有个互异节点,的次Lagrange插值多项式为
其中,插值型求积公式为
(4-4)
其中。可看出,仅由积分区间与插值节点确定,与被积函数的形式无关。求积公式(4-4)的截断误差为
(4-5)
求积公式
如其系数,则称此求积公式为插值型求积公式。
形如(4-3)的求积公式至少有次代数精度的充分必要条件是插值型的。
证明如果求积公式(4-3)是插值型的,由公式(4-5)可知,对于次数不超过的多项式,其余项等于零,因而这时求积公式至少具有次代数精度。
反之,如果求积公式(4-3)至少具有次代数精度,那么对于插值基函数应准确成立,并注意到,即有
所以求积公式(4-3)是插值型的。
求积公式的收敛性与稳定性
在求积公式(4-3)中,若
其中,则称求积公式(4-3)是收敛的。
实际使用任何求积公式时,除截断误差外,还有舍入误差,因此我们必须研究其数值稳定性。在求积公式(4-3)中,由于计算可能产生误差,实际得到,即,记
如果对任给正数,只要误差充分小就有
(4-6)
它表明求积公式(4-3)计算是稳定的,由此给出:
对任给,若存在,只要就有(4-6)成立,则称求积公式(4-3)是稳定的。
若求积公式(4-3)中系数,则此求积公式是稳定的;若有正有负,计算可能不稳定。
证明对任给,若取,对都有,则有
注意对任何代数精度的求积公式均有
可见时,有
(4-3)是稳定的。
若有正有负时,假设,且,有
它表明初始数据的误差可能会引起计算结果误差的增大,即计算可能不稳定。
2 Newton-Cotes公式
Cotes系数
被积函数在积分区间内变化平