第6节 插值.pptVIP

先构造 l0(x)。有上表知， x1 , x2, …, xn 为 l0(x) 的零点，设 由l0(x0)=1，得 同理可设 由li (xi)=1，得 于是， 所以我们得到 n 次Lagrange插值多项式： 容易验证， Ln(xi) = f (xi), i = 0, 1, 2,…, n. 6.2 Newton插值 Lagrange插值的优缺点： 优点：形式整齐、规范，理论上保证插值的存在唯一性。 缺点：计算量大、不具有承袭性。 一、差商及其性质 一阶差商：f (x)关于点x0，x1的一阶差商记为 f [x0, x1]， 二阶差商： f (x)关于点x0，x1， x2的二阶差商记为 f [x0, x1, x2], 一般地，k 阶差商 f [x0, x1,…, xn] 定义为： 差商的性质 性质1 k 阶差商 f [x0, x1,…, xk]可表成节点上函数值 f(x0), f(x1), …, f(xk) 的线性组合，即 例如，k = 2时， (6.3) 性质 2 各阶差商具有对称性, 即改变差商中节点的次序不会  改变差商的值。设i0, i1, …, ik为0, 1, …, k的任一排列, 则 由性质1知，任意改变节点的次序，只改变(5.17)式右端求 和的次序，故其值不变。例如,由定义知， 性质 3 若 f (x)为 n 次多项式，则一阶差商 f [x, xi]为n ?1次 多项式。 由定义 令x = xi, 则分子为0, 说明分子中含有因子x ? xi, 与分母约去。 性质 4 若 f (x)在 [a, b] 存在 n +1阶导数，xi∈ [a, b] , i = 0,1,…,n, 固定 x∈[a, b], 则 n+1 阶差商与导数 存在如下关系： 二、差商的计算 ★ 由差商的定义 一阶差商是由节点上函数值定义的，二阶差商是由一阶差商定义的，…，依此构造差商表： 三、 n 次Newton插值公式 给定n+1个插值点(xi, f(xi)), i = 0, 1, 2,…, n, xi互异， 类似地，有二阶至 n 阶差商的定义得 (x?x0)× (x?x0)(x?x1) × …… (x?x0)…(x?xn?1) 上述所有n +1个等式相加，得 其中， 插值误差为： 定义 给定 [a, b]上 n +1个节点：a = x0 x1 … xn = b 及节点上的函数值 f(xi) = yi, i = 0, 1, 2,…, n, 若 S(x) 满足： (1) S(xi) = yi, i = 0,1,…, n ; (2) S(x)在[xi, xi+1]上至多是一个三次多项式； (3) S(x)∈C2[a, b]. 则称 S(x) 为 f(x) 关于节点 x0, x1,…, xn 的三次样条插值函数，称 xi 为样条节点。 在[xi, xi+1]上构造一个三次多项式 共有 4n 个未知量，因此需要 4n 个条件。 由定义中的(1)知，有n+1个条件； 由定义中的(3)知，有3n ?3个条件，即 再附加2个边界条件，就可得到4n个条件。三次分段样条插 值函数可唯一确定。 M关系式：用节点处的二阶导数表示样条插值函数。 m关系式：用节点处的一阶导数表示样条插值函数。 引入记号： ∵ S(x) 为 [xi, xi+1] 上的三次多项式, ∴ S(x)为一次函数, 我们用节点上的二阶导数值 Mi 表示线性函数。设 其中， 对(6.8)积分两次, 得 将S(xi) = yi, S(xi+1) = yi+1,代入(6.9)得二元一次方程组,解得 二、m关系式 ——用一阶导数表示的样条插值函数 给定插值点 (xi, yi), 设