拉格朗日差值法python实现

聪明乖巧的小狮子

已于 2022-03-21 16:58:38 修改

阅读量3.9k

点赞数 11

分类专栏： python 计算方法文章标签： python 线性代数

于 2022-03-11 23:42:48 首次发布

本文链接：https://blog.csdn.net/m0_52295416/article/details/123436358

版权

python 同时被 2 个专栏收录

16 篇文章 2 订阅

订阅专栏

计算方法

6 篇文章 2 订阅

订阅专栏

$L a g r a n g e 插值法：$ $给定一组点 (x_{1},y_{1}),(x_{2},y_{2})...(x_{n},y_{n})$ $\ \ P_{i}(x)=\prod_{j=1,j\neq i}^{n}\frac{x-x_{j}}{x_{i}-x{j}}$ $拉格朗日插值多项式\ \ L_{n}(x)=\sum_{i=1}^{n}y_{i}P_{i}$
举个例子：假如我们只有三个点, 那么
$L_{2}(x)=\frac{(x-x_{1})(x-x_{2})}{(x_{0}-x_{1})(x_{0}-x_{2})}y_{0}+\frac{(x-x_{0})(x-x_{2})}{(x_{1}-x_{0})(x_{1}-x_{2})}y_{1}+\frac{(x-x_{0})(x-x_{1})}{(x_{2}-x_{0})(x_{2}-x_{1})}y_{2}$
python 实现
方法一：

from scipy.interpolate import lagrange
x = [1,3,5]
y = [2,10,1]
print(lagrange(x, y))  #拟合曲线
print(lagrange(x, y)(10))  #估值

结果：

        2
-2.125 x + 12.5 x - 8.375
-95.875

方法二：

x = [1, 3, 5]
y = [2, 10, 1]
def lagrange_Interpolation(x1):
    P = []
    L_n = 0
    for i in range(len(x)):
        numerator=1  #分子初始化
        denominator=1  #分母初始化
        for j in range(len(x)):
            if j!=i:
                numerator*=(x1-x[j])
                denominator*=(x[i]-x[j])

        P.append(numerator/denominator)

    for i in range(len(y)):
        L_n+=y[i]*P[i]
    return round(L_n,3)  #保留3位小数
print(lagrange_Interpolation(10))

结果：

-95.875

看拟合曲线图像：

import numpy
import matplotlib
from matplotlib import pyplot
x2=numpy.linspace(-10,10,20,endpoint=False) #从-10到10 20个点去拟合曲线
print(x2)
y2=[]
for i in range(len(x2)):
    y2.append(lagrange_Interpolation(x2[i]))
print(y2)

pyplot.plot(x2,y2)
pyplot.scatter(x,y, marker=".")
pyplot.show()