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

《地球物理计算方法》

第2章数值积分—Newton-cotes积分公式若此求积公式至少具有n阶的代数精度，则称此求积公式为n阶的Newton-Cotes公式．把区间[a，b]分为n等份，步长为h（h=(b－a)/n）,则n+1个点（取等分点）分别为：xk=a+kh;,由这n＋1个点构造的插值型求积公式为：*积分公式的构造：精度概念定义对于则此求积公式对于应成立等式条件满足。1、公式的导出*（2）插值公式方法一般多项式积分*既然柯特斯系数跟被积函数、积分区间无关，只与代替被积函数的多项式次数n有关，就可以在求积分前预先算出来。只要计算出柯特斯系数，就可以按下式计算积分（如n=4)不同n的科特斯系数表及代数精度基本性质 当时，其系数为正；当时，其系数有正、有负．1、2、3、对称性；1、梯形公式牛顿-科特斯公式当n＝1时（即2个点，1等份），有梯形公式（1次代数精度）：2、辛普生公式当n＝2时（即3个点，2等份），有辛普生公式（3次代数精度）：3.柯特斯公式其代数精度至少是4;当n＝4时（即5个点，4等份），有柯特斯公式（5次代数精度）例:分别用梯形公式、辛卜生公式和柯特斯公式计算定积分的近似值(计算结果取5位有效数字)(1)用梯形公