拉格朗日插值法理论误差分析.doc

该【拉格朗日插值法理论误差分析 】是由【雨林书屋】上传分享，beplayapp体育下载一共【7】页，该beplayapp体育下载可以免费在线阅读，需要了解更多关于【拉格朗日插值法理论误差分析 】的内容，可以使用beplayapp体育下载的站内搜索功能，选择自己适合的beplayapp体育下载，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此beplayapp体育下载到您的设备，方便您编辑和打印。浅析拉格朗日插值法目录:一、前言二、插值及多项式插值的介绍三、拉格朗日插值的理论及实验四、拉格朗日插值多项式的截断偏差及适用预计式五、参照文件一、前言插值在数学发展史上是个古老问题。插值是和拉格朗日(Lagrange)、牛顿Newton)、高斯(Gauss)等有名数学家的名字连在一同的。在科学研究和日常生活中,经常会碰到计算函数值等一类问题。插值法有很丰富的历史渊源,它最先根源人们对天体研究——有若干观察点(我们称为节点)计算随意时刻星球的地点(插值点和插值)。此刻,人们在诸如机械加工等工程技术和数据办理等科研都有很好的应用,最常有的应用就是气象预告。插值理论和方法能解决在实质中当很多函数表达式未知或形式复杂,怎样去结构近似表达式及求得在其余节点处的值的问题。二、插值及多项式插值1、插值问题的描绘设已知某函数关系yf(x)在某些失散点上的函数值:x0x1yy0y1LLxyn1xnn1yn插值问题:依据这些已知数据来结构函数yf(x)的一种简单的近似表达式,以便于计算点xxi,i0,1,L,n的函数值f(x),或计算函数的一阶、二阶导数值。2、插值的几何意义插值的几何意义如图1所示:图13、(x)是定义在区间a,b上的未知或复杂函数,但向来该函数在点ax0x1Lxnb处的函数值y0,y1,Lyn。找一个简单的函数,比如函数P(x),使之知足条件P(x)yi,i0,1,2,L,n,()往常把上述x0x1Lxn称为插值节点,把P(x)称为f(x)的插值多项式,条件()称为插值条件,并把求P(x)的过程称为插值法。:Pm(x)a0xma1xm1Lam1xam那么插值函数的结构就是要确定Pm(x)表达式中的m+1个系数L独立式,只需m=n便可证明插值函数多a0,a1,am1,am。因为插值条件包括n+1项式是独一存在。实质上,由n+1个插值条件可得naxn1Laxaya0x0001n10na0x1na1x1n1Lan1x1any1Ma0xnna1xnn1Lan1xnanyn这是一个对于a0,a1,Lan的n+1阶线性方程组,且其系数矩阵对应的队列式是线性代数中有名的范德蒙(Vandemonde)队列式。该队列式得值为niVn(x0,x1,Lxn)(xixj)i1j0因为ij时,xixj,所以Vn(x0,x1,Lxn)0。进而证了然上述线性方程组的阶是独一存在的。既知足插值条件的多项式独一存在。三、拉格朗日插值的理论及实验1、拉格朗日插值的理论拉格朗日(Lagrange)插值公式的基本思想是把Pn(x)的结构问题转变为n+1个插值基函数li(x)(i0,1,L,n)。第一我们利用节点直接结构以下多项式:ln(x)n1(x)(xxi)n'1(xi)此中n1(x)(xx0)(xx1)L(xxn),n'1(x)(xix0)L(xixi1)(xixi1)L(xixn)简单考证该多项式拥有性质li0,ji1,ji所以,n次多项式nLn(x)l0(x)y0l1(x)y1Lln(x)ynlk(x)ykk0必定拥有性质nLn(xi)lk(x)ykli(xi)yi,i0,1,L,n,k0既知足插值条件。我们称Ln(x)为拉格朗日插值多项式,li(x)称为拉格朗日插值及函数。一次拉格朗日插值多项式又叫做线性插值多项式。二次拉格朗日插值多项式又叫做抛物线插值多项式。2、拉格朗日插值实验经过学****掌握拉格朗日插值的理论,学致使用,使学到的知识运用到现实生活中,并运用计算机来解决我们在学****中碰到的一些问题。以下为运用MATLAB软件平台上计行拉格朗日插值问题:x024681012141618202224262830y例:已知在[0,30]内对应的节点x以及函数值y如表所示,利用拉格朗日插值多项式求在区间x=,x=,x=,x=。在已知数表函数的条件下,拉格朗日插值多项式可用来计算复杂函数或未知函数的函数值,为此我们第一编写以下利用拉格朗日插值多项式方法计算函数值的程序:functiony=lagrange(x0,y0,x)n=length(x0);m=length(x);fori=1:mz=x(i);s=;fork=1:np=;forj=1:nifj~=kp=p*(z-x0(j))/(x0(k)-x0(j));endends=p*y0(k)+s;endy(i)=s;end上述三重循环给出了拉格朗日插值计算多项式计算任何点x处的函数值的过程,,接下来我们在MATLAB平台长进行上述例子中的数值试验。在CommandWindow中输入的命令及结果以下所示:x=[0:2:30];y=[];lagrange(x,y,)ans=(x,y,)ans=(x,y,)ans=(x,y,)ans=,我们依据拉格朗日插值结果,利用plot命令画出未知函数的图像,命令程序以下:x0=[0:2:30];y0=lagrange(x,y,x0);plot(x0,y0)获得的未知函数图像为:四、拉格朗日插值多项式的截断偏差及适用预计式1、截断偏差在[a,b]区间上用Ln(x)近似未知或复杂函数f(x),其截断偏差是指RnxfxLnx()往常称Rnx为拉格朗日插值余额。注意到利用公式()预计截断偏差实质上特别困难。一是因为它要计算函数f(x)的高阶导数,当f(x)很复杂时,计算量很大,而当f(x)没有可用来计算的表达式时,导数没法正确计算;二是因为即便能获得高阶导数的分析式,但由于的详细地点不知道,所以要预计高阶导数在插值区间上的界一般是特别困难的事情。所以,公式()其实不适用。2、截断偏差的适用预计式既然公式()预计偏差时不适用,那么实质中怎样预计截断偏差呢?假定插值条件中包括n+2组数据f(xi)yi,i0,1,L,n,n1,那么利用n+1组数据我们能够结构一个n次拉格朗日插值多项式Ln(x),利用后n+1组数据我们能够结构另一个n次拉格朗日插值多项式*。利用公Ln(x)式()知,他们各自的插值余项为f(x)Ln(x)1f(n1)()(xx0)(xx1)L(xxn),(n1)!f(x)L*n(x)1f(n1)(*)(xx1)(xx2)L(xxn1),(n1)!两式相减得L*n(x)Ln(x)1fn1()(xx1)L(xxn)(xn1x0),(n1)!并可写成1f(n1)()(xx1)L(xxn)L*n(x)Ln(x).()(n1)!xn1x0注意到上式中利用fn1()fn1(*).该条件在好多状况下是建立的。利用式()可得Rn(x)f(x)Ln(x)Ln(x)L*n(x),x0xn1R*(x)f(x)L*(x)L*n(x)Ln(x),nnxn1x0)式()给出了用Ln(x)或L*n(x)作近似计算时的适用偏差预计式,它不需要计算高阶导数,也不用预计插值区间上高阶导数的界。总之,拉格朗日插值法的公式结构紧凑,在理论剖析中十分方便,但是在计算中,但插值点增添或减少时,所对应的基本多项式就得从头计算并且图像发生很大变化。像逐次线性插值法、牛顿插值法等都是在拉格朗日插值多项式的基础上延长出来的。我们依据实质中的详细问题,为减少插值偏差来选用相应的插值法来迅速的解决问题。五、参照文件数值计算原理李庆扬,关冶,,(第三版)颜庆津北京航空航天大学第一版社