两点三次 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+2xxkxk+1xk)(xxk+1xkxk+1)2yk+(1+2xxk+1xkxk+1)(xxkxk+1xk)2yk+1+(xxk)(xxk+1xkxk+1)2mk+(xxk+1)(xxkxk+1xk)2mk+1.

例题 1

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

f(ξ)=6f(ba)2(dc)+e(ba)(ba)3.

解答

作三次 Hermite 插值多项式:

H3(x)=c+e(xa)+(dc)e(ba)(ba)2(xa)2+f(ba)2(dc)+e(ba)(ba)3(xa)2(xb)

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(dc)e(ba)(ba)2(xa)+f(ba)2(dc)+e(ba)(ba)3[2(xa)(xb)+(xa)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

设函数 fC[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 ,则 gC[0,1]D3(0,1) 那么 {xk=0,xk+1=1yk=2,yk+1=0mk=0,mk+1=0 ,因此

H3(x)=2(1+x)(x1)2=4x36x2+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=0f(ξ)=24ξ+24.