计算方法(5-1).pptVIP

由于在每个小区间 上， 是用一次多项式去近似的，由插入法可知： 其中 记 为每个小区间 上用小梯形面积近似代替小曲边梯形面积，所产生的截断误差。 由积分中值定理可知，如果 在区间 上连续，则在积分区间 内存在一点 ，使得 成立。 所以 所以，在区间 上总的截断误差为： 又因为 在区间 上连续， 所以，在区间 内必存在一点 ，使得： 这就是复化梯形公式的截断误差公式。 假设 则： 5.3.2 复化Simpson公式 如果用分段二次插值函数近似被积函数，即在小区间上用Simpson公式计算积分近似值，就导出复化Simpson公式。 整理后得到： 这就是复化Simpson公式。 若 在区间 上有四阶连续导数，那么辛甫生公式的截断误差为： 因为f (4)(x) 连续，故存在??(a, b)，使得： 上式表明，步长h越小，截断误差越小。与复化梯形公式的分析相类似，可以证明，当n ? ?时，用复化Simpson公式所求得的近似值收敛于积分值，而且算法具有数值稳定性。 将区间 分成4等份,节点 ，这样就得到五个数据点： 通过这五个数据点，就可以做一个四次插值多项式 ，那么类似于辛甫生公式的推导，就可以得到： 5.3.3 复化柯特斯公式 为了提高求积精度将区间[a, b]分成 n 等分，分点为： 在每个小区间： 上，以分为4等份，内分点依次为： 在每个子区间上应用柯特斯公式，然后相加，便得到复化的柯特斯公式： 柯特斯公式的截断误差： 梯形公式： 具有1次代数精度。 辛甫生公式： 具有3次代数精度。 柯特斯公式： 具有5次代数精度。 当然，也可以由三个公式的截断误差，推出三个公式的所具有的代数精度。 辛甫生公式的插值节点只比梯形公式多一个，但其代数精度却比梯形公式高2，它们都是最为常用的数值积分公式，尤其是辛甫生公式，逻辑结构简单，且精度又比较高。 例 6: 根据函数表 k xk f (xk)=Sin xk/ xk k xk f (xk)=Sin xk/ xk 0 0 1.000000 5 0.625 0.936156 1 0.125 0.997398 6 0.75 0.9088522 0.25 0.989616 7 0.825 0.8771933 0.375 0.976727 8 1 0.8414714 0.5 0.958851 [解]：（1）由复化梯形公式,n=8, h=1/8： （2）由复化Simpson 公式,n=4,h=1/4： 与准确值I=0.9460831…比较，显然用复化Simpson公式计算精度较高。 事实上，由误差公式可知：RT (f )=O(h2), RS (f )=O(h4)，故当h比较小时，用复化Simpson公式计算误差较小。 由误差估计公式不仅可以计算所求近似值的误差，反之，亦可由给定的精度估计应取多大步长。 5.4逐次分半算法(变步长方法) 当实际应用计算机进行数值积分计算时，一般先提出一定的精度要求，然后根据误差估计式对于预先给定的精度给出步长h或n, 但由于误差估计式中要估计高阶导数，而这一点往往很困难。 因此实际计算时，常采用变步长方法：逐步缩小步长，每次将步长缩小一半，或者说逐次等分区间,反复利用复化求积公式，直到相邻两次计算结果相差不大为止或者满足给定精度为止。 5.4.1 梯形法的递推公式 因此计算梯形序列{T2m}可按： 4. 设将区间[a , b] n 等分，共有n+1个分点，如果将积分区间再等分一次，则分点增为2n+1个，将等分前后两个积分值联系起来加以考察： 注意到每个子区间 经过等分只增加了一个分点： 用复化梯形公