数值积分与数值微分.docVIP

第六章 数值积分与数值微分 第一节 值积分的基本概念 7.1.1求积公式与代数精确度 积分中值定理告诉我们,如果函数在区间上连续,则在积分区间内存在一点,使成立。由于的具体位置一般是未知的,因而难以准确地计算出。如果能够提供一种求的算法,相应地便得到一种数值求积方法。 若近似地 用积分区间端点处的函数值与的算术平均值替,便导出计算积分的梯形公式 (7.1.1) 若近似地用积分区间中点处的函数值代替,导出计算积分的中矩形公式 (7.1.2) 一般地,所谓数值求积方法是指，在积分区间［a,b］上适当地选取某些节点 ,然后用加权平均得到,这样构造出的求积公式为 (7.1.3) 或写为 (7.1.4) 其中,称为求积节点,权称为求积系数,它仅仅与节点有关,称为余项。 为了保证数值求积公式的精度,我们自然希望求积公式能够对尽可能多的函数f(x)都准确成立,这在数学上常用代数精确度这一概念来说明。 定义 如果某个求积公式对于次数不超过的一切多项式都准确成立,而对 某个次多项式并不准确成立,则称该求积公式的代数精确度为。 显然,梯形公式(7.1.1)与中矩形公式(7.1.2)均具有一次代数精确度。 一般地，欲使求积公式(7.1.3)具有次代数精确度,只要令它对于都准确成立即可,即要求 (7.1.5) (7.1.5)式由个方程组成，包含有个节点以及个待定的求积系数。如果我们事先选定并且取，求(7.1.5)可确定，从而使求积公式(7.1.3)至少具有次代数精确度，如果适当选择及，求解(7.1.5)可能使求积公式(7.1.3)具有次代数精确度，由此可知，构造数值求积公式实际上是求与的代数问题。 7.1.2 插值型的求积公式 设给定一组节点，且已知函数在这些节点上的值，作插值函数,我们取 （7.1.6） 作为积分的近似值，这样构造出的求积公式 (7.1.7) 称为插值型的,其中求积函数通过插值基函数积分得出 (7.1.8) 由插值余项定理可知,对于插值型的求积公 式(7.1.7),其余项 (7.1.9) 其中，与有关， 对于次数的多项式，其余项等于0，因而插值型求积公式(7.1.7)至少具有次代数精确度；反之，如果求积公式(7.1.7)至少具有次代数精确度，此时它对于插值基函数应准确成立，即注意到，因而式(7.1.8)成立，即(7.1.7)为插值型的。 综上所述，我们有下面结论： 定理 形如(7.1.7)的求积公式至少有次代数精确度的充分必要条件是它是插值型的。 第二节 牛顿-科兹公式 7.2.1 牛顿-科兹公式 将区间划分为等分，步长为，其分点以此分点为节点构造出的插值型求积公式 （7.2.1） 称为科兹公式(Newton-Cotes)公式，其中 （7.2.2） 称为科兹系数，令,则有 （7.2.3） 当时，由(7.2.3)式可得科兹系数为 相应的求积公式是下列梯形公式 （7.2.4） 当时，由(7.2.3)式算得科兹系数为 相应的求积公式也称为抛物线公式(或辛普森(Simpson)公式)。 （7.2.5） 当时，牛顿-科兹公式为 （7.2.6） 它也特别称为科兹公式，其中， 表7-1中列出柯特斯系数表开头的一部分，从而可以建立相应的求积公式。 表7-1 1 2 3 4 5 6 7