定积分近似计算以及误差估计.doc

定积分的近似计算方法 作者: 操乐青 指导老师: 邢抱花 摘要 关键词 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 几种近似计算定积分方法的比较分析及误差估计 例计算积分精确到0.001.(一) 利用矩形公式计算, 因为对于,有(如果12),所以按照公式 . 0. 如果取=10,则我们公式的余项的余数得,我们还必须加进由于在计算函数值实行四舍五入所产生的误差的界限相差于0.16,为了这个目的只要计算的值到四位小数精确到0.00005就够了.我们有 和6.9284 故计算结果为 。 方法(二) 按照梯形公式作同样的计算,在这种情况下,作公式 在这儿也试一试取=10,虽然此时仅可以证,纵坐标是 y和 故计算结果为 方法(三） 用公式做同样的计算 作公式 并且=5时有.实行计算到五位数字,精确到0.000005 . 将区间等分,节点为 (步长),)在每个小区间上采用梯形公式得