Hermite.pdf

2.3 埃尔米特 (Hermite) 插值
朱升峰
华东师范大学数学科学学院
2021 年 3 月
. . . . . . . . . . . . . . . . . . . .

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
重节点均差
根据均差性质，易证
定理 1
设 f(x) ∈ Cn
[a, b]，x0, x1, · · · , xn 为 [a, b] 上互异节点，则 f [x0, x1, · · · , xn] 是各
变量的多元连续函数。
定义重节点的一阶均差
f[x0, x0] := lim
x→x0
f [x0, x] = lim
x→x0
f(x) − f (x0)
x − x0
= f′
(x0)
定义重节点的二阶均差
f [x0, x0, x1] =
f [x0, x1] − f [x0, x0]
x1 − x0
当 x1 → x0 时，由性质 f[x0, x1, · · · , xn] = f(n)
(ξ)
n! 得
f [x0, x0, x0] = lim
x1→x0
x2→x0
f [x0, x1, x2] =
1
2
f′′
(x0) .
朱升峰 (ECNU) Hermite 插值 2021.03 2 / 10

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
泰勒（Taylor）插值多项式
定义重节点的 n 阶均差
f [x0, x0, · · · , x0] = lim
xi→x0
f [x0, x1, · · · , xn] =
1
n!
f(n)
(x0) .
Newton 均差插值多项式中若令 xi → x0(i = 1, 2, · · · , n) 则得泰勒插值多项式
Pn(x) = f (x0) + f′
(x0) (x − x0) + · · · +
f(n)
(x0)
n!
(x − x0)
n
.
满足插值条件
P(k)
n (x0) = f(k)
(x0) , k = 0, 1, · · · , n
余项为
Rn(x) =
f(n+1)
(ξ)
(n + 1)!
(x − x0)
n+1
, ξ ∈ (a, b).
Lagrange 插值 ←→ Newton 插值 ←→ Taylor 插值

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
Hermite 插值
n + 1个含函数值和导数值的插值条件可构造次数 ⩽ n 的 Hermite 插值多项式
三点三次 Hermite 插值多项式
两点三次 Hermite 插值多项式
三点三次 Hermite 插值多项式: 构造插值多项式 P = P(x) 满足插值条件
P(xi) = f(xi) (i = 0, 1, 2) 且P′
(x1) = f′
(x1)
f 的二次 Newton 插值多项式
P2(x) = f (x0) + f [x0, x1] (x − x0) + f [x0, x1, x2] (x − x0) (x − x1)
满足 P2(xi) = f(xi)(i = 0, 1, 2) 。因而
(P − P2)(xi) = 0 ⇒ P(x) − P2(x) = A(x − x0)(x − x1)(x − x2)

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
插值余项
由 P′
(x1) = f′
(x1) 得
A =
f′
(x1) − f [x0, x1] − (x1 − x0) f [x0, x1, x2]
(x1 − x0) (x1 − x2)
设
R(x) := f(x) − P(x) = k(x)(x − x0)(x − x1)2
(x − x2),
其中 k(x) 为待定函数。类似 Lagrange 插值余项的推导得
定理 2
若 f(x) ∈ C4
[a, b]，则对任何 x ∈ [a, b]，插值余项
R(x) =
f(4)
(ξ)
4!
(x − x0)(x − x1)2
(x − x2),
ξ ∈ (a, b) 且依赖于 x。

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
例：给定 f(x) = x3/2
, x0 = 1
4 , x1 = 1, x2 = 9
4 , 试求三次 Hermite 插值多项式，并
求出余项表达式。
解
f0 = f
(
1
4
)
=
1
8
, f1 = f(1) = 1, f2 = f
(
9
4
)
=
27
8
,
f[x0, x1] =
7
6
, f[x0, x1, x2] =
11
30
令
P(x) =
1
8
+
7
6
(
x −
1
4
)
+
11
30
(
x −
1
4
)
(x − 1)
+ A
(
x −
1
4
)
(x − 1)
(
x −
9
4
)
由条件P′
(1) = f′
(1) =
3
2
可得
A = −
14
225
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
所求 Hermite 插值多项式
P(x) =
1
8
+
7
6
(
x −
1
4
)
+
11
30
(
x −
1
4
)
(x − 1) −
14
225
(
x −
1
4
)
(x − 1)
(
x −
9
4
)
余项
R(x) = f(x) − P(x) =
f(4)
(ξ)
4!
(
x −
1
4
)
(x − 1)2
(
x −
9
4
)
.
=
1
4!
9
16
ξ−5/2
(
x −
1
4
)
(x − 1)2
(
x −
9
4
)
, ξ ∈
(
1
4
,
9
4
)
.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
两点三次 Hermite 插值多项式
已知插值节点 xk，xk+1 及节点上函数值和导数值求插值多项式 H3(x)，满足
H3 (xk) = yk, H3 (xk+1) = yk+1,
H′
3 (xk) = mk, H′
3 (xk+1) = mk+1.
插值基函数方法: 基函数 αk(x), αk+1(x), βk(x), βk+1(x)
H3(x) = αk(x)yk + αk+1(x)yk+1 + βk(x)mk + βk+1(x)mk+1 ,
满足插值条件
α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) = a′
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

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
由 αk (xk+1) = 0, α′
k (xk+1) = 0, 可设
αk(x) = (ax + b)(x − xk+1)2
又由
1 = αk (xk) = (axk + b)(xk − xk+1)2
0 = α′
k (xk) = 2(axk + b)(xk − xk+1) + a(xk − xk+1)2
解得
a = −
2
(xk − xk+1)3
, b =
(
1 +
2xk
xk − xk+1
)
1
(xk − xk+1)2
,
αk(x) =
(
1 + 2
x − xk
xk+1 − xk
) (
x − xk+1
xk − xk+1
)2
同理
αk+1(x) =
(
1 + 2
x − xk+1
xk − xk+1
) (
x − xk
xk+1 − xk
)2

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
求 βk(x): 可令
βk(x) = a (x − xk) (x − xk+1)2
由 β′
k(xk) = 1 得 a = 1/(xk − xk+1)2
，从而
βk(x) = (x − xk)
(
x − xk+1
xk − xk+1
)2
同理，
βk+1(x) = (x − xk+1)
(
x − xk
xk+1 − xk
)2
H3(x) =
(
1 + 2
x − xk
xk+1 − xk
) (
x − xk+1
xk − xk+1
)2
yk +
(
1 + 2
x − xk+1
xk − xk+1
) (
x − xk
xk+1 − xk
)2
yk+1
+ (x − xk)
(
x − xk+1
xk − xk+1
)2
mk + (x − xk+1)
(
x − xk
xk+1 − xk
)2
mk+1
余项 R3(x) = f(x) − H3(x)
R3(x) =
1
4!
f(4)
(ξ) (x − xk)2
(x − xk+1)2
, ξ ∈ (xk, xk+1)

Hermite.pdf

More Related Content

Similar to Hermite.pdf (20)

Hermite.pdf