第一节插值方法(5-6学时).pptVIP

1、5牛顿插值公式 2、差商及其性质 对给定函数 ，记 表示关于节点 的 阶差商， 3、差商形式的插值公式 4、差分形式的插值公式 1、6 埃尔米特插值 * 第三讲 教学内容： 差商、差分的概念与性质，Newton插值公式及其余项。 重点难点： 差商表、差分表，Newton插值公式的构造。 教学目标： 理解差商、差分的定义及其性质，掌握Newton插值公式 及其余项。 1、具有承袭性的插值公式 先考察线性插值的插值公式 线性插值公式可以写成如下形式： 把 看做是f(x)的零次多项式 令修正系数： 则（22）式可写成： 再修正 可以进一步得到拋物插值公式： 因为： 所以： 从而有： 一阶差商定义为： 二阶差商定义为： 零阶差商定义为： 由此可以看出，差商具有明显的承袭性， 因而可以从作为零阶差商的函数值 出发，逐步生成高阶差商。 一般地， 阶差商可以递推定义为： f(x0,x1,x2,x3) f(x1,x2,x3) f(x2,x3) f(x3) x3 f(x0,x1,x2) f(x1,x2) f(x2) x2 f(x0,x1) f(x1) x1 f(x0) x0 三阶差商 二阶差商 一阶差商 f(x) x 差商表 n 阶差商关于 是对称的，即与节点的排列次序无关。 定理1、1 f(x)关于 的 n 阶差商是在这些点上的函数值 的线性组合，即 证： 对一阶差商有： 命题成立 对二阶差商有： 命题成立 一般地，由数学归纳法可知： 根据上式，如果调换两个节点的顺序，只是意味着改变求和的次序，并不会改变其值，故差商与节点的排列次序无关。 即n 阶差商关于 是对称的 按差商定义有： ……………… 反复用后一个式子代入前面的式子 (25) 令 (26) (27) 则有： 因为： 所以： 故： 即： 这种差商形式的插值公式 称为牛顿插值公式 根据定理2，拉格朗日插值问题的解是唯一的。牛顿插值公式其实是拉格朗日插值公式的一种变形。 定理4： 在节点 所界定的范围 内存在一点ξ，使成立： 证： 由拉格郎日余项定理知：必有 因为： 使得： 所以 故： 例： 已知f(-2)=2，f(-1)=1，f(0)=2，f(0.5)=3，试选用合适的 插值节点利用二次插值多项式计算F(-0.5)的近似值，使之精度尽量可能高。 解： 根据牛顿插值余项公式 选用 作为插值节点。 通过计算，可以得到如下差商表： 2/3 2 3 0.5 1 2 0 1 -1 二阶差商 一阶差商 F(x) x 由牛顿插值公式 可以得到该问题的二次牛顿插值多项式 所以 例： 设 ,求差商 解： f(1)=7，f(2)=169，f(4)=16705 在实际应用Newton插值多项式时，经常遇到插值节点是等距的，即n+1个插值节点： 这里间距 h 为定数，称为步长，于是在差商中，分母部分将变得简单，计算量主要集中在分子（两节点处函数值的差）。 定义：设函数f(x)在等距节点 上的值为： 则称： 为 f(x)在 处一阶差分 称： 为 f(x)在 处二阶差分 称： 为 f(x)在 处n 阶差分 依据所给数据 可以逐步求出它的各阶差分，而生成如下形式的差分表： 例：给定 的函数表如下： 0 0.1 0.2 0.3 0.4 0.5 0.6 1 0.995 0.98007 0.955340.92106 0.87758 0.82534 在节点等距的情况下，差商可用差分来表示。 有 这样，对于等距节点的情况，将牛顿插值公式中的差商换成相应的差分。 令： 则： 于是，牛顿插值公式中的一般项 则 称此公式为函数插值的有限差分公式 例： 已知函数 y=sin x 的如下函数值表，利用插值法计算 Sin(0.42351) 的近似值。 0.56464 0.47943 0