§4.2Newton-Cotes公式.ppt

第二节 Newton-Cotes公式 * ? 2009, Henan Polytechnic University * §2 Newton-Cotes公式 第四章 数值积分与数值微分 * 4.2.1 牛顿—柯特斯(Newton-Cotes)求积公式 中,当所取节点是等距时称为牛顿-柯特斯公式。 在插值求积公式 其中 插值多项式 求积系数 公式的推导 　设将积分区间[a,b]n等分，求积节点为 那么， 令x=a+th，则t=(x-a)/h，且由　　　　 可知　　　　　． 所以 记： Cotes系数 且 则有： Newton-Cotes求积公式 Cotes系数性质 几种常用的Newton-Cotes求积公式 梯形公式，辛普生公式，Cotes公式 1. n=1时的梯形求积公式 按Newton-Cotes系数公式计算得 　　　　　　　 a b a b 故求积系数A0, A1为 记 求积公式为 -----梯形求积公式 容易验证梯形公式的代数精确度的次数为1. 在[a, b]上积分，可得 考虑梯形求积公式的误差估计 假定 用推广的积分中值定理，将过(a, f(a)), (b, f(b)) 点的线性插值的余项　 2. n=2时的Simpson（抛物线）求积公式 按Newton-Cotes系数公式可以计算出　　　　　　　　 a b 容易验证Simpson求积公式具有３次的代数精确度． 余项公式为： 所以 上述公式称为Simpson求积公式。 3. n=4时的Cotes求积公式 按Newton-Cotes系数公式可以计算出 由此可得Cotes求积公式： a b 余项公式为： n阶Newton-Cotes求积公式当n为偶数时代数精度为n+1 x 0 0.25 0.5 0.75 1 f(x) 1 0.9896158 0.958851 0.9088516 0.8414709 例：分别用梯形公式，辛普生公式和柯特斯公式计算 准确值为：0.9460831 解： 利用梯形公式可得： x 0 0.25 0.5 0.75 1 f(x) 1 0.9896158 0.958851 0.9088516 0.8414709 利用辛普生公式得： 利用柯特斯公式得： 例 : 用辛普森公式和柯特斯公式计算定积分 的近似值,并估计其误差(计算结果取5位小数) 解: 辛普森公式 由辛普森公式余项 知其误差为 由于 柯特斯公式 误差为 该定积分的准确值 这个例子告诉我们，对于同一个积分，当n≥2时，公式却是精确的，这是由于辛普森公式具有三次代数精度，柯特斯公式具有五次代数精度，它们对被积函数为三次多项式当然是精确成立的。 误差与区间长度有关，区间长度越长，误差越大。 利用积分区间可加性，将较长区间分成若干个小区间，在每个小区间上分别应用Newton-Cotes求积公式，再相加，即可得到复化的积分公式。 4.2.2 复化求积公式及其收敛性 常用复化求积公式 复化梯形公式 2. 复化辛普生公式 3. 复化柯特斯公式 1.复化梯形公式 复化梯形公式 复化梯形公式（单击播放） 2.复化辛普生（抛物线）公式 复化辛普生（抛物线）公式