三点插值公式的推导过程

参考原文链接https://blog.csdn.net/m0_54563743/article/details/116096015

第四章 插值法与曲线拟合

4.1、拉格朗日插值法

拉格朗日插值法可以找到一个多项式，其恰好在各个观测的点取到观测到的值.

就是给定些“点”，可以找到一个式子，把这些点穿起来

• 线性插值（两点插值）（两点确定一条直线）

• (x₀ , y₀)、(x₁ , y₁)

  • 用中学学过的直线方程的求解方法把这条直线表示出来。

  • 再将直线方程写为插值基函数×函数值的形式。
    $$L_1(x)=y_0\frac{x-x_1}{x_0-x_1}+y_1\frac{x-x_0}{x_1-x_0}$$

• 三点插值公式的推导（抛物线插值）：待定系数法

  • (x₁ , y₁),(x₂ , y₂),(x₃ , y₃)。这三点不在同一条直线上，则它们可以确定一条抛物线

  • 设抛物线的方程为：
    $$y=a{x}^2+bx+c$$

  • 将这三个点代入该方程中有：
    $$\{\begin{array}{c}a{x}_1^2+b{x}_1+c=y_1\\ a{x}_2^2+b{x}_2+c=y_2\\ a{x}_3^2+b{x}_3+c=y_3\end{array}$$

  • 其中x₁、x₂、x₃、y₁、y₂、y₃是已知量，将a,b,c当作未知量。写成矩阵乘法的形式(3×3 × 3×1==>3×1)，如下所示：
    $$\left[\begin{array}{ccc}{x}_1^2& x_1& 1\\ x_2^2& x_2& 1\\ x_3^2& x_3& 1\end{array}\right]\left[\begin{array}{c}a\\ b\\ c\end{array}\right]=\left[\begin{array}{c}{y}_1\\ y_2\\ y_3\end{array}\right]$$

  • 等式两边同时左乘A的逆后再同时左乘[x² , x , 1]，即可得到该抛物线的方程。

    • 在这里解释一下为什么左乘[x² , x , 1]，因为这样一乘就可以直接得出该抛物线方程的表达式：1×3 × 3×1 ==>式子

    因此该抛物线方程可以写成如下形式，(求逆的时候我用的方法是有伴随除以行列式)
    $$A^{-1}=\frac{A^{\ast }}{|A|}$$

    $$y=\left[\begin{array}{ccc}{x}^2& x& 1\end{array}\right]\ast \left[\begin{array}{ccc}{x}_2-x_3& -(x_1-x_3)& x_1-x_2\\ -(x_2^2-x_3^3)& (x_1^2-x_3^2)& -(x_1^2-x_2^2)\\ x_2{x}_3(x_2-x_3)& -x_1{x}_3(x_1-x_3)& x_1{x}_2(x_1-x_2)\end{array}\right]\ast \frac{1}{|A|}\ast \left[\begin{array}{c}{y}_1\\ y_2\\ y_3\end{array}\right]$$

  • 接下来就是计算这个乘法，首先用左边的行向量乘A^*,结果记为L(一行三列的向量)。

  • $$L=\left[\begin{array}{ccc}{x}^2(x_2-x_3)-x(x_2^2-x_3^2)+x_2{x}_3(x_2-x_3)& -x^2(x_1-x_3)+x(x_1^2-x_3^2)-x_1{x}_3(x_1-x_3)& x^2(x_1-x_2)-x(x_1^2-x_2^2)+x_1{x}_2(x_1-x_2)\end{array}\right]$$

    提取公因子后将L化简为如下对称的形式：
    $$L=\left[\begin{array}{ccc}(x_2-x_3)(x-x_2)(x-x_3)& (x_3-x_1)(x-x_3)(x-x_1)& (x_1-x_2)(x-x_1)(x-x_2)\end{array}\right]$$

  • 此时该抛物线的方程可以表示为如下形式:

  • $$y=L\ast \frac{1}{|A|}\ast \left[\begin{array}{c}{y}_1\\ y_2\\ y_3\end{array}\right]$$

  • 将系数提到最外面先算两个向量的乘积，此时方程形式为：

  • $$y=y_1\ast \frac{(x_2-x_3)(x-x_2)(x-x_3)}{|A|}+y_2\ast \frac{(x_3-x_1)(x-x_3)(x-x_1)}{|A|}+y_3\ast \frac{(x_1-x_2)(x-x_1)(x-x_2)}{|A|}$$

  • 求行列式的值时根据有y(x₁) = y₁，所以把x₁代入上式可以发现后两项的值均为零。即可得如下关系

  • $$y_1=y_1\ast \frac{(x_2-x_3)(x_1-x_2)(x_1-x_3)}{|A|}$$

  • 所以可以直接得到行列式的值为：

  • $$|A|=(x_2-x_3)(x_1-x_2)(x_1-x_3)$$

  • 不放心的话也很好验证，用线代里面行列式的解法算一下就可以了。

  • 此时就可以把|A|代入抛物线的方程中，发现每个分子都可以和分母约掉一项。

  • 所以最后就得到了该抛物线的方程为：

  • $$y=y_1\ast \frac{(x-x_2)(x-x_3)}{(x_1-x_2)(x_1-x_3)}+y_2\ast \frac{(x-x_3)(x-x_1)}{(x_2-x_1)(x_2-x_3)}+y_3\ast \frac{(x-x_1)(x-x_2)}{(x_3-x_1)(x_3-x_2)}$$

  • 其中在约分的时候可能出现一些正负号的问题。如第二项约分时(x₁-x₃)与（x₃-x₁)差个符号，这时可以调整分母的(x₁-x₂)变为(x₂-x₁)，保持形式的一致。这就是三点插值的公式的推导过程。