安庆师范学院数学与计算科学学院2014届毕业论文 第1页 共11页 定积分的近似计算方法与误差估计 作者: 操乐青 指导老师: 邢抱花

摘要 本文主要讨论了一元函数常见的数值积分方法,例如插值型求积公式，高斯求积公式等近

似计算方法,在用这些方法计算定积分时,会产生一些误差,为了减少误差, 可以利用复化求积公式、复化高斯公式等.本文围绕这些方法,系统介绍它们的计算公式以及截断误差,并用例题分析它们产生误差的大小、计算量等. 关键词 插值型积分 高斯积分 误差分析 近似计算

1引言 在计算定积分的值()baIfxdx时,常常根据微积分学基本定理求出)(xf的一个原函数)(xF,再用牛顿-莱布尼茨公式求得积分,()()()baIfxdxFbFa.但这种方法只限于解

决一小部分定积分的求值问题.当函数没有具体表达式,只是一些实验测得数据形成的表格或图形或者是()Fx无法用初等函数表示,例如,2bxaedx,2sinbaxdx等等,这就需要我们用一些近似方法来求积分值. 与数值积分一样,把积分区间细分,在每个小区间上,找到简单函数)(x来近似代替()fx,

且()baxdx的值容易求的.这样就把计算复杂的()bafxdx转化为求简单的积分值()baxdx.因此,定积分的近似计算实质上就是被积函数的近似计算问题. 2 定积分的近似计算——常见数值方法 2.1 矩形公式 根据定积分的定义，每一个积分和都可以看作是定积分的一个近似值，即

1()d()nbiiaifxxfx

在几何意义上，这是用一系列小矩形面积近似小曲边梯形的结果，所以把这个近似计算方法称为矩形法．不过，只有当积分区间被分割得很细时，矩形法才有一定的精确度．

针对不同i的取法，计算结果会有不同，常见的取法有：

（1）左端点法，即1ii

x,iabniixxfdxxf11)()(

（2）右端点法，即 iix,iniiabxxfdxxf1)()( 安庆师范学院数学与计算科学学院2014届毕业论文 第2页 共11页 （3）中点法，即12iiixx,iniiiabxxxfdxxf11)2()(

例1 用矩形公式近似计算积分 12 0d1xx（取100n）. 解 对1,0作n等分 bxinabaxxaxni10，由定义知：



niiniifnxfxdx111

1

0)(1)(

1

(1)左点法：在区间],[1iixx上取左端点，即取1iix，ni2,1 12 0

1d()1niiixfxx





0.78789399673078，

理论值 12 0d14xx，此时计算的相对误差 0.7878939967307840.0031784

（2）右点法：在区间],[1iixx上取右端点，即取iix，ni2,1 12 0

1d()1niiixfxx





0.78289399673078，

理论值 12 0d14xx，此时计算的相对误差 0.7828939967307840.0031884

（3）中点法：取 在区间1[,]iixx上取中点，即取12iiixx，ni2,1 12 0

1d()1niiixfxx





0.78540024673078，

理论值 12 0d14xx，此时计算的相对误差 60.7854002467307842.653104 安庆师范学院数学与计算科学学院2014届毕业论文 第3页 共11页 如果在分割的每个小区间上采用一次或二次多项式来近似代替被积函数，那么可以期望得到比矩形法效果好得多的近似计算公式．下面介绍的梯形法和抛物线法就是这一指导思想的产物．

2.2梯形公式 等分区间 bxinabaxxaxni10，nabx

相应函数值为 nyyy,,,10（nixfyii,,1,0),(）．

曲线)(xfy上相应的点为

nPPP,,,10（niyxP

iii,,1,0),,(

）

将曲线的每一段弧iiPP1用过点1iP，iP的弦iiPP1（线性函数）来代替，这使得每个],[1iixx

上的曲边梯形成为真正的梯形，其面积为

xyyii21，ni,,2,1．

于是各个小梯形面积之和就是曲边梯形面积的近似值， 11 11()d()22nnbiiiiaiiyyxfxxxyy





即 011 ()d()22bnnayybafxxyyn 称此式为梯形公式． 例2 用梯形公式近似计算定积分 12 0d1xx（取100n）.

解 10112 0d()122nnyyxbayyxn0.78539399673078， 理论值 12 0d14xx，此时计算的相对误差 60.7853939967307845.305104

很显然，这个误差要比简单的矩形左点法和右点法的计算误差小得多. 安庆师范学院数学与计算科学学院2014届毕业论文 第4页 共11页 2.3 抛物线公式 由梯形法求近似值，当)(xfy为凹曲线时，它就偏小；当)(xfy为凸曲线时，它就偏大．若每段改用与它凸性相接近的抛物线来近似时，就可减少上述缺点，这就是抛物线法．

将积分区间],[ba作n2等分，分点依次为

bxinabaxxaxni2102

，

nabx2

，

对应函数值为

nyyy210,,,（nixfyii2,,1,0),(），

曲线上相应点为 nPPP210,,,（niyxPiii2,,1,0),,(）．

现把区间],[20xx上的曲线段)(xfy用通过三点),(000yxP，),(111yxP，),(222yxP的抛物线 )(12xpxxy

来近似代替，然后求函数)(1xp从0x到2x的定积分： 20 1 ()dxxpxx20

2 ()dx

xxxx)()(2)(30220223032xxxxxx



]4)(2)()()[(62022022202002xxxxxxxxxx

由于2201xxx，代入上式整理后得 20

1 ()dx

xpxx)](4)()[(612122202002xxxxxx

xx

)4(621002yyyxx)4(6210yyynab

同样也有 42

2 ()dx

xpxx)4(6432yyynab

…… 222 ()dnnxnxpxx

)4(621222nnnyyynab



将这n个积分相加即得原来所要计算的定积分的近似值： 安庆师范学院数学与计算科学学院2014届毕业论文 第5页 共11页 222 22212 11()d()d(4)6iinnbxiiiiaxiibafxxpxxyyyn



，

即 0213212422 ()d[4()2()]6bnnnabafxxyyyyyyyyn

这就是抛物线公式，也称为辛卜生（Simpson）公式． 例3 用抛物线公式近似计算积分 12 0d1xx(取100n).

解 102132124222 0d[4()2()]16nnnxbayyyyyyyyxn =0.78539816339745， 理论值 12 0d14xx，此时计算的相对误差

160.7853981633974542.827104

2.4 几种近似计算定积分方法的比较分析及误差估计 例4 计算积分211ln2dxx，精确到0.001.

解 方法(一) 利用矩形公式计算, 因为对于xxf1)(,有320()2fxx(如果1所以按照公式 0)2(dxbaxab. 0

如果取n=10,则我们公式的余项的余数得31010.84101200R,我们还必须加进由于在计算函数值实行四舍五入所产生的误差的界限相差于0.16310,为了这个目的只要计算1x的值到四位小数精确到0.00005就够了.我们有

1232527292132152172192

1.051.151.251.351.551.651.751.851.95xxxxxxxxx 5128.05405.05714.06061.06897.07407.08.08696.09524.02192172152132927252321yy

yyyyyyy

Y的和计算6.9284 故计算结果为 69284.0109284.6。