[工学]第二章1 插值法2123.ppt

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 如果要表示 附近的函数值 ，也可使用牛顿插值 公式（3.6），但为了降低误差，插值点应按 的次序排列， 作变换 ，并利用公式均差与向后差分 关系公式（4.8）， 这时 得 * 称其为牛顿后插公式， （4.13） 其中 （4.12） 其余项 * 通常求开头部分插值点附近函数值时使用牛顿前插 公式，求插值节点末尾附近函数值时使用牛顿后插公式. 如果用相同节点进行插值，则向前向后两种公式只 是形式上差别，其计算结果是相同的. * 为使用牛顿插值公式，先构造差分表（表2-4）. 给出 在 处的函数值，试用4次等距节点插值公式计算 及 例5 的近似值并估计误差. 解 根据题意，插值条件为 由于 接近 ，所以应用牛顿向前插值公式计算 的近似值. * (注意：表中带下划线的数据为 点的各阶向前差分，双下划线为 点的各阶向后差分 .) * 取 则 用表2-4上半部的各阶向前差分，得 * 由余项公式（4.11）得误差估计 其中 * 于是 计算 应使用牛顿向后插值公式， 用差分表2-4中下半部的各阶向后差分，得 这里 * 其中 由余项公式（4.13）得误差估计 * * * * * * * * * * * * * * * * * * * * * * * * 由（2.17）,其截断误差 其中 于是 * （2） 用抛物插值计算，由公式（2.5）得 * 由（2.18）,截断误差限 其中 于是 这个结果与6位有效数字的正弦函数表完全一样， 这说明查表时用二次插值精度已相当高了. * * 2.3 均差与牛顿插值公式 2.3.1 插值多项式的逐次生成 利用插值基函数很容易得到Lagrange插值多项式，公 式结构紧凑，在理论分析中甚为方便，但当插值节点增减 时全部插值基函数 均要重新计算. 【Lagrange插值多项式的缺陷】 * 利用点斜式直线方程得 为了克服这一缺点，我们设计一种逐次生成插值多项式的方法： 对n=1,插值多项式 满足 它可看成零次插值 的修正： 其中 是函数 的差商. * 其中 对n=2，插值多项式 可表示为 这里 是函数 的“差分的差分”，称为“二阶差分”，也称“均差”. * （3.1） 其中 为待定系数， 确定 . 一般地，插值多项式 表示为如下便于计算的形式 可由 个插值条件 （3.2） * 称 为函数 关 于点 的一阶均差. 称为 的二阶均差. 定义2 2.3.2 均差及其性质 * （3.3） 一般地，称 为 的 阶均差 （均差也称为差商）. * 均差有如下的基本性质： （3.4） 这个性质可用归纳法证明. 1° 阶均差可表为函数值 的线 性组合， 这性质也表明均差与节点的排列次序无关，称为均差 的对称性. 即 * 3° 若 在 上存在 阶导数，且节点 （3.5） 这公式可直接用罗尔定理证明. (3.3’) 2° 由性质1°及（3.3）可得 即 则 阶均差与导数关系如下： * 均差计算可列均差表如下（表2-1）. * 2.3.3 Newton 插值多项式 根据均差定义，把 看成 上一点， 可得 * 只要把后一式代入前一式，就得到 其中 （3.6） * （3.7） 是由（2.10）定义的. 显然，由（3.6）确定的多项式 满足插值条件， 且次数不超过 ， 称 为牛顿（Newton）均差插值多项式. 系数 就是均差表2-1中加横线的各阶均差，它比拉格朗日插值计算量省，且便于程序设计. 其系数为 它就是形如（3.1）的多项式， * 但