Newton插值多项式拟合

Newton插值法

引子：在工程应用中，函数$y=f(x)$可能不是一个具体的解析表达式，而是通过实验、测量或者中间计算得到的一组数据；或者解析表达式较复杂，不便于计算和使用。因此需要一个简单函数$y=y(x)$来近似代替函数$y=f(x)$,使得

$$y(x_i)=f(x_i),i=0,1,2,..,n$$
称函数$y=y(x)$为函数$y=f(x)$在点$x_0,x_1,\dots ,x_n$处的插值函数。

插商

设函数$f(x)$在$[a,b]$上有定义，$[a,b]$上互异的节点$x_0,x_1,\dots ,x_n$的函数值为$f(x_0),f(x_1),f(x_2),\dots ,$称
$$f[x_0,x_1]=\frac{f(x_0)-f(x_1)}{x_0-x_1}$$
为$f(x)$在$x_0,x_1$处的一阶差商。
称$$f[x_0,x_1,x_2]=\frac{f[x_0,x_1]-f[x_1,x_2]}{x_0-x_2}$$为$f(x)$在$x_0,x_1,x_2处的二阶差商。

插值多项式

由差商的定义$f(x)=f(x_0)+f[x,x_0](x-x_0]$
$f[x,x_0]=f[x_0,x_1]+f[x,x_0,x_1](x-x_1)$
$f[x,x_0,x_1]=f[x_0,x_1,x_2]=f[x_0,x_1,x_2](x-x_2)$
$$\dots$$
$f[x,x_0,x_1,\dots ,x_{n-1}]=f[x_0,x_1,\dots ,x_{n-1},x_n]+f[x,x_0,x_1,\dots ,x_{n-1},x_n](x-x_n)$
依次将后一个灯饰代入前一个等式就有
$f(x)=f(x_0)+f[x_0,x_1](x-x_0)+f[x_0,x_1,x_2](x-x_0)(x-x_1)+\cdots +f[x_0,x_1,\dots ,x_{n-1},x_n](x-x_0)(x-x_1)\dots (x-x_{n-1})+f[x,x_0,x_1,\dots ,x_n](x-x_0)(x-x_1)\dots (x-x_n)$
记为$f(x)=N_n(x)+f[x,x_0,x_1,\dots ,x_{n-1},x_n]P_{n+1}(x)$
$$P_{n+1}(x)=\prod _{i=0}^n(x-x_i)$$
将上式中的x取为插值节点，由$N_n(x_i)=f(x_i),i=0,1,\dots n$,即$N_n(x)$是满足插值条件的$n$次多项式，称之为$Newton$插值多项式，其余项为
$E(x)=f(x)-N_n(x)=f[x,x_0,x_1,\dots ,x_n]p_{n+1}(x)$

差商表

插值节点	零阶差商	一阶差商	二阶差商	三阶差商	$\dots$
$x_0$	$f(x_0)$
$x_1$	$f(x_1)$	$f[x_0,x_1]$
$x_2$	$f(x_2)	$f[x_1,x_2]$	$f[x_0,x_1,x_2]$
$x_3$	$f(x_3)$	$f[x_2,x_3]$	$f[x_1,x_2,x_3]$	$f[x_0,x_1,x_2,x_3]$
$\dots$	$\dots$	$\dots$	$\dots$	$\dots$

例题

设有如下表的数据，求f(0.596)的近似值

$x_i$	0.40	0.55	0.65	0.80	0.90	1.05
$f(x_i)$	0.41075	0.57815	0.69675	0.88811	0.102652	1.25382

代码实现

X=input('请输入X的值');
Y=input('请输入Y的值');
x=input('请输入要拟合的节点的值');
n=length(X);
A=zeros(n,n+1);
A(:,1)=X;
A(:,2)=Y;
%计算差商表
for i=3:n+1
    A(i-1:n,i)=(A(i-2:n-1,i-1)-A(i-1:n,i-1))./(A(1:n+2-i,1)-A(i-1:n,1));
end

%计算插值多项式
B=A(:,2:n+1);
disp(A);
disp(B);
C=zeros(n);
C(1,1)=1;
for i=2:n
    C(i,i) = C(i-1,i-1)*(x-X(i-1));
end
disp(C);
f=0;
for j=1:n
    f=f+B(j,j)*C(j,j);
end
%展示计算结果
disp("f(0.596)的拟合值为：");
disp(f);

%运行结果
请输入X的值[0.40,0.55,0.65,0.80,0.90,1.05]
请输入Y的值[0.41075,0.57815,0.69675,0.88811,1.02652,1.25382]
请输入要拟合的节点的值0.596
    0.4000    0.4108         0         0         0         0         0
    0.5500    0.5782    1.1160         0         0         0         0
    0.6500    0.6967    1.1860    0.2800         0         0         0
    0.8000    0.8881    1.2757    0.3589    0.1973         0         0
    0.9000    1.0265    1.3841    0.4335    0.2130    0.0312         0
    1.0500    1.2538    1.5153    0.5249    0.2287    0.0314    0.0003

    0.4108         0         0         0         0         0
    0.5782    1.1160         0         0         0         0
    0.6967    1.1860    0.2800         0         0         0
    0.8881    1.2757    0.3589    0.1973         0         0
    1.0265    1.3841    0.4335    0.2130    0.0312         0
    1.2538    1.5153    0.5249    0.2287    0.0314    0.0003

    1.0000         0         0         0         0         0
         0    0.1960         0         0         0         0
         0         0    0.0090         0         0         0
         0         0         0   -0.0005         0         0
         0         0         0         0    0.0001         0
         0         0         0         0         0   -0.0000

  f(0.596)的拟合值为：
    0.6319