几种插值法的应用和比较

1、插值法的应用与比较 信科1302 万贤浩 132710381格朗日插值法在数值分析中，拉格朗日插值法是以法国十八世纪数学家约瑟夫路易斯拉格朗日命名的一种多项式插值方法.许多实际问题中都用函数来表示某种内在联系或规律，而不少函数都只能通过实验和观测来了解.如对实践中的某个物理量进行观测，在若干个不同的地方得到相应的观测值，拉格朗日插值法可以找到一个多项式，其恰好在各个观测的点取到观测到的值.这样的多项式称为拉格朗日（插值）多项式.数学上来说，拉格朗日插值法可以给出一个恰好穿过二维平面上若干个已知点的多项式函数.拉格朗日插值法最早被英国数学家爱德华华林于1779年发现，不久后由莱昂哈德欧拉再次发现

2、.1795年，拉格朗日在其著作师范学校数学基础教程中发表了这个插值方法，从此他的名字就和这个方法联系在一起.1.1拉格朗日插值多项式图1已知平面上四个点：(9, 5), (4, 2), (1, 2), (7, 9)，拉格朗日多项式：（黑色）穿过所有点.而每个基本多项式：, 以及各穿过对应的一点，并在其它的三个点的值上取零.对于给定的若个点,，对应于它们的次数不超过的拉格朗日多项式只有一个.如果计入次数更高的多项式，则有无穷个，因为所有与相差的多项式都满足条件.对某个多项式函数，已知有给定的个取值点：,其中对应着自变量的位置，而对应着函数在这个位置的取值.假设任意两个不同的都互不相同，那么应用拉

3、格朗日插值公式所得到的拉格朗日插值多项式为：，其中每个为拉格朗日基本多项式（或称插值基函数），其表达式为：，拉格朗日基本多项式的特点是在上取值为1，在其它的点， 上取值为0.例：设有某个多项式函数，已知它在三个点上的取值为： ， ， ，要求的值.首先写出每个拉格朗日基本多项式：；然后应用拉格朗日插值法，就可以得到的表达式（为函数的插值函数）：，此时数值就可以求出所需之值：.1.2插值多项式的存在性与唯一性 存在性对于给定的个点：拉格朗日插值法的思路是找到一个在一点取值为，而在其他点取值都是的多项式.这样，多项式在点取值为，而在其他点取值都是.而多项式就可以满足，在其它点取值为的多项式容易找到，

4、例如：，它在点取值为：.由于已经假定两两互不相同，因此上面的取值不等于.于是，将多项式除以这个取值，就得到一个满足“在取值为，而在其他点取值都是的多项式”：，这就是拉格朗日基本多项式.唯一性次数不超过的拉格朗日多项式至多只有一个，因为对任意两个次数不超过的拉格朗日多项式：和，它们的差在所有个点上取值都是,因此必然是多项式的倍数.因此，如果这个差不等于，次数就一定不小于.但是是两个次数不超过的多项式之差，它的次数也不超过，所以也就是说.这样就证明了唯一性.1.3性质拉格朗日插值法中用到的拉格朗日基本多项式（由某一组 确定）可以看做是由次数不超过的多项式所组成的线性空间：的一组基底.首先，如果存在

5、一组系数：使得，那么，一方面多项式是满足的拉格朗日插值多项式，另一方面是零多项式，所以取值永远是.所以，这证明了 是线性无关的.同时它一共包含个多项式，恰好等于的维数.所以 构成了 的一组基底.拉格朗日基本多项式作为基底的好处是所有的多项式都是齐次的（都是次多项式）.1.4优点与缺点拉格朗日插值法的公式结构整齐紧凑，在理论分析中十分方便，然而在计算中，当插值点增加或减少一个时，所对应的基本多项式就需要全部重新计算，于是整个公式都会变化，非常繁琐.这时可以用重心拉格朗日插值法或牛顿插值法来代替.此外，当插值点比较多的时候，拉格朗日插值多项式的次数可能会很高，因此具有数值不稳定的特点，也就是说尽管

6、在已知的几个点取到给定的数值，但在附近却会和“实际上”的值之间有很大的偏差.这类现象也被称为龙格现象，解决的办法是分段用较低次数的插值多项式.2 重心拉格朗日插值法重心拉格朗日插值法是拉格朗日插值法的一种改进.在拉格朗日插值法中，运用多项式，图（2）拉格朗日插值法的数值稳定性：如图(2)，用于模拟一个十分平稳的函数时，插值多项式的取值可能会突然出现一个大的偏差（图中的14至15中间）可以将拉格朗日基本多项式重新写为：，定义重心权，上面的表达式可以简化为：，于是拉格朗日插值多项式变为： ， （1）即所谓的重心拉格朗日插值公式（第一型）或改进拉格朗日插值公式.它的优点是当插值点的个数增加一个时，将

7、每个都除以，就可以得到新的重心权，计算复杂度为，比重新计算每个基本多项式所需要的复杂度降了一个量级.将以上的拉格朗日插值多项式用来对函数插值，可以得到：，因为是一个多项式.因此，将除以后可得到：， （2）这个公式被称为重心拉格朗日插值公式（第二型）或真正的重心拉格朗日插值公式.它继承了（1）式容易计算的特点，并且在代入值计算的时候不必计算多项式它的另一个优点是，结合切比雪夫节点进行插值的话，可以很好地模拟给定的函数，使得插值点个数趋于无穷时，最大偏差趋于零.同时，重心拉格朗日插值结合切比雪夫节点进行插值可以达到极佳的数值稳定性.第一型拉格朗日插值是向后稳定的，而第二型拉格朗日插值是向前稳定的，

8、并且勒贝格常数很小.3.分段线性插值对于分段线性插值，我们看一下下面的情况.3.1问题的重诉 已知,用分段线性插值法求插值，绘出插值结果图形，并观察插值误差.1.在-6，6中平均选取5个点作插值；2.在-6，6中平均选取11个点作插值；3.在-6，6中平均选取21个点作插值；4.在-6，6中平均选取41个点作插值.3.2问题的分析在数值计算中，已知数据通常是离散的，如果要得到这些离散点以外的其他点的函数值，就需要根据这些已知数据进行插值.而本题只提供了取样点和原函数.分析问题求解方法如下：（1）利用已知函数式计算取样点对应的函数值；将作为两个等长的已知向量，分别描述采样点和样本值.因此被插值函

9、数是一个单变量函数，可利用一维插值处理该数据插值问题.一维插值采用的方法通常有拉格朗日多项式插值（本题采用3次多项式插值），3次样条插值法和分段线性插值.（2）分别利用以上插值方法求插值.以0.5个单位为步长划分区间-6，6，并将每一点作为插值函数的取样点.再根据插值函数计算所选取样点的函数值.最后再利用所得函数值画出相应的函数图象，并与原函数的图象进行对比.3.3问题的假设 为了解决上述分析所提到的问题，本题可以作出如下假设：（1）假设原函数仅作为求解取样点对应的样点值的函数关系式.而其他各点的函数值都是未知量，叙用插值函数计算. （2）为了得到理想的对比函数图象，假设为已知的标准函数.可以

10、选取0.5个单位为步长划分区间-6，6，分别计算插值函数和标准函数在该区间的取样点的函数值.画出函数图象进行对比. 3.4分段线性插值原理给定区间, 将其分割成，已知函数在这些插值结点的函数值为；求一个分段函数,使其满足： (1) ，； (2) 在每个区间上, 是个一次函数.易知,是个折线函数, 在每个区间上,于是, 在上是连续的，但其一阶导数是不连续的.于是即可得到如下分段线性插值函数：,其中 3.5问题的求解在MATLAB中实现分段线性插值，最近点插值，3次多项式插值，3次样条插值的命令为interp1,其调用格式为: 1=interp1(,1，method)函数根据,的值，计算函数在1处

11、的值.,Y是两个等长的已知向量，分别描述采样点和样本值，1是一个向量或标量，描述欲插值点，1是一个与1等长的插值结果.method是插值方法，包括：linear：分段线性插值.它是把与插值点靠近的两个数据点用直线连接，然后在直线让选取对应插值点的数.nearest：近点插值法.根据已知两点间的插值点与这两点间的位置远近插值.当插值点距离前点远时,取前点的值,否则取后点的值.cubic：3次多项式插值.根据已知数据求出一个3次多项式，然后根据多项式进行插值.spline：3次样条插值.在每个分段（子区间）内构造一个3次多项式，使其插值函数除满足插值条件外，还要求个节点处具有光滑条件.再根据已知数据求出样条函数后,按照样条函数插值.运用Matlab工具软件编写代码，并分别画出图形如下：(一)在-6，6中平均选取5个点作插值：（二）在-6，6中平均选取11个点作插值：（三）在-6，6中平均选取21个点作插值：（四）在-6，6中平均选取41个点作插值3.6 分段插值方法的优劣性分析从以上对比函数图象可以看出，分段线性插值其总体光滑程度不够.在数学上，光滑程度的定量描述是函数(曲线) 的阶导数存在且连续，则称该曲线具有阶光滑性.一般情况下，阶数越高光滑程度越好.分段线性插值具有零阶光滑性，也就是不光滑.3次样条插值就是较低次数的多项式而达到较高阶光滑性的方法.总体上分段线性插值具有以下