《地球物理计算方法》第2章 数值积分——Newton-cotes积分公式.pptxVIP

《地球物理计算方法》

第2章数值积分——Newton-cotes积分公式

若此求积公式至少具有n阶的代数精度，则称此求积公式为n阶的Newton-Cotes公式．把区间[a，b]分为n等份，步长为h（h=(b－a)/n）,则n+1个点（取等分点）分别为：xk=a+kh;,由这n＋1个点构造的插值型求积公式为：

3积分公式的构造：精度概念定义对于则此求积公式对于应成立等式条件满足。

4（2）插值公式方法一般多项式积分

5既然柯特斯系数跟被积函数、积分区间无关，只与代替被积函数的多项式次数n有关，就可以在求积分前预先算出来。只要计算出柯特斯系数，就可以按下式计算积分（如n=4)

不同n的科特斯系数表及代数精度

基本性质 当时，其系数为正；当时，其系数有正、有负．1、2、3、对称性；

1、梯形公式牛顿-科特斯公式当n＝1时（即2个点，1等份），有梯形公式（1次代数精度）：2、辛普生公式当n＝2时（即3个点，2等份），有辛普生公式（3次代数精度）：

3.柯特斯公式当n＝4时（即5个点，4等份），有柯特斯公式（5次代数精度）

10低阶求积公式的代数精度所以，n为偶数时的积分公式更加常用。

11低阶求积公式的余项牛顿-柯特斯公式是一个插值型数值求积公式，当用插值多项式pn(x)代替f(x)积分时(代数精度为n)，其截断误差R(f)为：便可得出它的截断误差：其中