设f(x)f(x)在节点a≤x0,x1,⋯,xn≤ba≤x0,x1,⋯,xn≤b处的函数值为f0,f1,...,fnf0,f1,...,fn,设P(x)为f(x)P(x)为f(x)在区间[a,b][a,b]上的具有一阶导数的插值函数 
(1)若要求P(x)P(x)在[a,b][a,b]上具有一阶导数(一阶光滑度)

 
P(xi)=f(xi)=fiP′(xi)=f′(xi)=f′,i=0,1,...,nP(xi)=f(xi)=fiP′(xi)=f′(xi)=f′,i=0,1,...,n

P(x)P(x)可以是最高次数为2n+1次多项式,两个节点就可以用2×1+1=32×1+1=3次多项式作为插值函数。 
(2)同样,若要求P(x)在[a,b]P(x)在[a,b]上具有m阶导数(m阶光滑度),即P(x)P(x)在节点x0,x1,...,xnx0,x1,...,xn处必须满足:

 
P(xi)=f(xi)=fiP′(xi)=f′(xi)=f′i⋯P(m)(xi)=f(m)(xi)=f(m)i,i=0,1,...,nP(xi)=f(xi)=fiP′(xi)=f′(xi)=fi′⋯P(m)(xi)=f(m)(xi)=fi(m),i=0,1,...,n

定义:称满足(1)或(2)式的插值问题为Hermite插值,称满足(1)或(2)式的插值多项式P(x)P(x)为Hermite插值多项式,记为Hk(x)Hk(x),k为多项式次数。

三次Hermite插值

hermite 相关算法整理-LMLPHP 
hermite 相关算法整理-LMLPHP 
hermite 相关算法整理-LMLPHP 
hermite 相关算法整理-LMLPHP 
hermite 相关算法整理-LMLPHP 
hermite 相关算法整理-LMLPHP 
hermite 相关算法整理-LMLPHP 
hermite 相关算法整理-LMLPHP

三次Hermite插值余项

 
R3(x)=f(4)(ξ)4!(x−x0)2(x−x1)2,x0≤ξ≤x1R3(x)=f(4)(ξ)4!(x−x0)2(x−x1)2,x0≤ξ≤x1

当f(4)(x)在[x0,x1]f(4)(x)在[x0,x1]上存在时,上述公式成立。

例题

hermite 相关算法整理-LMLPHP 
hermite 相关算法整理-LMLPHP

05-22 23:03