两点三次 Hermite 插值

两点三次 Hermite 插值

基函数

采用基函数方法,令

H3(x)=αk(x)yk+αk+1(x)yk+1+βk(x)mk+βk+1(x)mk+1

其中 αk(x),αk+1(x),βk(x),βk+1(x) 是关于点 xk 及 xk+1 的三次 Hermite 插值基函数,它们应分别满足以下条件:

αk(xk)=1αk(xk+1)=0αk′(xk)=αk′(xk+1)=0αk+1(xk)=0αk+1(xk+1)=1αk+1′(xk)=αk+1′(xk+1)=0βk(xk)=βk(xk+1)=0βk′(xk)=1βk′(xk+1)=0βk+1(xk)=βk+1(xk+1)=0βk+1′(xk)=0βk+1′(xk+1)=1

解得基函数,那么有

H3(x)=(1+2x−xkxk+1−xk)(x−xk+1xk−xk+1)2yk+(1+2x−xk+1xk−xk+1)(x−xkxk+1−xk)2yk+1+(x−xk)(x−xk+1xk−xk+1)2mk+(x−xk+1)(x−xkxk+1−xk)2mk+1.

例题 1

若函数 f(x) 在 [a,b] 上三阶可导,且 f(a)=c,f(b)=d,f′(a)=e,f′(b)=f ,则 ∃ξ∈(a,b) ,使得

f‴(ξ)=6f(b−a)−2(d−c)+e(b−a)(b−a)3.

解答

作三次 Hermite 插值多项式:

H3(x)=c+e(x−a)+(d−c)−e(b−a)(b−a)2(x−a)2+f(b−a)−2(d−c)+e(b−a)(b−a)3(x−a)2(x−b)

则 H3(x) 与 f(x) 在 a,b 两点有相同的函数值与导数值,设 F(x)=f(x)−H3(x) ,则 F(a)=F(b)=0 ,由 Rolle 定理, ∃ξ1∈(a,b) ,使得 F′(ξ1)=0 ,
由于

F′(x)=f′(x)−H3′(x)=f′(x)+e+2(d−c)−e(b−a)(b−a)2(x−a)+f(b−a)−2(d−c)+e(b−a)(b−a)3[2(x−a)(x−b)+(x−a)2]

则 F′(a)=F′(b)=0 ,在 [a,ξ1] 和 [ξ1,b] 上应用 Rolle 定理,则 ∃ξ2∈(a,ξ1),ξ3∈(ξ1,b) ,使得

F″(ξ2)=F″(ξ3)=0

最后在 [ξ2,ξ3] 上应用 Rolle 定理,即 ∃ξ∈(ξ2,ξ3) ,使得 F‴(ξ)=0 ,即得证.

例题 2

设函数 f∈C[0,1]∩D3(0,1), f(0)=2, f(1)=1, f′(0)=0, f′(1)=4 证明: 存在一点 ξ∈(0,1) 使得 f‴(ξ)=24+24ξ.

解答

令 g(x)=f(x)−x4 ,则 g∈C[0,1]∩D3(0,1) 那么 {xk=0,xk+1=1yk=2,yk+1=0mk=0,mk+1=0 ,因此

H3(x)=2(1+x)(x−1)2=4x3−6x2+2

构造 F(x)=g(x)−H3(x) 有 F(0)=F(1)=0 ,由 Rolle 定理 ∃ξ1∈(0,1) 使得 F′(ξ1)=0 ,
又 F′(0)=F′(ξ1)=F′(1)=0 ,由 Rolle 定理得 ∃ξ2∈(0,ξ1),ξ3∈(ξ1,1) 使得 F″(ξi)=0,i=2,3 ,
再根据 Rolle 定理得 ∃ξ∈(ξ2,ξ3)⊂(0,1) 使得 F‴(ξ)=0 即

f‴(ξ)−4!ξ−24=0⇒f‴(ξ)=24ξ+24.