关于拉格朗日插值算法的定义可以查看维基百科拉格朗日插值法。
这里直接上代码,当前代码不是最优算法

class Algorithm:@staticmethoddef LagrangeInterpolation():grid_x = np.array([4, 5, 6])  # the defining pointsk = 3   # the number of defining points of Lagrange polyvalue = np.array([10, 5.25, 1]) # the corresponding value on each defining pointx = 18  # the point whose value we are interestedresult = 0  # later to save final resultfor j in range(k):result_l = 1for i in range(k):if i != j:result_l = result_l * (x - grid_x[i]) / (grid_x[j] - grid_x[i])result = result + value[j] * result_lreturn resultif __name__ == '__main__':main = Algorithm()result_ = main.LagrangeInterpolation()print(result_)
"""
result:
-11.0
"""

我们可以看到最后的结果为-11.0,该结果与维基百科上给出的例子结果一致。因此,我们成功给出了当x=18时,对应通过这三个点的曲线所对应的值,即我们成功地使用拉格朗日插值法完成了插值。

码字不易,如果大家觉得有用,请高抬贵手给一个赞让我上推荐让更多的人看到吧~

Python实现拉格朗日插值算法相关推荐

  1. python之拉格朗日插值算法完成计算

    一.提出问题 拉格朗日插值法,如何计算数值? 二.详解拉格朗日插值算法 插值法就是利用已知的点(x,y)建立合适的插值函数y=f(x),将未知点w带入插值函数 f(x)可以求出函数值 f(w).对于平 ...

  2. 函数插值的python实现——拉格朗日、牛顿插值

    函数插值的python实现--拉格朗日.牛顿插值 1. 拉格朗日(Larange)插值 实现代码 2. 牛顿(Larange)插值 2.1 牛顿插值多项式的基本形式 2.2 牛顿均差插值多项式 (1) ...

  3. 【数值分析】python实现拉格朗日插值

    基于python实现拉格朗日插值,可自定义节点数量n. 一.拉格朗日插值公式 二.python代码 import matplotlib.pyplot as plt import numpy as np ...

  4. Python用拉格朗日法进行数据横向插值

    在工作中,好多数据存在缺失值,对于Excel 中对数据采用列插值,张良均的<Python数据分析与挖掘实战>已有代码,但是我今天遇到数据需要横向插值. import pandas as p ...

  5. Python实现拉格朗日求解最优问题【椭圆内接长方体的最大体积】

    文章目录 一.问题描述 二.拉格朗日手工求解问题 三.Python 编程实现 四.参考

  6. 拉格朗日乘数法及python实现拉格朗日乘数法

    拉格朗日乘数法(Lagrange Multiplier Method)基本思想 作为一种优化算法,拉格朗日乘子法主要用于解决约束优化问题,它的基本思想就是通过引入拉格朗日乘子来将含有n个变量和k个约束 ...

  7. python缺失值拉格朗日插值_拉格朗日插值-python

    在数据库中,有些数据是异常值或者空值,这些值在分析的时候应该特殊处理,比如最简单的忽略掉或者通过算法推测它的值.其中拉格朗日插值就是通过其他已经知道的值,对x位置缺失的值插入的算法. 假定我们已经知道 ...

  8. Python实现拉格朗日插值法

    已知sinx的一组x,y对应关系,用拉格朗日插值法估计sin(0.3367)的值. x x0.32 0.34 0.36 y 0.314567 0.333487 0.352274 / / class I ...

  9. python实现拉格朗日插值公式_Python实现拉格朗日插值法

    已知sinx的一组x,y对应关系,用拉格朗日插值法估计sin(0.3367)的值. x x0.32 0.34 0.36 y 0.314567 0.333487 0.352274 / / class I ...

最新文章

  1. 枚举 ---- B. Power Sequence[Codeforces Round #666 (Div. 2)][暴力]
  2. “人工智能”浪潮退去,“AI四小龙”怎么了?
  3. 最速下降法的C语言实现
  4. java 验证时间格式 效率_有关Java中日期格式化的效率实录
  5. SpringSecurity 跨域问题 Cors
  6. 阿里云服务器linux系统上安装git版本控制系统
  7. cv2中函数名的规则
  8. 分布式文件存储FastDFS之环境搭建篇
  9. c语言超市收银台程序,超市收银电脑操作流程
  10. 塞班S60v3版平台手机证书权限内容大解析
  11. 等级保护三级安全建设
  12. Bartender 连接MySql数据库并对模板赋值打印标签二维码
  13. 常用编码:Shift_JIS, GBK,EUCKR,Big5,UTF8,CP1252
  14. ALFA缺陷检测软件外观检测自学习人工智能软件
  15. 《数据结构与算法自学与面试指南》01-01:图灵奖得主尼古拉斯·沃斯
  16. 异构计算-1-10x10:一种异质性和提高能量效率的通用架构方法
  17. AWS GWLB集成FortiGate防火墙
  18. 霜语服务器优势服,暴风城都被毛了!魔兽世界“毛人风”所在霜语服务器联盟已经凉凉...
  19. [CISCN2019 华北赛区 Day2 Web1]Hack World
  20. Android以太网框架情景分析之启动简介

热门文章

  1. 常用激活函数的总结与比较
  2. 用计算机随机模拟掷骰子的实验,云南师范大学概率论实验报告随机事件的模拟--模拟掷均匀硬币的随机试验...
  3. 如何用css实现div高度自适应占满屏幕
  4. 爬虫python不用虚拟环境吗_Python界虚拟环境第三把神器,你不知道就out啦!
  5. Linux+shell学习记录和思维导图
  6. 【luogu P4590】游园会(DP套DP)
  7. 基于AU人脸运动单元 识别假笑男孩的真实表情
  8. 自己动手DIY macos下的绘图软件Pencil之原生菜单
  9. 微信小程序及Java后台测试登录
  10. matlab怎么验证正规阵,正规矩阵的性质及判定