真.浅析拉格朗日插值法

真.浅析拉格朗日插值法

拉格朗日插值法用来干什么

mmmm…，这个不好解释
形象的来说,就是给你一个函数$f(x)$(不知道函数是什么的自己查)，这个函数你不知道是长什么样子的，然后告诉你这个函数的最高次数是$n$,再给你$n+1$数，以及这n+1个点每个的函数值也告诉你，让你求$f(k)$,k可以为任何数

从数学角度来看待:
$$\begin{gathered} 有一个函数f(x),最高次数项为n \\ 已知n+1个数x_0,x_1,x_2...x_n,和f(x_0),f(x_1)....f(x_n) \\ 求f(k) \end{gathered}$$
说实话，我敢称拉格朗日插值为数学读心术

拉格朗日插值怎么做

首先，我们先抛开拉格朗日插值法不谈，就只看题目，你的想法是什么?

一阶同学:?,二阶同学:?,三阶同学:?,长郡oj学生(周禹坤/姜元程/…):高斯消元!!!

没错，确实可以用高斯消元来解,问题是怎么解?
$$\begin{gathered} 我们知道这个函数的最高次数项为n,所以说这个函数可以写成如下形式: \\ f(x)=a_0+a_{1}x+a_2{x}^2+a_3{x}^3+...+a_n{x}^n \\ {a}_i是i次项的系数 \\ 因此可以写出高斯消元组 \\ {a}_0+a_1{x}_0+a_2{x}_0^2+a_3{x}_0^3+...+a_n{x}_0^n=f(x_0) \\ {a}_0+a_1{x}_1+a_2{x}_1^2+a_3{x}_1^3+...+a_n{x}_1^n=f(x_1) \\ . \\ . \\ . \\ {a}_0+a_1{x}_n+a_2{x}_n^2+a_3{x}_n^3+...+a_n{x}_n^n=f(x_n) \\ 把每个系数求出来之后，再将k带入进去，就可以得到f(k)了 \end{gathered}$$

我们想想这种方法的时间复杂度是多少?

$O(n^3+n)=O(n^3)$

可以看到复杂度特别的高，对于$n\le 10^3$的数据就过不了了

那么我们再来看拉格朗日插值

定理一:$n+1$个点，确定一个$n$次多项式

所以说我们只要自己构造一个n次函数$g(x)$，他可以穿过那$n+1$个点，那么他就是$f(x)$

那我们考虑如何构造。我们可以构造$n+1$个函数$f_0(x),f_1(x),f_2(x)...f_n(x)$

构造要求:$f_i(x),当x为x_i时,f_i(x)为1,当x为x_j(j不等于i)时f_i(x)为0,其他情况后面再说$

$先看当x为x_j(j不等于i)时f_i(x)为0这一条，我们知道,0乘任何数都等于0，可以考虑f_i(x)={\prod }_{j=0,j不等于i}^n(x-x_j),这样子就满足了这一条$

$那么f_i(x),当x为x_i时,f_i(x)为1,这条怎么和上条结合呢，哦，我们只需要在除以一个{\prod }_{j=0,j不等于i}^n(x_i-x_j)就可以了$

$所以总结下来，公式就是f_i(x)=\frac{{\prod }_{j=0,j不等于i}^n(x-x_j)}{{\prod }_{j=0,j不等于i}^n(x_i-x_j)}$

那么构造这些函数有什么用呢，我们看一下，如果将函数$f_i(x)$去乘上一个$f(x_i)$,你就会发现这个函数$f(x_i)f_i(x)$在$x=x_i$时，函数值为$f(x_i)$,$当x为x_j(j不等于i)时f(x_i)f_i(x)为0$

将每个函数加起来，得到的最终结果就是原函数了

$最终公式为f(x)={\sum }_{i=0}^{n}f(x_i)\frac{{\prod }_{j=0,j不等于i}^n(x-x_j)}{{\prod }_{j=0,j不等于i}^n(x_i-x_j)}$

注意：在题目中通常会说，过n+1个点$(x_i,y_i)$,这里的$y_i就是f(x_i)$