第7章数值微积分2概述.ppt

第7章 数值微积分 基础教学部数学教研室 彭 晓 华 立体化教学资源系列——数值分析 复习： 牛顿-柯特斯公式 其中 称为柯特斯系数， 当 相应的求积公式称为梯形公式 当 时相应的求积公式称为辛普森公式 时称为柯特斯公式 当 其中， 截断误差 牛顿-柯特斯公式的截断误差（即余项）为 当 时，梯形公式的截断误差为 如果 连续， 辛普森公式的截断误差为 柯特斯公式的截断误差为 为了提高求积公式的精度，通常可把积分区间分成若干 子区间，再在每个子区间上用低阶求积公式，这种方法 称为复化求积方法(Composite Numerical Integration). 复化求积方法的实质是利用分段低次插值多项式的积分， 的积分.本节讨论复化梯形公式、复化辛 7.3 复化求积公式 去逼近 普森公式和复化柯特斯公式. 将 区间 等分，步长 ，分点 ，则 等分为 个子区间，即 ，子区间 7.3.1 复化梯形公式 上利用梯形公式，则 从而 记 为复化梯形公式(Composite Trapzoidal Rule) （7.12） 在每个子区间 7.3.2 复化辛普森公式 上利用辛普森公式，则 其中， 是区间 的中点.从而 记 为复化辛普森公式(Composite Simpson’s Rule) （7.13） 在 7.3.3 复化柯特斯公式 上取 的柯特斯公式，则 其中， 因此 在 记 为复化柯特斯公式 （7.14） 7.3.4 复化求积公式的截断误差 ，则复化梯形公式的截断误差 （7.15） 证明 因为在 上梯形公式的截断误差为 所以 在 上连续，由连续函数性质知，在 中存在点 ，使 设 其中 因此 证毕. （7.17） （7.16） 类似于复化梯形公式截断误差的推导，可得复化 辛普森公式的截断误差 复化柯特斯公式的截断误差 【注】 复化梯形公式、复化辛普森公式和复化柯特斯公 式都是有效的求积计算公式. 实际上，只要 ，则可得到收敛性，即 （2） 稳定性.由于 和 所以由定理2可知复化梯形公式、复化辛普森公式和复化柯 特斯公式均是稳定的求积公式. 柯特斯公式的误差阶分别为 收敛性都是显然的. 的求积系数均为正数， （1） 收敛性.由截断误差公式（7.15）、（7.16）、 （7.17）可知，复化梯形公式、复化辛普森公式和复化 例2 计算积分 解 1） 将积分区间[0，1]八等分，分点及分点处的函 数值见表7-3，用复化梯形公式计算，得 . 0 1/8 1/4 3/8 1/2 4 3 3 3 3 5/8 3/8 7/8 1 2 2 2 2表7-3 例2数据 2） 将积分区间[0，1]四等分，用复化辛普森公式计算，得 两种复化方法都用到表7-3中九个点上的函数值，它 们的计算工作量基本上相同，但所得结果与积分真值 相比较，复化辛普森所得近似值 远比复化 要精确.因此，在实际计算时， 梯形所得近似值 较多地应用复化辛普森公式. 例3 分别用复化梯形公式和复化辛普森公式计算 使误差不超过 解 由复化梯形公式的截断误差，有 所以， 因此，复化梯形公式至少应取361个节点. 由复化辛普森公式的截断误差，有 所以， ， 因此，复化辛普森公式至少应取10个节点. 问各取多少个节点？ 龙贝格积分(Romberg Integration)以逐次对分区间和复化求积为出发点，逐次修正近似值的方法构造出一种新的求积算法，称为龙贝格算法. 7.4 龙贝格求积公式 7.4.1 梯形法的递推化 上取步长 ， 可取 等正整数.用复化梯形公式求积，可得 将每个子区间对分，即 区间 等分， ，每个子区间的新增分点为 ，则 在区间 即对 梯形公式 ， 由递推公式（7.18）可知，在已经算出 的基础上再计算 时，只要计算 这与直接利用复化梯形公式求 节省了一半. ，求得梯形法的递推计算，称为变步长 （7.18） 个区间上新增分点的函数值. 相比较，计算量几乎 例5 利用递推公式（7.18）重新计算积分 解 1） 首先对区间[0，1]使用梯形公式，得 2） 将区间二等分，新增分点 ，由递推公式 3） 再将各小区间二等分，新增分点 由递推公式（7.18）得 （7.18）得 0 1 2 3 4