插值型求积公式的比较.docxVIP

插值型求积公式的比较信科1304 魏佳铭问题引入在计算过程中要求I(f)=根据Newton-Leibniz公式=F(a)-F(b)可以计算得出但是在实际计算中存在原函数表达式复杂，求解困难；f(x)表达式未知，只有通过测量或实验得来的数据表等问题，使问题无法求解。二、方法引进为得出数值积分的算法，我们对其做离散化处理来简化问题，方便计算。若存在实数,且任取f(x)∈C[a,b]都有为一个数值求积公式。称为求积系数，称为求积节点，而称 R(f)=为求积余项评价一个求积公式的优劣可以用求积余项说明，通常用与求积余项有关的所谓代数精度来评价求积余项。三、插值型求积公式插值型求积公式借助多项式插值函数来构造求积公式。常用的插值型求积公式有Newton-Cotes求积公式及Gauss求积公式。Newton-Cotes求积公式n点的Newton-Cotes公式（等距节点求积公式）为了使插值型求积系数计算更简单，将求积节点取为[a,b]上的等距节点=令积分变量x=a+th作变换，当x∈[a,b]时，有t∈[0,n-1],于是有插值型求积公式的求积系数为记=则有得2点的Newton-Cotes公式（梯形公式）几何意义：用一条过两点的直线近似代替被积函数的曲线，从而用一个梯形的面积来近似代替一个曲边梯形的面积。梯形公式余项3点的Newton-Cotes公式（Simpson公式或抛物线公式）几何意义：用通过A,B,C三点的抛物线代替y=f(x)所得曲边梯形的面积。抛物线公式余项评价：n点Newton-Cotes公式，其特征是节点是等距的，这种特点使得求积公式便于构造，求积公式易于形成。但同时也限制了公式的精度。n是偶数时，代数精度为n+1，n是奇数时，代数精度为n。Newton-Cotes求积公式在求积系数不为负数时是数值稳定的。求积节点数n越大，对应的求积公式精度越高。但当n8时，有正有负，随n增大而增大，从而导致舍入误差增大。故n＞8时，Newton-Cotes公式是数值不稳定的。因而一般不用n＞8的Newton-Cotes公式来做定积分计算。从关于稳定性的分析可知，数值计算中应用低阶的Newton-Cotes公式。Gauss型求积公式设是n+1次正交多项式的n+1个零点，则插值型求积公式是Gauss型求积公式（1）Gauss型求积公式的构造[1]用待定系数法构造Gauss型求积公式[2]利用正交多项式构造Gauss型求积公式：1.以n+1次正交多项式的零点作为积分点2.用高斯点对f(x)作Lagrange插值多项式f(x)≈代入积分式=因此，求积系数为（i=0，1……n）（2）Gauss型求积公式的求积余项由Gauss型求积公式的代数精度为2n-1及积分中值定理有Gauss求积余项为R（f）=，m∈[a,b]（3）Gauss型求积公式的数值稳定性用n+1个点构造的插值型求积公式的代数精度不超过2n+1。 Gauss公式是插值型求积公式中代数精度最高的，高斯求积公式的求积系数全是正的，且是稳定的算法。（2）常用的Gauss型求积公式 [1]Gauss-Legendre求积公式求积余项R(f)= m∈（-1,1）Gauss-Legendre的优点：计算精度高；可计算无穷区间上的积分和奇异积分。Gauss-Legendre的缺点：需计算Gauss点和Gauss系数；增加节点时需重新计算。[2]Gauss-Chebyshev求积公式求积余项R(f)= m∈（-1,1）[3]Gauss-Laguerre求积公式求积余项R(f)= m∈（0，+∞）[4]Gauss-Hermite求积公式求积余项R(f)= m∈（-∞，+∞）评价：高斯型求积公式是插值型求积公式，其系数由高斯点唯一确定，且高斯型求积公式是代数精度最高的求积公式，具有代数精度高，且总是收敛、稳定的优点。但当求积点数增加时，前面的函数值不能在后面利用。四、例题比较用不同方法计算，并作比较用Newton-Cotes公式[1]运用梯形公式I≈[2]运用Simpson公式I≈=≈02）用Gauss公式令x=,[1]用2个节点的Gauss公式 I≈+=0.9460411[2]用3个节点的Gauss公式I≈0.5555556×+0.8888889×+0.5555556×=0.9460831算法比较：对Newton-Cotes公式，当n=1时只有1位有效数字，当n=2时有三位有效数字，而用Gauss公式用了3个函数值，得到了7位有效数字。五、归纳总结插值法是常见的求积分方法，解决了积分的计算问题，Newton-Cotes公式和Gauss型求积公式是两种不同的插值法。前者是等距节点下的