最新《数值积分》第二章-插值法.pptVIP

X -2 -1 0 1 2 3 -5 1 1 1 7 25 X 2.4 2.5 2.6 f(x) 0.0025 -0.0484 -0.0968 X y 2.4 0.0025 -0.0509 2.5 -0.0484 0.0025 2.6 -0.0968 -0.0484 X 1 2 3 y 2 4 12 3 x y f[0,1] f[0,1,2] f[0,1,2,3] 1 2 2 1 2 2 4 3 5 2 4 8 3 12 因此要定这4n个待定系数，还缺少两个条件，这两个条件常在插值区间[a,b]的边界点a,b处给出，称为边界条件。 虽然可利用方程组(1)和边界条件求出所有待定系数 从而得到三次样条函数在 的表达式 ，但这样做工作量太大，不便于应用，下面介绍一种简便求法。 设在节点 处 的二阶导数为 因为在子区间 上 是不高于三次的多项式，其二阶导 数 必是线性函数（或常数），于是有： 引入 差商的概念后可以像 确定前插公式中的系数 那样逐步地确定（1）的系数： 的 n 次插值多项式为： 故满足插值条件 称之为Newton 的基本插值多项式。 x 一阶均差 二阶均差 100 10 121 11 144 12 解: 利用此公式来计算非等距节点上的函数值。 基本插值公式计算 的近似值。 Newton基本插值多项式的截断误差为： 由表可见Newton基本插值多项式（6）中各系数依次为： 故用线形插值所求的近似值为： 用抛物插值所求的近似值为： 这里所求结果与前面例中计算一致。但计算量却明显减少。 可以利用公式（7）计算上例中的截断误差。 是我们要计算的，故不能准确的计算出 确值，只能作出一种估计，例如，当四阶差商变化不大时 最后，我们指出，在带有差商的差商型的插值公式中， 为了计算差商需要进行多次除法，因此当节点等距时，应 当用差分代替差商，若节点是可以随意取的，则自然应当 选为等距的。此外，利用差分作插值多项式，也要合理 地决定差分的阶数，以避免高阶差分的误差积累，而这种 积累有时是很严重的。请看下例。 例：求序列{0，0，0， ，0，0，0}的差分表 表中仅有一个微小误差 ，但在各阶差分中，随着差分