Newton插值法 | 差商 + Newton插值公式 + 插值余项

Newton插值法

Aitken逐次插值法虽然具有承袭性的特点，但其插值公式是递推型的，不便于进行理论分析。为此，可以把n次插值多项式改写成升幂的形式：
$$N_n(x)=c_0+c_1(x-x_0)+c_2(x-x)(x-x_1)+\cdots +c_n(x-x_0)(x-x_1)\cdots (x-x_{n-1})\text{\hspace{1em}(10)}$$
其中，$c_0,c_1,c_2,\cdots ,c_n$为待定系数。根据定理1，该多项式是惟一存在的。将插值条件
$$N_n(x_i)=f(x_i)(i=0,1,2,\cdots ,n)$$
代入（10）式中，即可以惟一确定出系数$c_0,c_1,c_2,\cdots ,c_n$，从而得到$N_n(x)$的升幂形式。为了方便计算，在具体计算这些系数之前，先引入差商的概念。

1. 差商的定义与性质

定义3：已知顺序排列的节点$x_0,x_1,x_2,\cdots ,x_{k-1},x_k,\cdots ,x_n$所对应的函数值为$f(x_0),f(x_2),\cdots$。

定义
$$f[x_0,x_k]=\frac{f(x_0)-f(x_k)}{x_0-x_k}(1\le k\le n)$$
函数$f(x)$在点$x_0,x_k$处的一阶差商；

定义
$$f[x_0;x_1,x_k]=\frac{f[x_0,x_1]-f[x_0,x_k]}{x_1-x_k}(2\le k\le n)$$
为函数$f(x)$在点$x_0,x_1,x_k$处的二阶差商；

定义
f[x0,x1;x2,xk]=f[x0;x1,x2]−f[x0;x1,xk]x2−xk(3≤k≤n) f[x_0,x_1;x_2,x_k]=\frac{f[x_0;x_1,x_2]-f[x_0;x_1,x_k]}{x_2-x_k} \quad(3\leq k\leq n) f[x0​,x1​;x2​,xk​]=x2​−xk​f[x0​;x