第二单元插值法.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)和边界条件求出所有待定系数 从而得到三次样条函数在 的表达式 ，但这样做工作量太大，不便于应用，下面介绍一种简便求法。 设在节点 处 的二阶导数为 因为在子区间 上 是不高于三次的多项式，其二阶导 数 必是线性函数（或常数），于是有： 在等距节点 下，还可以引入向后差分和中心差分。定义如下：y=f(x)在xk处以h为步长的一阶中心差分和m阶向后差分为 y=f(x)在xk点处以h为步长的一阶中心差分和m阶中心差分: 各阶向后差分与中心差分的计算均可类似如向前差分，构造向后差分与中心差分表来完成计算。 若将节点 排列顺序，则（1）可以写为： 利用向后差分也可用简化牛顿插值公式，导出类似牛顿向后插值公 式这是因为： 利用插值公式的条件 与导出（2）的方法类似可得一个向后差分表示的插值多项式： ( 其中 t0 ), 称之为牛顿向后插值公式简称向后插值公式， （4）适合于计算函数值表表尾 x 处的函数值。 例：已知函数表同上表，计算sin(0.58)的近似值，并估计截断误差. 解：因x=0.58位于表尾 附近，故用后插公式计算sin(0.58)的近似值。 插值余项公式为: 为了计算函数在 处的各阶向后差分，应构造向后差分表。但由向前向后差定义可见，对于同一函数表来说，构造出来的向前差分向后差分再数值上完全相同，因此上表中的数据依次给出了sinx 在 x=0.6 处的函数值和各阶向后差分值。 于是由后插公式（4）可得： 因三阶向后差分接近于常数，故用三次插值进行计算，且 由（5）可知截断误差为： 与自变量之差 的比值 设f(x)在一串互异点 上的值依