_数值积分与数值微分.PPTVIP

对于一般区间 上的积分，也可以利用教材中表格的数据 写出Gauss型求积公式。其原理与方法是：先作变量替换，令 则将 上积分化为 上的积分： 记 ，则上式化为： 利用表中数据，对于给定的n=1，2，3，4，可以写出Gauss型 公式： 即 代入（A）得： 其中系数 和节点 可查表得出，由变量替换公式 易见，由于求积公式（C）对变量t不高于2n+1 的多项式准确成立，从而求积公式（D）对自变量x的不高于2n+1 的多项式也准确成立，即（D）是Gauss型求积公式。 例：利用四点Gauss求积公式计算 的近似值。 解：由教材中表示Gauss型求积公式（D）得： 其中a=0，b=1， = -0 = -0 将上述各数据代入上公式中有： 优点：在此例计算过程中，只涉及到 四个点上的函数值，可 见Gauss型求积公式具有计算工作量小，所得近似值精确度 高的优点，是一种高精度的求积公式。 Gauss型求积公式的缺点是：当n改变大小时，系数和节点几 乎都在改变。虽然可以通过其他资料查到较大n的系数和节点， 但应用时却十分不便。同时，余项却涉及到高阶导数（被积 函数的），要利用它们来控制精度也十分困难。 为克服这些缺点，在实际计算中较多地采用复合求积 的方法。例如，先把积分区间 分成m个等长的小 区间 然后，在每个小区间上使用同 一低阶（如二点的，三点的，… ）高斯型求积公式算 出积分近似值，再相加即将积分 的近似值： 其中， 由查表可得。同时在实际计算中， 还常用相邻两次计算结果 和 的关系式 相当于相对误差）即算出 后，观察 是否 成立，以判定是否终止计算。请同学们据此编程计算。 来控制运算（当 时， 证明：以 为节点构造次数 的多项式H(x) 使满足条件： 由第二章Hermite插值知，H(x)为Hermite插值多项式。 由于Gauss公式具有2n+1次代数精度，则它对H(x)能准 确成立，即 其中 定理：对于Gauss公式（1），其余项为： 再注意到函数 在 上保号，应用积分中值定理即可得结论。 对比Newton-cotes公式，Gauss不但具有高精度，而且是数值稳 定的。Gauss公式之所以能够保证数值稳定性，是由于其求积系数 具有稳定性。 证明:由于 是n次多项式，因而 是2n次 多项式，故Gauss公式对 能准确成立。即 的求积系数 全是正的。即 定理：Gauss求积公式 注意到 从而上式右端实际上等于 从而有： 故定理得证。 ，由于实际问题中，通常 不一定提供准确的数据 ，而只是给出含有误差的数据 （如：由于计算机字长的限制而产生的舍入误差）。因而实际 求得的积分值为： 那么，原始数据的误差对积分值的影响能否加以控制呢？ 上面已经证明Gauss公式的求积系数 则 又由于Gauss公式对f(x)=1准确成立,即 下面讨论求积过程的数值稳定性问题。 利用公式 计算积分时 则 由此可以断定Gauss公式是稳定的。 ，其中 称为权函数。 考察积分 当 时即为普通积分 。 可以仿照处理普通积分的方法讨论带权的积分。如，求积分公式 如果对于任意次数不超过2n+1的多项式均能准确地成立，则称 之为Gauss型的。上述Gauss公式求积