Python实现拉格朗日插值算法
关于拉格朗日插值算法的定义可以查看维基百科拉格朗日插值法。
这里直接上代码,当前代码不是最优算法
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实现拉格朗日插值算法相关推荐
- python之拉格朗日插值算法完成计算
一.提出问题 拉格朗日插值法,如何计算数值? 二.详解拉格朗日插值算法 插值法就是利用已知的点(x,y)建立合适的插值函数y=f(x),将未知点w带入插值函数 f(x)可以求出函数值 f(w).对于平 ...
- 函数插值的python实现——拉格朗日、牛顿插值
函数插值的python实现--拉格朗日.牛顿插值 1. 拉格朗日(Larange)插值 实现代码 2. 牛顿(Larange)插值 2.1 牛顿插值多项式的基本形式 2.2 牛顿均差插值多项式 (1) ...
- 【数值分析】python实现拉格朗日插值
基于python实现拉格朗日插值,可自定义节点数量n. 一.拉格朗日插值公式 二.python代码 import matplotlib.pyplot as plt import numpy as np ...
- Python用拉格朗日法进行数据横向插值
在工作中,好多数据存在缺失值,对于Excel 中对数据采用列插值,张良均的<Python数据分析与挖掘实战>已有代码,但是我今天遇到数据需要横向插值. import pandas as p ...
- Python实现拉格朗日求解最优问题【椭圆内接长方体的最大体积】
文章目录 一.问题描述 二.拉格朗日手工求解问题 三.Python 编程实现 四.参考
- 拉格朗日乘数法及python实现拉格朗日乘数法
拉格朗日乘数法(Lagrange Multiplier Method)基本思想 作为一种优化算法,拉格朗日乘子法主要用于解决约束优化问题,它的基本思想就是通过引入拉格朗日乘子来将含有n个变量和k个约束 ...
- python缺失值拉格朗日插值_拉格朗日插值-python
在数据库中,有些数据是异常值或者空值,这些值在分析的时候应该特殊处理,比如最简单的忽略掉或者通过算法推测它的值.其中拉格朗日插值就是通过其他已经知道的值,对x位置缺失的值插入的算法. 假定我们已经知道 ...
- Python实现拉格朗日插值法
已知sinx的一组x,y对应关系,用拉格朗日插值法估计sin(0.3367)的值. x x0.32 0.34 0.36 y 0.314567 0.333487 0.352274 / / class I ...
- python实现拉格朗日插值公式_Python实现拉格朗日插值法
已知sinx的一组x,y对应关系,用拉格朗日插值法估计sin(0.3367)的值. x x0.32 0.34 0.36 y 0.314567 0.333487 0.352274 / / class I ...
最新文章
- 枚举 ---- B. Power Sequence[Codeforces Round #666 (Div. 2)][暴力]
- “人工智能”浪潮退去,“AI四小龙”怎么了?
- 最速下降法的C语言实现
- java 验证时间格式 效率_有关Java中日期格式化的效率实录
- SpringSecurity 跨域问题 Cors
- 阿里云服务器linux系统上安装git版本控制系统
- cv2中函数名的规则
- 分布式文件存储FastDFS之环境搭建篇
- c语言超市收银台程序,超市收银电脑操作流程
- 塞班S60v3版平台手机证书权限内容大解析
- 等级保护三级安全建设
- Bartender 连接MySql数据库并对模板赋值打印标签二维码
- 常用编码:Shift_JIS, GBK,EUCKR,Big5,UTF8,CP1252
- ALFA缺陷检测软件外观检测自学习人工智能软件
- 《数据结构与算法自学与面试指南》01-01:图灵奖得主尼古拉斯·沃斯
- 异构计算-1-10x10:一种异质性和提高能量效率的通用架构方法
- AWS GWLB集成FortiGate防火墙
- 霜语服务器优势服,暴风城都被毛了!魔兽世界“毛人风”所在霜语服务器联盟已经凉凉...
- [CISCN2019 华北赛区 Day2 Web1]Hack World
- Android以太网框架情景分析之启动简介
热门文章
- 常用激活函数的总结与比较
- 用计算机随机模拟掷骰子的实验,云南师范大学概率论实验报告随机事件的模拟--模拟掷均匀硬币的随机试验...
- 如何用css实现div高度自适应占满屏幕
- 爬虫python不用虚拟环境吗_Python界虚拟环境第三把神器,你不知道就out啦!
- Linux+shell学习记录和思维导图
- 【luogu P4590】游园会(DP套DP)
- 基于AU人脸运动单元 识别假笑男孩的真实表情
- 自己动手DIY macos下的绘图软件Pencil之原生菜单
- 微信小程序及Java后台测试登录
- matlab怎么验证正规阵,正规矩阵的性质及判定