数值分析-数值积分.pptVIP

三点高斯－勒让德求积公式的形式是其它的高斯－勒让德求积公式的形式根据常用高斯－勒让德求积公式的节点和系数表，自己可以写出来.常用高斯－勒让德求积公式的节点和系数表n节点数求积节点xk求积系数Ａk010212±0.5773503123±0.77459666920.00000000000,55555555560.888888888934±0.8611363116±0.33998104360.34785484510.652145154945±0.9061798459±0.53846931010.00000000000.23692688510.47862867050.5688888889例6利用两点Gauss-Legendre求积公式计算n=1解:因为为偶函数参见书p103-例6.5.1.本章介绍的几种求积方法各具特点：(1)梯形和抛物形求积公式是低精度的方法，但对于光滑性较差的被积函数有时效果比用高精度的方法还好，再加上公式简单，因而使用非常广泛.特别在计算机上，复化的梯形公式和抛物形公式便于采用逐次对分的方法，计算程序十分简单.(2)龙贝格求积方法，其算法简单，程序也便于实现，且当节点加密时，前面的计算结果直接参与后面的计算，因而大大减少了计算量.此方法的一个最大缺点是节点的增加是成倍的.(3)高斯型求积公式，该方法是最高代数精度的求积方法，但它的节点和求积系数都没有规则，当节点增加时，前面的计算结果不能被利用，只能重新计算.因此上机计算时，需要事先输入节点数和各种高斯型求积公式的节点与系数表.它的最大优点是适用于某些无穷区间上的广义积分的计算.n前四切记住11/21/221/64/61/631/83/83/81/847/9032/9012/9032/907/90519/28875/28850/28850/28875/28819/288641/840216/84027/840272/84027/840216/84041/840常用的柯特斯系数表考，理解，记住例题解：由梯形公式得由辛普森公式得例题分别用梯形公式、辛普森公式和柯特斯公式计算积分由柯特斯公式得积分的精确值四、复化求积公式高次插值有Runge现象，怎么办？可采用分段低次插值来解决高阶Newton-Cotes公式会出现数值不稳定。而低阶Newton-Cotes公式有时又不能满足精度要求，怎么办？可将积分区间分成若干小区间，在每个小区间上用低阶求积公式计算，然后求和。看看?复化梯形公式：在每个上用梯形公式：=Tn/*中值定理*/复化梯形公式积分法?复化Simpson公式：44444=Sn复化Simpson公式积分法?复化Cotes公式：=Cn理解例:利用数据表01/83/81/25/83/47/811/422.265492.460002.876403.200003.506853.764703.938464计算积分解：这个问题有明显的答案取n=8用复化梯形公式=3.138988494取n=4用辛卜生公式=3.141592502运算量基本相同Romberg方法对于不易直接用积分公式计算的原函数，通常用复合梯形求积公式或复合抛物线求积公式等方法，但这些方法精度不高，收敛的速度缓慢。为了提高收敛速度，减少计算量，人们寻求其他方法.Romberg方法也称为逐次分半加速法。它是在梯形公式、辛卜生公式和柯特斯公式之间的关系的基础上，构造出一种加速计算积分的方法。?作为一种外推算法,它在不增加计算量的前提下提高了误差的精度.龙贝格求积公式梯形法的递推化上节介绍的复化求积方法可提高求积精度，实际计算时若精度不够可将步长逐次分半.设将区间[a,b]分为n等分，共有n+1个分点，如果将求积区间再分一次，则分点增至2n+1个，我们将二分前后两个积分值联系起来加以考虑.并注意到每个子区间[xk,xk+1]经过二分只增加了一个分点设hn=(b-a)/n,xk=a+khn(k=0,1,?,n),在[xk,xk+1]上用梯形公式得在[xk,xk+