拉格朗日插值法

1. 拉格朗日插值法简介

在数值分析中,拉格朗日插值法是以法国十八世纪数学家约瑟夫·拉格朗日命名的一种多项式插值方法。许多实际问题中都用函数来表示某种内在联系或规律,而不少函数都只能通过实验和观测来了解。如对实践中的某个物理量进行观测,在若干个不同的地方得到相应的观测值,拉格朗日插值法可以找到一个多项式,其恰好在各个观测的点取到观测到的值。这样的多项式称为拉格朗日(插值)多项式。数学上来说,拉格朗日插值法可以给出一个恰好穿过二维平面上若干个已知点的多项式函数。拉格朗日插值法最早被英国数学家爱德华·华林于1779年发现,不久后(1783年)由莱昂哈德·欧拉再次发现。1795年,拉格朗日在其著作《师范学校数学基础教程》中发表了这个插值方法,从此他的名字就和这个方法联系在一起。
拉格朗日插值法_第1张图片

约瑟夫·拉格朗日伯爵(1736 -1813)

2.定义

对某个多项式函数,已知有给定的 k + 1 k + 1 k+1个取值点: ( x 0 , y 0 ) , … , ( x k , y k ) (x_{0},y_{0}),\ldots ,(x_{k},y_{k}) (x0,y0),,(xk,yk),其中 x j x_j xj对应着自变量的位置,而 y j y_{j} yj对应着函数在这个位置的取值。假设任意两个不同的 x j x_j xj都互不相同,那么应用拉格朗日插值公式所得到的拉格朗日插值多项式为:
L ( x ) : = ∑ j = 0 k y j ℓ j ( x ) L(x):=\sum _{ {j=0}}^{ {k}}y_{j}\ell _{j}(x) L(x):=j=0kyjj(x)
其中每个 ℓ j ( x ) \ell _{j}(x) j(x)为拉格朗日基本多项式(或称插值基函数),其表达式为:

ℓ j ( x ) : = ∏ i = 0 ,   i ≠ j k x − x i x j − x i = ( x − x 0 ) ( x j − x 0 ) ⋯ ( x − x j − 1 ) ( x j − x j − 1 ) ( x − x j + 1 ) ( x j − x j + 1 ) ⋯ ( x − x k ) ( x j − x k ) \ell _{j}(x):=\prod _{ {i=0,\,i\neq j}}^{ {k}}{\frac {x-x_{i}}{x_{j}-x_{i}}}={\frac {(x-x_{0})}{(x_{j}-x_{0})}}\cdots {\frac {(x-x_{ {j-1}})}{(x_{j}-x_{ {j-1}})}}{\frac {(x-x_{ {j+1}})}{(x_{j}-x_{ {j+1}})}}\cdots {\frac {(x-x_{ {k}})}{(x_{j}-x_{ {k}})}} j(x):=i=0,i=jkxjxixxi=(xjx0)(xx0)(xjxj1)(xxj1)(xjxj+1)(xxj+1)(xjxk)(xxk)

拉格朗日基本多项式 ℓ j ( x ) \ell _{j}(x) j(x)的特点是在 x j x_{j} xj上取值为1,在其它的点 x i ,   i ≠ j x_{i},\,i\neq j xi,i=j上取值为0。

即:对于给定的 n + 1 n+1 n+1个点 ( x 0 , y 0 ) , ( x 1 , y 1 ) , … , ( x n , y n ) (x_{0},y_{0}),(x_{1},y_{1}),\ldots ,(x_{n},y_{n}) (x0,y0),(x1,y1),,(xn,yn),对应于它们的次数不超过 n n n的拉格朗日多项式 L L L只有一个。

3. 证明

3.1 存在性证明

对于给定的 k + 1 k+1 k+1个点: ( x 0 , y 0 ) , … , ( x k , y k ) (x_{0},y_{0}),\ldots ,(x_{k},y_{k}) (x0,y0),,(xk,yk),拉格朗日插值法的思路是找到一个在一点 x j x_{j} xj取值为1,而在其他点取值都是0的多项式 ℓ j ( x ) \ell _{j}(x) j(x)。这样,多项式 y j ℓ j ( x ) y_{j}\ell _{j}(x) yjj(x)在点 x j x_{j} xj取值为 y j y_{j} yj,而在其他点取值都是0。而多项式 L ( x ) : = ∑ j = 0 k y j ℓ j ( x ) L(x):=\sum _{ {j=0}}^{ {k}}y_{j}\ell _{j}(x) L(x):=j=0ky

你可能感兴趣的:(#,密码学与安全技术,拉格朗日插值法)