一、拉格朗日的基本思想:

二、线性插值

三、多个点

四、代码实现

def lagrange(xx,y):l=len(y)l_n = 0for k in range(l):xxx=xx.copy()x_k = xxx[k]xxx.pop(k)l_k = 1for i in range(len(xxx)):l_k *= (x - xxx[i]) / (x_k -xxx[i])l_n += y[k] * l_kreturn 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 = 0for k in range(l):xxx=xx.copy()x_k = xxx[k]xxx.pop(k)l_k = 1for i in range(len(xxx)):l_k *= (x - xxx[i]) / (x_k -xxx[i])l_n += y[k] * l_kreturn 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代码实现相关推荐

  1. 缺失值处理 - 拉格朗日插值法 - Python代码

    目录 缺失值处理 拉格朗日差值法的理论基础 拉格朗日插值法代码实现 其他数据预处理方法 缺失值处理 处理缺失值常用的办法可分为三类:删除记录.数据插补.不处理. 其中常见的数据插补法有: 如果通过删除 ...

  2. 解读 拉格朗日插值法python,保你学明白

    文章目录 概述拉格朗日插值法 什么是插值法 拉格朗日插值法的原理 拉格朗日公式 拉格朗日插值法的代码实现 Python 进行拉格朗日插值的主要知识点 Polyfit 函数 Polyval 函数 Lin ...

  3. 拉格朗日插值python代码_Python实现的拉格朗日插值法示例

    本文实例讲述了Python实现的拉格朗日插值法.分享给大家供大家参考,具体如下: 拉格朗日插值简单介绍 拉格朗日插值法是以法国十八世纪数学家约瑟夫·拉格朗日命名的一种多项式插值方法. 许多实际问题中都 ...

  4. 拉格朗日插值法--python

    数据分析 数据清洗:缺失值处理.1删除记录 2数据插补 3不处理 数据 链接:https://pan.baidu.com/s/1jiIOoselsqVQR4P_EaS3pA 提取码:t970 常见插补 ...

  5. 拉格朗日插值法+python实现

    拉格朗日插值法 1.数学原理 1.1 线性插值 1.2二次插值 1.3 n次拉格朗日插值多项式 2.Python实现 本文将介绍拉格朗日插值法得数学原理,并用Python实现. 1.数学原理 首先我们 ...

  6. 拉格朗日插值法python实现

    目录 1.原理 2.涉及的Python库 3.例子 1.原理 对某个多项式函数有已知的k+1个点,假设任意两个不同的都互不相同,那么应用拉格朗日插值公式所得到的拉格朗日插值多项式为: L(x)=∑j= ...

  7. python数值积分_数值积分 python代码实现

    老规矩,数学原理什么的就不写了. 直接贴代码和实例演示,以下代码基于python和numpy. 在这里,我将用代码实现复化梯形算法.复化 Simpson 算法.Romberg 积分算法和三点 Gaus ...

  8. 拉格朗日插值法《python数据分析与挖掘实践》

    数据清洗 数据清洗主要是删除原始数据集中地无关数据.重复数据.平滑噪声数据,筛选掉与挖掘主题无关的数据,处理缺失值.异常值等. 缺失值处理的方法可分为3类:删除记录.数据插补和不处理. 删除含有缺失值 ...

  9. python拉格朗日插值法_Python实现的拉格朗日插值法示例

    本文实例讲述了Python实现的拉格朗日插值法.分享给大家供大家参考,具体如下: 拉格朗日插值简单介绍 拉格朗日插值法是以法国十八世纪数学家约瑟夫·拉格朗日命名的一种多项式插值方法. 许多实际问题中都 ...

最新文章

  1. offsetTop,clientX,clientTop,clientWidth,offsetWidth 坐标,一次弄明白
  2. 经典C语言程序100例之二
  3. 像Excel一样使用python进行数据分析(2)
  4. centOS 5.5 安装 ORACLE 11G RAC 问题汇总
  5. github上面如何编辑README
  6. 浅谈web应用的负载均衡、集群、高可用(HA)解决方案
  7. matlab repmat_利用MATLAB实现KNN
  8. 计算机视觉农作物检测,基于计算机视觉的农作物病害检测系统的研究
  9. Oracle 12C 新特性之扩展数据类型(extended data type)
  10. #研发中间件介绍#异步消息可靠推送Notify
  11. VMProtect修复导入表的插件
  12. 数值积分方法的c语言实现,编程实现数值积分的几种--方法 c语言.doc
  13. (转) Csrss进程剖析
  14. 联系人备份--vcf
  15. 阿铭Linux_网站维护学习笔记201903026
  16. 计算机一级用写在简历里面吗,写简历要小心,简历填写最忌讳的是什么,你知道吗?...
  17. 软件工程之美学习笔记二十五 24 | 技术债务:是继续修修补补凑合着用,还是推翻重来?
  18. php 网状结构,数据库模型-数据结构-网状模型
  19. 手机您的浏览器不允许第三方_手机壁纸:我的世界不允许你的消失,不管结局是否完美...
  20. java判断是不是质数(素数)超短方法

热门文章

  1. 网桥 - Linux实现交换机
  2. 归一法的计算方法讲解_第七讲 归一法
  3. javaMail学习(四)——使用javaMail给网易账户发简单邮件
  4. 软件测试别再被“薪资陷阱”困扰了,这份攻略带你轻松查薪资
  5. 软间隔和支持向量机回归
  6. 2021橙子百度域名URL批量采集软件【急速】
  7. 更精确SSR的交点检测
  8. 玩转基因组浏览器之自定义IGV的参考基因组
  9. 整理了一些可以发外链的网站,分享几十个发外链的网站
  10. oracle sql stddev,STDDEV和VARIANCE这两个函数是怎么用