拉格朗日插值代码python实现(不掉包)

今天我们来讲一下,使用拉格朗日插值公式进行插值,通过python实现
那么拉格朗日插值公式是什么样的呢?
百度百科定义如下:

当然如果你没有看懂的话,可以再去查询一些相关资料
那么它的python代码实现如下:


import numpy as np
import os
x=[0,2,4,5,6,7]
y=[2,7,10.5,12.8,14.6,18.6]def inter(x,e):result=0for i in  range(len(x)):# print(i)x_i=x[i]y_i=y[i]x_c=x.copy()#  y_c=y.copy()x_c.pop(i)# y_c.pop(i)#     print(x_c)X=np.array(x_c)Y=y_iX_f=X-x_iX_e=X-e#print(np.cumprod(X_f))#print(np.cumprod(X_e))#print(np.cumprod(X_f)[-1])#print(np.cumprod(X_e)[-1])result=result+np.cumprod(X_e)[-1]/np.cumprod(X_f)[-1]*Y#print("***")#print(np.cumprod(X_e)[-1]/np.cumprod(X_f)[-1]*Y)return result
result=inter(x,2.5)
print(result)os.system("pause")

拉格朗日插值代码python实现(不掉包)相关推荐

  1. matlab编写拉格朗日插值代码函数

    要求:根据拉格朗日多项式插值法原理,设计算法流程并且编写拉格朗日插值代码函数. 代码如下: function[y]=lagrange(x0,y0,x)   %建立一个函数名为lagrange的函数,输 ...

  2. [数值计算-14]:拉格朗日插值与Python代码实现

    作者主页(文火冰糖的硅基工坊):https://blog.csdn.net/HiWangWenBing 本文网址:https://blog.csdn.net/HiWangWenBing/article ...

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

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

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

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

  5. 拉格朗日插值和牛顿插值的龙格现象

    文章目录 一.实验目的 二.实验设备信息 三.实验内容 (一)拉格朗日插值多项式 (二)牛顿插值多项式 四.实验步骤 (一)拉格朗日插值函数实现 (二)牛顿插值函数实现 (三)观察拉格朗日插值和牛顿插 ...

  6. Matlab实现线性插值、抛物插值、牛顿插值、拉格朗日插值、分段抛物插值、分段线性插值

    目录 线性插值 原理 流程图 代码 抛物插值 原理 流程图 代码 拉格朗日插值 代码 牛顿插值 原理 代码 分段线性插值 代码 线性插值 原理 流程图 单个点的线性插值代码 X=[0.2 0.4]; ...

  7. 【插值】牛顿插值、拉格朗日插值、三次样条插值的Python代码实现

    插值简介 插值即根据有限的离散点绘制出穿过所有样本点的曲线,从直观上想象似乎画一条穿过n个特定点的曲线有无数种画法,但从数学意义上来说我们希望画出的曲线能够尽量平滑,震荡幅度尽量小能够在非样本点上符合 ...

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

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

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

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

最新文章

  1. 第三节 线程传参详解、detach()大坑、成员函数做线程函数
  2. matlab fft simulink,在Simulink中复制MATLAB的FFT?
  3. Spring在Java领域的统治地位:86% 的Java开发者依赖它!
  4. ubuntu 修改卷标
  5. 水泵怎么做_泳池设备日常怎么维护和保养?
  6. NOIP模拟测试15「建造城市city(插板法)·轰炸·石头剪刀布」
  7. Oracle/PLSQL While Loop
  8. python之变量的私密处理
  9. python免费教程视频-微软推出 Python 免费在线教程视频
  10. QT 图片与base64互转
  11. JS浏览器对象-Location对象
  12. redis基础知识——菜鸟教程
  13. 超市管理系统java_java实现超市管理系统
  14. vue使用百度编辑器ueditor,ueditor1.5.0下载
  15. 度分秒转度数和度数转度分秒
  16. Google离开我们快十年了
  17. visual assist x太卡了_河西区专业达格XD7500租赁 X-RAY -哪家好-老友网
  18. ppt 计算机图标不见了,显示桌面图标不见了如何恢复.ppt
  19. 一个例子学会使用Jetpack Compose Modifier
  20. js将数字转换成中文大写

热门文章

  1. GraphicsLab Project之HDR渲染
  2. DataGrip深入实战教程
  3. 计算机毕业设计太难不会做怎么办?
  4. 大话微服务:Spring Cloud gateway+OAuth2 实现单点登录和权限控制(二) OAuth2.0 四种模式的通俗理解
  5. Oracle SGA PGA UGA
  6. 超越光速:时间旅行不是梦
  7. plotly.js 常见图形使用 常见图形操作 折线图 热力图 轮廓图 泡泡图 图点击事件
  8. WebGl之.wrl生成.stl文件
  9. 如何编写更高层次嵌入式C代码?
  10. 名词解释第八讲:公钥