高斯求积公式.pptVIP

高斯求积公式 引言 求积公式 高斯求积公式的系数和余项 举例 引言 n+1个节点的插值求积公式 的代数精确度不低于n求积公式,能不能在区间[a,b]上适当选择n个节点x1,x2,……,xn,使插值求积公式的代数精度高于n？ 答案是肯定的，适当选择节点，可使公式的精度最高达到2n+1，这就是所要介绍的高斯求积公式。 为考虑一般性,设求积公式为 　 定理: 若f(2n)(x)在[a,b]上连续，则高斯求积公式的余项为 其中??(a,b),w(x)=(x-x1)(x-x2)…..(x-xn)。 高斯求积公式的系数Ak恒为正,故高斯求积公式是稳定的. Guass求积公式有多种,他们的Guass点xk, Guass系数Ak都有表可以查询. 常用的高斯求积公式 　　1.Gauss - Legendre 求积公式 　 　　　　　　　　　　　　　　　 (1) 　其中高斯点为Legendre多项式的零点 Ln(x)= 对于一般有限区间[a,b]，用线性变换x=(a+b)/2+(b-a)t/2使它变成为[-1,1]。 n xk(n) Ak(n) Rn 1 0 2 2 -0.5773503 1 +0.5773503 1 3 -0.7745967 5/9=0.5555556 +0.7745967 5/9=0.5555556 0 8/9=0.8888889 4 -0.8611363 0.3478548 -0.3399810 0.6521452 +0.3399810 0.6521452 +0.8611363 0.3478548 5 -0.9061799 0.2369269 -0.5384693 0.4786287 0 0.5688889 +0.5384693 0.4786287 +0.9061799 0.2369269 例题利用高斯求积公式计算 ［解］令x=1/2 (1+t), 则 用高斯-Legendre求积公式计算.取n=5 积分精确值为 　　　I=ln2=0 　　由此可见，高斯公式精确度是很高的 2.Gauss - Chebyshev 求积公式 　 　　　　　　　　　　　　　　　 (2) 　 其中高斯点为Chebyshev 多项式Tn(x)的零点 Tn(x)=cos(narccos(x)) 　3.Gauss - Laguerre 求积公式 　 　　　　　　　　　　　　　　　 (3) 　 4 .Gauss - Hermite 求积公式 (4) 例题:分别用不同方法计算如下积分,并做比较 令I= 各种做法比较如下： 一、Newton-Cotes公式 当n=1时，即用梯形公式，I=0.9270354 当n=2时, 即用Simpson公式,I=0.9461359 当n=3时,I=0.9461090 当n=4时,I=0.9460830 当n=5时,I=0.9460831 二:用复化梯形公式 令h=1/8=0.125 三：用复化抛物线 令h=1/8=0.125 四、 Romberg公式 K Tn Sn Cn Rn 0 0.9207355 1 0.9397933 0.9461459 2 0.9445135 0.9460869 0.9400830 3 0.9456906 0.9460833 0.9460831 0.9460831 五、Gauss公式