数值分析课件 第4章5-6节.ppt

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.1） 对于列表函数 运用插值原理，可以建立插值多项式 作为它的近似. 由于多项式的求导比较容易，我们取 的值作为 的近似值，这样建立的数值公式 （6.3） 统称插值型的求导公式. 即使 与 的值相差不多， 与导数的真值 仍然可能差别很大. 导数的近似值 因而在使用求导公式(6.3)时应特别注意误差的分析. 依据插值余项定理，求导公式(6.3)的余项为 式中 （6.3） 但如果限定求某个节点 上的导数值，那么第二项中 由于 是 的未知函数，所以对随意给出的点 ， 误差是无法预估的. 因式 变为零，这时余项公式为 （6.4） 下面仅考察节点处的导数值并假定所给节点是等距的. 1. 两点公式 设已给出两个节点 上的函数值 对上式两端求导，记 ，有 做线性插值 于是有下列求导公式： 利用余项公式(6.4)知，带余项的两点公式是 （6.4） 2. 三点公式 设已给出三个节点 上的函数值， 做二次插值 令 上式可表示为 两端对 求导，有 （6.5） 式中撇号（′）表示对变量 求导数. 分别取 得到三种三点公式： 带余项的三点求导公式为 （6.6） 其中的公式(6.6)是中点公式， 它比其余两个三点公式少用 了一个函数值. 用插值多项式 作为 的近似函数，还可以建立 高阶数值微分公式： 例如，将式(6.5)再对 求导一次，有 于是有 而带余项的二阶三点公式如下： （6.7） 微分是积分的逆运算，因此可利用数值积分的方法来 计算数值微分. 设 是一个充分光滑的函数， （6.8） 对上式右边积分采用不同的求积公式就可得到不同的数值 微分公式. 则有 例如，用中矩形公式(1.2)，则得 从而得到中点微分公式 若对(6.8)右端积分用辛普森求积公式，则有 （1.2） （6.8） 略去上式余项，并记 的近似值为 则得到 辛普森数值微分公式 这是关于 的 个方程组， 已知， （6.9） 若 则可得 这是关于 的三对角方程组，且系数矩阵为严格 对角占优的，可用追赶法求解(见第5章5.4节). 如果端点导数值不知道，那么对(6.9)中第1个和第 个方程可分别用 及 的中点微分公式近似， 然后求 即为 的近似值. 即取 （6.9） 例8 给定 的一张数据表(表4-9左部)， 并给定 及 的值(见表4-9). 解 解之得 利用辛普森数值微分公式求 在 上的一阶导数. 结果见表4-9. 根据(6.9)有 三次样条函数 与 ，不但函数值很接近，而且 导数值也很接近，并有 （6.10）