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

《数 值 方 法》 课程设计 内容提要： 实际问题中常常需要计算积分，有些数值方法，如微分方程和积分方程的求解，也都和积分计算相联系。依据微积分基本定理，对于积分只要找到被积函数f(x)的原函数F(x)。F’(x)=f(x)，便有牛顿―莱布尼兹公式来求得积分．但实际运用这种方法往往有困难，因为大量的被积函数，找不到用初等函数表示的原函数；另外当f(x)是由实验测量或数值计算给出的一张数据表时，牛顿―莱布尼兹公式也不能直接运用，因此，实际往往通过其他的一些方法如利用复化梯形公式、龙贝格公式、牛顿―莱布尼兹公式等来求解。 复化梯形求积法通常把积分区间等分成若干个子区间，在每一个子区间上用低阶的求积公式，对所有的子区间求和即得整个区间 [a,b] 上的积分公式。 龙贝格算法是在区间逐次分半过程中，对梯形公式的近似值进行加权平均获得准确程度较高的积分近似值的一种方法，它具有公式简练、计算结果准确、使用方便、稳定性好等优点，因此在等距情形宜采用龙贝格求积公式。龙贝格算法也称为逐次分半加速法，它是在梯形公式、辛普森公式和牛顿―莱布尼兹公式之间的关系的基础上．构造出一种加速计算积分的方法。作为一种外推算法，它在不增加计算量的前提下提高了误差的精度。在等距基点的情况下，用计算机计算积分值通常都采用把区间逐次分半的方法进行。这样，前一次分割得到的函数值在分半以后仍可被利用．并且易于编程。 目 录 一、程序设计主要内容 ……………………………………………………3 1. 程序设计目的…………………………………………………………3 2.程序设计背景 …………………………………………………………3 3．程序设计要求…………………………………………………………3 4．程序设计意义…………………………………………………………4 5．所做工作：……………………………………………………………4 二、理论分析…………………………………………………………………4 1、问题分析………………………………………………………………4 2、理论依据及求解对策…………………………………………………4 三、方法详解……………………………………………………………………5 1、推导 …………………………………………………………………5 2．分析及程序框图 ……………………………………………………5 四、问题解决…………………………………………………………………9 五、结果分析 ………………………………………………………………13 六、体会 ……………………………………………………………………14 七、参考文献 …………………………………………………………………14 一．程序设计主要内容 1. 程序设计目的： 1）学会用数值积分避开求f(x)的原函F(x)的繁琐步骤，并可以有效的控制结果，使其在要求的误差范围之内。 2）在某些求积函数中，用数值积分求解一些原函数F(x)不能用初等函数表示成有限形式 。 3）熟练掌握用复化梯形法和龙贝格法求解积分。 4）编程实现复化梯形法的递推算法。 5）编程实现龙贝格(Romberg)积分法。 2.程序设计背景： 对于较大积分区间、复杂被积函数、较高精度要求的数值积分问题，需要较多的求积节点。如果采用高阶插值型求积公式，当被积函数f(x)不是多项式函数时，求积过程可能不稳定，因此这时只能采用复化求积。而龙贝格求积是对逐次分半梯形公式求积公式加速的一种外推方法，收敛速度较复化求积更快，也是一种实用的数值积分方法。 在一元函数的积分学中,我们已经熟知,若函数f(x)在区间［a, b］上连续且其原?函数为F(x) ,则可用牛顿―莱布尼兹公式来求定积分。牛顿―莱布尼兹公式虽然在理论上或在解决实际问题中都起了很大的作用,?但它并不能完全解决定积分的计算问题。因为定积分的计算常常会碰到以下三种情况： (1)被积函数f(x)的原函数F(x)不易找到。许多很简单的函?数,例 等,其原函数都不能用初等函数表示成有限形式。 (2)被积函数f(x)没有具体的解析表达式。其函数关系由表格或图形表示,无法求出原函数。 (3)尽管f(x)的原函数能表示成有限形式但其表达式相当复杂。 另外，许多实际问题中的被积函数往往是列表函数或其他形式的非连续函数，对这类函数的定积分，也不能用不定积分方法求解。由于以上原因，数值积分的理论与方法一直是计算数学研究的基本课题。 3．程序设计要求： 1）用复化梯形法的递推算法计算不同积分，要求误差不超过10-10 ,并输出积分区间的分割数。 2）用Romberg 积分法计算不同积分，要求误差不超过10-10 ，并与Simpson 方法比较计算量。 4．程序设计意义： 数值计算的广泛性和深入性已成为现代