无穷限广义积分的数值计算[文献综述] .pdfVIP

毕业论文文献综述

信息与计算科学

无穷限广义积分的数值计算

一．前言部分

[1]

定积分的数值近似称为数值求积.它起源于古代用铺贴小方块近似计算不规则图形或

曲边形的面积．在近似积分中，主要从定义积分的黎曼和出发，用被积函数在积分区间上有

限个点上值的加权和来近似计算积分．

我们一般使用牛顿-科茨求积公式，梯形公式及其复合公式，辛普森公式及其复合公式，

Gauss求积公式，切比雪夫求积法，三次样条函数求积法，自适应积分法等方法来进行数值

求积．

在讨论积分时有两个最基本的限制：积分区间的有穷性和被积函数的有界性．但在很多

实际问题中往往需要突破这些限制，考虑无穷区间上的“积分”．

根据函数的变化率，利用定积分我们可以计算函数在指定区间上的增量，利用变限定积

分可以把握函数变化区间上增量的变化，为了把握函数在无穷区间上增量的变化，我们还需

[2]

要引进并讨论无穷限积分．

比如现在人类要发射人造地球卫星或发射完成星际航行的飞行器，就要摆脱地球

强大的引力，那如何离开地球呢?

地球上的物体要脱离地球引力成为环绕太阳运动的人造行星，需要的最小速度是

第二宇宙速度．第二宇宙速度为11．2公里/秒，是第一宇宙速度的2倍．地面物体获

得这样的速度即能沿一条抛物线轨道脱离地球．

我们可以运用无穷限广义积分解决第二宇宙速度问题．

在黎曼积分的定义中，被积函数和积分区间都是有界的．若被积函数或积分区间无界，





a，bfa,b

则称为广义积分．对无界区间，如，如果对任何有限的，在区间上可积，

并且下列极限存在且为有限数，则广义积分的定义为

b





fxdxlimfxdx

．

aba

对无界的积分区间，可以使用有限区间上的标准求积程序计算广义积分，具体方法如下：

•用有限的积分区间代替无限的积分区间．选择积分范围时要注意所截掉的部分应是极

小的，另外应对这一部分在整个积分中所占的份额作出估计．同时这个有限区间也不应太大，

以免在利用自适应求积程序时，陷入无休止的积分函数调用之中．

•通过适当的变换将无界区间变成有界区间．典型的变换包括，xlnt或者

t

x．但是在变换的时候一定要注意不要引入新的奇异点或产生其它问题．



1t

还有一种方法就是采用专门计算无界区间积分的求积公式，比如说高斯-拉盖尔

（Gauss-laguerre）或者高斯-艾尔米特求积公式．

一般采用变量替换，无穷区间的截断，无穷区间上的高斯求积公式，极限过程等方法去

解决无穷限广义积分的数值计算．

二．主题部分

2．1数值积分的一般方法



许多定积分都无法用解析方法求出．对于那些并不知道函数fx的表达式只能通过实



验得到fx在一系列点上的值的积分问题也只能用数值方法．[3]

2．1．1梯形法则[4]



把以曲线fx为曲边的曲边梯形分解成小曲边梯形以后，估计小曲边梯形面积的一个

方法是用左矩形或右矩形面积代