第三章多项式插值方法 教学目的及要求: 要求掌握基本的定理及各种插值方法。 插值方法是数学分析中很古老的一个分支它有悠久的历史等距结点内插公 式是由我国隋朝数学家刘焯(公元544610年)首先提出的而不等距结点内插公 式是由唐朝数学家张遂(公元683-727年)提出的这比西欧学者相应结果早 千年 插值方法在数值分析的许多分支(例如,数值积分,数值微分,微分方程数值 解曲线曲面拟合函数值近似计算等等)均有应用下面仅以近似计算函数值为例 来说明 设已知某个函数关系y=f(x)的列表函数值 yo 而x≠x(=01…n)间应该如何估值j=f(x)对于函数关系y=f(x),我们所知 道仅仅上述的表列值它们常常是间接求得的例如是由实验(观测)得来的或者是 从级数或微分方程求得的 我们可以使用插值方法估计y插值方法的目的是寻求简单的连续函数o(x) 使它在n+1个点x0,x,…,xn处取给定值o(x,)=y1=f(x)(=0,1…,n),而在别处 希望它也能近似地代表函数f(x)因为(x)已是有解析表达式的简单函数所以 它在x=x处的值可以按表达式精确地计算出来这样我们就可以将o()看成 y=f(x)的近似值了
第三章 多项式插值方法 教学目的及要求： 要求掌握基本的定理及各种插值方法。 插值方法是数学分析中很古老的一个分支.它有悠久的历史.等距结点内插公 式是由我国隋朝数学家刘焯(公元 544—610 年)首先提出的;而不等距结点内插公 式是由唐朝数学家张遂(公元 683—727 年) 提出的.这比西欧学者相应结果早一 千年. 插值方法在数值分析的许多分支(例如, 数值积分, 数值微分, 微分方程数值 解,曲线曲面拟合,函数值近似计算,等等)均有应用.下面仅以近似计算函数值为例 来说明 设已知某个函数关系 y = f (x) 的列表函数值 n n y y y y x x x x 0 1 0 1 而 x x (i n) i = 0,1, 问应该如何估值 y = f ( x). 对于函数关系 y = f (x),我们所知 道仅仅上述的表列值,它们常常是间接求得的.例如是由实验(观测)得来的,或者是 从级数或微分方程求得的. 我们可以使用插值方法估计 y. 插值方法的目的是寻求简单的连续函数 (x), 使它在 n+1 个点 n x , x , , x 0 1 处取给定值 (x ) y f (x )(i 0,1, ,n) i = i = i = ,而在别处 希望它也能近似地代表函数 f (x).因为 (x) 已是有解析表达式的简单函数,所以 它在 x = x 处的值可以按表达式精确地计算出来.这样我们就可以将 (x) 看成 y = f ( x). 的近似值了
给定点x,x1…x为插值结点称函数(x)为函数f(x)的关于x0,x1…,xn的 插值函数称y=f(x)为被插函数 严格的说插值方法一词只用于x落在给定点x,x1…,xn之间的情形所以也 称它为内插法如果x落在给定点x0,x1…,xn之外并且仍以插值函数o(x)在x处 近似地代替f(x),则一般称这种近似计算函数的方法为外插法 本章我只研究多项式插值,亦即(x)是x的多项式的情形这不仅仅因为多 项式是最简单的函数,而且因为在许多场合,函数f(x)容易用多项式近似地表 示出来此外用多项式作插值函数可满意地解决一系列有应用价值的重要问题 特别是数值积分与数值微分的问题 本章讲不涉及三角插值法其实只要理解了代数多项式插值方法的实质读者就不 难自行导出关于三角多项式插值方法的一系列相应与代数多项式插值方法的理 论结果 §1 Lagrange插值公式 设y=f(x)是实变量x得单值函数,且已知f(x)在给定的n1个互异点 xn处的值y0,y1…yn,即 y=f(x)i=0,…,n 插值的基本间题是寻求多项式p(x),使得 p(x)=y,i=0.…n 设p(x)是一个m次多项式 x=ao+a,x+a, +a x 则插值问题是,如何确定p(x)中的系数an,a12…,an,使得(11)式得以满足所以该 问题等价于求解下述的线性方程组
给定点 n x , x , , x 0 1 为插值结点.称函数 (x) 为函数 f (x) 的关于 n x , x , , x 0 1 的 插值函数.称 y = f (x) 为被插函数. 严格的说,插值方法一词只用于 x 落在给定点 n x , x , , x 0 1 之间的情形,所以也 称它为内插法.如果 x 落在给定点 n x , x , , x 0 1 之外,并且仍以插值函数 (x) 在 x 处 近似地代替 f ( x).,则一般称这种近似计算函数的方法为外插法. 本章我只研究多项式插值，亦即 (x) 是 x 的多项式的情形.这不仅仅因为多 项式是最简单的函数，而且因为在许多场合，函数 f (x) 容易用多项式近似地表 示出来.此外,用多项式作插值函数可满意地解决一系列有应用价值的重要问题. 特别是数值积分与数值微分的问题. 本章讲不涉及三角插值法.其实,只要理解了代数多项式插值方法的实质读者就不 难自行导出关于三角多项式插值方法的一系列相应与代数多项式插值方法的理 论结果 §1. Lagrange 插值公式 设 y = f (x) 是实变量 x 得单值函数，且已知 f (x) 在给定的 n+1 个互异点 n x , x , , x 0 1 处的值 n y , y , , y 0 1 ,即 y f (x ), i 0, ,n. i = i = 插值的基本问题是,寻求多项式 p(x),使得 p(x ) y , i 0, n. (1.1) i = i = 设 p(x) 是一个 m 次多项式 ( ) , 0 2 = 0 + 1 + 2 + + m m p x a a x a x am x a 则插值问题是，如何确定 p(x) 中的系数 a a am , , , 0 1 ,使得(1.1)式得以满足.所以该 问题等价于求解下述的线性方程组:
a0 +axo +axo +.+anxo =yo, a0 +a,x,+axi+.+amx=y, ao+a,+ 上述的线性方程组的系数矩阵为 它是一个(n+1)×(m+1)矩阵 当m>n时,A的列数大于行数不难证明矩阵A的秩数为n+1.因为A的前 n+1列所组成的行列式为(称为 Vandermonde行列式) s ,)def 我们有 (rn…xn,x,)=∏I(x-x 为证(1.3)考虑n次多项式 显然x0,x1…xn均为它的零点且它的x”系数恰为W(x0n…x)即 W(x0…xn-1x)=W(x0…xn-)x-x0)(x-xn) 从而有下述递推关系式 W(x0,…xn-1,xn)=(xn-x0)(x,n-xn)(x0…xn)
(1.2) , , , 2 0 1 2 1 1 2 0 1 1 2 1 0 0 2 0 1 0 2 0 + + + + = + + + + = + + + + = n m n n m n m m m m a a x a x a x y a a x a x a x y a a x a x a x y 上述的线性方程组的系数矩阵为 = m n m m n n x x x x x x x x x A 1 0 2 2 1 1 2 0 0 1 1 1 它是一个(n+1)×（m+1）矩阵. 当 m>n 时，A 的列数大于行数.不难证明矩阵 A 的秩数为 n+1.因为 A 的前 n+1 列所组成的行列式为(称为 Vandermonde 行列式) ( ) m n m m n n n n x x x x x x x x x W x x x def 1 0 2 2 1 1 2 0 0 0 1 1 1 1 , . , − 我们有 ( , . , ) ( ) (1.3) 0 1 − = − j i n n j i W x x x x x 为证(1.3),考虑 n 次多项式 ( ) n n n n n n n n n x x x x x x x x x x x x W x x x 2 1 2 1 1 1 2 1 1 0 2 0 0 0 1 1 1 1 1 , . , − − − − = 显然 0 1 1 , , , n− x x x 均为它的零点,且它的 n x 系数恰为 ( ) 0 1 , . n− W x x 即 ( ) ( )( ) ( ) 0 1 0 1 0 1 , . , , . n− = n− − − n− W x x x W x x x x x x 从而有下述递推关系式 ( ) ( ) ( ) ( ) 0 1 0 1 0 1 , . , , . n− n = n − n − n− n− W x x x x x x x W x x
运用它即可证明(1.3)试式 根据(1.3)并注意到诸x0,x12…,xn互异,从而线性方程组(12)的系数矩阵的秩 数为n+1它表明(1.2)解是不唯一的,即插值问题(1.1)解不唯 当m<n时,矩阵A的行数大于列数按照(1.3)式线性方程组(12)的每m+1个 方程组成的方程组均有唯一一组解a0,a12…an但一般说来,如此求出的各组 a,a12…,an未必相同 即此时(1.2)可能是矛盾方程组 鉴于以上情形,看来取m=n是最为适宜的现在我们重提多项式插值问题:给 定n+1个互异点x02x1,…,xn,对任意一组数y,y1…yn,是否存在唯一的 p(x)∈P(x),使得如下插值条件被满足 p(x) 0, (14) 该问题的答案是肯定的今采用构造性方法把所要求的多项式p(x)求出来 试设想如果可求出具有如下性质的特殊的插值多项式l(x)∈P(=0,…m) j≠i,j=0, 1,j (15) 则多项式 p(x)=∑y1(x) (1.6) 必为满足(14)的多项式但(15)中上面的等式指出x0,…x中除x外均为l(x)的 零点因此1(x)=c(x-x0)…(x-x)x-x1)(x-xn) 其中c为常数,但(1.5)中下面的等式指出 (x-xo)-(x-x_(x-x)-(x-x) 所以
运用它即可证明(1.3)式 根据(1.3),并注意到诸 n x , x , , x 0 1 互异，从而线性方程组(1.2)的系数矩阵的秩 数为 n+1 .它表明(1.2)的解是不唯一的，即插值问题(1.1)的解不唯一。 当 m<n 时,矩阵 A 的行数大于列数.按照(1.3)式,线性方程组(1.2)的每 m+1 个 方程组成的方程组均有唯一一组解 a a am , , , 0 1 .但一般说来，如此求出的各组 a a am , , , 0 1 未必相同. 即此时(1.2)可能是矛盾方程组. 鉴于以上情形，看来取 m=n 是最为适宜的.现在我们重提多项式插值问题: 给 定 n+1 个互异点 n x , x , , x 0 1 ,对任意一组数 n y , y , , y 0 1 ，是否存在唯一的 p(x) P(x) ，使得如下插值条件被满足 p(x ) y , i 0, n. (1.4) i = i = 该问题的答案是肯定的 .今采用构造性方法把所要求的多项式 p(x) 求出来。 试设想,如果可求出具有如下性质的特殊的插值多项式 l (x) P (i 0, n). i n = ( ) (1.5) 1, 0, , 0, , = = = j i j i j n l x i 则多项式 ( ) ( ) (1.6) 0 = = n i i i p x y l x 必为满足(1.4)的多项式.但(1.5)中上面的等式,指出 n x , , x 0 中除 i x 外,均为 l (x) i 的 零点,因此 ( ) ( ) ( )( ) ( ) i i i n l x = c x − x x − x x − x x − x 0 −1 +1 其中 c 为常数，但(1.5)中下面的等式指出 ( ) ( )( ) ( ) . 1 i 0 i i 1 i i 1 i n x x x x x x x x c − − − − = − + 所以
(x-x)…(x-x)(x-xn)…(x-xn) (x2-x0)(x2-x1)(x2-x)…( 记n{x)=(x-x0)-(x-x,),则l(x)又可表示为更简洁的形式 (x-x, ww'(x 总之n次多项式 ()=∑ Jw(x) (1.7) 满足插值条件(1.4) 若q(x)∈P也满足插值条件(14),则mx)=q(x)-p(x)∈P必以x2…x为 零点即(x1)=0,=0,…n这样一来n次多项式小(x)竟然有n+1个不同的零点 是故q(x)=p(x)所以由(1)表示的n次多项式(严格地说,是次数不超过n的多 项式)是P中满足插值条件的唯一多项式它常常称作为 Lagrange插值多项式 并记为 Ln(x)=∑y 台(x-x,hn(x) 按前述推理可知 Lagrange插值多项式Ln(x)也可视为是从下面的行列式方程 中解出来的 ln 0 9) y (请读者自行补证)由(1.9)式表示的公式便于推广到一般形式的插值问题由于篇 幅所限此处不能祥述 由(1.1)所示的条件称为插值条件点组x,x1,…xn称为插值结点上面所得到
( ) ( ) ( )( ) ( ) ( ) ( )( ) ( ) . 0 1 1 0 1 1 i i i i i i n i i n i x x x x x x x x x x x x x x x x l x − − − − − − − − = − + − + 记 ( ) ( ) ( ) n w x = x − x x − x 0 ，则 l (x) i 又可表示为更简洁的形式 ( ) ( ) (x x )w (x) w x l x i i − = 总之 n 次多项式 ( ) ( ) ( ) ( ) (1.7) 0 = − = n i i i x x w x w x p x y 满足插值条件（1.4） 若 ( ) Pn q x 也满足插值条件（1.4）,则 ( ) ( ) ( ) Pn x = q x − p x 必以 n x , , x 0 为 零点,即 (xi ) = 0, i = 0, n .这样一来,n 次多项式 (x) 竟然有 n+1 个不同的零点. 是故 q(x) p(x). 所以由(1.7)表示的 n 次多项式（严格地说，是次数不超过 n 的多 项式）是 Pn 中满足插值条件的唯一多项式.它常常称作为 Lagrange 插值多项式, 并记为 ( ) ( ) ( ) ( ) = − = n i i n i x x w x w x L x y 0 按前述推理可知 Lagrange 插值多项式 L (x) n 也可视为是从下面的行列式方程 中解出来的: ( ) 0 (1.9) 1 1 1 1 2 1 2 1 1 1 0 2 0 0 0 2 = n n n n n n n n n y x x x y x x x y x x x L x x x x (请读者自行补证).由(1.9)式表示的公式便于推广到一般形式的插值问题由于篇 幅所限,此处不能祥述. 由(1.1)所示的条件称为插值条件,点组 n x , x , , x 0 1 称为插值结点.上面所得到
的结果可以从几何上解释为有且仅有一条n次代数曲线,通过平面上事先给定 的n+1个点(x,y)1=01…n,其中x≠x(≠ Lagrange插值公式(l.8)具有结构清晰、紧凑的特点,因而适合于作理论分 析和应用 例1已知f(-1)=2,f()=1,f(2)=1。求f(x)的 Lagrange插值多项式 解依公式有(xo=-1x1=1,x2=2) 1()=(x0-x)x0-x)6(2-3x+2 x l(x)= x1-x0x1-x2 l2(x) Xo XI 从而 P2(x)=2(x2-3x+2)-3(x2-x-2)+2(x2-1)=(x2-3x+8) 例2设∫(x)=e则f(-)=0.36787,()=271828./(2)=738906 依 Lagrange插值公式,有 e'≈p2(x)=1.16519x2+17520x+0.37788 §2. Newton插值公式 Lagrange插值公式的缺点是,当插值结点的个数有所变动时(例如,为了提 高精度有时需增加插值节点的个数) Lagrange因子l(x=01,…,n)就要虽之发 生变化从而整个公式的结构也要发生变化这在计算实践中是不方便的为了克 服它的上述缺点在这一节中我们引进了 Newton形的插值公式 显然n+1个结点x0,x1…,x上的n次 Lagrange插值多项式也可以写成下列 形式 p、(x)=a+a1(x-x0)+…+an(x-x0x-x1)(x-xn)
的结果可以从几何上解释为,有且仅有一条 n 次代数曲线，通过平面上事先给定 的 n+1 个点 (x , y ),i 0,1, n, i i = .其中 x x (i j) i j Lagrange 插值公式（1.8）具有结构清晰、紧凑的特点，因而适合于作理论分 析和应用. 例 1 已知 f (−1) = 2, f (1) =1, f (2) =1 。求 f (x) 的 Lagrange 插值多项式 解 依公式有 ( 1, 1, 2) x0 = − x1 = x2 = ( ) ( )( ) ( )( ) ( ) ( ) ( )( ) ( )( ) ( ) ( ) ( )( ) ( )( ) ( 1). 3 1 2 , 2 1 3 2 , 6 1 2 2 0 2 1 0 1 2 2 1 0 1 2 0 2 1 2 0 1 0 2 1 2 0 = − − − − − = = − − − − − − − = = − + − − − − = x x x x x x x x x l x x x x x x x x x x x l x x x x x x x x x x x l x 从而 ( ) ( ) ( ) ( 3 8). 6 1 2 3 2 3 2 2 1 6 1 ( ) 2 2 2 2 p2 x = x − x + − x − x − + x − = x − x + 例 2 设 ( ) x f x = e 则 f (−1) = 0.367 87, f (1) = 2.718 28, f (2) = 7.389 06 依 Lagrange 插值公式，有 ( ) 1.16519 1.175 20 0.377 88 2 e p2 x = x + x + x §2. Newton 插值公式 Lagrange 插值公式的缺点是，当插值结点的个数有所变动时（例如，为了提 高精度,有时需增加插值节点的个数）Lagrange 因子 l (x)(i n) i = 0,1, , 就要虽之发 生变化,从而整个公式的结构也要发生变化,这在计算实践中是不方便的.为了克 服它的上述缺点,在这一节中我们引进了 Newton 形的插值公式 显然,n+1 个结点 n x , x , , x 0 1 上的 n 次 Lagrange 插值多项式也可以写成下列 形式： ( ) ( ) ( )( ) ( ) n = 0 + 1 − 0 + + n − 0 − 1 − n−1 p x a a x x a x x x x x x
下面,来确定上式中的a0,a12…an 令Pn(x)表示n个结点x02x1,…xm上的(n-1)次 Lagrange插值多项式,由于 pn(x)=pn(x)=y,(=0.…n-1), 所以 Pn p, c(x-xoXx-x1)( 此处c为常数,由条件p(xn)=yn可以定出 c=.=x x, -x).(x, -x-d 又因P-1(xn)=∑yl(x) 故又有 Xxn-x1)…(x (x2-x0)…(x1-x1)Xx1-x1)…(x-xn) 引进记号 f(xx…x)=c=∑y{n(x-x 得p(x)与P(x)之间的关系 P P 同理 Pa-1(x)=pn-2(x)+f(xo, x,, , x,-1Xx-xo,)-(x-xn-2) 继续下去,最终得到
下面，来确定上式中的 , , . a0 a1 an 令 p (x) n−1 表示 n 个结点 0 1 1 , , , n− x x x 上的（n-1）次 Lagrange 插值多项式，由于 ( ) ( ) , ( 0, 1) pn xi = pn−1 xi = yi i = n − , 所以 ( ) ( ) ( )( ) ( ) n − n−1 = − 0 − 1 − n−1 p x p x c x x x x x x 此处 c 为常数，由条件 ( ) n n p x = y 可以定出 ( ) ( )( ) ( ) 0 1 1 1 − − − − − − = n n n n n n n x x x x x x y p x c 又因 ( ) ( ) − = − = 1 0 1 n i n n i i n p x y l x 故又有 ( )( ) ( ) ( ) ( )( ) ( ) ( ) 1 0 0 0 1 1 0 1 1 − = = − + − = − − − − − + − − − = n i n l i l i i l i i i i i i n i n n n n n y x x x x x x x x x x y x x x x x x y c 引进记号 ( ) ( ) 1 0 0 0 1 , , , − = = = = − n i n l i l n i i l f x x x c y x x 得 p (x) n 与 p (x) n−1 之间的关系 ( ) ( ) ( )( )( ) ( ) 1 0 1 0 1 1 , , , n = n− + n − − − n− p x p x f x x x x x x x x x 同理 ( ) ( ) ( )( )( ) ( ) 1 2 0 1 1 0 1 2 , , , n− = n− + n− − − − n− p x p x f x x x x x x x x x 继续下去，最终得到
pn、(x)=f(x0)+f(xn,x;)x-x0)+…+∫(xn,x,…,xn)x-x0)(x-x)…(x-xn)(2) 公式(22)就是 Newton型插值公式系数f(x)f(x,x)…,f(x0,x1,…,xn)由(21) 式确定 Newton插值公式的导数很不好记因此有必要另寻方法确定它们为此我们引 进差商的概念并指出 Newton插值公式中各导数f(x0,x1,…,x)=1…,n)即是 f(x)的i阶差商设已知不同的自变量xx1…,x上的函数值f(x)(=1…,n)称 ≠ x 为f(x)的一阶差商或均差)一阶差商的一阶差商 fl, x,)-5, x) ≠ x -x 叫做f(x)的二阶差商一般说来我们称(n-1)阶差商的一阶差商 为函数∫(x)的n阶差商 差商有以下诸性质 1.若F(x)=cf(x),c为常数,则 F(xn,xn1…,x)=cf(xn,xn-1…,x 2.若F(x)=f(x)+g(x),则 xo=flax 3.若∫(x)=xm,m为自然数则 (xn,xm…,x)={1 诸x的m-n次的齐次函数,n<m 4.差商∫(xn,xn1…,x)是x,x1…xn的对称函数,亦即当任意调换
( ) ( ) ( , )( ) ( , , , )( )( ) ( ) (2.2) n = 0 + 0 1 − 0 + + 0 1 n − 0 − 1 − n−1 p x f x f x x x x f x x x x x x x x x 公式（2.2）就是 Newton 型插值公式.系数 ( ) ( ) ( ) n f x , f x , x , , f x , x , , x 0 0 1 0 1 由（2.1） 式确 定(1) . Newton插值公式的导数很不好记,因此有必要另寻方法确定它们.为此我们引 进差商的概念,并指出 Newton 插值公式中各导数 f (x x x )(i n) i , , , 1, , 0 1 = 即是 f (x) 的 i 阶差商.设已知不同的自变量 n x , x , , x 0 1 上的函数值 f (x )(i n) i =1, , 称 ( ) ( ) ( ) (i j) x x f x f x f x x i j i j i j − − , = 为 f (x) 的一阶差商(或均差). 一阶差商的一阶差商 ( ) ( ) ( ) (i k) x x f x x f x x f x x x i k i j j k i j k − − = , , , , 叫做 f (x) 的二阶差商.一般说来我们称(n-1)阶差商的一阶差商 ( ) ( ) ( ) 0 1 1 1 2 0 1 0 , , , , , , , , , x x f x x x f x x x f x x x n n n n n n n − − = − − − − 为函数 f (x) 的 n 阶差商 差商有以下诸性质 1．若 F(x) = cf (x)，c 为常数，则 ( ) ( ) 1 0 1 0 F x , x , , x cf x , x , , x n n− = n n− 2．若 F(x) = f (x)+ g(x) ，则 ( ) ( ) ( ) 1 0 1 0 1 0 F x , x , , x f x , x , , x g x , x , , x n n− = n n− + n n− 3．若 ( ) m f x = x ，m 为自然数,则 ( ) − = − = . 1, , 0, , , , , 1 0 x m n n m n m n m f x x x i n n 诸 的 次的齐次函数， 4．差 商 ( ) 1 0 f x , x , , x n n− 是 n x , x , , x 0 1 的对称函数 , 亦 即 当 任 意 调 换
x02x1…,x的位置时差商的值不变例如 ∫(x0,x…xn)=∫(xn,xn1…x)=f(xn,x0…,xn1 5.差商可以表示成两行列式之商 注0规定,当n=0时,{∏(x-x1)=1 f(xo, xu x,) f(x0)f(x1)…f(xn)|x 性质1和性质2由定义可以直接推出。现在我们证明性质3。x"的一阶差商 可根据定义直接计算出来 fo =x1+x1x0++xm- 如所见,它是x12x0的m-1次齐次函数。 相继作出各阶差商并依完全归纳法,可证实下列公式 f(x,xmx)=∑xxx, +Fn-1+…+l=m-n 此处求和运算遍及所有可能的形如xnxm1x8的xn2xn1x0的m-n次齐次项。 这样便证明了性质3。 再来证明性质4。作出相继的各阶各差商之后,读者不难看出它们是由形如 f(x)/∏(x-x)的(m+1个项的和表示出来的。由完全归纳法,易求得 ∫(x0,x1xn)可由(2.1)式的右端表出。使用前面的记号
n x , x , , x 0 1 的位置时,差商的值不变,例如 ( ) ( ) ( ) 0 1 1 0 0 1 , , , , , , , , n = n n− = n n− f x x x f x x x f x x x 5. 差商可以表示成两行列式之商: 注 (1) 规定，当 n=0 时， ( ) − = n l i l i l x x 0 =1 f (x0 , x1 ,...xn ) = ( ) ( ) ... ( ) ... ... ... ... ... ... 1 1 ... 1 0 1 1 1 1 1 0 0 1 n n n n n n f x f x f x x x x x x x − − − : n n n n n n n n n x x x x x x x x x ... ... ... ... ... ... ... 1 1 ... 1 0 1 1 1 1 1 0 0 1 − − − . 性质 1 和性质 2 由定义可以直接推出。现在我们证明性质 3。 m x 的一阶差商 可根据定义直接计算出来： ( , ) ... . 1 0 0 2 1 1 1 1 0 1 0 1 0 − − − = + + + − − = m m m m m x x x x x x x x f x x 如所见，它是 1 0 x , x 的 m−1 次齐次函数。 相继作出各阶差商并依完全归纳法，可证实下列公式： ( , ,... ) ... , 0 1 1 0 0 1 n r n r r n n f x x x =x x x − ... . rn + rn−1 + + r0 = m − n 此处求和运算遍及所有可能的形如 1 0 1 0 ... r r n r n x x x n n− − 的 1 0 x , x ,...x n n− 的 m− n 次齐次项。 这样便证明了性质 3。 再来证明性质 4。作出相继的各阶各差商之后，读者不难看出它们是由形如 = − n l i l i i l f x x x 0 ( ) ( ) 的 (n +1) 个项的和表示出来的。由完全归纳法，易求得 ( , ,... ) 0 1 n f x x x 可由（ 2.1 ）式的右端表出。使用前面的记号
o(x)=(x-x0)(x-x1).(x-xn)也可将它写成 f(,) 如此便证明了性质4 最后,用完全归纳法同样可以证明性质5。 为了作数值计算,常利用形式如下的差商表 f(x) 一阶差商 二阶差商 阶差商 f(x0) f(x,) f(x0,x1) f( x2 f(x1,x2) )|f(x,x,x2,x) f(x3) f(x2,x3) 由性质4得知 Newton插值公式(2.2)中的系数 f(x),∫(x0,x)…∫(x0,x1…,x)恰标出)。因此,当已知y2=f(x1)(i=0,1…,n)时 利用差商表可以很容易地算出f(x)的各阶差商的值,而不必去记忆公式(2.1)。 因为在(n+1)个不同的点x0,x1…x上取给定值的次数不超过n的多项式 是唯一的,所以次数相同的 Newton插值多项式与 Lagrange插值多项式是恒等 的,它们的差异仅是书写形式不同而已。但是,这种差异却为计算实践带来了很 大的方便。实际上,对于 Newton插值公式来说,当需要增加一个插值结点时, 只需在原插值多项式的后面再添加一个新项就可以了。 例1已知列表函数: 3 5 2 4 求这个函数的插值多项式
( ) ( ) ( )...( ), 0 1 n ω x = x − x x − x x − x 也可将它写成 . '( ) ( ) ( , ... ) 0 0 1, = = n i i i n x f x f x x x ω 如此便证明了性质 4。 最后，用完全归纳法同样可以证明性质 5。 为了作数值计算，常利用形式如下的差商表： f (x) 一阶差商 二阶差商 三阶差商 3 2 1 0 x x x x ( ) ( ) ( ) ( ) 3 2 1 0 f x f x f x f x ( , ) ( , ) ( , ) 2 3 1 2 0 1 f x x f x x f x x ( , , ) ( , , ) 1 2 3 0 1 2 f x x x f x x x ( , , , ) 0 1 2 3 f x x x x 由性质 4 得 知 Newton 插 值 公 式 (2.2) 中 的 系 数 ( ), ( , ), , ( , , , ) 0 0 1 0 1 n f x f x x f x x x 恰标出）。因此，当已知 y f (x )(i 0,1, ,n) i = i = 时 利用差商表可以很容易地算出 f (x) 的各阶差商的值，而不必去记忆公式 (2.1) 。 因为在 (n +1) 个不同的点 n x , x , , x 0 1 上取给定值的次数不超过ｎ的多项式 是唯一的，所以次数相同的 Newton 插值多项式与 Lagrange 插值多项式是恒等 的，它们的差异仅是书写形式不同而已。但是，这种差异却为计算实践带来了很 大的方便。实际上，对于 Newton 插值公式来说，当需要增加一个插值结点时， 只需在原插值多项式的后面再添加一个新项就可以了。 例１ 已知列表函数： x ２ ３ ５ ６ y ５ ２ ３ ４ 求这个函数的插值多项式
