计算方法讲义课件 七 数值积分.docVIP

第七章 数值积分 如果函数f(x)在区间[a,b]上连续，且原函数为F(x)，则可用牛顿―莱布尼兹公式：来求得定积分。然而很多函数无法用牛顿―莱布尼兹公式求积分。 一个简单被积函数，例如 QUOTE ，其不定积分可能很复杂，见下面的MATLAB实例： syms a b c x int(sqrt(a+b*x+c*x*x),x) ans=1/4*(2*c*x+b)/c*(a+b*x+c*x^2)^(1/2)+1/2/c^(1/2)*log((1/2*b+c*x)/c^(1/2)+(a+b*x+c*x^2)^(1/2))*a-1/8/c^(3/2)*log((1/2*b+c*x)/c^(1/2)+(a+b*x+c*x^2)^(1/2))*b^2 所以有必要研究简单、高效的计算定积分的方法（即数值积分方法）。数值积分的基本思想是构造一个简单函数Pn(x)来近似代替被积分函数f (x)，然后通过求得的近似值。 7.1 插值型求积公式 设，插值型求积公式就是构造插值多项式Pn(x)，使。 构造以a，b为结点的线性插值多项式，称为梯形公式。 以a, ，b为三个插值节点，构造二次插值多项式，则可以推出，，，。由此得公式：，称为辛卜生（Sinpson）求积公式。 根据经典拉格朗日插值公式，代入求定积分则有，引入记号，，?k为求积系数，xk为求积节点。注意：一积分结果为函数值的一个代数和，二是。 如果积分区间比较大，直接使用上述求积公式精度难以保证。可对f (x)用分段抛物插值。通常采取的办法是复化求积方法： （1）等分求积区间，比如取步长，分[a, b]为n等分，分点为，k = 0, 1, 2,…, n。 （2）在区间 [xk, xk+1]上使用以上求积公式求得Ik。 （3）取和值作为整个区间上的积分值。 将梯型公式和辛卜生公式应用于各子区间上得到子区间的定积分，再将子区间的定积分加起来得到整个区间的定积分近似值，相关公式称为复化梯型公式和复化辛卜生公式。相对于复化梯型公式，复化辛卜生公式是一种精度较高的求积公式。例如对于复化梯型公式，令，则 xk 0 1/8 1/4 3/8 1/2 5/8 3/4 7/8 1 f (xk) 4 3.93846 3.7647 3.50685 3.2 2.8764 2.46 2.26549 2 例 利用数据表计算积分 ()。解：取n = 8用复化梯形公式： 7.2 变步长梯形方法 使用复化求积公式须给出合适的步长，步长太大精度难保证，步长太小会增加计算量，事先给出一个合适的步长是十分困难的。递推公式避免了老节点的重复计算，使计算量减少了一半。 变步长积分法思想是将区间逐次对分，比较前后两次计算结果，若满足精度要求就停止，否则再次对分，直到到达精度要求为止。设将区间[a, b] n等分，共有n+1个分点，按复化梯形公式计算Tn，需要计算n+1个f (x)的值。T2n 的全部分点中有n+1个是原有的点。小区间[xk, xk+1]经过二分增加分点后，用复化梯形公式得积分为：，因此有 例 计算 QUOTE 。 解：根据梯形公式和复化梯形公式，，于是有 n 1 2 4 8 16 32 Tn 0.9397 0.9445 0.9456 0.9459 0.9461 0.9461 7.3 求积公式的误差 Pn (x)是f (x)的n次插值多项式，当本身就是次数不超过n的多项式时，求积公式是精确的。由于，若f (xk)的舍入误差小于? ，则。所以舍入误差对数值积分的影响不大。 应用插值多项式余项定理，对于插值多项式次数为1的情况有：，可以证明梯形公式的截断误差为：（注意：这里用到了积分中值定理：设在区间[a,b]上连续，在[a,b]区间上可积且不变号，则在[a,b]区间上至少有一个满足）。 将[a, b]区间n等分，取考虑复化梯形积分公式的误差，这个误差是n个等分区间段上得误差之和，即（注意：这里用到了介值定理，即对于连续函数和自然数n，存在使）。 7.4 收敛条件及收敛加速 梯形法简单，但精度低，收敛的速度慢。如何提高收敛速度？ 设I是精确积分值，根据复化梯形公式的余项表达式可知：，。假定, 则有。整理得：。可见只要二分前后Tn与T2n 相当接近，就可以保证T2n的误差很小。T2n的误差大致等于，用误差值作为T2n的补偿，可期望所得到的，可能是更好的结果。 也可以这样考虑，将所有Tn看做构成一个函数T，变量是h2。当h趋近0时，T(h2)接近I*，即T(0)=I*，连接点和得一直线，其方程为：，延伸该直线与Y轴相交（见下图），。这就是一种迭代的加速。 图 一种迭代加速 7.5 高斯型求积公式 在插值型求积公式中，插值节点是事先固定的，有时还进一步限定是等距的。是否可以在[a, b]上