9512.net
甜梦文库
当前位置:首页 >> >>

Hermite插值公式



§5.5 Hermite插值公式 插值公式
Newton插值和Lagrange插值虽然构造比较简单,但都存 在插值曲线在节点处有尖点,不光滑,插值多项式在节 点处不可导等缺点. 为了保证插值多项式 Pn (x)能更好地逼近 f (x),对 Pn (x)增 加一些约束条件,例如要求 Pn (x) 在某些结点处与 f (x) 相切,即具有相同的导数值.

一、Hermite插值问题
求一个次数不大于n+r+1的代数多项式 H (x) ,满足:

H ( xi ) = f ( xi ), i = 0,1,2,L, n ------(1) H ′( xi ) = f ′( xi ), i = 0,1,2,L, r (r ≤ n)
1

称以上的插值问题为Hermite插值问题.

注意:
式(1)包含n+r+2个条件,所以能够确定次数不大于 n+r+1的代数多项式H (x) .

二、Hermite插值公式推导 令

H ( x) = ∑ hk ( x) f ( xk ) + ∑ hk ( x) f ′( xk ) ------(2)
k =0 k =0

n

r

其中,hk ( x) (k = 0,1,L, n)和hk ( x) (k = 0,1,L, r ) 都是 n+r+1次待定多项式,并且它们满足以下条件:.
2

1 hk ( xi ) = 0

i=k i≠k

i, k = 0,1,L, n

------(3)

′ hk ( xi ) = 0, k = 0,1,L, n; i = 0,1,L, r

1 hk ( xi ) = 0 ′

i=k i≠k

i, k = 0,1,L, r

------(4)

hk ( xi ) = 0, k = 0,1,L, r; i = 0,1,L, n
显然满足条件(3),(4)的多项式(2)的次数不大于n+r+1次, 且满足插值条件(1).

1.求解 hk ( x) (k = 0,1,L n)
由条件(3)知 xi 零点 .

(i = 0,1,L, r ; i ≠ k ) 是 hk (x) 的二重
3

且由条件(3)知 xi 的零点 .

(i = r + 1, r + 2,L, n; i ≠ k ) 是 hk (x)

(1)当0 ≤ k ≤ r时, hk ( x)具有如下形式:

hk ( x) = ( Ax + B)( x x0 ) 2 L( x xk 1 ) 2 ( x xk +1 ) 2 L( x xr ) 2 ( x xr +1 )L ( x xn ) = ( Ax + B)∏ ( x xi ) 2 ∏ ( x xi )
i =0 i≠k i = r +1 r n

------(5)

其中,A,B是待定系数

由条件(3)知hk ( xk ) = 1, hk ( xk ) = 0

4



( Axk + B)∏ ( xk xi ) 2 ∏ ( xk xi ) = 1
i =0 i≠k i = r +1

r

n

A∏ ( xk xi ) 2 ∏ ( xk xi ) +
i =0 i≠k i = r +1

r

n

2( Axk + B )∑ ( xk x j )∏ ( xk xi )
j =0 i =0 i≠ j i≠k n r n

r

r

2

i = r +1

∏ (x

n

k

xi ) +

( Axk + B ) ∑ ∏ ( xk xi ) 2 ∏ ( xk xi ) = 0
j = r +1 i =0 i≠k i = r +1 i≠ j

由上述两式解得:
5

A=

n 1 1 2∑ +∑ j =0 xk x j j = r +1 xk x j

r

( xk xi ) 2 ∏ ( xk xi ) ∏
i =0 i≠k i =r +1

r

n

B=

1 Axk ( xk xi ) 2 ∏ ( xk xi ) ∏
i =0 i≠k i = r +1 r n

将A,B代入式(5),得
′ ′ hk ( x) = {1 ( x xk )[lkn ( xk ) + lkr ( xk )]}lkn ( x)lkr ( x) ------(6) k = 0,1,L, r
6

其中,

x xi lkn ( x) = ∏ i =0 xk xi
n i≠k r

x xi lkr ( x) = ∏ i =0 xk xi
i≠k n

1 ′ lkn ( xk ) = ∏ i =0 xk xi
i≠k
r

1 ′ lkr ( xk ) = ∏ i =0 xk xi
i≠k
7

(2)当r + 1 ≤ k ≤ n时, hk ( x)具有如下形式:

hk ( x) = C ∏ ( x xi )
i =0

r

2

i = r +1 i≠k

∏ (x x )
i

n

------(7)

由条件(3)知hk ( xk ) = 1

C=

1

∏ (x
i =0

r

k

xi )

2

i =r +1 i≠k

∏ (x

n

k

xi )

将C代入式(7),得

wr ( x) hk ( x) = lkn ( x), wr ( xk )

k = r + 1, r + 2,L, n -(8)
8

其中,

wr ( x) = ∏ ( x xi )

r

wr ( xk ) = ∏ ( xk xi )
i =0
n

i =0 r

x xi lkn ( x) = ∏ i =0 xk xi
i≠k

综合(1)(2)得到 hk ( x) ( k = 0,1,L n) 即式(6),(8)

2.求解 hk ( x) (k = 0,1,L n)
由条件(4)知 xi 零点 .

(i = 0,1,L, r ; i ≠ k ) 是 hk (x) 的二重
9

且由条件(4)知 xi 的零点 .

(i = k , r + 1, r + 2,L, n) 是
n r

hk (x)

当0 ≤ k ≤ r时, hk ( x)具有如下形式:

hk ( x) = D∏ ( x xi )∏ ( x xi )
i =0

-----(9)

′ 由条件(4)知hk ( xk ) = 1
D= 1

i =0 i≠k

∑∏ (x
j =0 i =0 i≠ j

n

n

k

xi )∏ ( xk xi ) + ∑ ∏ ( xk xi )∏ ( xk xi )
i =0 i≠k j =0 i =0 j ≠k i =0 i≠k i≠ j

r

r

n

r

将D代入式(9),得

hk ( x) = ( x xk )lkn ( x)lkr ( x), k = 0,1,L, r -----(10)
10

其中,

x xi lkn ( x) = ∏ i =0 xk xi
n i≠k r

x xi lkr ( x) = ∏ i =0 xk xi
i≠k

由式(2)(6)(8)(10)所表示的多项式称为Hermite插值多项式 其中由式(6)(8)(10)所表示的多项式称为Hermite插值基函数

定理1 Hermite插值问题式(1)的解H ( x)存在而且唯一
证明:
11

存在性已由上面推导,下证唯一性. 反证法,设插值问题式(1)有两个不同的解 H1 ( x), H 2 ( x) 令 G ( x ) = H1 ( x ) H 2 ( x )

并且其为次数不大于n + r + 1的多项式, 且满足
G ( xi ) = 0, i = 0,1,L, n G′( xi ) = 0, i = 0,1,L, r

于是G ( x)必含有因式( x xi ) 2 (i = 0,1,L, r ) 和( x xi )(i = r + 1, r + 2,L, n)

故G ( x)的次数至少为n + r + 2, 矛盾.证毕
12

定理2 ( Hermite插值余项定理) Hermite插值公式的余项为 f ( n+r +2 ) (ξ ) f ( x) H ( x) = wn ( x) wr ( x) (n + r + 2)! 其中,ξ是插值区间(a, b)内的某一点.
证明:

引进辅助函数 wn (t ) wr (t ) F (t ) = f (t ) H (t ) [ f ( x) H ( x)] wn ( x) wr ( x)

由条件(1)知

F ( x) = F ( x0 ) = F ( x1 ) = L = F ( xn ) = 0

F ′( x0 ) = F ′( x1 ) = L = F ′( xr ) = 0
13

即F (t ) = 0有n r + 1个单根x, xr +1 , xr +2 ,L, xn 和r + 1个二重根x0 , x1 ,L, xr .

由Rolle定理知,F ′(t )在(a, b)内至少有n + r + 2个零点 F ′′(t )至少有n + r + 1个零点,依此类推可知:

F ( n+r +2 ) (t )在(a, b)内至少有一个零点ξ,因此
f
( n+ r + 2 )

(n + r + 2)! (ξ ) [ f ( x) H ( x)] = 0 wn ( x) wr ( x)

f ( n+r +2 ) (ξ ) 即得 f ( x) H ( x) = wn ( x) wr ( x) (n + r + 2)!
14

若r = n, 则相应的Hermite插值多项式为
H ( x) = ∑ hk ( x) f ( xk ) + ∑ hk ( x) f ′( xk )
k =0 k =0 n n

-----(11)

其中
2 ′ ( xk )]lkn ( x), k = 0,1,L, n hk ( x) = [1 2( x xk )lkn 2 hk ( x) = ( x xk )lkn ( x), k = 0,1,L, n

余项公式为: f ( 2 n+2 ) (ξ ) 2 f ( x) H ( x) = wn ( x), ξ ∈ (a, b) (2n + 2)!

特别当r = n = 1时,插值条件为:
15

H ( xi ) = f ( xi ), H ′( xi ) = f ′( xi ), i = 0,1

由此得三次Hermite插值多项式:
H ( x) = h0 ( x) f ( x0 ) + h1 ( x) f ( x1 ) + h0 ( x) f ′( x0 ) + h1 ( x) f ′( x1 ) --(12)

x x0 x x1 2 h0 ( x) = (1 + 2 )( ) x1 x0 x0 x1 x x1 x x0 2 h1 ( x) = (1 + 2 )( ) x0 x1 x1 x0

16 多项式(12)常用作分段低次插值,称为分段三次Hermite插值.

x x1 2 h0 ( x) = ( x x0 )( ) x0 x1 x x0 2 h1 ( x) = ( x x1 )( ) x1 x0

例1.

已知f ( x )在节点1,处的函数值为f (1) = 2 , f ( 2 ) = 3 2 f ( x )在节点1,处的导数值为f ′(1) = 0 , f ′( 2 ) = 1 2

求f ( x )的两点三次插值多项式, 及f ( x )在x = 1.5 ,1.7处的函数值.

解:

x0 = 1, x1 = 2

y0 = 2 , y1 = 3

′ ′ y 0 = 0 , y1 = 1

′ ′ H 3 ( x) = y0 h0 ( x) + y1h1 ( x) + y0 h0 ( x) + y1h1 ( x)

x x0 x x1 x x1 x x0 = y0 1 + 2 x x + y1 1 + 2 x x x x x1 x0 0 1 0 0 1 1

2

2

x x0 x x1 ′ ′ + y0 ( x x0 ) + y1 ( x x1 ) x x x x 1 0 1 0
2

2

17

2 H 3 ( x ) = 2 (1 + 2( x 1)) ( x 2 )2 + 3(1 2( x 2 )) ( x 1)

( x 2 ) ( x 1)2 = 3 x 3 + 13 x 2 17 x + 9

f ( 1.5) ≈ H 3 ( 1.5) = 2.625 f (1.7 ) ≈ H 3 (1.7 ) = 2.931

作为多项式插值,三次已是较高的次数,次数再高就有 可能发生Runge现象,因此,对有n+1节点的插值问题, 我们可以使用分段两点三次Hermite插值
18



更多相关文章:
数值分析实验四Hermite插值法
数值分析实验报告 专业:计算机科学与技术 班级:14 汉(2) 学号:20141501069 姓名:于童 指导教师:马季骕老师 实验项目 Hermite 插值法 学会 Hermite 插值法,并应用...
分段三次Hermite插值
要想得到在插值区间上光滑的分段线性插值多项式,可采用分段埃尔 米特(Hermite)插值,这里我们考虑在整个 , 上用分段三次埃尔米特 插值多项式来逼近()。 一般的将...
分段Hermite插值
分段Hermite插值_IT/计算机_专业资料。6.6 分段埃尔米特插值及其 MATLAB 程序分段...1≤ x ≤1 H n ,3 ( x) 的误差公式和误差限. 解在 MATLAB 工作窗口...
分段三次Hermite插值
分段三次 Hermite 插值 1. 目的意义: 可以得到在插值区间上光滑的分段插值多项式 2. 数学模型(数学公式) : ? H1 ( x), x ? [ x0 , x1 ] ? H ( x...
埃尔米特(Hermite)插值
要求编程构造埃尔米特插值多项式 H2n+1(t),满足 H2n+1(ti)=si, H'2n+1(ti)=vi,对所有 i=0, 1, ..., n 成立,并据此计算 m 个给定时刻的里程和...
牛顿形式的埃尔米特插值多项式
牛顿形式的埃尔米特插值多项式 【摘要】 :在了解了插值法以后,陆续的又接触和学习到多项式插值、拉格朗日 插值、牛顿插值多项式等,但在有些实际问题中,仍需要其它...
三次埃尔米特插值
x 2 ? 1 输入 6 个节点的分段线性插值 4 输入 6 个节点的分段三次埃尔米特插值 请输入插值函数:x^2+1 请选择插值法(1 为分段线性插值,2 为分段三次...
常见插值法
插值 光滑 节点处导数相同 Hermite 插值法 导数不 易得到 plot(xi,yi_Lagarange,':g'); %画出Hermite插值 由边界条件推 导各节点导数 三次样条 插值法 ...
Hermite插值matlab编程
Hermite插值matlab编程_物理_自然科学_专业资料。数学建模 !!!1.Hermite 插值 function...Matlab的Hermite插值法 1页 1下载券 实习:Matlab作业hermit... 11页 2下...
插值多项式汇总
插值多项式汇总_数学_自然科学_专业资料。请求高手帮忙给个 hermite 插值法和拉格朗日插值法的 matlab 程序, 要求实现 x=[-5.0000 -3.8889 -2.7778 -1.6667 -0...
更多相关标签:

All rights reserved Powered by 甜梦文库 9512.net

copyright ©right 2010-2021。
甜梦文库内容来自网络,如有侵犯请联系客服。zhit325@126.com|网站地图