当前位置:文档之家› 插值计算与插值多项式

插值计算与插值多项式


公式,求 7
p2(x) = +
(x–x1)(x–x2) (x0–x1)(x0–x2) (x–x0)(x–x1) (x2–x0)(x2–x1)
y0 + y2
(x–x0)(x–x2) (x1–x0)(x1–x2)
y1
x0=1, x1=4, x2=9
y0=1, y1=2, y2=3
(7–4)(7–9)
(7–1)(7–9)
i=0, 1, 2
解: 用待定系数法, 将各节点值依次代入所求多项式, 得
解上述方程, 将求出的a0, a1, a2 代入 p(x) = a0 + a1x + a2x2 即得所求二次多项式
6.2.2 拉格朗日插值多项式
❖ 我们看到,两个插值点可求出一次插值多项式
p1(x),而三个插值点可求出二次插值多项式p2(x) 。 当 插 值 点 增 加 到 n+1 个 时 , 我 们 可 以 利 用
6.1 插值法的数学描述
设函数y=f(x) 在区间[a, b]上连续, x0 , x1 , , xn 是 [a, b]上取定的n+1个互异节点,且在这些点处的函数值 为已知 f (x0 ), f (x1), , f (xn ) ,即 yi f (xi ) 若存在一个 f(x)的近似函数 (x),满足
x1
(给定的三个点在一条直线上)
例6.6 已知f (x)的观测数据
x 0124
f (x) 1 9 23 3
构造Lagrange插值多项式
解 四个点可构造三次Lagrange插值多项式:基函数为
l0
(x)

(x (0
1)( x 1)(0

2)( x 2)(0

4) 4)


1 8
1 3
3 1
2
f (1.5) p(1.5) 1.25
(2) 抛物插值
抛物插值又称二次插值,它也是常用的代数插值之
一。设已知f(x)在三个互异点x0,x1,x2的函数值y0 ,y1,y2,要构造次数不超过二次的多项式
P(x) a2 x2 a1x a0
使满足二次插值条件:
P(xi ) yi (i 0,1,2)
为了构造满足插值条件 p(xi ) f (xi ) (i=0,1,2,…,n ) 的便于使用的插值多项式P(x),先考察几种简单情形, 然后再推广到一般形式。 6.2.1 线性插值与抛物插值 (1)线性插值
线性插值是代数插值的最简单形式。假设给定了函数 f(x)在两个互异的点 x0 ,x1 的值,y0 f (x0 ), y1 f (x1)
y1
的线性组合得到,其系数分别为 y0,y1
称 l0 (x),l1(x)为节点 x0 , x1的线性插值基函数
线性插值基函数 l0 (x),l1(x) 满足下述条件
xi
x0
x1
l0 (x)
1
0
l1 ( x)
0
1
并且他们都是一次函数。
注意他们的特点对下面的推广很重要 于是线性插值函数可以表示为与基函数的线性组合
x3

7 8
x2

7 4
x
1
l1 ( x)

(x 0)(x 2)(x 4) (1 0)(1 2)(1 4)

1 3
x3

2x2

8 3
x
l2
(x)

(x (2

0)( x 0)(2
1)( x 1)( 2
即:
lk
(xi
)

ki

1 0
(i k) (i k)
由条件 lk (xi ) 0 ( i k )知, x0 , x1, , xk1, xk1, , xn 都是n次 lk (x) 的零点,故可设
lk (x) Ak (x x0 )( x x1 ) (x xk1 )( x xk1 ) (x xn )
,现要求用线性函数 p(x) ax b 近似地代替f(x)。选
择参数a和b, 使 p(xi ) f (xi )(i 0,1) 。称这样的线性函数 P(x)为f(x)的线性插值函数 。
线性插值
线性插值多项式
由直线两点式可知,通过A,B的直线方程为
p(x)

y0

y1 x1

y0 x0
l0 (x) c(x x1)( x x2 )
1
再由另一条件 l0 (x0 ) 1
确定系数 c (x0 x1)(x0 x2 )
从而导出
l0
(x)

(x (x0

x1)( x x2 ) x1)( x0 x2 )
P(x)的参数 a0 , a1, a2 直接由插值条件决定,
jk
j0
jk
称 lk (x) 为关于基点 xi 的n次插值基函数(i=0,1,…,n)
以n+1个n次基本插值多项式lk (x)(k 0,1, , n) 为基础,就能直接写出满足插值条件
P(xi ) f (xi ) (i 0,1,2, , n) 的n次代数插值多项式。
P(x) l0 (x) y0 l1(x) y1 ln (x) yn 事实上,由于每个插值基函数 lk (x)(k 0,1, , n) 都是n次值多项式,所以他们的线性组合
方程组的解唯一。
类似地可以构造出满足条件:l1(x1) 1, l1(x0 ) 0, l1(x2 ) 0
的插值多项式
l1 (x)

(x ( x1

x0 )( x x2 ) x0 )( x1 x2 )
及满足条件:l2 (x2 ) 1, l2 (x0 ) 0, l2 (x1) 0 的插值多项式
(x
x0 )

p1 ( x)
它也可变形为
l0 (x)

x x1 x0 x1
,
l1 (x)

x x0 x1 x0
显然有:

l0 (x)

x x1 x0 x1
l1 ( x)

x x0 x1 x0
可以看出
L1 ( x)

x x1 x0 x1
y0

x x0 x1 x0
p(xi ) f (xi ) (i 1,2, , n)
则称 p(x) 为f(x)的一个插值函数, f(x)为被插函数, 点 xi为插值节点, R(x)= f (x) p(x) 称为插值余项, 区间 [a, b]称为插值区间, 插值点在插值区间内的称为内插, 否则称外插
插值的几何意义
6.2 拉格朗日(Lagrange)插值
)
y0

(x (x1

x0 x0
)( )(
x x2) x1 x2 )
y1

(x (x2

x0 x0
)( )(
x x1) x2 x1
)
y2
p(x) (x 1)(x 2) 1 (x 0)(x 2) 2 (x 0)(x 1) 3 (0 1)(0 2) (1 0)(1 2) (2 0)(2 1)
这就是二次插值问题。其几何意义是用经过3个点
(x0, y0 ), (x1, y1), (x2, y2 ) 的抛物线 y P(x) 近似代替曲线
y f (x) , 如下图所示。因此也称之为抛物插值。
抛物插值函数
y
y0
y1
y=L2(x) y1
y=f(x)
O
x0
x1
x2
x
因过三点的二次曲线为抛物线,故称为抛物插值。
x x0 x1 y y0 y1

xn1 xn yn1 yn
• 插值问题:根据这些已知数据来构造函数 y f (x) 的一种
简单的近似表达式,以便于计算点 x xi ,i 0,1,L , n 的函 数值 f (x),或计算函数的一阶、二阶导数值。
y=p(x) y=f(x)
简单的说,插值的目的就是根据给定的数据表,寻 找一个解析形式的函数p(x),近似代替f(x)
y 115 p(115) 10.714
例6.2 已知y=f(x)的函数表
X1 3
y12
求线性插值多项式, 并计算x=1.5 的值
解: 由线性插值多项式公式得
p(x)
x x1 x0 x1
y0

x x0 x1 x0
y1
x 3 1 x 1 2 1 (x 1)
p2(7) =
(1–4)(1–9)
*1 + (4–1)(4–9)
*2
(7–1)(7–4)
+ (9–1)(9–4) * 3
= 2.7
例6.4 已知函数y=f(x)在节点上满足
x x0 x1 x2
y y0 y1 y2 求二次多项式 p(x) = a0 + a1x + a2x2
使之满足 p(xi) = yi

x1 x1
)( )(
x x2) x0 x2
)
y0

(x (x1

x0 x0
)( )(
x x2) x1 x2 )
y1

(x (x2

x0 x0
)( )(
x x1) x2 x1
)
y2
容易看出,P(x)满足条件 P(xi ) yi (i 0,1,2)
例6.3 已知x=1, 4, 9 的平方根值, 用抛物插值
相关主题