高斯求积公式.ppt

数值分析例:分别用不同方法计算如下积分,并做比较各种做法比较如下：1、用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.9460830I准=0.9460831数值分析第23页，共39页，星期日，2025年，2月5日数值分析2:用复化梯形公式令h=1/8=0.1253：用复化辛卜生公式令h=1/8=0.125I准=0.9460831数值分析第24页，共39页，星期日，2025年，2月5日数值分析4、用Romberg公式KTnSnCnRn00.920735510.93979330.946145920.94451350.94608690.940083030.94569060.94608330.94608310.9460831I准=0.9460831数值分析第25页，共39页，星期日，2025年，2月5日数值分析5、用Gauss公式解：令x=(t+1)/2,I准=0.9460831（2）用3个节点的Gauss公式（1）用2个节点的Gauss公式数值分析第26页，共39页，星期日，2025年，2月5日数值分析算法比较此例题的精确值为0.9460831...由例题的各种算法可知：对Newton-cotes公式，当n=1时只有1位有效数字，当n=2时有3位有效数字，当n=5时有7位有效数字。对复化梯形公式有2位有效数字，对复化辛卜生公式有6位有效数字。用复合梯形公式，对积分区间[0，1]二分了11次用2049个函数值，才可得到7位准确数字。用Romberg公式对区间二分3次，用了9个函数值，得到同样的结果。用Gauss公式仅用了3个函数值，就得到结果。数值分析第27页，共39页，星期日，2025年，2月5日数值分析2.Gauss-Chebyshev公式常用的高斯求积公式数值分析第28页，共39页，星期日，2025年，2月5日关于高斯求积公式第1页，共39页，星期日，2025年，2月5日数值分析考虑更一般形式的数值积分问题定义：若求积公式对一切不高于m次的多项式p(x)都等号成立，即R(p)=0;而对于某个m+1次多项式等号不成立，则称此求积公式的代数精度为m.一、构造高斯型求积公式的基本原理和方法

数值分析第2页，共39页，星期日，2025年，2月5日数值分析定理1：设节点x0,x1…,xn∈[a,b]，则求积公式的代数精度最高为2n+1次。分别取f(x)=1,x，x2，...xr代入公式，并让其成为等式，得：A0+A1+……+An=∫ab1dx.=b-ax0A0+x1A1+……+xnAn=∫abxdx.=（b2-a2)/2......x0rA0+x1rA1+……+xnrAn=∫abxrdxr=(br+1-ar+1)(r+1)数值分析第3页，共39页，星期日，2025年，2月5日数值分析事实上,取2n+2次多项式g(x)=(x-x0)2(x-x1)2….(x-xn)2代入求积公式,这里x0,x1…,xn是节点，有左?右,故等式不成立,求积公式的代数精度最高为2n+1次。