在中,拉格朗日插值法是以十八世纪命名的一种方法。许多实际问题中都用来表示某种内在联系或规律,而不少函数都只能通过实验和观测来了解。如对实践中的某个量进行观测,在若干个不同的地方得到相应的观测值,拉格朗日插值法可以找到一个,其恰好在各个观测的点取到观测到的值。这样的多项式称为拉格朗日(插值)多项式。上来说,拉格朗日插值法可以给出一个恰好穿过二维上若干个已知点的多项式函数。拉格朗日插值法最早被数学家于1779年发现,不久后(1783年)由再次发现。1795年,拉格朗日在其著作《师范学校数学基础教程》中发表了这个插值方法,从此他的名字就和这个方法联系在一起。
对于给定的若n+1个点,对应于它们的次数不超过n的拉格朗日多项式只有一个。如果计入次数更高的多项式,则有无穷个,因为所有与相差的多项式都满足条件。例子:
以及 各穿过对应的一点,并在其它的三个点的 x值上取零。
定义
对某个多项式,已知有给定的k + 1个取值点:
其中对应着的位置,而对应着函数在这个位置的取值。
假设任意两个不同的xj都互不相同,那么应用拉格朗日插值公式所得到的拉格朗日插值多项式为:
其中每个为拉格朗日基本多项式(或称插值基函数),其表达式为:
拉格朗日基本多项式的特点是在上取值为1,在其它的点上取值为0。
范例
假设有某个二次多项式函数,已知它在三个点上的取值为:
要求的值。
首先写出每个拉格朗日基本多项式:
然后应用拉格朗日插值法,就可以得到的表达式(为函数的插值函数):
此时代入数值就可以求出所需之值:。
证明
存在性
对于给定的k+1个点:,拉格朗日插值法的思路是找到一个在一点取值为1,而在其他点取值都是0的多项式。这样,多项式在点取值为,而在其他点取值都是0。而多项式就可以满足
在其它点取值为0的多项式容易找到,例如:
它在点取值为:。由于已经假定两两互不相同,因此上面的取值不等于0。于是,将多项式除以这个取值,就得到一个满足“在取值为1,而在其他点取值都是0的多项式”:
这就是拉格朗日基本多项式。
唯一性
次数不超过k的拉格朗日多项式至多只有一个,因为对任意两个次数不超过k的拉格朗日多项式:和,它们的差在所有k+1个点上取值都是0,因此必然是多项式的倍数。因此,如果这个差不等于0,次数就一定不小于k+1。但是是两个次数不超过k的多项式之差,它的次数也不超过k。所以,也就是说。这样就证明了唯一性。
几何性质
拉格朗日插值法中用到的拉格朗日基本多项式(由某一组确定)可以看做是由次数不超过n的多项式所组成的:的一组。首先,如果存在一组:使得,
- ,
那么,一方面多项式P是满足的拉格朗日插值多项式,另一方面P是零多项式,所以取值永远是0。所以
- 。
这证明了是的。同时它一共包含n+1个多项式,恰好等于的维数。所以构成了的一组基底。
拉格朗日基本多项式作为基底的好处是所有的多项式都是齐次的(都是n次多项式)。
优点与缺点
拉格朗日插值法的公式结构整齐紧凑,在理论分析中十分方便,然而在计算中,当插值点增加或减少一个时,所对应的基本多项式就需要全部重新计算,于是整个公式都会变化,非常繁琐。这时可以用重心拉格朗日插值法或来代替。此外,当插值点比较多的时候,拉格朗日插值多项式的次数可能会很高,因此具有的特点,也就是说尽管在已知的几个点取到给定的数值,但在附近却会和“实际上”的值之间有很大的偏差(如右下图)。这类现象也被称为,解决的办法是分段用较低次数的插值多项式。
重心拉格朗日插值法
重心拉格朗日插值法是拉格朗日插值法的一种改进。在拉格朗日插值法中,运用多项式
可以将拉格朗日基本多项式重新写为:
定义重心权
上面的表达式可以简化为:
于是拉格朗日插值多项式变为:
即所谓的重心拉格朗日插值公式(第一型)或改进拉格朗日插值公式。它的优点是当插值点的个数增加一个时,将每个都除以,就可以得到新的重心权,计算复杂度为,比重新计算每个基本多项式所需要的复杂度降了一个量级。
将以上的拉格朗日插值多项式用来对函数插值,可以得到:
因为是一个多项式。
因此,将除以后可得到:
这个公式被称为重心拉格朗日插值公式(第二型)或真正的重心拉格朗日插值公式。它继承了(1)式容易计算的特点,并且在代入x值计算的时候不必计算多项式。它的另一个优点是,结合进行插值的话,可以很好地模拟给定的函数,使得插值点个数趋于无穷时,最大偏差趋于零。同时,重心拉格朗日插值结合进行插值可以达到极佳的数值稳定性。第一型拉格朗日插值是的,而第二型拉格朗日插值是向前稳定的,并且很小。
参考来源
- E. Waring. Problems Concerning Interpolations. Philosophical Transactions of the Royal Society of London. 1779, 69: 59–67.
- (英文)E. Meijering. A chronology of interpolation: From ancient astronomy to modern signal and image processing,. Proceedings of the IEEE: 323.
- (英文)Julius Orion Smith III. . Center for Computer Research in Music and Acoustics (CCRMA), Stanford University.
- 冯有前,《数值分析》,第63页
- 李庆扬,《数值分析》第4版,第31页
- 冯有前,《数值分析》,第64页
- ^ Jean-Paul Berrut, Lloyd N. Trefethen. . Review. 2004, 46 (3): 501–517.:.
- 王兆清,李淑萍,唐炳涛. 一维重心型插值:公式、算法和应用. 山东建筑大学学报. 2007, 22 (5): 447–453.
- NICHOLAS J. HIGHAM. . IMA Journal of Numerical Analysis. 2004, 24 (4): 547–556.
- (中文)李庆扬,王能超,易大义. 《数值分析》第4版. 清华大学出版社. 2001. .
- (中文)冯有前. 《数值分析》. 清华大学出版社. 2001. .
- (中文). 太原理工大学.