定积分近似计算方法(共17页)

1、精选优质文档-倾情为你奉上 定积分的近似计算方法 摘要 本文主要讨论了一元函数常见的数值积分方法,例如插值型求积公式、龙贝格求积公式、高斯求积公式等近似计算方法,在用这些方法计算定积分时,会产生一些误差,为了减少误差, 可以利用复化求积公式、复化高斯公式等.本文围绕这些方法,系统介绍它们的计算公式以及截断误差,并用例题分析它们产生误差的大小、计算量等.关键词 插值型积分 龙贝格积分 高斯积分 误差分析 近似计算 1引言 在计算定积分的值时,常常根据微积分学基本定理求出的一个原函数,再用牛顿-莱布尼茨公式求的积分,.但在实际应用中,这种方法只限于解决一小部分定积分的求值问题.当函数没有具体表达式

2、,只是一些实验测得数据形成的表格或图形或者是无法用初等函数表示,例如,等等,这就需要我们用一些近似方法求的积分值.与数值积分一样,把积分区间细分,在每个小区间上,找到简单函数来近似代替,且的值容易求的.这样就把计算复杂的转化为求简单的积分值.因此,定积分的近似计算实质上是就是被积函数的近似计算问题.2常见数值方法 2.1牛顿-科茨数值方法牛顿-科茨求积公式是求积节点等距离分布的插值型求积公式.利用插值多项式来构造数值积分公式是最常用、最基本的方法,具体做法是： 给定区间上一组节点,以及节点处函数,作的次拉格朗日多项式 , 其中 ,将插值公式专心-专注-专业 . 其中 ,依赖于变量, 上式积分得

3、 若记 . (1) , (2)则有 (3)称式(3)为插值求型公式,其中. 与无关,叫求积系数, 为求积节点, 为求积公式余项,其中求积系数由(1)决定. 2.1.1梯形求积公式1梯形公式当插值节点分别选取区间端点时,由式(3)分别求出求积系数 , .从而的求积公式 . (4)称求积公式(4)为梯形求积公式,简称梯形公式. 2梯形公式截断误差: . (5) 3梯形求积公式的代数精度:1当时,式(5)中 .精确成立.2.1.2 辛普森求积公式 1辛普森求积公式 当选取节点为时,由式(1)求下列求积系数 , . . 从而求积公式 . （6）称式(6)为抛物线积分公式或辛普森积分公式.2抛物线求积公

4、式误差估计 定理1.若在上有四阶连续导数,则抛物线公式(6)的余项为: . (7)3抛物线公式的代数精度为3. 易验证,当时,式(6)精确成立,而当时,式(6)不能精确成立. 2.1.3 牛顿-科茨公式 1牛顿-科茨公式 在等距离节点下,其中. .作为变量替换,那么由求积公式(1),得系数: (8) 则 (9) 于是差值求积公式为: (10)称公式(10)为牛顿-科茨求积公式,其中称为科茨系数.显然,科茨系数与被积函数及积分区间无关,它指依赖于,且为多项式积分.因此,只要给出,就能看出,并写出相应地牛顿-科茨公式.2牛顿-科茨公式的截断误差与代数精度. 当与情况分析牛顿-科茨公式的截断误差为

5、牛顿-科茨公式的截断误差还可以写成 为奇数） (11)其中,且不依赖于,对为任何并不超过次多项式,均有,因而,即精确成立,也就是说,牛顿-科茨公式的代数精度至少为,牛顿-科茨公式在为偶数时,至少具有次代数精度,在为奇数情况时,至少具有次代数精度.2.1.4复化梯形求积公式将区间等分,节点为 (步长),)在每个小区间上采用梯形公式(4)得 (12) 称式(12)为复化梯形公式.复化梯形公式余项为 (13)2.1.5复化辛普森求积公式在每个小区间上,辛普森公式(6)得 （14） 记 (15)式中,为的中点,即.式(15)称为复化辛普森公式,其余项为 , 故 (16) 为复化辛普森的截断误差. 2.

6、1.6复化科茨求积公式 将区间等分, ,为正整数,在每个子区间上用科茨求积公式得到复化求积公式: (17)其中 , 其截断误差为. 2.1.7 变步长复化求积方法复化求积公式虽然计算简单,也达到了提高精度的目的,但为了满足精度要求必须顾及误差,利用误差公式往往很困难,因为误差表达式中含有未知函数的导数,而估计各阶导数的最大值不太容易.我们可以采取把积分的区间细分的办法,在计算积分时将步长逐步折半,利用前后两次结果进行误差估计,如此继续,直到相邻两次结果相差不大,取最小的步长算出的结果为积分值,这种方法称为变步长积分法.以复化梯形公式为例,把区间分成等分,设复化梯形公式的近似值为,原积分值为,由

7、复化梯形公式误差公式(14)知: 再把区间分成等分,得近似值,则 假定在上变化不大,既有.由上式得 . 于是 (18)式(18)表明若用作为的近似值,其截断误差约为 (19) 2.2 龙贝格求积公式龙贝格积分法的基本思想是采用复化梯形求积方法不断折半步长过程中,在积分结果中加入时候误差估计值进行补偿,使积分计算的收敛性加速,就可以加工出精度较高的积分结果.由式(19), 的误差大致为,因此,可用这个误差值作为的一种补偿,加到上,则可得到积分准确值,比的更好近似值. (20)式(20)左端时 记 恰好为上应用辛普生公式(16)的结果.在每个小区间应用辛普生公式: 代入式(20)的左端得 从而复化

8、辛普森公式与复化梯形公式公式有以下关系式 (21)类似也可以推证,在辛普森序列基础上,利用以下关系式 (22)可以造出收敛速度更快的科茨序列将此推行下去,在科茨序列基础上,通过 (23)构造出收敛速度比科茨序列更快的龙贝格序列.以上这种通过逐步构造龙贝格序列的积分近似值法就称为龙贝格积分法.2.3高斯求积公式 由定理知,插值型求积公式的代数精度与求积节点的个数有关,具有个节点的插值型求积公式至少具有次代数精度.不仅如此,代数精度与节点的选取有关,在构造牛顿-科茨求积公式时,为了简化处理过程,限定用等分节点作为求积节点,这样做,虽然公式确实得到简化,但同时也限制了公式的代数精度.设积分本段讨论如

9、下求积公式 (24) 对任意积分区间,通过变 可以转换到区间上,这时 此时,求积公式写为 若一组节点使插值型求积公式(24)具有次代数精度,则称此组节点为高斯点,并称相应求积公式(24)为高斯求积公式. 2.3.1 高斯求积公式的余项 其中 ,且不依赖于. 2.3.2 复化高斯求积公式 复化高斯求积公式的基本思想是:将积分区间分成个等长小区间,然后在低阶()高斯求积公式算出近似值,最后将他们相加的积分的近似值,即 (25)其中,与可由书中表中查出.3 应用3.1插值型积分的应用例1 用牛顿-科茨公式()计算积分. 解 时 例2 利用复化梯形求积公式计算积分 解 设,分点个数为1,2,4,5时,

10、求出相应积分, 列表如下:=1的计算结果见表1-1所列10.50.00.51.00.80.45=2的表格如下20.250.000.250.501.000.0.800. =4时计算结果如下表40.1250.000.1250.250.3750.501.000.0.0.0.800.= 5时计算结果如下50.10.00.10.20.30.40.51.00.0.0.917430.0.80. 例3 利用复化求积公式,问积分区间为多少等分才能得证有5位有效数字？ 解 由式（14）知 有,当时,在,所以 由于的准确值具有一位整数,所以要使近似值具有5位有效数字,必须满足 取对数有 .即将区间19等分可满足给定

11、的精度要求. 例4 利用复化抛物线求积公式计算 .解 设,取=1,2, 3时,公式当=1,2,3时结果如下表所示当=1时10.251.00.0.800.当=2时20.1251.00.0.0.0.800.当=3时30.83331.00.0.0.0.90.852070.80.4636 例5 用复化梯形公式,辛普森公式和科茨公式计算积分的近似值. 解按精度要求确定分多少等分,即确定步长,要使,只需令,则 所以只要取=4即可, 当时,在每个子区间上用式(25）,或(14),或(17),结果3.2 龙贝格积分公式应用 例6 用龙贝格算法计算积分的近似值,要求误差小于. 解 步骤如下: 得 计算由此得 .

12、 (3)算出从而 (4)计算从而得到: , (5)再计算 从而得到: , , 所以 3.3高斯求积公式的应用例7 用两点复化高斯求积公式计算要求允许误差 解 在本算法中取时,其中 =2时, =, =3时, =. 3.4 几种方法的比较分析例8 计算积分，精确到0.001.(1)利用矩形公式计算, 因为对于,有(如果1<<2),所以按照公式 . 0<<. 如果取=10,则我们公式的余项的余数得,我们还必须加进由于在计算函数值实行四舍五入所产生的误差的界限相差于0.16,为了这个目的只要计算的值到四位小数精确到0.00005就够了.我们有 和6.9284 (2) 按照梯形公

13、式作同样的计算,在这种情况下,作公式 在这儿也试一试取=10,虽然此时仅可以证,纵坐标是 和 (3) 用辛普森公式做同样的计算 作公式 并且=5时有.实行计算到五位数字,精确到0. . 由此可见，用辛普森公式计算得到的值误差最小，计算量相对一般；而用矩形公式计算得到的值误差较大，计算量也比较大；用梯形公式计算的值误差比用矩形公式得到的值要误差小，计算量也是如此.所以我们计算定积分时用辛普森公式往往得到的值误差小，而对没有要求误差大小的，则可以选择辛普森或者是梯形公式，因为这两种方法计算量相对较小. 结 束 语 本文只讨论了一些一维数值积分方法及其它们的应用,误差分析等有关内容.其中最常用的方法

14、是插值型积分以及复化方法、龙贝格积分方法和高斯积分方法,并讨论了相关求积方法的代数精度和误差分析,并给出了一些例题,分析各种方法的近似值,得出误差分析最小的近似方法.由于篇幅有限,对于高维数值积分方法本文便不再讨论. 参考文献 1 华东师范大学数学系，数学分析（第一版）M，北京:高等教育出版社,2001.2 李庆阳,关治,白峰杉，数值计算原理（第二版）M，北京: 清华大学出版社, 2008.3 肖筱南，现代数值计算方法（第一版）M，北京: 北京大学出版社, 1999.4 菲赫金格尔茨，微积分学教程（第三版）M，北京: 高等教育出版社, 2005. 5 裴礼文，数学分析中的典型问题与方法（第一版

15、）M ，北京: 北京大学出版社,2004.6 李桂成，计算方法（第三版）M，北京: 高等教育出版社,2010.7 ,.Calculating Skillfully the Curve Integral and Surface Integral Type 2 by Symmetry, ,The Approximate Numerical Method of the Definite IntegralAbstract This paper mainly discusses common numerical methods of unary function, such as approximate calculation method of interpolation integral, Lebesgue integral and Gauss integration. With these methods in calculating the integral, it will produce some error. In order to reduce the error, we can use after the formula for product and after the Gauss formula. This paper focus on these met