计算方法牛顿-柯特斯求积公式与复合求积公式.ppt

《计算方法牛顿-柯特斯求积公式与复合求积公式.ppt》由会员分享，可在线阅读，更多相关《计算方法牛顿-柯特斯求积公式与复合求积公式.ppt（48页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第7次牛顿-柯特斯求积公式与复合求积公式计算方法计算方法(Numerical Analysis)1.牛顿柯特斯求积公式2.牛顿-科特斯求积公式的例子3.复合求积公式4.复合求积公式的例子附录：复合梯形公式与复合辛普生公式算法实现与流程图牛顿柯特斯求积公式采用等距节点的插值型求积公式4.2牛顿柯特斯求积公式是插值基函数。有关系式定义：在插值求积公式中,当所取节点时称为牛顿-柯特斯(Newton-Cotes)公式，其中：是等距利用等步长的特点计算积分系数Ak求积节点为：因此：将积分区间a,b划分为n等分,步长a=x0 x2xkx1xixn=b 可以推出：a=x0 x2xkx1xixn=b 作变量替

2、换并注意得：代入插值求积公式(4.1)有称为牛顿-柯特斯求积公式,Ck称为柯特斯系数引进记号(柯特斯系数)则1.将区间a,b分为n等分，则n+1个柯特斯系数之和为1证：由于插值型积分公式的系数Ak之和等于（b-a）柯特斯系数的性质由关系：得：2.Ck是不依赖于积分区间a,b以及被积函数f(x)的常数，只要给出n，就可以算出柯特斯系数。例如，当n=1时似曾相识当n=2时，由 P104表4-1给出了n从18的柯特斯系数。当n=8时，出现了负系数，从而影响稳定性和收敛性，因此，实用的只是低阶公式。似曾相识Newton-Cotes公式柯特斯系数列表：当n=8的时候，出现负值，不稳定对n=6,7,8的情

3、况，见教材。nCk11/21/221/62/31/631/83/83/81/847/9016/452/1516/457/90519/28825/9625/144 25/14425/9619/288几个重要的低阶求积公式在牛顿-柯特斯求积公式中n=1,2,4时，就分别得到下面的梯形公式、辛卜生公式和柯特斯公式。定理4.2（梯形公式的误差）设f(x)在a,b上具有连续的二阶导数，则梯形公式的误差（余项）为当b-a1时，误差较大；b-a2时，误差较大；b-a4时，误差较大；b-a4时，误差较小总结：Newton-Cotes公式给出了等距节点的插值型求积公式的统一计算公式。定义：在插值求积公式中,当所

4、取节点时称为牛顿-柯特斯公式：是等距k=0,nn=1，梯形公式；n=2,辛普生公式；n=4，牛顿-柯特斯公式.Home牛顿-柯特斯求积公式例题例4.11分别用梯形公式、辛卜生公式和柯特斯公式计算定积分的近似值.(1)用梯形公式计算0.51(2)用辛卜生公式0.510.75误差(3)用柯特斯公式计算(n=4,4等份，5个节点)，系数为 0.510.750.8750.625积分的准确值为可见，三个求积公式的精度逐渐提高。梯形公式辛卜生公式柯特斯公式精度0.4267770.430930.43096例4.12用辛卜生公式和柯特斯公式计算定积分的近似值,并估计其误差(计算结果取5位小数)解:辛卜生公式由

5、于f(x)是3阶多项式，所以辛卜生公式余项132解:柯特斯公式 知其误差为该定积分的准确值，此例说明，对于同一个积分，当n2时，两个公式都是精确的。原因：辛卜生公式具有3次代数精度，柯特斯公式具有5次代数精度，它们对被积函数为3次多项式当然是精确成立的。复合求积公式4.3复合求积公式一般地，应用牛顿-柯特斯求积公式求积分的近似解的时候，随着求积节点数的增多，对应公式的精度也会相应提高。但n8时，开始出现负值的柯特斯系数。因此，可能导致舍入误差增大，且往往难以估计。不能单纯用增加求积节点数（例如，8个节点，确定一个7次多项式来近似被积函数）的方法来提高计算精度。新想法：将积分区间分成若干个小区间

6、，在每个小区间上采用低阶求积公式（低阶多项式），然后把所有小区间上的计算结果整合起来，得到整个区间上的求积公式。此即复合求积公式的基本思想。4.3.1复合梯形公式及其误差上应用梯形公式，得：将积分区间a,b划分为n等分,步长为求积节点为abxkxk+1在每个小区间y=f(x)abx1x2x3x4x5x6使用复合梯形公式计算xy然后将Ik累加求和,用作为所求积分I的近似值。得到复合梯形公式：(4.5)设f”(x)在a,b上连续,根据上次课关于梯形公式误差的讨论，知在上梯形公式的余项为：复合梯形公式的误差在a,b上的余项为(*)因为在a,b上连续，根据连续函数的介值定理知，存在，使因此，注意到h=

7、(b-a)/n，得余项 意义：当h0的时候，复合梯形求积公式收敛于f(x)在a,b上的定积分值，收敛阶为O(h2)。可以调整h的值，使得计算的结果满足预先定义的精度。有界4.3.2复合辛卜生公式及其误差合并xk和xk+1项将积分区间a,b划分为n等分,记子区间的中点为,在每个小区间上应用辛卜生公式，则有y=f(x)x1x2x3xyx0类似于复合梯形公式余项的讨论，设在a,b上连续，则复合辛卜生公式(4.6)的求积余项为：得到复合辛卜生公式：(4.6)复合辛卜生公式的误差O(h4)有界若将每个子区间4等分:xkxk+1x内分点记为:复合柯特斯公式:在区间，柯特斯求积公式为：复合柯特斯公式的误差：

8、O(h6)有界类似于复合梯形公式余项的讨论，设在a,b上连续，则复合辛卜生公式(4.6)的求积余项为：得复合柯特斯公式:复合求积公式的余项表明，只要被积函数f(x)所涉及的各阶导数在a,b上连续，那么复合梯形公式、复化辛卜生公式与复化柯特斯公式所得近似值的余项和步长的关系依次为因此当h0(即n)时,都收敛于积分真值，且收敛阶一个比一个高，收敛速度一个比一个快。复合求积公式的例子例4.13依次用n=8的复合梯形公式、n=4的复合辛卜生公式计算定积分解1:由复合梯形公式计算，将区间0,18等分，即n=8,h=0.125 01使用计算器计算f(x)在各个节点上的值（积分准确值I=0.9460831）

9、解2：由复合辛卜生公式计算，将区间0,14等分，即n=4，h=0.25.01评：这两种方法计算量基本相同，然而精度却差别较大，同积分的准确值比较，却有六位有效数字。而利用复化辛卜生法计算结果积分准确值：I=0.9460831利用复合梯形法计算结果只有两位有效数字。例4.14用复合梯形公式计算定积分解:,则,又区间长度b-a=1，对复合梯形公式有余项即,n213，取n=213，即将区间0,1分为213等份时，用复合梯形公式计算误差不超过 问区间0,1应分多少等份，才能使误差不超过Home同学们计算作业作业P135:3，4题P136:6题Home复合梯形求积算法实现（1）复合梯形公式计算步骤确定步长h=(b-a)/N(N为等分数)对k=1,2,N，计算T=T+f(a+kh)T=h f(a)+2T+f(b)/2附录：复合梯形公式与复合辛普生公式算法实现与流程图（2 2）复复合合梯梯形形公公式式的的流流程程图图（1）复合辛卜生公式计算步骤确定步长h=(b-a)/N,S1=f(a+h/2),S2=0(N为等分数)对k=1,2,N-1，计算S1=S1+f(a+kh+h/2),S2=S2+f(a+kh)S=h f(a)+4S1+2S2+f(b)/6复合辛卜生求积算法实现（2 2）复复合合辛辛卜卜生生公公式式流流程程图图