(完整word版)几种插值法的应用和比较.docVIP

插值法的应用与比较 信科1302 万贤浩1格朗日插值法 在数值分析中，拉格朗日插值法是以法国十八世纪数学家约瑟夫•路易斯•拉格朗日命 名的一种多项式插值方法.许多实际问题中都用函数来表示某种内在联系或规律，而不少函数 都只能通过实验和观测来了解.如对实践中的某个物理量进行观测，在若干个不同的地方得到 相应的观测值，拉格朗日插值法可以找到一个多项式，其恰好在各个观测的点取到观测到的 值.这样的多项式称为拉格朗日(插值)多项式 .数学上来说，拉格朗日插值法可以给出一个 恰好穿过二维平面上若干个已知点的多项式函数 .拉格朗日插值法最早被英国数学家爱德 华•华林于1779年发现，不久后由莱昂哈德•欧拉再次发现.1795年，拉格朗日在其著作 《师 范学校数学基础教程》中发表了这个插值方法，从此他的名字就和这个方法联系在一起 • 1.1拉格朗日插值多项式 图1 已知平面上四个点：(-9, 5), (-4, 2), (-1, -2), (7, 9)，拉格朗日多项式： L(x)(黑 色)穿过所有点.而每个基本多项式： y0l0(x), y1l1 (x), y2l2(x)以及y丨(x)各穿过对应的 一点，并在其它的三个点的 x值上取零. 对于给定的若n 1个点(X。，y°) , (x「yj , (x^yj，对应于它们的次数不超过 n的拉 格朗日多项式L只有一个.如果计入次数更高的多项式，则有无穷个，因为所有与 L相差 (x x°)(x xj (x xn)的多项式都满足条件. 对某个多项式函数，已知有给定的 k 1个取值点： 假设任意两个不同的 Xj都互不相同，那么应用拉格朗日插值公式所得到的拉格朗日插值多项 式为： k L(x) yjlj(x)， j o 其中每个lj(x)为拉格朗日基本多项式(或称插值基函数)，其表达式为:  拉格朗日基本多项式lj X的特点是在xj上取值为1,在其它的点xi，i j上取值为o. 例：设有某个多项式函数 f，已知它在三个点上的取值为: ? f (4) 10， ? f (5) 5.25， ? f(6) 1， 要求f(18)的值. 首先写出每个拉格朗日基本多项式:  f(4)lo(x) f(5)h(x) f(6)l2(x) (x 4)(x 6) 1 (x 4)(x 5) (5 4)(5 6) (6 4)(6 5) -(x2 28x 136)， 4此时数值18就可以求出所需之值： f(18) p(18) 11. 1.2插值多项式的存在性与唯一性 存在性 对于给定的k 1个点：(x0, y0), (xk, yk)拉格朗日插值法的思路是找到一个在一点 Xj取值 为1,而在其他点取值都是 0的多项式lj(x).这样，多项式yjlj(x)在点Xj取值为yj. k 而在其他点取值都是 0 .而多项式L x yjlj (x)就可以满足 j 0 k L(x) yilj(x) 0 0 i 0 在其它点取值为0的多项式容易找到，例如:  它在点 Xj 取值为：(Xj X°) (Xj Xj 1) (Xj 上面的取值不等于0 .于是，将多项式除以这个取值，就得到一个满足“在 Xj取值为1，而在 其他点取值都是0的多项式”： (x Xj 1) (X Xj 1) (X Xk) X X 1 (X X。) 1 j Xj X (Xj X。) (Xj Xj 1) (Xj Xj 1) (Xj Xk) 这就是拉格朗日基本多项式 . 唯一性 次数不超过k的拉格朗日多项式至多只有一个， 因为对任意两个次数不超过 k的拉格朗日多项 式：P1和P2，它们的差P1 P2在所有k 1个点上取值都是0,因此必然是多项式 (x x°)(x X1) (x Xk)的倍数.因此，如果这个差 P1 P2不等于0，次数就一定不小于 k 1 .但是p1 p2是两个次数不超过k的多项式之差，它的次 数也不超过k，所以 P1 P2 0也就是说P1 P2.这样就证明了唯一性 1.3性质 拉格朗日插值法中用到的拉格朗日基本多项式 l0 ,l1, , l n (由某一组x0 x1 xn 确定)可以看做是由次数不超过 n的多项式所组成的线性空间： nX的一组基底.首先，如 果存在一组系数：0,1, , n使得， P 01 0 1l1 n1 n 0， 那么，一方面多项式 p是满足P(Xo) o,P(Xi) i, ,P(Xn) n的拉格朗日插值多 项式，另一方面 p是零多项式，所以取值永远是 0.所以 0 1 n 0， 这证明了 1o,1i, ,1n是线性无关的•同时它一共包含n 1个多项式，恰好等于 nX的维数. 所以10,11, ,1n构成了 n X 的一组基底. 拉格朗日基本多项式作为基底的好处是所有的多项式都是齐次的(都是 n次多项式)• 1.4优点与