拉格朗日插值法(图文详解)

Posted

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了拉格朗日插值法(图文详解)相关的知识,希望对你有一定的参考价值。

数值分析中,拉格朗日插值法是以法国十八世纪数学家约瑟夫·拉格朗日命名的一种多项式插值方法。许多实际问题中都用函数来表示某种内在联系或规律,而不少函数都只能通过实验和观测来了解。如对实践中的某个物理量进行观测,在若干个不同的地方得到相应的观测值,拉格朗日插值法可以找到一个多项式,其恰好在各个观测的点取到观测到的值。这样的多项式称为拉格朗日(插值)多项式数学上来说,拉格朗日插值法可以给出一个恰好穿过二维平面上若干个已知点的多项式函数。拉格朗日插值法最早被英国数学家爱德华·华林于1779年发现[1],不久后(1783年)由莱昂哈德·欧拉再次发现。1795年,拉格朗日在其著作《师范学校数学基础教程》中发表了这个插值方法,从此他的名字就和这个方法联系在一起[2]

对于给定的若n+1个点技术分享,对应于它们的次数不超过n的拉格朗日多项式技术分享只有一个。如果计入次数更高的多项式,则有无穷个,因为所有与技术分享相差技术分享的多项式都满足条件。例子:

技术分享

技术分享
已知平面上四个点:(?9, 5),(?4, 2)(?1, ?2)(7, 9),拉格朗日多项式:Lx(黑色)穿过所有点。而每个基本多项式:y0?0(x)y1?1(x)y2?2(x)以及y3?3(x)各穿过对应的一点,并在其它的三个点的x值上取零。

定义

对某个多项式函数,已知有给定的k + 1个取值点:

技术分享

其中技术分享对应着自变量的位置,而技术分享对应着函数在这个位置的取值。

假设任意两个不同的xj都互不相同,那么应用拉格朗日插值公式所得到的拉格朗日插值多项式为:

技术分享

其中每个技术分享拉格朗日基本多项式(或称插值基函数),其表达式为:

技术分享[3]

拉格朗日基本多项式技术分享的特点是在技术分享上取值为1,在其它的点技术分享上取值为0

范例

假设有某个二次多项式函数技术分享,已知它在三个点上的取值为:

  • 技术分享
  • 技术分享
  • 技术分享

要求技术分享的值。

首先写出每个拉格朗日基本多项式:

技术分享
技术分享
技术分享

然后应用拉格朗日插值法,就可以得到技术分享的表达式(技术分享为函数技术分享的插值函数):

技术分享
技术分享
技术分享

此时代入数值技术分享就可以求出所需之值:技术分享

证明

存在性

对于给定的k+1个点:技术分享,拉格朗日插值法的思路是找到一个在一点技术分享取值为1,而在其他点取值都是0的多项式技术分享。这样,多项式技术分享在点技术分享取值为技术分享,而在其他点取值都是0。而多项式技术分享就可以满足

技术分享

在其它点取值为0的多项式容易找到,例如:

技术分享

它在点技术分享取值为:技术分享。由于已经假定技术分享两两互不相同,因此上面的取值不等于0。于是,将多项式除以这个取值,就得到一个满足“在技术分享取值为1,而在其他点取值都是0的多项式”:

技术分享

这就是拉格朗日基本多项式。

唯一性

次数不超过k的拉格朗日多项式至多只有一个,因为对任意两个次数不超过k的拉格朗日多项式:技术分享技术分享,它们的差技术分享在所有k+1个点上取值都是0,因此必然是多项式技术分享的倍数。因此,如果这个差技术分享不等于0,次数就一定不小于k+1。但是技术分享是两个次数不超过k的多项式之差,它的次数也不超过k。所以技术分享,也就是说技术分享。这样就证明了唯一性[4]

几何性质

拉格朗日插值法中用到的拉格朗日基本多项式技术分享(由某一组技术分享确定)可以看做是由次数不超过n的多项式所组成的线性空间技术分享的一组基底。首先,如果存在一组系数技术分享使得,

技术分享

那么,一方面多项式P是满足技术分享的拉格朗日插值多项式,另一方面P是零多项式,所以取值永远是0。所以

技术分享

这证明了技术分享线性无关的。同时它一共包含n+1个多项式,恰好等于技术分享的维数。所以技术分享构成了技术分享的一组基底。

拉格朗日基本多项式作为基底的好处是所有的多项式都是齐次的(都是n次多项式)。

优点与缺点

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

重心拉格朗日插值法

重心拉格朗日插值法是拉格朗日插值法的一种改进。在拉格朗日插值法中,运用多项式

技术分享
技术分享
技术分享
拉格朗日插值法的数值稳定性:如图,用于模拟一个十分平稳的函数时,插值多项式的取值可能会突然出现一个大的偏差(图中的14至15中间)

可以将拉格朗日基本多项式重新写为:

技术分享

定义重心权[7][8]

技术分享

上面的表达式可以简化为:

技术分享

于是拉格朗日插值多项式变为:

技术分享

即所谓的重心拉格朗日插值公式(第一型)或改进拉格朗日插值公式。它的优点是当插值点的个数增加一个时,将每个技术分享都除以技术分享,就可以得到新的重心权技术分享,计算复杂度为技术分享,比重新计算每个基本多项式所需要的复杂度技术分享降了一个量级。

将以上的拉格朗日插值多项式用来对函数技术分享插值,可以得到:

技术分享

因为技术分享是一个多项式。

因此,将技术分享除以技术分享后可得到:

技术分享[7]

这个公式被称为重心拉格朗日插值公式(第二型)或真正的重心拉格朗日插值公式。它继承了(1)式容易计算的特点,并且在代入x值计算技术分享的时候不必计算多项式技术分享[7]。它的另一个优点是,结合切比雪夫节点进行插值的话,可以很好地模拟给定的函数,使得插值点个数趋于无穷时,最大偏差趋于零[7]。同时,重心拉格朗日插值结合切比雪夫节点进行插值可以达到极佳的数值稳定性。第一型拉格朗日插值是向后稳定的,而第二型拉格朗日插值是向前稳定的,并且勒贝格常数很小[9]

参考来源

  1. ^ E. Waring. Problems Concerning Interpolations. Philosophical Transactions of the Royal Society of London. 1779, 69: 59–67.
  2. ^ (英文)E. Meijering. A chronology of interpolation: From ancient astronomy to modern signal and image processing,. Proceedings of the IEEE: 323.
  3. ^ (英文)Julius Orion Smith III. Lagrange_Interpolation. Center for Computer Research in Music and Acoustics (CCRMA), Stanford University.
  4. ^ 冯有前,《数值分析》,第63页
  5. ^ 李庆扬,《数值分析》第4版,第31页
  6. ^ 冯有前,《数值分析》,第64页
  7. 7.0 7.1 7.2 7.3 Jean-Paul Berrut, Lloyd N. Trefethen. Barycentric Lagrange InterpolationSIAM Review. 2004, 46 (3): 501–517.doi:10.1137/S0036144502417715.
  8. ^ 王兆清,李淑萍,唐炳涛. 一维重心型插值:公式、算法和应用. 山东建筑大学学报. 2007, 22 (5): 447–453.
  9. ^ NICHOLAS J. HIGHAM. The numerical stability of barycentric Lagrange Interpolation. IMA Journal of Numerical Analysis. 2004, 24 (4): 547–556.

以上是关于拉格朗日插值法(图文详解)的主要内容,如果未能解决你的问题,请参考以下文章

拉格朗日插值法求2的平方根

拉格朗日插值公式

拉格朗日插值法理论误差怎么得的

拉格朗日插值公式

拉格朗日插值法

拉格朗日插值法