族谱网 头条 人物百科

拉格朗日插值法

2020-10-16
出处:族谱网
作者:阿族小谱
浏览:561
转发:0
评论:0
定义对某个多项式函数,已知有给定的k+1个取值点:其中xj{displaystylex_{j}}对应着自变量的位置,而yj{displaystyley_{j}}对应着函数在

定义

对某个多项式函数,已知有给定的k + 1个取值点:

其中xj{\displaystyle x_{j}}对应着自变量的位置,而yj{\displaystyle y_{j}}对应着函数在这个位置的取值。

假设任意两个不同的xj都互不相同,那么应用拉格朗日插值公式所得到的拉格朗日插值多项式为:

其中每个ℓ ℓ -->j(x){\displaystyle \ell _{j}(x)}为拉格朗日基本多项式(或称插值基函数),其表达式为:

拉格朗日基本多项式ℓ ℓ -->j(x){\displaystyle \ell _{j}(x)}的特点是在xj{\displaystyle x_{j}}上取值为1,在其它的点xi,i≠ ≠ -->j{\displaystyle x_{i},\,i\neq j}上取值为0。

范例

假设有某个二次多项式函数f{\displaystyle f},已知它在三个点上的取值为:

f(4)=10{\displaystyle f(4)=10}

f(5)=5.25{\displaystyle f(5)=5.25}

f(6)=1{\displaystyle f(6)=1}

要求f(18){\displaystyle f(18)}的值。

首先写出每个拉格朗日基本多项式:

然后应用拉格朗日插值法,就可以得到p{\displaystyle p}的表达式(p{\displaystyle p}为函数f{\displaystyle f}的插值函数):

此时代入数值 18{\displaystyle \ 18}就可以求出所需之值: f(18)=p(18)=− − -->11{\displaystyle \ f(18)=p(18)=-11}。

证明

存在性

对于给定的k+1个点:(x0,y0),… … -->,(xk,yk){\displaystyle (x_{0},y_{0}),\ldots ,(x_{k},y_{k})},拉格朗日插值法的思路是找到一个在一点xj{\displaystyle x_{j}}取值为1,而在其他点取值都是0的多项式ℓ ℓ -->j(x){\displaystyle \ell _{j}(x)}。这样,多项式yjℓ ℓ -->j(x){\displaystyle y_{j}\ell _{j}(x)}在点xj{\displaystyle x_{j}}取值为yj{\displaystyle y_{j}},而在其他点取值都是0。而多项式L(x):=∑ ∑ -->j=0kyjℓ ℓ -->j(x){\displaystyle L(x):=\sum _{j=0}^{k}y_{j}\ell _{j}(x)}就可以满足

在其它点取值为0的多项式容易找到,例如:

它在点xj{\displaystyle x_{j}}取值为:(xj− − -->x0)⋯ ⋯ -->(xj− − -->xj− − -->1)(xj− − -->xj+1)⋯ ⋯ -->(xj− − -->xk){\displaystyle (x_{j}-x_{0})\cdots (x_{j}-x_{j-1})(x_{j}-x_{j+1})\cdots (x_{j}-x_{k})}。由于已经假定xi{\displaystyle x_{i}}两两互不相同,因此上面的取值不等于0。于是,将多项式除以这个取值,就得到一个满足“在xj{\displaystyle x_{j}}取值为1,而在其他点取值都是0的多项式”:

这就是拉格朗日基本多项式。

唯一性

次数不超过k的拉格朗日多项式至多只有一个,因为对任意两个次数不超过k的拉格朗日多项式:P1{\displaystyle P_{1}}和P2{\displaystyle P_{2}},它们的差P1− − -->P2{\displaystyle P_{1}-P_{2}}在所有k+1个点上取值都是0,因此必然是多项式(x− − -->x0)(x− − -->x1)⋯ ⋯ -->(x− − -->xk){\displaystyle (x-x_{0})(x-x_{1})\cdots (x-x_{k})}的倍数。因此,如果这个差P1− − -->P2{\displaystyle P_{1}-P_{2}}不等于0,次数就一定不小于k+1。但是P1− − -->P2{\displaystyle P_{1}-P_{2}}是两个次数不超过k的多项式之差,它的次数也不超过k。所以P1− − -->P2=0{\displaystyle P_{1}-P_{2}=0},也就是说P1=P2{\displaystyle P_{1}=P_{2}}。这样就证明了唯一性。

几何性质

拉格朗日插值法中用到的拉格朗日基本多项式ℓ ℓ -->0,ℓ ℓ -->1,⋯ ⋯ -->,ℓ ℓ -->n{\displaystyle \ell _{0},\ell _{1},\cdots ,\ell _{n}}(由某一组x0<x1<xn{\displaystyle x_{0}确定)可以看做是由次数不超过n的多项式所组成的线性空间:Kn[X]{\displaystyle \mathbb {K}系数{n}[X]}的一组基底。首先,如果存在一组系数:λ λ -->0,λ λ -->1,⋯ ⋯ -->,λ λ -->n{\displaystyle \lambda _{0},\lambda _{1},\cdots ,\lambda _{n}}使得,

那么,一方面多项式P是满足P(x0)=λ λ -->0,P(x1)=λ λ -->1,⋯ ⋯ -->,P(xn)=λ λ -->n{\displaystyle P(x_{0})=\lambda _{0},P(x_{1})=\lambda _{1},\cdots ,P(x_{n})=\lambda _{n}}的拉格朗日插值多项式,另一方面P是零多项式,所以取值永远是0。所以

这证明了ℓ ℓ -->0,ℓ ℓ -->1,⋯ ⋯ -->,ℓ ℓ -->n{\displaystyle \ell _{0},\ell _{1},\cdots ,\ell _{n}}是线性无关的。同时它一共包含n+1个多项式,恰好等于Kn[X]{\displaystyle \mathbb {K} _{n}[X]}的维数。所以ℓ ℓ -->0,ℓ ℓ -->1,⋯ ⋯ -->,ℓ ℓ -->n{\displaystyle \ell _{0},\ell _{1},\cdots ,\ell _{n}}构成了Kn[X]{\displaystyle \mathbb {K} _{n}[X]}的一组基底。

拉格朗日基本多项式作为基底的好处是所有的多项式都是齐次的(都是n次多项式)。

优点与缺点

拉格朗日插值法的公式结构整齐紧凑,在理论分析中十分方便,然而在计算中,当插值点增加或减少一个时,所对应的基本多项式就需要全部重新计算,于是整个公式都会变化,非常繁琐。这时可以用重心拉格朗日插值法或牛顿插值法来代替。此外,当插值点比较多的时候,拉格朗日插值多项式的次数可能会很高,因此具有数值不稳定的特点,也就是说尽管在已知的几个点取到给定的数值,但在附近却会和“实际上”的值之间有很大的偏差(如右下图)。这类现象也被称为龙格现象,解决的办法是分段用较低次数的插值多项式。

重心拉格朗日插值法

重心拉格朗日插值法是拉格朗日插值法的一种改进。在拉格朗日插值法中,运用多项式

拉格朗日插值法

拉格朗日插值法的数值稳定性:如图,用于模拟一个十分平稳的函数时,插值多项式的取值可能会突然出现一个大的偏差(图中的14至15中间)

可以将拉格朗日基本多项式重新写为:

定义重心权

上面的表达式可以简化为:

于是拉格朗日插值多项式变为:

即所谓的重心拉格朗日插值公式(第一型)或改进拉格朗日插值公式。它的优点是当插值点的个数增加一个时,将每个wj{\displaystyle w_{j}}都除以(xj− − -->xk+1){\displaystyle (x_{j}-x_{k+1})},就可以得到新的重心权wk+1{\displaystyle w_{k+1}},计算复杂度为O(n){\displaystyle {\mathcal {O}}(n)},比重新计算每个基本多项式所需要的复杂度O(n2){\displaystyle {\mathcal {O}}(n^{2})}降了一个量级。

将以上的拉格朗日插值多项式用来对函数g(x)≡ ≡ -->1{\displaystyle g(x)\equiv 1}插值,可以得到:

因为g(x)≡ ≡ -->1{\displaystyle g(x)\equiv 1}是一个多项式。

因此,将L(x){\displaystyle L(x)}除以g(x){\displaystyle g(x)}后可得到:

这个公式被称为重心拉格朗日插值公式(第二型)或真正的重心拉格朗日插值公式。它继承了(1)式容易计算的特点,并且在代入x值计算L(x){\displaystyle L(x)}的时候不必计算多项式ℓ ℓ -->(x){\displaystyle \ell (x)}。它的另一个优点是切比雪夫比雪夫节点进行插值的话,可以很好地模拟给定的函数,使得插值点个数趋于无穷时,最大偏差趋于零。同时,重心拉格朗日插值结合切比雪夫节点进行插值可以达到极佳的数值稳定性。第一型拉格朗日插值是向后稳定的,而第二型拉格朗日插值是向前稳定的,并且勒贝格常数很小。

参考文献

来源

李庆扬、王能超、易大义. 《数值分析》 第4版. 清华大学出版社. 2001. ISBN 7-302-04561-5 (中文(简体)‎). 

冯有前. 《数值分析》. 清华大学出版社. 2001. ISBN 7-810-82495-3 (中文(简体)‎). 

拉格朗日插值多项式. 太原理工大学 (中文(简体)‎). 

参见

约瑟夫·拉格朗日

多项式插值

样条插值

伯恩斯坦多项式

牛顿多项式

龙格现象


免责声明:以上内容版权归原作者所有,如有侵犯您的原创版权请告知,我们将尽快删除相关内容。感谢每一位辛勤著写的作者,感谢每一位的分享。

——— 没有了 ———
编辑:阿族小谱

相关资料

展开
发表评论
写好了,提交
{{item.label}}
{{commentTotal}}条评论
{{item.userName}}
发布时间:{{item.time}}
{{item.content}}
回复
举报
点击加载更多
打赏作者
“感谢您的打赏,我会更努力的创作”
— 请选择您要打赏的金额 —
{{item.label}}
{{item.label}}
打赏成功!
“感谢您的打赏,我会更努力的创作”
返回

更多文章

更多精彩文章
打赏
私信

推荐阅读

· 拉格朗日中值定理
内容拉格朗日中值定理的几何意义文字叙述如果函数f(x){\displaystylef(x)}满足:在闭区间[a,b]{\displaystyle[a,b]}上连续;在开区间(a,b){\displaystyle(a,b)}内可导;那么在(a,b){\displaystyle(a,b)}内至少有一点ξξ-->(a<b){\displaystyle\xi(a,使等式成立。逻辑语言的叙述若函数f(x){\displaystylef(x)}满足:f(x)∈∈-->C[a,b]{\displaystylef(x)\inC[a,b]};f(x)∈∈-->D(a,b){\displaystylef(x)\inD(a,b)};则∃∃-->ξξ-->∈∈-->(a,b),s.t{\displaystyle\exists\xi\in(a,b),s.t}证明令g(x)=...
· 洛朗·拉福格
生平拉福格1986年入读巴黎高等师范学院。他在巴黎第十一大学奥赛数学实验室算术及代数几何小组受GérardLaumon指导完成博士论文,也在那里做初级研究员。现时在法国国家科学研究中心任高级研究员至今。2000年起他担任法国高等科学研究所数学教授。2002年在中国北京举行的第24届国际数学家大会上,他获颁菲尔兹奖,表彰他对数论和代数几何贡献突出。他推广同是菲尔兹奖得主乌克兰数学家弗拉基米尔·德林费尔德的方法,证明朗兰兹猜想在特征为正的代数曲线函数域的GLn{\displaystyleGL_{n}}上成立。2003年11月18日,他成为法兰西科学院院士。对教育的关注2004年,他开始关注法国教育体系,接触教育关注组织“拯救语文教育”(Sauverleslettres)。他和阿兰·孔涅和其他科学家共同署名发表文章《为科学和技术的未来所需的基础知识:如何...
· 欧拉-拉格朗日方程
第一方程设f=f(x,y,z){\displaystylef=f(x,\y,\z)},以及fy,fz{\displaystylef_{y},\f_{z}}在[a,b]××-->R2{\displaystyle[a,\b]\times\mathbb{R}^{2}}中连续,并设泛函若y∈∈-->C1[a,b]{\displaystyley\inC^{1}[a,\b]}使得泛函J(y){\displaystyleJ(y)}取得局部平稳值,则对于所有的x∈∈-->(a,b){\displaystylex\in(a,\b)},推广到多维的情况,记若y→→-->′(x)∈∈-->(C1[a,b])n{\displaystyle{\vec{y}}"(x)\in(C^{1}[a,b])^{n}}使得泛函J(y→→-->)=∫∫-->...
· 插值
定义给定n{\displaystylen}个离散数据点(称为节点)(xk,yk){\displaystyle(x_{k},y_{k})},k=1,2,...,n{\displaystylek=1,2,...,n}。对于x,(x≠≠-->xk,k=1,2,...n){\displaystylex,(x\neqx_{k},k=1,2,...n)},求x{\displaystylex}所对应的y{\displaystyley}的值称为内插。f(x){\displaystylef(x)}为定义在区间[a,b]{\displaystyle[a,b]}上的函数。x1,x2,x3...xn{\displaystylex_{1},x_{2},x_{3}...x_{n}}为[a,b]{\displaystyle[a,b]}上n个互不相同的点,G{\displaystyleG}为给定的某一函数类。若G{...
· 拉格朗日点
位置五个拉格朗日点之定义及位置如下:L1在M1和M2两个大天体的连线上,且在它们之间。例如一个围绕太阳旋转的物体,它距太阳的距离越近,它的轨道周期就越短。但是这忽略地球的万有引力对其产生的拉力的影响。如果这个物体在地球与太阳之间,地球引力的影响会减弱太阳对这物体的拉力,因此增加这个物体的轨道周期。物体距地球越近,这种影响就越大。在L1点,物体的轨道周期恰好等于地球的轨道周期。太阳及日光层探测仪(SOHO)即在日-地系统的L1点上运行。L2在两个大天体的连线上,且在较小的天体一侧。日地系统的L2在地球远离太阳的一侧。一般来讲,一个物体距太阳的距离越远,它的轨道周期通常就越长,但L2点上的物体还受到地球的引力,所以轨道周期变地与地球的相等。日地系统的L2通常用于放置空间天文台。因为L2上的物体可以保持背向太阳和地球的方位,易于保护和校准。威尔金森微波各向异性探测器已经在日-地系统的L2点上运行...

关于我们

关注族谱网 微信公众号,每日及时查看相关推荐,订阅互动等。

APP下载

下载族谱APP 微信公众号,每日及时查看
扫一扫添加客服微信