基于Python实现拉格朗日插值法

 更新时间:2022年12月12日 09:02:37   作者:泡泡怡  
拉格朗日插值法是以法国十八世纪数学家约瑟夫·拉格朗日命名的一种多项式插值方法。本文将利用Python语言实现这一插值法,需要的可以参考一下

一、拉格朗日的基本思想

拉格朗日插值法是以法国十八世纪数学家约瑟夫·拉格朗日命名的一种多项式插值方法。

许多实际问题中都用函数来表示某种内在联系或规律,而不少函数都只能通过实验和观测来了解。在若干个不同的地方得到相应的观测值,拉格朗日插值法可以找到一个简单函数,其恰好在各个现测的点取到观测到的值,这个函数可以是代数多项式,三角多项式等。

二、线性插值

三、多个点

四、代码实现

def lagrange(xx,y):
    l=len(y)
    l_n = 0
    for k in range(l):
        xxx=xx.copy()
        x_k = xxx[k]
        xxx.pop(k)
        l_k = 1
        for i in range(len(xxx)):
            l_k *= (x - xxx[i]) / (x_k -xxx[i])
        l_n += y[k] * l_k
    return expand(l_n) 

五、完整代码

import matplotlib.pyplot as plt
import numpy as np 
import pandas as pd
from sympy import expand
from sympy.abc import x
xx=[]
for i in range(7):
    xx.append(data['x'][i])
y=[]
for j in range(7):
    y.append(data['y'][j])
def lagrange(xx,y):
    l=len(y)
    l_n = 0
    for k in range(l):
        xxx=xx.copy()
        x_k = xxx[k]
        xxx.pop(k)
        l_k = 1
        for i in range(len(xxx)):
            l_k *= (x - xxx[i]) / (x_k -xxx[i])
        l_n += y[k] * l_k
    return expand(l_n)  
lagrange_interpolation_polynomial = lagrange(xx, y)
print("拉格朗日插值多项式为:",lagrange_interpolation_polynomial)
x2=np.linspace(-1,4,100)
y1=[]
for i in range(len(x2)):
    y1.append(lagrange_interpolation_polynomial.subs(x,x2[i]))
print(y1)
#绘制散点图,逼近函数
plt.figure(figsize=(8,4))
plt.scatter(xx,y,c='red')
plt.plot(x2,y1,'-')
plt.show()

六、pop()函数

七、结果展示

到此这篇关于基于Python实现拉格朗日插值法的文章就介绍到这了,更多相关Python拉格朗日插值法内容请搜索脚本之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持脚本之家!

相关文章

  • python时间序列数据相减的实现

    python时间序列数据相减的实现

    本文主要介绍了python时间序列数据相减的实现,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧
    2023-04-04
  • Python unittest生成测试报告过程解析

    Python unittest生成测试报告过程解析

    这篇文章主要介绍了Python unittest生成测试报告过程解析,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下
    2020-09-09
  • Django 全局的static和templates的使用详解

    Django 全局的static和templates的使用详解

    这篇文章主要介绍了Django 全局的static和templates的使用详解,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下
    2019-07-07
  • python中pandas.DataFrame对行与列求和及添加新行与列示例

    python中pandas.DataFrame对行与列求和及添加新行与列示例

    pandas是python环境下最有名的数据统计包,而DataFrame翻译为数据框,是一种数据组织方式,这篇文章主要给大家介绍了python中pandas.DataFrame对行与列求和及添加新行与列的方法,文中给出了详细的示例代码,需要的朋友可以参考借鉴,下面来一起看看吧。
    2017-03-03
  • Python中的数据类dataclass解读

    Python中的数据类dataclass解读

    这篇文章主要介绍了Python中的数据类dataclass使用,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教
    2024-01-01
  • 利用Pytorch实现获取特征图的方法详解

    利用Pytorch实现获取特征图的方法详解

    这篇文章主要为大家详细介绍了如何利用Pytorch实现获取特征图,包括提取单个特征图和提取多个特征图,文中的示例代码讲解详细,需要的可以参考一下
    2022-10-10
  • Pandas数据处理库画图与文件读取使用示例

    Pandas数据处理库画图与文件读取使用示例

    这篇文章主要为大家介绍了Pandas数据处理库画图与文件读取使用示例,有需要的朋友可以借鉴参考下,希望能够有所帮助,祝大家多多进步,早日升职加薪
    2022-10-10
  • Pycharm Terminal 与Project interpreter 安装包不同步问题解决

    Pycharm Terminal 与Project interpreter 安装

    本文主要介绍了Pycharm Terminal 与Project interpreter 安装包不同步问题解决,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧
    2023-02-02
  • python绘制y关于x的线性回归线性方程图像实例

    python绘制y关于x的线性回归线性方程图像实例

    这篇文章主要为大家介绍了python绘制y关于x的线性回归线性方程图像实例,有需要的朋友可以借鉴参考下,希望能够有所帮助,祝大家多多进步,早日升职加薪
    2023-10-10
  • ubuntu安装jupyter并设置远程访问的实现

    ubuntu安装jupyter并设置远程访问的实现

    Jupyter Notebook是Ipython的升级版,而Ipython可以说是一个加强版的交互式Shell,本文主要介绍了ubuntu安装jupyter并设置远程访问的实现,感兴趣的可以了解一下
    2022-03-03

最新评论