1-插值法-上课.pptVIP

容易验证，Newton插值满足插值条件： Nn(xi) = f (xi), i = 0, 1, 2,…, n. ★ 关于Lagrange插值和Newton插值的几点说明 由插值的唯一性质，Ln(x) = Nn(x)。因此，他们的误差也相同，即当 f(x)∈Cn+1[a, b]时，有 故得差商的性质4 * 2. 牛顿插值的误差不要求函数的高阶导数存在，所以更具有一般性。它对 f(x)是由离散点给出的函数情形或 f(x)的导数不存在的情形均适用。 3. 引入记号: f [x0] = f (x0) , t0(x) = 1, t1(x) = x ? x0, t2(x) = (x ? x0)(x ? x1), …, tn(x) = (x ? x0)(x ? x1) … (x ? xn?1)， 于是 n 次Newton插值公式可表为 称 t0(x), t1(x), t2(x), …, tn(x) 为Newton插值的基函数，而且满 * 足如下关系： ti(x) = ti?1(x)(x ?xi?1), i =1,2,…, n； ti(xj) = 0， j i，  ti(xj) ≠0，j ≥i。 4. Newton插值具有承袭性质，即 5. Newton插值公式的计算量 乘：1+2+…+ (n?1)+ n = n(n+1)/2 除：n + (n?1)+…+ 2+1 = n(n+1)/2 第 i 个节点以后非零 * 例 1.7 给定四个插值点(?2,17), (0,1), (1,2), (2,19), 计算  N2(0.9), N3(0.9)。 解 x0= ? 2，x1=0，x2=1，x3=2，由例1.5知， f [x0, x1] = ? 8， f [x0, x1, x2] = 3，f [x0, x1, x2, x3] = 5/4， 所以， N2(0.9) = 17 ? 8(0.9+2)+3(0.9+2)×0.9 = 1.63; N3(0.9) = N2(0.9)+1.25×0.9(0.9+2)(0.9 ? 1)= 1.30375. * ★ Newton插值的算法 用牛顿插值公式，首先要计算各阶差商值，然后再计算插值。牛顿插值由承袭性质容易实现，关键是实现差商。 1. 输入初始数据:节点数 n, 插值点序列{x，f(xi)}, i =0,1,2,…, n, 及要计算的函数点u； 2. 形成差商表 g [k], k =1,2, …, n ;（ g [k] = f [x0,…,xk]） 3. 形成插值基函数及插值 t =1, newton = f(x0) 对 i =1,2, …, n t = (u?xi?1) ·t ; (由tk = (u?xk?1)tk ? 1 形成(u?x0) … (u?xi?1) ) newton = newton + t·g[i] 4. 输出牛顿插值 Nn(u)=newton。 * 牛顿插值公式中只用到差商表中对角线上的差商，即 f [x0], f [x0, x1], f [x0, x1, x2], …, f [x0, x1,…, xn].可以分别用一维数组 g[i] 来存放这些差商，i = 0, 1, 2, … , n. 形成差商具体步骤： (1) 对g[i]初始化，即 g[i] = f (i), i = 0, 1, 2, …, n. (2) 除了g[0](即 f (x0)) 以外，其余函数值在计算一阶差商后 不再使用，因此可用来存放一阶差商，即 g[ j] = f [xj ?1, xj], j = n, n ?1,…, 2, 1 (3) 类似地，计算二阶差商后，除g[ 1] = f [x0, x1]外, 其余一 阶差商也不再使用, 可用g[ j]存放二阶差商 f [xj ?2, xj ?1, xj], 即 g[ j] = f [xj ?2,