11. SciPy拉格朗日插值
拉格朗日插值法是以法国十八世纪数学家约瑟夫·路易斯·拉格朗日命名的一种多项式插值方法。许多实际问题中都用函数来表示某种内在联系或规律,而不少函数都只能通过实验和观测来了解。如对实践中的某个物理量进行观测,在若干个不同的地方得到相应的观测值,拉格朗日插值法可以找到一个多项式,其恰好在各个观测的点取到观测到的值。这样的多项式称为拉格朗日(插值)多项式。
11.1 拉格朗日多项式
一般下,若已知函数$\;y = f(x)$在[a, b]上有定以,且有$\;n+1$个互异的值$\;x_0,x_1,\cdots,x_n$处的函数值为$\;y_0,y_1,\cdots,y_n$,则可以考虑构造一个过这$\;n+1 $个点的、次数不超过$n$的多项式$\;y = P_n(x)$,使其满足:
$$
P_n(x_k) = y_k, k \in[0, n]
$$
函数$\;y = P_n(x)$称为$\;f(x)$的插值函数,$x_0,x_1,\cdots,x_n$为插值点。
设集合$D_n= ${$0,1,\cdots,n$},对于任意$\;k \in D_n$都有$\;P_k(x)$ ,$\;B_k = ${$\;i\;\vert\; i \neq k\;,\;i \in D_n$}得到拉格朗日基本多项式(插值基函数)公式:
$$
P_k(x) = \prod_{i \in B_k} \frac{x-x_i}{x_k - x_i}
$$
应用拉格朗日插值基函数得到拉格朗日插值多项式:
$$
L_n(x)=\sum_{j = 0}^{n}y_jP_j(x)
$$