,单击此处编辑母版标题样式,单击此处编辑母版文本样式,*,牛顿,-,科特斯公式,等距节点的插值型求积公式称为,牛顿-科特斯公式,：,取等距节点：,x,i,=,a,+,i h,，,i=,1,2,n,令,x,=,a,+,t h,得：,插值型求积公式,其中,牛顿-科特斯公式 等距节点的插值型求积公式称为牛顿-科特斯公,1,牛顿,-,科特斯公式（续）,注：Cotes 系数,仅取决于,n,和,i,，,可通过查表得到。与被积函数,f,(,x,),及积分区间,a,b,均无关。,科特斯(,Cotes,),系数,牛顿-科特斯公式,：,牛顿-科特斯公式（续）注：Cotes 系数仅取决于 n 和,2,几个常见公式,n,=1:,代数精度,=,1,梯形求积公式,n,=2:,代数精度,=,3,抛物线求积公式,Simpson,求积公式,n,=4:,科特斯(Cotes)求积公式,T,S,C,几个常见公式n=1:代数精度=1梯形求积公式n=,3,科特斯系数表,科特斯系数表,4,系数特点和稳定性,科特斯系数具有以下特点：,(1),(2),(3),当,n,8,时，出现负数，,稳定性得不到保证,。而且当,n,较大时，由于Runge现象，,收敛性也无法保证,。,故,一般,不采用,高阶的牛顿-科特斯求积公式,。,当,n,7,时，牛顿,-,科特斯公式是稳定的。,系数特点和稳定性 科特斯系数具有以下特点：(1)(2)(,5,牛顿,-,科特斯公式的代数精度,定理,当,n,为偶数时，牛顿科特斯公式至少有,n,+1,阶代数精度。,证,：,只要证明当,n,为偶数时，公式对,f,(,x,),x,n,+1,精确成立。,由插值型求积公式的误差公式得,作变量代换,x,=,a,+,t h,，并将,x,i,=,a,+,i h,代入得,再作变量代换,t,=,n,-,s,，得,又,n,偶数,牛顿-科特斯公式的代数精度定理当 n 为偶数时，牛顿科特斯,6,余项,梯形公式的余项,中值定理,Simpson公式的余项,三次Hermite插值,余项 梯形公式的余项中值定理 Simpson公式的余项三次H,7,余项的一般形式,定理,(1)若,n,为偶数，,f,(,x,),C,n,+2,a,b,，则存在,(,a,b,),使得,设 ，则有,(2)若,n,为奇数，,f,(,x,),C,n,+1,a,b,，则存在,(,a,b,),使得,余项的一般形式定理(1)若 n 为偶数，f(x)C,8,举例（一）,例：,分别用梯形公式和simpson公式计算积分,解：,a,0,b,1,f,(,x,)=,e,-x,，,由,simpson,公式可,得,由,梯形公式可,得,与精确值,0.6321,相比,得误差分别为,0.0518,和,0.0002,。,举例（一）例：分别用梯形公式和simpson公式计算积分,9,复合求积公式,提高积分计算精度的常用两种方法,用,复合公式,用,非等距节点,复合求积公式：,将积分区间分割成多个小区间，然后在每个,小区间,上使用低次牛顿科特斯求积公式。,将,a,b,分成,n,等分,x,i,x,i,+1,，其中节点,(,i=,0,1,n,),复合求积公式 提高积分计算精度的常用两种方法 用 复合公式,10,复合梯形公式,复合梯形公式:,T,n,余项：,，,(,a,b,),复合梯形公式 复合梯形公式:Tn 余项：，(a,b),11,复合simpson公式,复合simpson公式:,S,n,余项：,，,(,a,b,),4,4,4,4,4,复合simpson公式 复合simpson公式:Sn 余项：,12,复合科特斯公式,复合cotes公式:,C,n,余项：,，,(,a,b,),复合科特斯公式 复合cotes公式:Cn 余项：，(a,13,举例（二）,解：,例：,设 ，利用下表中的数据分别用复合梯形公式和复合simpson公式计算积分,x,i,0,1/8,2/8,3/8,4/8,5/8,6/8,7/8,1.0,f,(,x,i,),1,0.997,0.990,0.977,0.954,0.936,0.909,0.877,0.841,举例（二）解：例：设,14,h,很小时的误差,i,(,x,i,x,i+,1,),(,h,0),定积分定义,即,同理,h 很小时的误差i(xi,xi+1)(h,15,收敛速度与误差估计,定义,若一个积分公式的误差满足 且,C,0，,则称该公式是,p,阶收敛,的。,例：,计算,解：,其中,=,3.1,38988494,其中,=,3.141592,502,运算量基本相同,收敛速度与误差估计定义 若一个积分公式的误差满足,16,Q:,给定精度,，如何取,n,?,例如：要求 ，如何判断,n,=？,?,上例中若要求 ，则,即：取,n,=409,Q:给定精度，如何取 n?例如：要求,17,通常采取将区间,不断对分,的方法，即取,n,=2,k,上例中2,k,409,k,=9,时，,T,512,=,3.141592,02,注意到区间再次对分时,可用来判断迭代,是否停止。,Q:,给定精度,，如何取,n,?,通常采取将区间不断对分的方法，即取 n=2k上例中2k,18,2.3,龙贝格算法,梯形法的递推化,龙贝格算法,理查森外推加速法,2.3 龙贝格算法梯形法的递推化,19,1,梯形法的递推化,方法思路:,复化求积方法可提高求积精度，实际计算,时可以将步长逐次分半。,在每个子区间x,k,x,k+1,经过二分只增加了一,个分点x,k+1/2,=1/2(x,k,+x,k+1,),用复化梯形公式,求得该子区间上的积分值为,1 梯形法的递推化方法思路:复化求积方法可提高求积精度，,20,注意，这里h=(a+b)/n代表二分前的步长。将每个子区间上的积分值相加得,从而可导出下列递推公式,1 梯形法的递推化,注意，这里h=(a+b)/n代表二分前的步长。将每个子区间上,21,龙贝格算法,龙贝格积分法是在计算梯形和序列的基础上应用了线性外推的加速方法，由此构成的一种具有超线性收敛的自动积分法,龙贝格算法龙贝格积分法是在计算梯形和序列的基础上应用了线,22,基本思想,根据复化梯形公式的余项表达式可知,基本思想根据复化梯形公式的余项表达式可知,23,将上式移项整理，可得,可以做这样的补偿,基本思想,将上式移项整理，可得可以做这样的补偿基本思想,24,同理,由此得到,同理,基本思想,同理由此得到同理基本思想,25,由此法，可得如下三角形数表,梯形,辛卜生,柯特斯,龙贝格,T,0,T,3,T,2,T,1,S,0,S,2,S,1,C,0,C,1,D,0,基本思想,由此法，可得如下三角形数表梯形辛卜生柯特斯龙贝格T0T3T2,26,样条插值积分,用三次样条插值函数,S,(,x,),近似被积函数,f,(,x,),，从而得到样条插值积分公式。,(,i=,0,1,n,),将,a,b,分,n,等分，,设,S,(,x,i,),m,i,，则,S,(,x,),在,x,i,x,i,+1,上为满足以下条件的三次多项式：,，,由三次 Hermite 插值多项式公式,（,P,.46）,可得,样条插值积分 用三次样条插值函数 S(x)近似被积函数 f,27,样条插值积分（续）,于是有,由于,S,(,x,),在,x,i,x,i,+1,上为三次多项式，所以simpson公式精确成立，即,于是得积分公式,U,n,样条插值积分（续）于是有 由于 S(x)在 xi,28,样条插值积分（续）,U,n,T,n,余项：,只增加计算两端点的导数，计算精度即由,O,(,h,2,),提高到,O,(,h,4,),样条插值积分（续）UnTn余项：只增加计算两端点的导数，计算,29,