2012数值分析 -lec22--外推技术+高斯型求积公式.pptVIP

Richardson外推算法 高斯(Gauss)求积公式 构造高斯型求积公式的基本原理和方法 (1) 用待定系数法构造高斯求积公式 （2）利用正交多项式构造高斯求积公式 常用的高斯求积公式 一般区间的Gauss - Legendre 求积公式 如果积分区间是[a,b]，用线性变换 例 利用高斯求积公式计算 解: 令x=1/2 (1+t), 则 用高斯-Legendre求积公式计算.取n=4 积分精确值为 　　　I=ln2=0 　　由此可见，高斯公式精确度是很高的. 高斯型求积公式的截断误差和稳定性分析 高斯型求积公式的截断误差和稳定性分析 作业 教材第177页习题：12、15、17、19、20 2.Gauss-Chebyshev公式 常用的高斯求积公式 数值分析 3.Gauss-Laguerre公式 数值分析 4.Gauss-Hermite公式 数值分析 数值分析 数值分析 已知Hermite插值误差是 因为对2n+1次多项式求积公式准确成立，即 代入上式 即有 数值分析 以下将证明高斯形求积公式的求积系数恒正 数值分析 数值分析 朱立永 北京航空航天大学 数学与系统科学学院 Email: numerical_analysis@ Password:beihang 答疑时间：星期三下午2:00－5:00 答疑地点：主216 第二十一讲外推技巧与高斯型求积公式 第六章 数值积分 外推算法的基本思路 由事后误差估计式， 则， 这启发我们，可以用低阶的公式组合后成为一个高阶的公式。 类似地， 理查逊外推法是数值方法中常用的一种加速收敛技术 整理后得到： 龙贝格(Romberg)方法 龙贝格(Romberg)算法是将理查逊(Richardson)外推法应用于数值积分，由低精度求积公式推出高精度求积公式的算法。 由R外推可得龙贝格积分法（逐次分半加倍法或梯形公式外推法）： 的计算误差为 。 前面介绍的 n+1个节点的 Newton -Cotes求积公式， 其特征是节点是等距的。这种特点使得求积公式便于 构造，复化求积公式易于形成。但同时也限制了公式 的精度。 n是偶数时，代数精度为n+1， n是奇数时， 代数精度为n 。 我们知道 n+1个节点的插值型求积公式的代数精 确度不低于n 。设想：能不能在区间[a,b]上适当选择 n+1个节点 x 0x1,x2,……,xn ,使插值求积公式的代数精 度高于n？ 答案是肯定的，适当选择节点，可使公式的精度 最高达到2n+1，这就是本节所要介绍的高斯求积公式。 数值分析 考虑更一般形式的数值积分问题 定义：若求积公式 对一切不高于m次的多项式p(x)都等号成立，即R(p)=0;而对于某个m+1次多项式等号不成立，则称此求积公式的代数精度为m. 数值分析 定理1：设节点x0, x1…,xn∈[a,b]，则求积公式 的代数精度最高为2n+1次。 分别取 f(x)=1, x，x2，...xr 代入公式，并让其成为 等式，得： A0 + A1 + …… + An =∫ab1dx.= b-a x0 A0 + x1 A1+ …… +xn An =∫abxdx.= （b2-a 2)/2 ...... x0 rA0 + x1 rA1+ …… +xn rAn =∫abxr dxr =(br+1-a r+1) (r+1) 数值分析 事实上,取 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, 这样导出求积公式的代数精