3.4 Gauss求积公式讲解.ppt

3.4.3 Gauss求积公式的余项与稳定性 定理 3.4.2 设 　,则Guass公式(3.4.1)的余项是 (3.4.10) 证 由Gauss点 构造次数不超过2n+1的Hermite插值 多项式H(x),满足条件 由于Gauss公式具有2n+1次代数精度它对于H(x)能准确成立,即 特别地， 对于两点Gauss-Legendre求积公式有 由Hermite插值多项式的插值余项有 再考虑到 在[a,b]上保号,应用积分中值定理得(3.4.10),定理得证. 对比Newton-Cotes求积公式,Gauss求积公式不但具有高精度,而且是数 值稳定的.Gauss公式的稳定性之所以能够得到保证,是由于它的求积系数 具有非负性. 引理3.4.2 Gauss求积公式(3.4.1)中的系数 全部为正. 对于两点Gauss-Chebyshev求积公式有 在实际计算积分的近似值 时， 不能精确地 取到，一般只能是近似值，设 实际求 得的积分值为 定理 3.4.3 对于函数值的变化所引起的求积公式的误差有 证 对于以Gauss点 为节点的插值基函数 是2n次多项式,故Gauss公式(3.4.1)对于它能准确, 即有 由于上式左端大于零,所以有 证毕. (3.4.11) 因此，（3.4.11）成立，定理成立。 由定理3.4.3可知，数据误差对于求积公式计算值的影响是可以控制的，即 Gauss求积公式在数值计算中是稳定的。 证 由于求积系数 因此有 在Gauss求积公式中，取f(x)=1，此时求积公式精确成立，即得 例1：选取常数a，求使积分公式 的代数精度尽量高，并问其代数精度为几次？ 解：取f(x)=1.则对上求积公式，左端 f(x)=x. 则 ，左端= 综合举例: 令左端=右端 左端= 再取 故，当取 时，求积公式具有3次代数精度。 例2 若用复化梯形公式计算积分 问积分区间要多少等分 才能保证有6位有效数字？ 解：由复化梯形公式截断误差式知 由于该积分有一位整数，所以要求使近似积分有6位有效数字，只需取n满足： （有效数字定义见第一章，表示不 超过某一位数字的一半） 即 即 因此，至少要将[0，1]区间212等分。若将同一问题改为复化 Simpson公式，则由复化Simpson公式截断误差式同样可得： 由此可见，Simpson公式复化型比复化梯形公式计算量少得多。 的近似值，要求误差 例3 用Romberg求积法计算 解：此时积分限为a=0,b=1.而 (本例主要说明Romberg过程） ①． 3.4 Gauss求积公式 3.4.3 Gauss求积公式的余项与稳定性 3.4.2 常用Gauss求积公式 3.4.1 Gauss求积公式的基本理论 3.4 Gauss求积公式 学习目标： 掌握高斯求积公式的用法。会用高斯?勒让德求积公式。 前面我们限定把积分区间的等分点作为求积节点,从而构造 出一类特殊的插值求积公式，即Newton—Cotes公式。这种做法 虽然简化了算法，但却降低了所得公式的代数精度。例如：在 构造形如 的两点积分公式时，如果限定求积节点 ，那么 所 得插 值型求积公式 ，其代数精度仅为1。 但是，如果我们对公式（*）中系数 和节点 都不加