几种定积分的数值计算方法摘要:本文归纳了定积分近似计算中的几种常用方法，并着重分析了各种数值方法的计算思想，结合实例，对其优劣性作了简要说明•关键词:数值方法;矩形法;梯形法;抛物线法;类矩形;类梯形Several Numerical Methods for Solvi ng Defi nite In tegrals Abstract:Several com mon methods for solvi ng defi nite in tegrals are summarized inthispaper. Mean time, the idea for each method is emphatically an alyzed. Afterwards, anu merical example is illustrated to show that the adva ntages and disadva ntages of these methods.Keywords: Numerical methods, Rectangle method, Trapezoidal method, Parabolic method, Class recta ngle, Class trapezoid1.引言在科学研究和实际生产中，经常遇到求积分的计算问题 ，由积分学知识可知，若函数f(x)在区间［a,b ］连续且原函数为F(x),则可用牛顿-莱布尼茨公式baf(x) F(b) F(a)求得积分.这个公式不论在理论上还是在解决实际问题中都起到了很大的作用 •在科学研究和实际生产中，经常遇到求积分的计算问题，由积分学知识可知 若函数f(x)在区间［a,b ］连续且原函数为F(x),则可用牛顿-莱布尼茨公式ba f(x) F(b) F(a)求得积分.这个公式不论在理论上还是在解决实际问题中都起到了很大的作用 •另外，对于求导数也有一系列的求导公式和求导法则.但是,在实际问题中遇到求积分的计算，经常会 有这样的情况：(1) 函数f(x)的原函数无法用初等函数给出.例如积分等，从而无法用牛顿-莱布尼茨公式计算出积分(2) 函数f(x)使用表格形式或图形给出，因而无法直接用积分公式或导数公式。
(3) 函数f(x)的原函数或导数值虽然能够求出，但形式过于复杂，不便使用. 由此可见,利用原函数求积分或利用求导法则求导数有它的局限性，所以就有了求解数值积分的很多方法，目前有牛顿一柯特斯公式法，矩形法梯形法，抛物线法，随机投点法，平均 值法,高斯型求积法，龙贝格积分法，李查逊外推算法等等，本文对其中部分方法作一个比较.2•几何意义上的数值算法s 在几何上表示以［a,b ］为底,以曲线y f(x)为曲边的曲边梯形的面积 A,因此，计算 s 的近似值也就是A 的近似值，如图1所示.沿着积分区间［a,b ］,可以把大的曲边梯形分割 成许多小的曲边梯形面积之和.常采用均匀分割,假设［a,b ］上等分n 的小区间b a1 2o e x dx . 1sin xdxX i-1 X i h,x。
a,X n b,其中h ——表示小区间的长度•n2.1矩形法矩形法就是用小矩形面积近似代替各个小曲边梯形面积 ，从面积得到S 的近似值.若取2.2梯形法梯形法则用小直边梯形的面积近似代替小曲边梯形面积，见图2,从而得到S 的近似值, 即A匚 f ⑻f(b) n1f (x i ).n 2 i i图2分割曲边梯形近似积分2.3抛物线法抛物线法以抛物线为曲边梯形的曲边，曲边梯形的面积近似代替小曲边梯形的面积，如 图3所示.图3抛物线积分yk£ E小区间左端点的函数值为小矩形的高 ，如图1中所示，则A(x)图1 分割曲边矩形近似积分*X o,X i,X2对应的曲线上的点P O,R,P2可以唯一地确定一条抛物线y ax2 bx c ,这条抛物线将作将代替从X o至X2的曲线段，此时积分可以转化为对抛物线积分，而抛物线的积分可以利用牛顿一莱布尼玆公式.第1、2个小区边梯形的面积：X2/ 2hA ( ax bx c)dx [f(x0) 4f(xJ f(x2)]X0 3上面利用了条件F0,R,P2是抛物线上的点以及等式X2 X o 2X1 .同理可证：A2 尹Q EQ f(X4)]hA n/2 -[ f (X n 2) 4f(X ni) f (X n )]n/ 2 n/2 1所以,S A1 A2 A n/2 普{[f(a) f (b)] 4 f (X2i 1) 2 f^)}i 1 i 13•概率意义上的数值算法概率算法是定积分问题数值求解的一类常用方法，其设计思想简单，易于实现•尽管算法要耗费较多计算时间，但是往往能得到问题的近似解，并且近似程度能随计算时间的增加而不断提高•概率算法可用于计算定积分的近似值•3.1平均值法b考虑定积分I f (x)dx的近似计算，其中f (x)在a,b内可积,用平均值法计算该积a分,首先随机产生n个独立的随机变量，且服从在a,b上均匀分布，即i(i 1,2, n);其次，计算I的近似值「,rn i1f(i).由中心极限定理知，若i (i 1,2, n)相互独立、同分布，且数学期望及标准差在，则当n 充分大时，随机变量Y -——-渐近服从正态分布N(0,1),即对任意的t 0,i这表明，用平均值法计算定积分的收敛速度较慢，在概率意义下的误差阶仅为0( 1).n3.2 “类矩形” Monte-Carlo 方法1由于平均值法计算定积分的收敛速度较慢，且在概率意义下的误差阶仅为 0(丄)，就Un 有对平均值法的改进，“类矩形” Monte-Carlo 方法,改进过程为：先将积分区间a,b n 等 分,随机产生n 个相互独立且服从0,1上均匀分布的随机变量序列｛ i ｝, (i 1,2, n)；然后 由这n个随机点类似于矩形公式构造计算公式 ，即作变换b a /ia( i ni 1), i1,2, n将｛ i ｝映射到子区间i 1 {a(b a), an 丄(b a)} n a,b ,i 1,2, ,n最后，计算1的近似值~,~ 口nf( i )•n i 1F 面用两个命题证明“类矩阵”方法的可行性.命题 1 设 f(x) C 1 a, b ,记 M maxf(x), x 0a,b ,有x a,b(b a)22证明：由Lagrange 中值定理得上式两边在a,b 积分，得由f (x)得连续性，得f(x)dxf (x) f(x °) f ( )(x X 。
)(介于x 与X 。
之间)f (x)dxf(x °)(b a) f ( )(x X o )dx1 O(-). n3.3 “类梯形” Monte-Carlo 方法再给出平均值法的另一种改进.首先将a,b n 等分,再在每个子区间上随机产生 2n 个 相互独立且服从［0,1］上均匀分布的随机变量序列，并两两分组，得{ 2i 1, 2i }, (i 1,2,3, ,n); 做变换b a . 2i 1a(2i 2 2i 1)2n 2ia (2i 1 2i )2n)(x x 0) dx M x x 0 dx2 Xo(a b)x o12(a2 b 2) AM.命题 2 设 f (x) C 1 a,b ,hM max f (x), M ix a ,bmax(i 1)h,aihf (x) ,i1.2.证明：1 〜bb a n a f(x)dxf( i )an i 1f (x)dxa (i 1)hf( i )a ih a (i 1)ha ih a(i1](艸f (x) dx hf ( i )hf( i )Xi1,2, ,n于是M(b a)2 2n〜与I 如上，则〜与I 的误差满足 Innna ih将2i 1, 2i 分别映射到子区间公式构造“类梯形”公式近类似a ih a (i 1)hf (X)dX.F 面证明“类梯形”方法可行性的两个命题:所以 P(X i ) f(X i ),i 1,2.令将X 看成a,b 上的一个定点，构造辅助函数(t) f(t) P(t) k(x)(t St X 2)由于(xj(X 2)(x) 0,由Rolle 中值定理，'(t)在a,b 内至少有两个零点，对'(t)再用Rolle 中值定理，知''(t)在a,b 内至少有一个零点，即存在 a,b ,使''()f''( ) 2k (x) 0,所以k(x)—° .将它代入(1)式,并两段同时从a 到b 积分,2i 1 [a (b a), a n [a 字(b2na),a字(b a)]2n 丄(b a)],in1,2,3, ,n然后在每个等分子区间上 [ai 1(b a),a n 丄(b a)］利用2i 1, 2i 两点类似于梯形nS i jf( n2i 1)f( 2i )]~ ~ b 最后计算I 的近似值1,1 ban f ( 2i 1)2f( 2i ) 命题3设f xC a, b ，记 M=max f'' xx a,b，则 X i ,X 2X| X 2，有bf x dxaf x 1 f x 2b12证明: 过 X i , f x i , X 2, f X 2两点的直线方程为 P(x) f x 1U —(X xjX 2 X-!R(x) f(x)P(x) k(x)(x X i )(x X 2)(1)12bb a f x dxa2bf''()f x 1 f x 2(x x 1)(x x 2)dx a2M b (7 (x〜aX)(Xaxj(x x 2)dxX i )(X x 2)dxx?(x N )(Xbx 2)dx (x x 1)(x x 2)dxX 2L(X i ，x ?)X i (x axj(xx 2)dxx ? (x x 1)(x x 2)dxX 1b(X X 2x 1)(x x 2)dx不妨设a x 1 x 2a b 、、2L x1(b a )(X 22 ) (X 2X 1)a b 、、2L X 2 (b a)( X 12) (X 2X 1)X 1 1— (3a b)b ,则将L(x i ,x ?)分别对求偏导数0 0解得唯一驻点: 3b)x ?，得；(a旦)(b416a) , L(a,b) (b a)3 6故当a x 1 x 2b 时,bf x dxax 1 f x 2b 结论成立. 命题4设f (x)2C a,b ,■^ydx ，该积分精确值为3.1416.下面分别给出本文所涉及计算方法对它x的计算结果:4.1用三种基于几何意义的算法：矩形算法,梯形法，抛物线法作比较，结果如表1:梯形3.1416496i.7e 5b—a, M max x a,b f(X ),M i xa max ’aih(x),i 1,2, nI 与I 如上，则I的误差满足:I证明: ba f(x)dxf (2i 1 2i)a ih由命题3,得a ih 于是 4.例题a(i 1)h f(X )dX/ f( 2i 门f(-2i )i 1a ih(i 1)h f(x)dxa (i 1) hf(x)dxa (i 1)hf (2i 1)hi 1121 or ). nf( 2i 1) f( 2i)hf( 2i )(b a)3M 12n 2M-h 3,i 1,2, ,n12对于积分 01抛物线 3.14160128 4.0e 114.2用平均值法，及其改进“类矩形” Monte-Carlo方法，“类梯形” Monte-Carlo方法计算结果如表2:表2概率意义算法的比较节点数算法近似值误差平均值法 3.13849728 3.1e 3类矩形法 3.14169037.0e 5104类梯形法 3.14160029 2.0e 95 .结语本文介绍的几种求积公式各有特点：梯形求积公式和抛物线法求积公式是低精度公式但对于光滑性较差的被积函数有时比用高精度方法能得到更好的效果，尤其是梯形求积公式.当被积函数为周期函数时，效果更为突出.由表1分析，一般情形下，三种基于几何的算法中矩形算法的误差最大，梯形法次之，抛物线法最高.抛物线法的积分精度远远高于另外两种方法，特别是在积分区间分割份数较小的情况下，仍然保持较高的近似程度.“类矩形” Monte-Carlo方法；“类梯形” Monte-Carlo方法是平均值法的改进提高了平均值法的精确度.通过表2可以看出，直接用平均值法计算定积分,104节点的计算已经很可观了，但计算结果只有2位有效数字，而选取同样的节点数，计算量几乎不变，类矩阵法就达到了4位有效数字，类梯形法则达到了8位有效数字，恰好与上述定理中误差阶的估计是一致的，从而也验证了“类矩形” Monte-Carlo方法和”类梯形” Monte-Carlo方法的高效性.从表2中也可以看出随着节点数的增大，积分精度会不断提高，当然计算复杂度就会增加.参考文献[1] 费祥历刘奋,马铭福•高等数学（第2版上册）[M].:石油大学出版社，2008: 211-287 .[2] 徐萃薇，孙绳武•计算方法引论（第三版）•:高等教育出版社,2007.[3] 王晓东计算机算法分析与设计[M].:电子工业出版社,2001:197-228.[4] 徐钟济•蒙特卡罗方法[M].:上海科学技术出版社,1985.⑸张平文，李铁军，数值分析[M].:北京大学出版社,2007.⑹阮宗利.计算一元定积分的若干数值算法及其比较[J].中国石油大学学报（科技教育）,2010:182-184.[7] 朱长青.计算方法及其应用.科学出版社,2006.[8] 张威,刘志军，李艳红.数值分析与科学计算.:清华大学出版社,200,5[9] 明万元，郑华盛.求解数值积分的两类新的Monte-Carlo 方法[J].南昌航空大学学报,2010,40（10）:180-186.[10] 刘长虹，关永亮等.蒙特卡洛在数值积分上的应用[J].上海工程技术大学学报,2010,24（1）:43-46.[11] 王岩.Monte-Carlo方法应用研究[J].云南大学学报（自然科学版）,2006,28（SI）:23-26.。