Ch2 插值法
/* Interpolation */
当精确函数 y = f(x) 非常复杂或未知时,在
一系列节点 x
0
… x
n
处测得函数值 y
0
= f(x
0
),
… y
n
= f(x
n
) ,由此构造一个简单易算的近似
函数 g(x) f(x) ,满足条件 g(x
i
) = f(x
i
) (i =
0, … n) 。这里的 g(x) 称为 f(x) 的插值函数。
最常用的插值函数是 … ?
多项式
x
0
x
1
x
2
x
3
x
4
x
g(x) f(x)
§1 拉格朗日多项式 /* Lagrange Polynomial */
niyxP
iin
,...,0,)(
求 n 次多项式
使得
n
nn
xaxaaxP
10
)(
条件:无重合节点,即
ji
xx
ji
n = 1
已知 x
0
, x
1
;
y
0
,
y
1
,
求
xaaxP
101
)(
使得
111001
)(,)( yxPyxP
可见 P
1
(x) 是过 ( x
0
, y
0
) 和 ( x
1
, y
1
) 两点的直
线。
)()(
0
01
01
01
xx
xx
yy
yxP
10
1
xx
xx
01
0
xx
xx
= y
0
+ y
1
l
0
(x) l
1
(x)
1
0
)(
i
ii
yxl
称为拉氏基函数 /* Lagrange Basis */ ,
满足条件 l
i
(x
j
)=
ij
/* Kronecker Delta */
§1 Lagrange Polynomial
n 1
希望找到 l
i
(x) , i = 0, …, n 使得
l
i
(x
j
)=
ij
;然后
令
n
i
i
in
y
xlxP
0
)
()(
,则显然有 P
n
(x
i
) = y
i
。
l
i
(x)
每个 l
i
有 n 个根 x
0
…
x
i
… x
n
n
j
j i
jiniii
xxCxxxxxxCxl
0
0
)())...()...(()(
j i
j
i
iii
xx
Cxl
)(
1
1)(
n
j
ij
ji
j
i
xx
xx
xl
0
)(
)(
)(
n
i
iin
yxlxL
0
)()(
Lagrange
Polynomial
与 有关,而与 无关
节点
f
§1 Lagrange Polynomial
定理
( 唯一性 ) 满足
的 n 阶插值多项式是唯一存在的。
niyxP
ii
,...,0,)(
证明:
反证:若不唯一,则除了 L
n
(x) 外还有另一 n 阶多
项式 P
n
(x) 满足 P
n
(x
i
) = y
i
。
考察
则 Q
n
的阶数
,)()()( xLxPxQ
nnn
n
而 Q
n
有 个不同
的根
n + 1 x
0
… x
n
注:若不将多项式次数限制为 n ,则插值多项式不唯一。
例如
也是一个插值多项式,其中
可以是任意多项式。
n
i
in
xxxpxLxP
0
)()()()(
)( xp