基于Python实现拉格朗日插值法
作者：泡泡怡 发布时间：2022-03-18 18:52:18
标签：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（）函数
七、结果展示
来源：https://blog.csdn.net/m0_72662900/article/details/128278422
0
投稿猜你喜欢
前言MySQL支持单机事务的良好表现毋庸置疑，那么在分布式系统中，涉及多个节点，MySQL又是如何实现分布式事务的呢？比如开发一个业务系统，- 用语言实现 好处： 1、可以减少对数据库的访问。 2、可移植性好。 坏处： 1、操作起来考虑的东西较多，修改一处就要修改别一处。也就是说是相
- LBound函数：返回指定数组维的最小可用下标；UBound函数：返回指定数组维的最小可用上标；LBound(arrayname, [,di
- 现在网上有很多python2写的爬虫抓取网页图片的实例，但不适用新手（新手都使用python3环境，不兼容python2），所以我用Pyth
- 本文实例讲述了Python基于list的append和pop方法实现堆栈与队列功能。分享给大家供大家参考，具体如下：#coding=utf8
- 栗子：计算斐波那契数列（任一个数都是前两个数之和的数字序列）Python2.7实现代码如下：<strong><span s
- 本文讲述了Python检测网络延迟的代码。分享给大家供大家参考，具体如下：#!/usr/bin/env python # coding: u
- 目录1、设定答题卡模板2、读取答题卡图像并对图像进行灰度化处理3、高斯模糊图像去噪点4、使用大津法二值分割图像5、使用开运算去噪点6、使用c
- 为了实现这样的小工具，我们先设想有下面这些功能。1、可以自定义整理某一个路径下面的所有需要被整理的文件。2、默认情况下，使用文件后缀作为同一
- 前言你知道a**3是什么意思吗？你知道怎么打印模块路径吗？你知道reduce（）函数怎么用和列表怎么拼接字符串吗？然而你可能全都知道但是还是
- 相对于C++的继承编写，Python更简洁，而且效率也是很高的，下面编写一个简单Python的继承例子。#!/usr/bin/python&
- Python函数参数类型检查有一个很经典的笑话：三个月之前，只有我和上帝知道这代码是干什么的。现在，只有上帝知道了。在Python中，不知道
- 本文实例为大家分享了Javascript实现单选框效果的具体代码，供大家参考，具体内容如下描述: 点击每一个li 将li的内容赋值给 div
- 信号和槽机制是 QT 的核心机制，要精通 QT 编程就必须对信号和槽有所了解。信号和槽是一种高级接口，应用于对象之间的通信，它是 QT 的核
- 自上一篇文章 Z Story : Using Django with GAE Python 后台抓取多个网站的页面全文 后，大体的进度如下：
- 项目中遇到一个需求，要把中国式的显示阿拉伯数字的方式改为欧式的，即每三位显示，中间用逗号隔开，比如12345678改成12，345，678的
- 以下是Python基础学习内容的学习笔记的全部内容，非常的详细，如果你对Python语言感兴趣，并且针对性的系统学习一下基础语言知识，下面的
- NMS 算法在目标检测，目标定位领域有较广泛的应用。算法原理非极大值抑制算法（Non-maximum suppression, NMS）的本
- json和dictpython中的dict类型要转换为json格式的数据需要用到json库：import json <json>
- 这次，我们再来用Ursina引擎来做一个太阳系行星模拟器吧！想要了解Ursina 3D引擎的基本使用方法的话，查看我的另一篇文章：详解Pyt
