关于埃尔米特插值多项式逼近介绍

关于埃尔米特插值多项式逼近介绍,第1张

关于埃尔米特插值多项式逼近介绍

[拼音]:Ai’ermite chazhi duoxiangshi bijin

[外文]:approximation by Hermite interpolation polynomials

埃尔米特插值是一种常见的插值方法。假设在区间[αb]上给定了n个互不相同的点x1,x2,…,xn以及一张数表

(*)

m=α1+α2+…+αn。早在 1878年C.埃尔米特就证明:存在惟一的次数不高于m-1的代数多项式Hn(x),使得

Hn(x)为表(*)的以 为结点组的埃尔米特插值多项式。如果定义在[αb]上的函数 ƒ(x)在xk(k=1,2,…,n)处有αk-1阶导数,并取,则称相应的Hn(x)为ƒ(x)的以为结点组的(α1,α2,…,αn)阶埃尔米特插值多项式。作为特殊情况,若诸αk都为1,则Hn(x)就是ƒ(x)的拉格朗日插值多项式;若n=1,则Hn(x)为ƒ(x)的α1-1阶泰勒多项式。最使人们注意的是诸αk都为2的情况,这时Hn(x)为次数不高于2n-1的代数多项式。如果写

Hn(x)可表示为

在这种情况下,常取,而给以适当的限制。这个想法大致起源于拉格朗日插值多项式的研究。为了改善插值多项式的逼近度,需对其导数作一定的要求。

为了简单,考虑定义区间为[-1,1]的情况。L.费耶尔首先让,称

为函数ƒ(x)的埃尔米特-费耶尔插值多项式。如果取切比雪夫多项式Tn(x)=cos(n arc cos x) 的零点全体为结点组, 则有绝对常数с,使得对于[-1,1] 上的任一连续函数ƒ(x)都有

式中-1≤x ≤1,ω(ƒ,δ)为ƒ(x)的连续性模。然而,用ƒn(ƒx)逼近ƒ(x)有其饱和性,逼近阶最多为1/n。若

关于[-1,1]上的x均匀成立,则ƒ(x)是个常数。但是对于其他结点组,会有较大的差异。例如,取勒让德多项式

的零点全体为结点组时,对于[-1,1]上的连续函数ƒ(x),相应的ƒn(ƒx)仅可能在(-1,1)中内闭一致收敛于ƒ(x),为了使n→∞时,Fn(ƒx)在[-1,1]上一致收敛于ƒ(x),充分必要条件是

这种在区间端点发生奇异的情况并不很稀有,它促使人们去改变端点的插值情况。P.图兰首先提出在区间端点xon=1,xn+1n=-1处取值与函数取值相同的要求。从而构造了拟埃尔米特-费耶尔插值多项式Q2n+1(ƒ,x),即假定结点组是取在开区间(-1,1)中的,而2n+1次代数多项式Q2n+1(ƒ,x)满足条件

这时,如取为Xn(x)的零点全体,则

当然也可以考虑仅在一端插值的情况。然而,倘若将端点作为结点,又会发生剧烈的变化。例如,取

则以 为结点组的埃尔米特-费耶尔插值多项式序列Fn+2(ƒ,x),对于 ƒ(x)=x2这样好的函数,也会在(-1,1)中处处发散。而取

为结点组时,相应的Fn+2(ƒ,x)对于连续函数ƒ(x)却有逼近阶

埃尔米特插值多项式可以从各方面扩充。例如,可以在某些结点处放弃对某些阶导数的要求,这就是所谓伯克霍夫插值。其中常见的是(0,2)插值,也即对于给定的结点组以及数组,,要确定一个次数不高于2n-1的代数多项式使得

,(k=1,2,…,n)。当取αkn=ƒ(γkn)时,考虑S2n-1(ƒx)对ƒ(x)的逼近,也可以考虑埃尔米特插值多项式对函数及其导数的同时逼近。例如,取

为结点,对于[-1,1]上的可微函数,考虑

ƒ(x)及ƒ'(x)的同时逼近。此时有

至于对于无限区间或周期函数的情形,自然也可作类似的讨论,只是在周期的情形,有时插值三角多项式却未必存在。

至于ƒ(x)的(α1,α2,…,αn)阶埃尔米特插值多项式Hn(x)对ƒ(x)的逼近,如果ƒ(x)在[αb]上有m阶导数,则在[αb]中有与x有关的点ξ使得

式中。

欢迎分享,转载请注明来源:内存溢出

原文地址: http://outofmemory.cn/bake/4700114.html

(0)
打赏 微信扫一扫 微信扫一扫 支付宝扫一扫 支付宝扫一扫
上一篇 2022-11-07
下一篇 2022-11-07

发表评论

登录后才能评论

评论列表(0条)

保存