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

定积分的近似计算方法与误差估计 作者 : 操乐青 指导老师 : 邢抱花 摘要  本文主要讨论了一元函数常见的数值积分方法  , 例如插值型求积公式，高斯求积公式等近 似计算方法  , 在用这些方法计算定积分时  , 会产生一些误差  , 为了减少误差  , 可以利用复化求积公 式、复化高斯公式等  . 本文围绕这些方法  , 系统介绍它们的计算公式以及截断误差  , 并用例题分析它 们产生误差的大小、计算量等  . 关键词  插值型积分  高斯积分  误差分析 近似计算 引言 在计算定积分的值 I b f ( x) 的一个原函数 f (x)dx 时 ,常常根据微积分学基本定理求出 a F (x) ,再用牛顿 -莱布尼茨公式求得积分 , I b F (b) F (a) . 但这种方法只限于解 f (x)dx a 决一小部分定积分的求值问题 .当函数没有具体表达式 ,只是一些实验测得数据形成的表格或 图形或者是 F (x) 无法用初等函数表示 b ex 2 b sin x2 dx 等等 ,这就需要我们用一些近 ,例如 , a dx , a 似方法来求积分值 . 与数值积分一样 , 把积分区间细分 ,在每个小区间上 ,找到简单函数 ( x) 来近似代替 f ( x) , b b b ( x)dx . 且 ( x)dx 的值容易求的 .这样就把计算复杂的 f ( x) dx 转化为求简单的积分值 a a a 因此 ,定积分的近似计算实质上就是被积函数的近似计算问题 . 定积分的近似计算——常见数值方法 2.1 矩形公式 根据定积分的定义，每一个积分和都可以看作是定积分的一个近似值，即 b n f (x)dx f ( i ) xi a 1 i 在几何意义上，这是用一系列小矩形面积近似小曲边梯形的结果，所以把这个近似计算方法称为矩形法．不过，只有当积分区间被分割得很细时，矩形法才有一定的精确度． 针对不同 i 的取法，计算结果会有不同，常见的取法有 ： （ 1）左端点法，即 xi 1 , a n f ( x) dx f ( xi 1 ) xi i b i 1 xi , a n （2）右端点法，即 f ( x)dx f ( xi ) xi i b i 1 xi 1 xi a n xi 1xi （ 3）中点法，即 i 2 , f (x)dx f ( 2 ) xi b i 1 例 1 用矩形公式近似计算积分 1 dx （取 n 100 ）. 0 1 x2 解 对 0,1 作 n 等分 x0 a x1 xi a b a i xn b ，由定义知： n 1 dx n f ( i ) xi 1 n f ( i ) 0 1 x i 1 n 1 1 (1) 左点法：在区间 [ xi 1 , xi ] 上取左端点，即取 i xi 1 ， i 1,2 n 1 dx n f ( i ) xi 0.78789399673078 ， 0 1 x2 i 1 理论值 1 dx ，此时计算的相对误差 x2 0 1 4 0.78789399673078 4 0.003178 4 （ 2）右点法：在区间 [ xi 1 , xi ] 上取右端点，即取 i xi ， i 1,2 n 1 dx n f ( i ) xi 0.78289399673078 ， 0 1 x2 i 1 理论值 1 dx ，此时计算的相对误差 0 1 x2 4 0.78289399673078 4 0.003188 4 （ 3）中点法：取 在区间 [ xi 1 , xi ] 上取中点，即取 i xi 1 xi ， i 1,2 n 2 1 dx n f ( i ) xi 0.78540024673078 ， 0 1 x2 i 1 理论值 1 dx ，此时计算的相对误差 0 1 x2 4 0.78540024673078 4 2.653 10 6 4 如果在分割的每个小区间上采用一次或二次多项式来近似代替被积函数，那么可以期望得到比矩形法效果好得多的近似计算公式．下面介绍的梯形法和抛物线法就是这一指导思想的产物． 2.2 梯形公式 等分区间 x0 a x1 xi b a i xn b ， b a a x n n 相应函数值为 y0 , y1 , , yn （ yi f (xi ), i 0,1, , n）． 曲线 y f ( x) 上相应的点为 P0 , P1, , Pn （ Pi ( xi , yi ), i 0,1, , n ） 将曲线的每一段弧 Pi 1 Pi 用过点 Pi 1 ，Pi 的弦 Pi 1 Pi（线性函数） 来代替， 这使得每个 [ xi 1