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

浅析拉格朗日插值法

目录:

一、 引言

二、 插值及多项式插值的介绍 三、 拉格朗日插值的理论及实验

四、 拉格朗日插值多项式的截断误差及实用估计式 五、 参考文献

一、引言

插值在数学发展史上是个古老问题。插值是和拉格朗日(Lagrange)、牛顿(Newton)、高斯(Gauss)等著名数学家的名字连在一起的。在科学研究和日常生活中,常常会遇到计算函数值等一类问题。插值法有很丰富的历史渊源,它最初来源人们对天体研究——有若干观测点(我们称为节点)计算任意时刻星球的位置(插值点和插值)。现在,人们在诸如机械加工等工程技术和数据处理等科研都有很好的应用,最常见的应用就是气象预报。插值理论和方法能解决在实际中当许多函数表达式未知或形式复杂,如何去构造近似表达式及求得在其他节点处的值的问题。

二、插值及多项式插值

1、插值问题的描述

设已知某函数关系y?f(x)在某些离散点上的函数值:

xx0x1yy0y1????xn?1yn?1xnyn

插值问题:根据这些已知数据来构造函数y?f(x)的一种简单的近似表达式,以便于计算点x?xi,i?0,1,?,n的函数值f(x),或计算函数的一阶、二阶导数值。

2、插值的几何意义

插值的几何意义如图1所示:

图1 3、多项式插值 3.1 基本概念

假设y?f(x)是定义在区间??a,b??上的未知或复杂函数,但一直该函数在点a?x0?x1???xn?b处的函数值y0,y1,?yn。找一个简单的函数,例如函数

P(x),使之满足条件

P(x)?iy,? (3.1) i0,1?,2, n ,通常把上述x0?x1???xn 称为插值节点,把P(x)称为f(x)的插值多项式,条件(3.1)称为插值条件,并把求P(x)的过程称为插值法。 3.2 插值多项式的存在性和唯一性 如果插值函数是如下m次的多项式:

mm?1P??am?1x?am m(x)?a0x?a1x那么插值函数的构造就是要确定Pm(x)表达式中的m+1个系数

a0,a1,?am?1a,m。由于插值条件包含n+1独立式,只要m=n就可证明插值函数多项式是唯一存在。

实际上,由n+1个插值条件可得

nn?1?a0x0?a1x0??an?1x0?an?y0?nn?1?a0x1?a1x1??an?1x1?an?y1?

??nn?1??a0xn?a1xn??an?1xn?an?yn 这是一个关于a0,a1,?an的n+1阶线性方程组,且其系数矩阵对应的行列式是线性代数中著名的范德蒙(Vandemonde)行列式。该行列式得值为 Vn(x0,x1,?xn)???(xi?xj)

i?1j?0ni因为i?j时,xi?xj,所以Vn(x0,x1,?xn)?0。从而证明了上述线性方程组的阶是唯一存在的。既满足插值条件的多项式唯一存在。

三、 拉格朗日插值的理论及实验 1、拉格朗日插值的理论

拉格朗日(Lagrange)插值公式的基本思想是把Pn(x)的构造问题转化为n+1个插值基函数li(x)(i?0,1,?,n)。首先我们利用节点直接构造如下多项式:

ln(x)?其中

?n?1(x) '(x?xi)?n?1(x)i ?n?1(x)?(x?x0)(x?x1)?(x?xn),

' ?n?1(x)?(xi?x0)?(xi?xi?1)(xi?xi?1)?(xi?xn)

容易验证该多项式具有性质

?0,j?ili??

?1,j?i

因此,n次多项式

Ln(x)?l0(x)y0?l1(x)y1??ln(x)yn??lk(x)yk

k?0n一定具有性质

联系客服:779662525#qq.com(#替换为@) 苏ICP备20003344号-4