拉格朗日插值 python scipy

拉格朗日插值数学原理：

根据数学知识，对于平面上已知的n个点（无两点在一条直线上）可以找到一个 n-1 次多项式：

$$y=a_0 +a_1x+ a_2x^2+ ... + a_{n-1}x^{n-1}$$

为了进行根据新的x， 求出对应的 y值，需要求出上式中的系数 $a_0, a_1, a_2.... a_{n-1}$

因为n个点$(x_1,y_1), (x_2,y_2) ...(x_n,y_n)$在以上多项式上，代入每个点，得

$$y_1 = a_0 + a_1x_1 + a_2x_1^2 + ... + a_{n-1}x_1^{n-1}$$
$$y_2 = a_0 + a_1x_2 + a_2x_2^2 + ... + a_{n-1}x_2^{n-1}$$
$$......$$
$$y_n = a_0 + a_1x_n + a_2x_n^2 + ... + a_{n-1}x_n^{n-1}$$

以上方程组中，n个方程，n个未知数$a_0, a_1 ...a_{n-1}$,所以方程可解，可以利用线性代数中的行列式求解。
解出的拉格朗日插值多项式为：

$$L(x) = y_1\frac{(x-x_2)(x-x_3)...(x-x_n)}{(x_1-x_2)(x_1-x_3)...(x1-x_n)} +y_2\frac{(x-x_1)(x-x_3)...(x-x_n)}{(x_2-x_1)(x_2-x_3)...(x2-x_n)} + ......$$
$$+ y_n\frac{(x-x_1)(x-x_2)...(x-x_{n-1})}{(x_n-x_2)(x_n-x_3)...(x_n-x_{n-1})}$$
$$=\sum_{i=0}^n y_i \prod_{j=0,j\neq i} \frac{x-x_i}{x_i -x_j}$$

实例：
假设一个2次多项行式$f$, 按照前面的介绍，取三个点为
$f(3)=10, f(6) = 8, f(9)=4$

$l_0(3)= \frac{(x-6)(x-9)}{(3-6)(3-9)}$
$l_1(6)= \frac{(x-3)(x-9)}{(6-3)(6-9)}$
$l_2(9)= \frac{(x-3)(x-6)}{(9-3)(9-6)}$

$L(x)=f(3)l_0(3)+f(6)l_1(6)+f(9)l_2(9)$
$=10 \frac{(x-6)(x-9)}{(3-6)(3-9)}+8 \frac{(x-3)(x-9)}{(6-3)(6-9)}+4 \frac{(x-3)(x-6)}{(9-3)(9-6)}$
$=\frac{-x^2+3x+90}{9}$

插值$f(10)=\frac{20}{9}$

python 实现

$$y=a_0 +a_1x+ a_2x^2$$

from scipy.interpolate  import lagrange
x = [3, 6, 9]
y = [10, 8, 4]
lagrange(x,y)
#poly1d([ -0.11111111,   0.33333333,  10.        ])

以上 $lagrange(x,y)$ 的输出值 $poly1d([ -0.11111111, 0.33333333, 10. ])$ 值的是多项式的三个系数
即： $a_0= -0.11111 , a_1 = 0.3333333, a_2=10.$

如果要进行插值操作，可以：

lagrange(x, y)(10)
# 2.222222