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过程） ①． ③． ④． ⑤． ⑥． ⑦． ⑧． ⑨． ⑩． ②． 如此继续算得： 由于 这个实例表明,Romberg求积法计算过程不便于手工计算，但由于计算程序具有规律性，不必存储求积系数和节点，而且精度较高，因此适合于在电子计算机上进行计算。 用Romberg方法计算时，是把区间逐次分半的，因此有时称该法 为逐次分半加速法。 例4 构造三个节点的Gauss-Legendre求积公式，并给出余项估计式。 解：由于三次Legendre多项式为： 其三个零点分别为： 令它对 准确成立 则三点Gauss-Legendre求积公式为： 余项为： （节点数） 例如，若要计算 的近似值，则由上积分公式得： 上述积分准确值为： 若利用三点Simpson求积公式。则 可见在节点数目相同的情况下，Gauss求积公式的精度是相当高的。 例5 给出计算积分 的两点计算 公式， 使得对f(