数值计算方法插值法讲解

又因为lk1(xk1) 1，故a(xk1 xk )(xk1 xk1) 1，得：
a

(xk 1

xk
1 )( xk 1

xk1) ，从而lk1(x)

(x (xk 1

xk xk
)(x xk1) ， )(xk1 xk1)
lk
(x)

(x ( xk

问题的提出
插值问题的数学提法：已知函数y f (x)在n 1个 点x0 , x1, , xn上的函数值yi f (xi ), (i 0,1, , n), 求一 个多项式y P(x)，使其满足P(xi ) yi , (i 0,1, , n). 即要求该多项式的函数曲线要经过y f (x)上已知的
平面上两点 xk , yk , xk1, yk1 ，求一条直线过该已
知两点。
线性插值
插值函数和插值基函数
由直线的点斜式公式可知：
P1(x)

yk

yk 1 xk 1

yk xk
(x
xk )，把此式按照
yk和yk1写成两项:P 1(x)
x xk1 xk xk 1
由于li (xk ) 0,k i，故li (x)有因子：
(x x0 ) (x xi1)(x xi1) (x xn )，因其已经是n 次多项式，故而仅相差一个常数因子。令：
li (x) a(x x0 ) (x xi1)(x xi1) (x xn )
插值法的概念
已知函数在n+1个点x0 ,x1 ,…,xn 上的函数值 yi=f(xi ), (i=0,1,…,n) ，求一个简单函数y=P(x),使其满 足: P(xi )=yi ,(i=0,1,…,n) 。即要求该简单函数的 曲线要经过y=f(x)上已知的这个n+1个 点： (x0 ,y0 ),(x1 ,y1 ),…,(xn ,yn ),同时在其它 x∈[a,b]上要估计误差: R(x) = f(x) - P(x)
由li (xi ) 1，可以定出a，进而得到：
li
(x)

(x x0 ) (xi x0 )
(x xi1)(x xi1) (xi xi1)(xi xi1)
(x xn ) (xi xn )
插值基函数
过n 1个不同的点分别决定n 1个n次插值函数 l0 (x),l1(x), ,ln (x)。每个插值基本多项式li (x)满足： （1）li (x)是n次多项式；

(x (20
10)(x 15) 10)(20 15)

1 50

x
10
故P2 (x)

y0l0
(
x
)
xyii
y l lg xi 11
(
x
)

y2l2 (x)

1 50
x 20
x 15
1.1761 x 10 x 20 1.3010 x 10 x 15
求一个次数不超过二次的多项式P2 (x)，使其满足
P2 (xk1)

yk1，P2 (xk )

yk，P2 (xk1)

yk
。
1
其几何意义为：已知平面上的三个点：
xk1, yk1 ， xk , yk , xk1, yk1 ，求一个二次
抛物线，使得该抛物线经过这三点。
二次插值基本多项式
从而Pn (x)是一个次数不超过n的多项式，且满足
Pn (xi ) yi，i 0,1, , n
以下就是拉格朗日基本多项式：
拉格朗日型n次插值多项式
已知函数y f (x)在n+1个不同的点x0, x1, , xn上 的函数值分别为y0, y1, , yn，求一个次数不超过n次 的多项式Pn (x)，使其满足
为什么要插值
1.在工程技术和科学研究中，有时对一个函数 f(x)只能通过实验或观测的手段得到它在某个 区间[a,b]上的有限个不同点上的函数值，也就 是只知道一张函数表，却没有明确的表达式。
2.虽然函数有明确的表达式，但由于形式复杂， 不便于计算和使用，所以人们往往希望做出一 个既能反映函数的特性，又便于计算的简单函 数P(x)去近似替代f(x)。
将(xi , f (xi )) 依次代入Pn (x ) ao a1x a2 x2 an xn 中
得到线性方程组
方程组的系数行列式是范德蒙（Vandermonde） 行列式：
当 xi 互异时，
，所以
线性方程组的解存在且惟一。即这样的插值
多项式Pn (x) ao a1x a2 x2 an xn 存 在而且惟一。（这同时也说明了不同形式的
y
。
f(x)
P(x)
y0
y1
y2
x0
x1
x2
当 n1 时,求一次多项式
yn-1 xn-1
yn
xn
x
插值法的分类
一，拉格朗日插值法 二，牛顿插值法 三，埃尔米特插值法 四，分段多项式插值法 五，样条插值法
一，拉格朗日插值法
流程：线性插值（一次插值）→二次插值→ n次拉格朗日插值法的方程组法证明→用中国
P2 (xi ) yi , i k 1, k, k 1
二次插值基本多项式
因为lk1(xk ) 0, lk1(xk1) 0，故lk1(x)有因子
(
x

xk
)(
x

lxkk(1x1()x)
xk 1
)，而其已经是一个二次多项式，仅
相差一个常数倍，可设lk1(x) a(x xk )(x xk1)，
有三个插值结点xk
1
,
xk
,
xk
，构造三个插值
1
基本多项式，要求满足： lxkk(1x1()x)
（1）基本多项式为二次多项式；
（2）它们的函数值满足下表：
xk 1
lk 1 ( x)
1
lk (x)
0
lk 1 ( x)
0
xk
xk 1
0
0
1
0
0
1
拉格朗日型二次插值多项式
由前述，拉格朗日型二次插值多项式： P2 (x) yk1lk1(x) yklk (x) yk1lk1(x)，P2 (x)是 三个二次插值多项式的线性组合，因为它是次数 不超过二次的多项式，且满足：
1
从而，P1(x) yklk (x) yk1lk1(x), 此形式称之为 拉格朗日型插值多项式。其中，插值基函数与
yk、yk
1无关，而由插值结点xk、xk
决定
1
例子
例1：已知lg10 1 , lg 20 1.3010，利用插值一次 多项式求 lg12的近似值。 解：f (x) lg x，f (x) lg x，f (10) 1，f (20) 1.3010 设x0 10，x1 20，y0 1，y1 1.3010， 则插值基本多项式为：
xk xk
1 1
)( )(
x xk1) xk xk1
)
，lk
1
(
x)

(x ( xk 1

xk1)(x xk ) xk1)(xk1 xk )
例子
例2：已知
xi xyii lg xi
10
yi lg xi 1
15
20
1.1761 1.3010
利用此三值的二次插值多项式求 lg12的近似值
n 1个点 x0, y0 , x1, y1 , , xn , yn ,同时在其他x a,b
上要估计误差R(x) f (x) P(x)
n1
个点
x0 , y0 , x1, y1 , , xn, yn ,
同时在其它
插值问题 上要估计误差
R(x) f (x) P(x)
其中P(x)为f(x)的插值函数，x0 ,x1 ,…,xn 称为插值节 点，包含插值节点的区间[a,b] 称为插值区间，求插值 函数P(x)的方法称为插值法。若P(x)是次数不超过n的 代数多项式，就称P(x)为插值多项式，相应的插值法 称为多项式插值。若P(x)是分段的多项式，就是分段 插值。若P(x)是三角多项式，就称三角插值。
插值多项式本质上是同一个）

当
{xi
}n i0
互异时，
，即行列试不为0，
由克莱默法则得知，线性方程组的解存在且惟一。即f(x)关
于
{的xi }nin次0插值多项式
存Pn (x) ao a1x a2 x2 an xn存在而且惟一。
注：（这同时也说明了不同形式的插值多项式本质n上是同
y
f(x)
y0 x0
y1 x1
y2 x2
P1(x)
x
拉格朗日插值公式
线性插值（一次插值）
已知函数f (x)在区间 xk , xk1的端点上的函数值
yk f (xk ), yk1 f (xk1)，求一个一次函数y P1(x) 使得yk P1(xk ), yk1 P1(xk1)。其几何意义是已知
（2）li (xi ) 1,而在其它n个li (xk ) 0,k i。
n次拉格朗日型插值多项式Pn(x)
Pn (x)是n 1个n次插值基本多项式l0 (x),l1(x), ,ln (x) 的线性组合，相应的组合系数是y0 , y1, , yn。即
n
Pn (x) y0l0 (x) y1l1(x) ynln (x) yklk (x) k 0
解：设x0 10，x1 15，x2 20，则