计算方法笔记.docVIP

拉格朗日插值法 已知平面上四个点：(?9, 5), (?4, 2), (?1, ?2), (7, 9)，拉格朗日多项式：L（x）（黑色）穿过所有点。而每个基本多项式：y0?0(x), y1?1(x),y2?2(x)以及y3?3(x)各穿过对应的一点，并在其它的三个点的x值上取零。 在数值分析中，拉格朗日插值法是以法国十八世纪数学家约瑟夫·拉格朗日命名的一种多项式插值方法。许多实际问题中都用函数来表示某种内在联系或规律，而不少函数都只能通过实验和观测来了解。如对实践中的某个物理量进行观测，在若干个不同的地方得到相应的观测值，拉格朗日插值法可以找到一个多项式，其恰好在各个观测的点取到观测到的值。这样的多项式称为拉格朗日（插值）多项式。数学上来说，拉格朗日插值法可以给出一个恰好穿过二维平面上若干个已知点的多项式函数。拉格朗日插值法最早被英国数学家爱德华·华林于1779年发现[1]，不久后（1783年）由莱昂哈德·欧拉再次发现。1795年，拉格朗日在其著作《师范学校数学基础教程》中发表了这个插值方法，从此他的名字就和这个方法联系在一起[2]。 对于给定的若n+1个点，对应于它们的次数不超过n的拉格朗日多项式只有一个。如果计入次数更高的多项式，则有无穷个，因为所有与相差的多项式都满足条件。 对某个多项式函数，已知有给定的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]。这类现象也被称为龙格现象，解决的办法是分段用较低次数的插值多项式。 1