拉格朗日插值法

线性插值法

线性插值法是指使用链接两个已知量的直线来确定在这两个已知量之间的一个未知量的值的方法。
假设已知坐标(x0,y0)(x1,y1),要得到[x0,x1]区间内某一位置x在直线上的值,根据图中所示,得到两点式直线方程:

yy0y1y0=xx0x1x0y=y0+y1y0x1x0(xx0)

拉格朗日插值法

线性插值近似法:

p(x)==f(x0)+f(x1)f(x0)x1x0(xx0)

函数上两点之间的近似随着所近似的函数的二阶导数的增大而逐渐变差。从直观上来看也是这样:函数的曲率越大,简单线性插值近似的误差也越大。

拉格朗日插值法

对某个多项式,已知有k+1个取值点:
x0,y0,,(xk,yk)
其中xj对应自变量的位置,yj对应取值。
对于任意两个不同的xj,拉格朗日插值多项式为:

L(x)=Σkj=0yilj(x)

其中,每个lj(x)为拉格朗日基本多项式(或称插值奇函数),其表达式为:
拉格朗日插值法
lj,j=0,1,...,n之间线性无关。
优点与缺点
拉格朗日插值法的公式结构整齐紧凑,在理论分析中十分方便,然而在计算中,当插值点增加或减少一个时,所对应的基本多项式就需要全部重新计算,于是整个公式都会变化,非常繁琐。这时可以用重心拉格朗日插值法或牛顿插值法来代替。此外,当插值点比较多的时候,拉格朗日插值多项式的次数可能会很高,因此具有数值不稳定的特点,也就是说尽管在已知的几个点取到给定的数值,但在附近却会和“实际上”的值之间有很大的偏差。这类现象也被称为龙格现象,解决的办法是分段用较低次数的插值多项式。

重心拉格朗日插值法

设多项式l(x)=(xx0)(xx1)(xxk)
可将拉格朗日基本多项式重写为:

lj(x)=l(x)xxj1Πki=0,ij(xjxi)

定义重心权:
wj=1Πki=0,ij(xjxi)

则上式简化为:
lj(x)=l(x)wjxxj

拉格朗日多项式变为:
L(x)=l(x)Σkj=0wjxxjyj

上式即为重心拉格朗日插值多项式(第一型)或改进拉格朗日插值公式。
将以上拉格朗日插值多项式用来对函数g(x)1插值,得到:
x,g(x)=l(x)Σkj=0wjxxj

因为g(x)1是个多项式
因此将Lx除以g(x)后得:
L(x)=Σkj=0wjxxjyjΣkj=0wjxxj

这个公式被称为重心拉格朗日插值公式(第二型)或真正的重心拉格朗日插值公式。它继承了(1)式容易计算的特点,并且在代入x值计算L(x)的时候不必计算多项式l(x), 第一型拉格朗日插值是向后稳定的,而第二型拉格朗日插值是向前稳定的。