高斯积分法讲义重点.pptVIP

4.5 高斯求积公式 4.5.2 高斯-勒让德求积公式 4.5.3 高斯-切比雪夫求积公式 4.6 数 值 微 分 4.6.2 插值型的求导公式 4.6.3 利用数值积分求导 4.6.4 三次样条求导 4.6.5 数值微分的外推算法 代入中点公式得 从截断误差的角度看，步长越小，计算结果越准确. 其中 且 （6.2） 再考察舍入误差. 按中点公式，当 很小时，因 与 很接 近，直接相减会造成有效数字的严重损失. 因此，从舍入误差的角度来看，步长是不宜太小的. 例如，用中点公式求 在 处的一阶导数 取4位数字计算. 结果见表4-8 导数的准确值 . 从表4-8中看到 的逼近效果最好，如果进一步 缩小步长，则逼近效果反而越差. 则计算 的舍入误差上界为 这是因为当 及 分别有差入误差 及 若令 它表明 越小，舍入误差 越大，故它是病态的. 用中点公式 6.1 计算 的误差上界为 要使误差 最小，步长 不宜太大，也不宜太小. 其最优步长应为 对于列表函数 运用插值原理，可以建立插值多项式 作为它的近似. 由于多项式的求导比较容易，我们取 的值作为 的近似值，这样建立的数值公式 （6.3） 统称插值型的求导公式. 即使 与 的值相差不多， 与导数的真值 仍然可能差别很大. 导数的近似值 因而在使用求导公式 6.3 时应特别注意误差的分析. 依据插值余项定理，求导公式 6.3 的余项为 式中 但如果限定求某个节点 上的导数值，那么第二项中 由于 是 的未知函数，所以对随意给出的点 ， 误差是无法预估的. 因式 变为零，这时余项公式为 （6.4） 下面仅考察节点处的导数值并假定所给节点是等距的. 1. 两点公式 设已给出两个节点 上的函数值 对上式两端求导，记 ，有 做线性插值 于是有下列求导公式： 利用余项公式 6.4 知，带余项的两点公式是 2. 三点公式 设已给出三个节点 上的函数值， 做二次插值 令 上式可表示为 两端对 求导，有 （6.5） 式中撇号（′）表示对变量 求导数. 分别取 得到三种三点公式： 带余项的三点求导公式为 （6.6） 其中的公式 6.6 是中点公式. 它比其余两个三点公式少用 了一个函数值. 用插值多项式 作为 的近似函数，还可以建立 高阶数值微分公式： 例如，将式 6.5 再对 求导一次，有 于是有 而带余项的二阶三点公式如下： （6.7） 微分是积分的逆运算，因此可利用数值积分的方法来 计算数值微分. 设 是一个充分光滑的函数， （6.8） 对上式右边积分采用不同的求积公式就可得到不同的数值 微分公式. 则有 例如，用中矩形公式 1.2 ，则得 从而得到中点微分公式 若对 6.8 右端积分用辛普森求积公式，则有 略去上式余项，并记 的近似值为 则得到 辛普森数值微分公式 这是关于 的 个方程组， 已知， （6.9） 若 则可得 这是关于 的三对角方程组，且系数矩阵为严格 对角占优的，可用追赶法求解 见第5章5.4节 . 如果端点导数值不知道，那么对 6.9 中第1个和第 个方程可分别用 及 的中点微分公式近似， 然后求 即为 的近似值. 即取 例8 给定 的一张数据表 表4-9左部 ， 并给定 及 的值 见表4-9 . 解 解之得 利用辛普森数值微分公式求 在 上的一阶导数. 结果见表4-9. 根据 6.9 有 4.5.1 一般理论 求积公式 含有 个待定参数 当 为等距节点时得到的插值求积公式其代数精度至少为 次. 如果适当选取 有可能使求积公式 具有 次代数精度，这类求积公式称为高斯 Gauss 求积公式. 为具有一般性，研究带权积分 这里 为权函数， 类似 1.3 ，求积公式为 （5.1） 为不依赖于 的求积系数. 使（5.1）具有 次代数精度. 为求积节点， 可适当选取 定义4 如果求积公式 5.1 具有 次代数精度， 则称其节点 为高斯点，相应公式 5.1 称为高斯求积公式. 根据定义要使 5.1 具有 次代数精度，只要对 （5.2） 当给定权函数 ，求出右端积分，则可由 5.2 解得 令 5.1 精确成立， 即 例5 （5.3） 解 令公式 5.3 对于 准确成立， 试构造下列积分的高斯求积公式： 得 （5.4） 由于 利用 5.4 的第1式，可将第2式化为 同样地，利用第2式化第3式，利用第3式化第4式，分别得 从上面三个式子消去 有 进一步整理得 由此解出 从而 这样，形如 5.3 的高斯公式是 由于非线性方程组 5.2 较复杂，通常 就很难求解. 故一般不通过解方