类似地可以证明N八(x)=f(x1)(=02,n) 由插值的唯一性知:N(x)≡Ln(x),因此他们的余式也相等 (n+1 即:O(x)f[x,x0,x1,…xn 0(x n l) 故有差商与导数的关系 (n+1) f[x,x0,x1…xn]= n 其中,ξ介于x,x,x1x的最大值与最小值之间
其中， 介于 的最大值与最小值之间。 故有差商与导数的关系 即： 由插值的唯一性知： 因此他们的余式也相等 类似地可以证明 n n n n n n n i i x x x x n f f x x x x x n f x f x x x x N x L x N x f x i n , , ,... 1)! ( ) [ , , ,... ] ( ) 1)! ( ) ( ) [ , , ,... ] ( ) ( ), ( ) ( ) ( 0,1,2,... ) 0 1 ( 1) 0 1 ( 1) 0 1 + = + = = = + +
重点插商 为使用方便,我们规定 flx,, xo,x, = lim f[x+h,x,xo,xir.,x h->0 x+h,xo,x, x,+fLx,lo,xy,,xn h→0 x 为重点插商
重点插商 为重点插商。 为使用方便，我们规定 [ , , ,..., ] [ , , ,..., ] [ , , ,..., ] [ , , , ,... ] [ , , , ,..., ] 0 1 0 1 0 1 0 0 1 0 0 1 lim lim n n n h n h n f x x x x dx d h f x h x x x f x x x x f x x x x x f x h x x x x = + + = = + → →
Newton插值计算 插商表1 x.f(x)一阶插商二阶插商三阶插商单元号 f(x0) f(0) f(x)f[xo,x, f(1) x Do.x. X.Xxx (2) x3 f(x3)fTxo,x,] 51xo 1 g] f[xo, x, 22, x, A3) x,f(x)/】x几xxx)
Newton插值计算 插商表1 一阶插商 二阶插商 三阶插商 单元号 F(0) F(1) F(2) F(3) … … … …… ……… F(n) ( ) k f x ( )0 f x ( )1 f x k x 0 x 1 x 2 x 3 x ( )2 f x ( )3 f x [ , ] 0 1 f x x [ , ] 0 2 f x x [ , ] 0 3 f x x [ , , ] 0 1 2 f x x x [ , , ] 0 1 3 f x x x [ , , , ] 0 1 2 3 f x x x x n x ( ) n f x [ , ] 0 n f x x [ , , ] 0 1 n f x x x [ , , , ] 0 1 2 n f x x x x
插商表2 k f(x)|阶差商 二阶差商 三阶差商 n阶差商 单元号 xo f(ro) F(0) x, f(x,)I f[xo, X,] F(1) f(x)x1,x2]f[x,x1,x2] F(2) x (x f[x1,x2,x3] F(3) x,I f(em) fIrm-I, n]I fIx,-2,m-Is x,]f[x, n n-34n-24n-1n
插商表2 k x ( ) k f x 一阶差商 二阶差商 三阶差商 n 阶差商 单元号 0 x ( )0 f x F(0) 1 x ( )1 f x 0 1 f x x [ , ] F(1) 2 x ( )2 f x 1 2 f x x [ , ] 0 1 2 f x x x [ , , ] F(2) 3 x ( )3 f x 2 3 f x x [ , ] 1 2 3 f x x x [ , , ] F(3) n x ( ) n f x 1 [ , ] n n f x x − 2 1 [ , , ] n n n f x x x − − 3 2 1 [ , , , ] n n n n f x x x x − − − 0 1 [ , , , ] n f x x x F(n)
求M(x) ■插商表1计算简单,好实现,但数值不稳 定 ■插商表2在计算杋上稳定性好,但算法复 杂 计算M(×)常采用秦九韶程序(取n=4)
求Nn (x) ◼ 插商表1计算简单，好实现，但数值不稳 定。 ◼ 插商表2在计算机上稳定性好，但算法复 杂。 ◼ 计算Nn (x)常采用秦九韶程序（取n=4）