数值计算程序设计数值积分.docVIP

《数 值 方 法》 课程设计 内容提要： 实际问题中常常需要计算积分，有些数值方法，如微分方程和积分方程的求解，也都和积分计算相联系。依据微积分基本定理，对于积分只要找到被积函数f(x)的原函数F(x)。F’(x)=f(x)，便有牛顿―莱布尼兹公式来求得积分．但实际运用这种方法往往有困难，因为大量的被积函数，找不到用初等函数表示的原函数；另外当f(x)是由实验测量或数值计算给出的一张数据表时，牛顿―莱布尼兹公式也不能直接运用，因此，实际往往通过其他的一些方法如利用复化梯形公式、龙贝格公式、牛顿―莱布尼兹公式等来求解。 复化梯形求积法通常把积分区间等分成若干个子区间，在每一个子区间上用低阶的求积公式，对所有的子区间求和即得整个区间 [a,b] 上的积分公式。 龙贝格算法是在区间逐次分半过程中，对梯形公式的近似值进行加权平均获得准确程度较高的积分近似值的一种方法，它具有公式简练、计算结果准确、使用方便、稳定性好等优点，因此在等距情形宜采用龙贝格求积公式龙贝格算法设计目的设计目的Romberg)积分法。 2.程序设计背景： 对于较大积分区间、复杂被积函数、较高精度要求的数值积分问题，需要较多的求积节点。如果采用高阶插值型求积公式，当被积函数f(x)不是多项式函数时，求积过程可能不稳定，因此这时只能采用复化求积。而龙贝格求积是对逐次分半梯形公式求积公式加速的一种外推方法，收敛速度较复化求积更快，也是一种实用的数值积分方法。 在一元函数的积分学中,我们已经熟知,若函数f(x)在区间［a, b］上连续且其原?函数为F(x) ,则可用牛顿―莱布尼兹公式来求定积分。牛顿―莱布尼兹公式虽然在理论上或在解决实际问题中都起了很大的作用,?但它并不能完全解决定积分的计算问题。因为定积分的计算常常会碰到以下三种情况： (1)被积函数f(x)的原函数F(x)不易找到。许多很简单的函?数,例 等,其原函数都不能用初等函数表示成有限形式。 (2)被积函数f(x)没有具体的解析表达式。其函数关系由表格或图形表示,无法求出原函数。 (3)尽管f(x)的原函数能表示成有限形式但其表达式相当复杂。 另外，许多实际问题中的被积函数往往是列表函数或其他形式的非连续函数，对这类函数的定积分，也不能用方法求解。由于以上原因，数值积分的理论与方法一直是计算数学研究的基本课题。10-10 ,并输出积分区间的分割数。 2）用Romberg 积分法计算不同积分，要求误差不超过10-10 ，并与Simpson 方法比较计算量。 4．程序设计意义： 数值计算的广泛性和深入性已成为现代科技发展的重要特征，数值计算和理 论研究、科学实验并列为科学研究的三大支柱。用数学方法、将计算机作为工具、提高解决各种实际问题的能力是现代社会的需要，也是知识更新的必要环节。数值积分方法作为数值计算方法中的一个重要部分，如果将其与c++等语言中程序设计相结合。用正确的数值计算方法编写相应的c++程序。我们就能很好地解决一些复杂的实际问题，做到既有速度又有精度。 5．所做工作： 用VC++编程实现复化梯形法的递推算法和龙贝格(Romberg)积分法，并控制了结果的精度，比较辛普森和龙贝格积分的计算量，使其在要求的误差范围之内。 二、理论分析： 1、问题分析： 因原问题要求误差不超过10-10 ，需把积分区间[a,b]分成若干个小区间，使被积函数f(x)在每一个小区间内足够平滑，然后用低阶数值积分公式求f(x)在每一个小区间内的积分，再对所有小区间的积分结果求和，就得到了f(x)在[a,b]上的积分。 2、理论依据及求解对策： 1）把积分区间[a,b]等分为m个小区间，令步长h=，求积节i=0，,1,2…m，等距节点复化梯形公式为 它的余项函数为 2）用Tm表示把［a,b］作m等分并按复合梯形公式求积的结果,将每一小段再对分,令新的小段的长h′=h/2,则T2m与Tm之间有如下关系: 其中 三、方法详解： 1、推导： 1）用梯形公式求f(x)在每个小区间上积分余项之和： 2）等距节点复化梯形公式的余项为每个小子区间上积分余项之和： 3）设中最小的1项为,最大的1项为。 由连续函数介值定理可知，在之间至少存在1点，使。 又因为 所以 所以 2．分析及程序框图： 1）x为1时，函数值为1 2)龙贝格中根据T1、T2和m计算下一个数并返回 3）龙贝格积分法求解过程 4）用复化梯形法的递推