chapter05插值与逼近.pptVIP

重节点差商 设x0,x1,…,xn为[a,b]上相异节点， f?C1[a,b].则 进一步，设f?C2[a,b]. f(x)在x0, x0, x0处的二阶差商 定理 设f?Cn[a,b], [a,b]上节点x0≤x1 ≤ … ≤ xn。 差商f[x0, x1, x2,…, xn,x]的导数 设f?C2[a,b]. 定理 Newton形式的Hermite插值多项式 相应的余项为： 而插值余项为 即,插值条件的第一式得到满足。 差商表的构造 z f(z) f [zi, zj] 例 解：构造差商表 xi f(xi) 一阶差商 二阶差商 三阶差商 1 0 1 0 1 2 0.693147 2 0.693147 0.5 例 设 f (x) = x7?x4+3x+1， 求 f [20,21], f [x,20,21,…,26] 和 f [x,20,21,…,27]. 解： 显然， f (7) (x) = 7!, f (8) (x) = 0, 由性质3得 Newton插值 给定两个插值点(x0, f(x0)), (x1, f(x1)), x0≠x1, 设 线性插值 N1(x) = a0 + a1(x ? x0) ………… 直线的点斜式 易知， 给定三个互异插值点(xi, f (xi))(i=0,1,2), 设 二次插值 N2(x) = N1(x)+a(x ? x0) (x ? x1) 易知， Newtom插值公式 1 2 … … … … n?1 1 + (x ? x0) ? 2 + … … + (x ? x0)…(x ? xn?1) ? n?1 Nn(x) Rn(x) ai = f [ x0, …, xi ] 设 插值多项式的余项（截断误差）为： 于是，满足插值条件:Nn(xi)=f(xi),i=0,1,2,…,n的Newton插值多项式为： 0 xi f (xi) 一阶差商 二阶差商 三阶差商 … n 阶差商 i x0 f (x0) 1 x1 f (x1) 2 x2 f (x2) 3 x3 f (x3) … … n xn f (xn) f [x0, x1] f [x1, x2] f [x2, x3] … f [xn ?1, xn] f [x0, x1, x2] f [x1, x2, x3] … f [xn ?2, xn ?1, xn] f [x0, x1, x2,x3] … f [xn ?3,…, xn] … … … f [x0, x1, …, xn] 例 已知函数 y=f(x)的部分函数值如下表 xi -2 0 1 2 f(xi) 17 1 2 19 求f(x)的Newton插值多项式N2(x), N3(x),并求N2(0.9), N3(0.9) 解： 构造差商表 3 1 -8 19 2 1 17 2 1 0 -2 3 2 1 0 三阶差商 二阶差商 一阶差商 f(xi) xi i 构造差商表 5/4 8 17 3 1 -8 19 2 1 17 2 1 0 -2 3 2 1 0 三阶差商 二阶差商 一阶差商 f(xi) xi i ? Newton插值的算法 用牛顿插值公式，首先要计算各阶差商值，然后再计算插值。 牛顿插值由承袭性质容易实现，关键是实现差商。 输入初始数据:节点数 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 ; newton = newton + t·g[i] 输出牛顿插值 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]以外，其余函数值在计算一阶差商后不再使用，因此可用来存放一阶差商，即 g[ j] = f [xj ?1, xj], j = n, n ?1,…, 2, 1 (3) 类似地，计算二阶差