高斯求积公式.ppt

数值分析二、高斯型求积公式的截断误差和稳定性分析数值分析第31页,共39页，星期六，2024年，5月数值分析已知Hermite插值误差是因为对2n+1次多项式求积公式准确成立，即代入上式即有数值分析第32页,共39页，星期六，2024年，5月数值分析以下将证明高斯形求积公式的求积系数恒正数值分析第33页,共39页，星期六，2024年，5月数值分析数值分析第34页,共39页，星期六，2024年，5月数值分析将积分区间[a,b]n等分，在每个小子区间上使用一个节点数较少的Gauss型求积公式,然后把它们加起来，就得到整个区间上Gauss型求积公式的复化形式。复化Gauss求积公式的基本思想：下面用Gauss-Legender求积公式推导复化Gauss型求积公式.将积分区间[a,b]n等分，三、复化Gauss求积公式数值分析第35页,共39页，星期六，2024年，5月数值分析数值分析第36页,共39页，星期六，2024年，5月数值分析例如,用2点的Gauss-Legender求积公式复合,由表9-4,取n=1,得Aj=1,xj=0.5773502692代入到上式中,得2点的复化Gauss-Legender求积公式再将上式应用Gauss-Legender求积公式就得到了复化Gauss型求积公式.数值分析第37页,共39页，星期六，2024年，5月数值分析数值分析二版习题P276----11(1)，16，三版习题P250-----11(1)，16，第38页,共39页，星期六，2024年，5月感谢大家观看第39页,共39页，星期六，2024年，5月关于高斯求积公式数值分析考虑更一般形式的数值积分问题定义：若求积公式对一切不高于m次的多项式p(x)都等号成立，即R(p)=0;而对于某个m+1次多项式等号不成立，则称此求积公式的代数精度为m.一、构造高斯型求积公式的基本原理和方法

数值分析第2页,共39页，星期六，2024年，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页，星期六，2024年，5月数值分析事实上,取2n+2次多项式g(x)=(x-x0)2(x-x1)2….(x-xn)2代入求积公式,这里x0,x1…,xn是节点，有左?右,故等式不成立,求积公式的代数精度最高为2n+1次。证毕.上式共有r+1个等式，2n+2个待定系数(变元),要想如上方程组有唯一解，应有方程的个数等于变元的个数,即r+1=2n+2,这样导出求积公式的代数精度至少是2n+1,下面证明代数精度只能是2n+1.数值分析第4页,共39页，星期六，2024年，5月数值分析定义:使求积公式达到最高代数精度2n+1的求积公式称为Guass求积公式。Guass求积公式的节点xk称为Guass点,系数Ak称为Guass系数.因为Guass求积公式也是插值型求积公式,故有结论:n+1个节点的插值型求积公式的代数精度d满足：n?d?2n+1。数值分析第5页,共39页，星期六，2024年，5月数值分析例：选择系数与节点，使求积公式（1）成为Gauss公式。