第三章数值积分.docVIP

第三章数值积分 第三章数值积分 PAGE 第三章数值积分 第三章 数值积分 在一元函数的积分学中，定积分 （） 的计算，主要是使用牛顿－莱布尼兹公式 （） 虽然上述公式在理论和实践中起着重要作用，但它并不能完全解决积分的计算问题，在许多情况下，上述公式是无法适用的.常见的有下列几种情况： （1）被积函数的原函数不是初等函数，如，等函数. （2）被积函数的具体表达式未知，只知道函数在某些点处的函数值. （3）的原函数虽然知道,但过于复杂，在实际应用中计算量较大，使用不便. 针对以上几种情况，就要考虑使用数值方法进行近似求解. 下面介绍几个计算定积分近似值的公式：矩形公式、梯形公式、辛普森(Simpson)公式、牛顿-柯特斯公式。 由积分中值定理可知，对于定义在区间上的连续函数，总存在一个使得 （） 由此受到启发，只要设法给出的一个近似值，便可得到一个计算定积分近似值的方法. 例如，用或，近似替代，则可以分别得到矩形公式 （） 或 （） 或 （） 进一步，若用近似替代，则可以得到下面梯形公式 ； （） 若用近似替代，则可以得到著名的辛普森(Simpson)公式 . （） 定理 若在上连续，则梯形公式的余项为 （） 定理 若在连续，则辛普森公式的余项为 . （） 练习： 分别用梯形公式、辛普森公式计算积分，并估计各种方法计算结果的误差. 一般地，若用；（，）近似替代，则得到定积分的近似计算公式 若令，则有 . （） （）通常称为机械求积公式，其中一般为常数，称为求积系数.为积分区间上固定点，称为求积节点.对于一个确定的求积公式，计算所得的值与被积函数有关.为方便起见，用来表示公（）式右端的值，即 显然，对于一个机械求积公式，通常需要考虑下列两个问题： （1） 确定求积系数和求积节点， （2） 讨论求积公式的误差. 而对于不同的求积公式，如何判断哪一个数值效果更好呢？对于节点相同的求积公式，通常用代数精度的高低作为判定这些公式优劣的一个标准.一般说来，代数精度越高，求出的积分近似值精度越好.下面给出代数精度的定义 定义 若数值积分公式 对任意不高于次的代数多项式都准确地成立，而对于却不能准确地成立，则称该数值积分公式的代数精度为次. 例 求积公式 具有3次代数精度. 事实上，当时，左端＝， 而右端＝，所以左端＝右端； 当时，左端＝， 而右端＝，所以左端＝右端； 当时，左端＝， 而右端＝，所以左端＝右端； 当时，左端＝， 而右端＝，所以左端＝右端； 当时，左端＝， 而右端＝右端. 所以该求积公式有3次代数精度. 同样可以证明公式（）、（）代数精度均为1次. 第二节 牛顿—柯特斯公式 一、 构造数值积分公式的基本方法 构造形如的数值积分公式的方法很多，常用的方法之一是利用被积函数的插值多项式来构造，即被积函数用其次拉格朗日型插值多项式的积分近似替代， 即 （） 其中， 为节点的插值基函数。 若记 （） 可以得到形如（）的公式，我们称由系数（）确定的数值积分公式（）为插值型数值求积公式. 积分的真值与其数值积分公式给出的近似值之差称为该数值积分公式的截断误差（或余项），记作.那么数值积分公式（）的误差为 = == 其中, . 显然，含有个节点的插值型数值积分公式的代数精度至少是次. 若将积分区间等分并取分点 , 相应的插值型积分公式为： （） 称上式为牛顿—柯特斯(Newton-Cotes)公式，其中称为柯特斯系数. 只要给出区间的等分数，就能算出，从而由公式（）可以得到相应的牛顿—柯特斯公式. 例如, 当时有 相应的牛顿—柯特斯公式为 （） 这就是前面所