五种插值法的对比研究Word文档格式.docx

1、最后使用 MATLAB 软件的编程实现，绘制出不同插值法下的函数曲线，从几何上再次进行对 比，得出结论。通过此次论文的写作，我对于插值法有了更深的理解和认知，对于 今后插值法的选择也会更加容易权衡把握。关键词： 插值法；对比；插值函数；多项式Abstract: Interpolation is one of the most basic methods in numerical analysis.There are many functions in practical problems,some give no expression,someonly supply discrete data

2、. So we only find it again from the adjacent number next to find its correct value and according to a certain relationship to the adjacent number corrected.The correct relationship is an interpolation in fact.In practical applications,the effect of approximation is also different when different inte

3、rpolation functions are used.We have contacted five basic interpolation methods,such as Lagrange interpolation,Newton interpolation, piecewise linear interpolation, piecewise three Hermite interpolation and spline interpolation function.Firstly,this paper introduces the basic concepts and ideas to s

4、olve problems of five kinds of interpolation methods.And then through the comparative analysis of the advantagesand disadvantagesof different interpolation methods in the process of solving typical problems.Finally,using MATLAB software programming,draw different interpolation method of function cur

5、ve,from geometry again contrast,draw conclusions.Through the writing of this paper,I have a deeper understanding and recognition of the interpolation method,and it will be easier to balance and select which interpolation methods to use in the future.Key Words: Interpolation method comparison interpo

6、lation function polynomial1 导言1.1 选题背景插值方法最早来源于生产实践，作为一种数学方法，其经历了漫长的历史考验 与证实。早在数千多年前，我们的祖先就凭借插值方法，利用已知的少部分日月五 星运行规律的观测值获得了相对较完整的运行规律。在一千多年前的隋唐时期，中 国的贤能之士就将插值技术应用到了制定历法的过程中。而到公元六世纪时，隋朝 的刘焯又把等距节点的二次插值应用于天文计算中。在 16-19 世纪，多项式插值被 用来解决航海学和天文学的一些重要问题。十七世纪时，牛顿（ Newto n）和格雷格里（Gregory）建立了等距结点上的一般插值公式，后来拉格朗日（

7、Lagrange）建立出了非等距结点插值公式。在微积分产生并且广泛应用之后，插值的基本理论和结 果随之有了进一步的完善，之后其应用也越来越广泛，尤其是在计算机普遍使用之 后，插值法在各领域中的地位也越来越重要，与此同时自身也得到了发展。经典的插值方法是基于泰勒插值（Taylor）和拉格朗日插值的，其实 Taylor插 值与拉格朗日插值的联系十分密切，即拉格朗日插值的极限形式可以视为 Taylor插值，反之， Taylor 插值的离散化形式就是拉格朗日插值。我们在建立拉格朗日插值 多项式时很是简单方便，但一旦节点增加，就不能再使用原来的多项式计算，需要 重新建立新的多项式，这无疑使计算变得繁琐起

8、来，而 Newt on （牛顿）插值就克服了这一问题。此外根据实际问题，插值法的应用在很多情况下都需要尽量满足插值 函数与原函数相差无异的前提，即要求在节点上插值函数与被插值函数的函数值和 导数值都是相等的，也就是另一种插值法， Hermite （埃尔米特）插值法。事实上，我们把Taylor插值和拉格朗日插值进行联系融合就能总结出 Hermite （埃尔米特）插值，这也推广了前两种插值法。现在，插值技术的应用在很多领域得到了普及，当我们需要认识某一事物的本 质时，常根据其观测点，利用插值技术对特定问题进行深入拓展和解决，以加深对 该事物的认识。多项式插值是函数插值中最常用的一种形式。在一般的插

9、值问题中，插值条件 可以唯一地确定一个次数不超过 n的插值多项式。从几何上可以解释为：可以从多 项式曲线中找出一些不超过n次的点通过平面上n 1个不同的点。插值多项式有两 种常用的表达式形式，一种是拉格朗日插值多项式，另一种是牛顿插值多项式，此 外拉格朗日插值公式与牛顿插值公式永远相等。此外，在进行高阶次插值时常常出现不稳定的情况，而采用样条插值和分段线 性插值法就可以防止这类情况的发生。分段线性插值或分段三次埃尔米特插值等此 种分段低次插值法可以使逼近效果加强，但却整体光滑而不收敛。为此，引入了更 理想化的三次样条插值法。1.2研究的目的和意义在数值分析中，对于插值函数的学习是必不可少的，因

10、为它能辅助我们把模糊 的数据准确化，把想当然的数据变得无懈可击。但是对于五种插值函数，他们具有 不同的优势和适用范围，五种方法对同一问题的处理的结果一定不同，这时对于方 法的选择显得至关重要。因此我们对于他们差异化的了解与认知是必不可少的。通过此篇论文的对比研究，我希望不但可以给数值分析领域中的学习者一些帮 助和启示甚至让他们在求知的路上少些磕绊，也能推动一些运用到插值函数知识的 社会工作领域的工作者的职业进步。2五种插值法2.1拉格朗日插值拉格朗日是n次多项式插值，解题方法是先构造插值基函数再求 n次插值多项式。对Lagrange n次插值多项式，首先要选取n 1个插值点x。/,xn上的n次

11、插值基函数，（X X0）.（X X）（x Xi 1）.（X Xn）li （X）=（Xi X）.（Xi Xi_1）（Xi Xi 十）（XiXn）（i =0,1,2,n）有了这n 1个n次插值基函数，就能很容易的写出n次Lagrange插值n多项式了，其具体的表达式为Ln（x） =f （人儿（x）。i =0拉格朗日插值原理：表1插值数值表XiX。X1x2 .Xnf （x）f （X0）f（xjf（X2）f（Xn）Lagrange插值的方法是：对于给定的n个插值节点乂0,人, xn和对应的函数值 y0, y1, y2,yn，我们利用n次Lagrange插值多项式，可以对插值区间上任意的x对 应的函数值

12、y利用下式Ln（x）来求解。表1中的n次Lagrange插值多项式Ln（x）的数学表达式为：中，（x） =（x x0）（x Xj.（x Xn）。2.2牛顿插值牛顿插值也是n次多项式插值，提出了构造插值多项式的另一种方法。它具有 继承性和易变化节点的特点。牛顿插值原理：Newton插值的方法：由表1构造的牛顿插值多项式为：N（X）= f （Xo） （X-Xo）fXo,Xi （X-XoXx-XjflXo,% .（x-Xo） . X（XnjL）fX0,X1, .）Cn,用上式插值时，首先要计算各阶差商，而各阶差商的计算可以归纳为一阶差商的逐次计算，一般的 fXo,XL，Xn = fXo，XlXk2X

13、k fX0,Xi,.,X2Xk Xk余项为：Rn（x）二 f （x） -N（x）二 fx,x0,xL，,xn （x），其中 ，（x） =（x xo）（x XL）.（X xn）2.3分段线性插值分段线性插值的意义在于克服拉格朗日插值法的非收敛性。其实分段线性插值就是利用每两个相邻的插值基点做线性插值，就可以得到分段线性插值函数：y（xli （x）fi li i（x）fi i，x Xi,，（i =0,1,2,n）设分段线性插值函数为y（x），则具有以下性质：1y（x）可以分段表示并且yi（x）在每个小区间4,Xi上都是线性函数;2yi （x）二 f （xj = fi， （i =0,1,2,n）；3

14、y（x）在整个区间a,b上连续3 0特点：插值函数的序列具有一致的收敛性，弥补了高阶拉格朗日插值方法的不足，可是存在插值精度低、基点处不光滑的缺陷，其中增加插值点可以提高插值精 度。几何上，分段线性插值是通过顺次连接各插值点形成线段，从而逼近原始曲线， 这也是计算机绘图的基本原理。2.4分段三次Hermite插值对于函数f（x），有时我们不仅知道它在一些点处的函数值，而且还能知道它在 这些点的导数值。当在这些点上的插值函数P（x）的函数值和导数值同时满足与f（x） 的函数值和导数值相等的要求时，此时的问题就是 Hermite插值问题或带有导数的插值问题。假定已知函数f（x）在插值区间p,q上的

15、n 1个互不相同的节点Xi （i =0,1，,n）处满足f（x）二fi及f（X）二fi （i =0,1,2,.,n），如果函数G（x）的存在满足下列条件：1G（x）在每个小区间上的多项式次数为 3；2G（x） C1a,b；3G（xJ = f（Xi），Gg = f（x），i=（0,1，,n）5就称G（x）是f （x）在n 1个节点xi上的分段三次埃尔米特插值多项式。所以，G（x）二hky。）hkk1（x）仏的丫宀二（1 2 乂-兀）（x-Xk 彳）2 yk .（1 . 2 x-x-1）（ x Xk）2 yk 1Xk 卅Xk Xkxk 卅 Xk Xk 卅 xk 卅xk（X-Xk）（x xk1）2y

16、k （X-Xk1）（ X x ）2yk 1Xk Xk* Xk 岀Xk2.5样条插值函数2.5.1样条插值的相关概念分段低次插值函数，虽然有收敛性，但平整度差。因此，早期的制图工程师在制图时首先会在样点处固定弹性木条，其他各处任意成形，这样就能画出一条曲线， 定义样条曲线。事实上，该曲线是由分段三次曲线并接而成，在连接点也就是样点 上必须要二阶连续可导，从数学角度加以归纳得到数学样条这个概念。利用样条插值方法得到的插值曲线光滑性好，但却不收敛。由此我们可以引用 三次样条函数以达到插值函数的收敛性且光滑度也更好了。2.5.2三次样条插值函数对于给定区间p, q上p = X。： x : .： Xn二

17、q这n 1个节点和在这些点上的函数值 f （xj =yj（i =0,1,.,n），若函数 g（x）满足：1在每个子区间Xi,Xi（i =1,2,., n）上，多项式g（x）的次数不超过3；2g（x）,g （x）,g （x）在p,q上连续；3满足g（xj =yi（i =0,1,.,n）的插值条件。则g（x）是函数f （x）关于n个节点xi处的三次样条插值函数。3五种插值法的对比研究3.1五种插值法的解题分析比较例1已知表2x11/2y-2 e-4 e请写出在以上3个节点处的牛顿插值（一次和二次）以及拉格朗日插值。解：（1）拉格朗日型插值多项式 构造过（0,1）（1,eJ的一次插值基函数X 为lo

18、（X）二 1 （X-1）Xo Xili（x）=L=xxi xo则一次插值多项式为:i（x）二 yl（x） yi（x） 一（x-1） xe，构造过x0, x2的二次插值基函数因此二次插值多项式为:2（x）二 yolo（x） yili（x） y2l2（x）1i= 2（x-i）（x ） 2x（x ）e -4x（x-i）e22（2）牛顿型插值多项式构造牛顿一次插值函数:因为fx0,Xi（Xi） f（X0）=e_iXi -Xo所以构造牛顿二次插值函数:1（X）二 f （Xo） （X Xo）fXo,Xi =1 x（e，-1）于是2（x） = f（Xo） （X-Xo）fXo，Xi （x-XoXx-XJflX

19、oMK=i x（e，-i） x（x-i）（2 2e，-4e）综上，由拉格朗日公式n（x） =7 yjlj（x），牛顿公式j=0n（X）二 f（Xo）（X -X）f X,Xi . （X Xo）.（X Xn）fXo,Xi,.,Xn及例题可以看出：（i）拉格朗日插值法优势：公式的结构整齐紧密，对于理论研究分析非常方便;缺点：当增加或减少一个插值点的计算，将需要重新计算相应的插值基函数, 然后插值多项式的公式代入结果也会改变，大大增加了计算量，解题十分繁琐。此 外，当插值点很多时，拉格朗日多项式的插值次数也会很高，使计算结果的值变得 动荡。换言之，即使在已知的几个点处得到正确的结果，但在附近的点处“事

20、实上” 的值和得到的结果之间的会有较大的差距。（2）牛顿插值法牛顿插值法的公式是另一种 n次插值多项式的构造形式，然而它却克服 了拉格朗日插值多项式的缺陷，它的一个显著优势就是每当增加一个插值节点，只 要在原牛顿插值公式中增加一项就可形成高一次的插值公式。此外，如果在实际应 用中遇到等距分布的插值节点，牛顿插值公式就能得到进一步的简化，从而得到等 距节点的插值公式，这样为缩短实际运算时间做出了很大的贡献。这种插值仅仅要求插值多项式在插值节点处与被插函数有相等的函数值， 而这种插值多项式却不能全面反映被插值函数的性态。然而在许多实际问题中，不 仅要求插值函数与被插值函数在所有节点处有相同的函数值

21、，它也需要在一个或全 部节点上插值多项式与被插函数有相同的低阶甚至高阶的导数值。对于这些情况， 拉格朗日插值和牛顿插值都不能满足。1 1例2 过0,1两点并且满足f（0） =1f （0） = , f=2, f=，构造一个三次埃尔2 2米特插值多项式。利用公式有ho（x） =（1 2Xo）（Xl ）2 =（1 2x）（x-1）2Xi -Xo X治h1（x）=（1 2）（）2 =（3-2x）x2Xo X1 x -xox 一 论 2 2Ho（x） =（x Xo）（ L） =x（x 1）Xo X1Hx） =（x_xj（ X X。）2 =x2（x -1）H （x） = yoho（x） yh（x） yoH

22、o（x） y；Hi（x）2 2 1 2 1 2=（1 2x）（x -1）2 2x2（3-2x） - x（x -1）2 x2（x-1）-x3 3x2 -x 1由这个例题2可以看出：对于埃尔米特插值，我们不仅已知函数在某些点处的 函数值，而且插值函数在这些点处的导数与被插函数相同。因此，（-优点：关于插值函数和被插函数的贴合程度，埃尔米特插值比多项 式的好。（2）缺点：埃尔米特插值只有在被插值函数在插值节点处的函数值和导 数值已知时才可以使用，而这在实际问题中是无法实现的，因为在一般情况下我们 是不可能也没必要知道函数在插值节点处的导数值。因此成为能否运用埃尔米特插 值的一个重要因素就是：我们知不

23、知道插值函数在节点处的导数值。例3对于函数f（x）=21 25x取等距节点x i（i =0,1,.,10），建立插值多项式1o（x），并探究它与f（x）的误10差。 根据题意知道多项式的次数为10,代入拉格朗日插值多项式的公式有10（X）八 f（XJIj（X）其中Xi - -1 i, i = 0,1,.,10 （x - X0 ）.（ X - x9 ）（ X -为0 ）（x Xo）.（Xi X9）（Xi X10）计算结果如下表所示:表310（Xi ）10 （Xi）f （Xi ） 21 + 25Xi1 +25Xj-1.000.03846-0.400.200000.19999-0.900.04706

24、1.57872-0.300.307690.23535-0.800.05882-0.200.50000-0.700.07547-0.22620-0.100.800000.84340-0.600.100000.001.00000-0.500.137930.25376对于0,1区间上的值可以由对称性得到，根据结果可以看出， 10（x）在原点附近能较好的逼近f（x），而在其余点处：10（x）与f（x）的差异较大，越靠近端点，逼近 效果就越不好。由例题3可以不难发现，在高次插值中拉格朗日插值多项式存在较大缺陷，因而为了弥补这种不足我们一般利用分段线性插值的方法。1 2例4给定函数y - 2，1込x1取等距节点人=T i（i = 0,1,.,10），作分段1+25x 10线性插值函数（x），并计算（0.9）的值解:首先计算出-1,0区间上的函数值表:表4-1-0.8-0.6-0.4-0.2对于区间0,1上的函数值可由对称性得到其次，构造各点的插值基函数:0 1EXE0.8、5（x0.8） 0.8 ex 兰1故得到分段线性插值函数 （x）（x） =0.03846（1。（x） l10（x） 0.05882（11（x） l9（x）