定积分的近似计算以及类误差估计.docVIP

安庆师范学院数学与计算科学学院2014届毕业论文 第 第 PAGE 5页 共11页 定积分的近似计算方法与误差估计 作者: 操乐青 指导老师: 邢抱花 摘要 本文主要讨论了一元函数常见的数值积分方法,例如插值型求积公式，高斯求积公式等近似计算方法,在用这些方法计算定积分时,会产生一些误差,为了减少误差, 可以利用复化求积公式、复化高斯公式等.本文围绕这些方法,系统介绍它们的计算公式以及截断误差,并用例题分析它们产生误差的大小、计算量等. 关键词 插值型积分 高斯积分 误差分析 近似计算 1引言 在计算定积分的值时,常常根据微积分学基本定理求出的一个原函数,再用牛顿-莱布尼茨公式求得积分,.但这种方法只限于解决一小部分定积分的求值问题.当函数没有具体表达式,只是一些实验测得数据形成的表格或图形或者是无法用初等函数表示,例如,,等等,这就需要我们用一些近似方法来求积分值. 与数值积分一样,把积分区间细分,在每个小区间上,找到简单函数来近似代替,且的值容易求的.这样就把计算复杂的转化为求简单的积分值.因此,定积分的近似计算实质上就是被积函数的近似计算问题. 2 定积分的近似计算——常见数值方法 2.1 矩形公式 根据定积分的定义，每一个积分和都可以看作是定积分的一个近似值，即 在几何意义上，这是用一系列小矩形面积近似小曲边梯形的结果，所以把这个近似计算方法称为矩形法．不过，只有当积分区间被分割得很细时，矩形法才有一定的精确度． 针对不同的取法，计算结果会有不同，常见的取法有： （1）左端点法，即, （2）右端点法，即 , （3）中点法，即, 例1 用矩形公式近似计算积分（取）. 解 对作n等分 ，由定义知： (1)左点法：在区间上取左端点，即取， 0.78789399673078， 理论值，此时计算的相对误差 （2）右点法：在区间上取右端点，即取， 0.78289399673078， 理论值，此时计算的相对误差 （3）中点法：取 在区间上取中点，即取， 0.78540024673078， 理论值，此时计算的相对误差 如果在分割的每个小区间上采用一次或二次多项式来近似代替被积函数，那么可以期望得到比矩形法效果好得多的近似计算公式．下面介绍的梯形法和抛物线法就是这一指导思想的产物． 2.2梯形公式 等分区间 ， 相应函数值为 （）． 曲线上相应的点为 （） 将曲线的每一段弧用过点，的弦（线性函数）来代替，这使得每个上的曲边梯形成为真正的梯形，其面积为 ，． 于是各个小梯形面积之和就是曲边梯形面积的近似值， 即 称此式为梯形公式． 例2 用梯形公式近似计算定积分 （取）. 解 0.78539399673078， 理论值，此时计算的相对误差 很显然，这个误差要比简单的矩形左点法和右点法的计算误差小得多. 2.3 抛物线公式 由梯形法求近似值，当为凹曲线时，它就偏小；当为凸曲线时，它就偏大．若每段改用与它凸性相接近的抛物线来近似时，就可减少上述缺点，这就是抛物线法． 将积分区间作等分，分点依次为 ，， 对应函数值为 （）， 曲线上相应点为 （）． 现把区间上的曲线段用通过三点，，的抛物线 来近似代替，然后求函数从到的定积分： 由于，代入上式整理后得 同样也有 …… 将这个积分相加即得原来所要计算的定积分的近似值： ， 即 这就是抛物线公式，也称为辛卜生（Simpson）公式． 例3 用抛物线公式近似计算积分 (取). 解 =0.78539816339745， 理论值，此时计算的相对误差 2.4 几种近似计算定积分方法的比较分析及误差估计 例4 计算积分，精确到0.001. 解 方法(一) 利用矩形公式计算, 因为对于,有(如果12),所以按照公式 . 0. 如果取=10,则我们公式的余项的余数得,我们还必须加进由于在计算函数值实行四舍五入所产生的误差的界限相差于0.16,为了这个目的只要计算的值到四位小数精确到0.00005就够了.我们有