无穷限广义积分的数值计算开题报告.docVIP

开题报告 无穷限广义积分的数值计算 一、选题的背景、意义 1．1 数值求积的定义 定积分的数值近似称为数值求积.[1]它起源于古代用铺贴小方块近似计算不规则图形或曲边形的面积．在近似积分中，主要从定义积分的黎曼和出发，用被积函数在积分区间上有限个点上值的加权和来近似计算积分． 1．2 无穷限广义积分的定义[2] 设函数定义在无穷区间上，且在任何有限区间上可积．如果存在极限 . 则称此极限为函数在上的无穷限反常积分(简称无穷积分)，记作 , 并称收敛．如果极限不存在，为了方便起见，亦称发散． 1．3 无穷限广义积分的发展 在讨论积分时有两个最基本的限制：积分区间的有穷性和有穷性和被积函数的有界性．但在很多实际问题中往往需要突破这些限制，考虑无穷区间上的“积分”． 比如现在人类要发射人造地球卫星或发射完成星际航行的飞行器，就要摆脱地球强大的引力，那如何离开地球呢 地球上的物体要脱离地球引力成为环绕太阳运动的人造行星，需要的最小速度是第二宇宙速度．第二宇宙速度为11．2公里/秒，是第一宇宙速度的2倍．地面物体获得这样的速度即能沿一条抛物线轨道脱离地球． 二、研究的基本内容与拟解决的主要问题 2．1数值积分的一般方法 许多定积分都无法用解析方法求出．对于那些并不知道函数的表达式只能通过实验得到在一系列点上的值的积分问题也只能用数值方法[3]． 2．1．1梯形法则[4] 把以曲线为曲边的曲边梯形分解成小曲边梯形以后，估计小曲边梯形面积的一个方法是用左矩形或右矩形面积代替小曲边梯形面积；但是这时误差会比较大．事实上，这种方法相当于用一系列的水平线逼近曲线．我们可以把这些水平线看成是函数的零次插值多项式． 一个更好的方法就是用一条折线逼近曲线；事实上，我们让小矩形的上边连续倾斜直到最好地拟合曲线．得到相应的求积公式是 , 对所有（即次数最多是1次的全体多项式）公式精确成立．此外，它的误差项是 , 其中．通过多项式逼近中的误差积分，再利用积分中值定理，可以确定梯形法则的误差项． 2．1．2复合梯形法则 如果划分区间为： . 那么在每个子区间上可应用梯形法则．这时结点未必是等距的．这样，我们得到复合梯形法则 . 对等间距及结点，复合梯形法则具有形式 ， 其中求和符号上的两撇表示求和式中的第一项和最后一项都被减半.复合梯形法则的误差项是 , 其中．对于每个子区间上的误差项求和并利用以下事实：在内存在一点使得，其中以及，即平均值，这样便得到总误差项． 2．1．3辛普森法则[5] 对任意区间的类似计算可得到熟悉的辛普森法则： . 从它的推导过程可知，对于所有次数的多项式辛普森法则是精确成立的．出乎意料的是， 对于所有次数的多项式它也精确成立． 与辛普森法则联系在一起的误差项是： ， 其中． 2．1．4 Gauss公式[6] 设有计算 的求积公式 , 其中求积节点，求积系数. 如果其代数精度为，则称为求积公式为Gauss-Legendre公式（简称Gauss公式），称相应的求积节点为Gauss点． 由代数精度的定义知，式为Gauss公式的充分必要条件是求积节点和求积系数满足下列方程组： . Gauss积分不但具有高精度，而且是稳定的，其原因是由于它的求积系数具有非负性． Gauss公式的求积系数全是正的． 高斯求积公式，[7]它不但具有最高的代数精度，而且收敛性和稳定性都有保证．因此是高精度的求积公式，高斯公式的主要缺点是节点和系数无规律，所以不便编程实现，在实际应用中，可以把低阶高斯公式进行复化． 2．2 无穷积分的敛散性判别 无穷积分的基本问题就是敛散性的判别问题，是求解无穷积分近似值的一个先决条件．由定义知道，无穷积分收敛与否，取决于函数在时是否存在极限．因此可由函数极限的柯西准则导出无穷积分的柯西准则． 无穷积分收敛的充要条件是：任给，存在，只要、，便有 ．