scipy笔记:scipy.interpolate.interp1d

2023-12-06 08:52

本文主要是介绍scipy笔记:scipy.interpolate.interp1d,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

1 主要使用方法

class scipy.interpolate.interp1d(x, y, kind='linear', axis=-1, copy=True, bounds_error=None, fill_value=nan, assume_sorted=False)

2 主要函数

x一维实数值数组,代表插值的自变量
y

N维实数值数组,其中沿着插值轴的 y 长度必须等于 x 的长度

默认的插值轴是 y 的最后一个轴

kind

str or int

指定插值类型的字符串或表示样条插值器阶数的整数

指定的插值类型有:

'linear'线性插值。这是最简单的插值形式,其中新点的值是通过直线连接相邻数据点来估计的
'nearest'最近邻插值。新点的值取自最接近的数据点的值
'nearest-up'类似于最近邻插值,但在处理半整数(如 0.5、1.5 等)时有所不同。'nearest-up' 会向上取整,而 'nearest' 则会向下取整
'zero'零阶插值。在相邻数据点之间,新点的值设置为左侧数据点的值
'slinear'次线性插值。这是线性插值的另一种形式
'quadratic'二次样条插值。新点的值基于二次多项式,这些多项式在每对相邻数据点间定义
'cubic'三次样条插值。它使用三次多项式在数据点之间进行插值,这种方法通常提供了较平滑的曲线。
'previous' 和 'next'这两种方法分别返回给定点之前和之后最近的数据点的值。'previous' 返回左侧的数据点值,而 'next' 返回右侧的数据点值
axisy 数组中对应于 x 坐标值的轴。默认值为 -1。
copy 若为 True,则该类会对 x 和 y 进行内部复制。若为 False,则使用 x 和 y 的引用。默认为 True
bounds_error

若为 True,在 x 范围外进行插值尝试时(需要外推)会引发 ValueError。

若为 False,则超出范围的值会被赋予 fill_value。

默认情况下,除非指定 fill_value="extrapolate",否则会引发错误

fill_value

如果为 ndarray(或浮点数),则在数据范围外的请求点将使用此值填充。

如果未提供,则默认为 NaN。

如果为两元素元组,则第一个元素用于 x_new < x[0],第二个元素用于 x_new > x[-1]。

使用两元素元组或 ndarray 需要 bounds_error=False

assume_sorted

如果为 False,x 的值可以任意排序,并且首先进行排序。

如果为 True,则 x 必须是单调递增的数组

3 举例

3.1 导入库&准备数据

# 导入库
import numpy as np
import matplotlib.pyplot as plt
from scipy.interpolate import *# 数据部分
x=np.arange(0,10)
y=np.exp(-x/3.0)

3.2 创建interp1d

x_new=np.arange(0,9,0.1)
y_new=f(x_new)
y_new
'''
array([1.        , 0.97165313, 0.94330626, 0.91495939, 0.88661252,0.85826566, 0.82991879, 0.80157192, 0.77322505, 0.74487818,0.71653131, 0.69621989, 0.67590847, 0.65559705, 0.63528563,0.61497421, 0.5946628 , 0.57435138, 0.55403996, 0.53372854,0.51341712, 0.49886335, 0.48430958, 0.46975582, 0.45520205,0.44064828, 0.42609451, 0.41154074, 0.39698698, 0.38243321,0.36787944, 0.35745121, 0.34702298, 0.33659475, 0.32616652,0.31573829, 0.30531006, 0.29488183, 0.2844536 , 0.27402537,0.26359714, 0.25612498, 0.24865283, 0.24118068, 0.23370852,0.22623637, 0.21876422, 0.21129206, 0.20381991, 0.19634776,0.1888756 , 0.18352157, 0.17816754, 0.17281351, 0.16745947,0.16210544, 0.15675141, 0.15139738, 0.14604335, 0.14068932,0.13533528, 0.13149895, 0.12766262, 0.12382629, 0.11998996,0.11615363, 0.11231729, 0.10848096, 0.10464463, 0.1008083 ,0.09697197, 0.09422312, 0.09147426, 0.08872541, 0.08597656,0.08322771, 0.08047886, 0.07773001, 0.07498115, 0.0722323 ,0.06948345, 0.06751381, 0.06554417, 0.06357454, 0.0616049 ,0.05963526, 0.05766562, 0.05569598, 0.05372634, 0.05175671])
'''
plt.plot(x,y,'o',x_new,y_new,'-')

3.3 不同kind不同结果

f_c=interp1d(x,y,kind='previous')
y_new=f_c(x_new)
plt.plot(x,y,'o',x_new,y_new,'-')

这篇关于scipy笔记:scipy.interpolate.interp1d的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/461238

相关文章

【学习笔记】 陈强-机器学习-Python-Ch15 人工神经网络(1)sklearn

系列文章目录 监督学习:参数方法 【学习笔记】 陈强-机器学习-Python-Ch4 线性回归 【学习笔记】 陈强-机器学习-Python-Ch5 逻辑回归 【课后题练习】 陈强-机器学习-Python-Ch5 逻辑回归(SAheart.csv) 【学习笔记】 陈强-机器学习-Python-Ch6 多项逻辑回归 【学习笔记 及 课后题练习】 陈强-机器学习-Python-Ch7 判别分析 【学

系统架构师考试学习笔记第三篇——架构设计高级知识(20)通信系统架构设计理论与实践

本章知识考点:         第20课时主要学习通信系统架构设计的理论和工作中的实践。根据新版考试大纲,本课时知识点会涉及案例分析题(25分),而在历年考试中,案例题对该部分内容的考查并不多,虽在综合知识选择题目中经常考查,但分值也不高。本课时内容侧重于对知识点的记忆和理解,按照以往的出题规律,通信系统架构设计基础知识点多来源于教材内的基础网络设备、网络架构和教材外最新时事热点技术。本课时知识

论文阅读笔记: Segment Anything

文章目录 Segment Anything摘要引言任务模型数据引擎数据集负责任的人工智能 Segment Anything Model图像编码器提示编码器mask解码器解决歧义损失和训练 Segment Anything 论文地址: https://arxiv.org/abs/2304.02643 代码地址:https://github.com/facebookresear

数学建模笔记—— 非线性规划

数学建模笔记—— 非线性规划 非线性规划1. 模型原理1.1 非线性规划的标准型1.2 非线性规划求解的Matlab函数 2. 典型例题3. matlab代码求解3.1 例1 一个简单示例3.2 例2 选址问题1. 第一问 线性规划2. 第二问 非线性规划 非线性规划 非线性规划是一种求解目标函数或约束条件中有一个或几个非线性函数的最优化问题的方法。运筹学的一个重要分支。2

【C++学习笔记 20】C++中的智能指针

智能指针的功能 在上一篇笔记提到了在栈和堆上创建变量的区别,使用new关键字创建变量时,需要搭配delete关键字销毁变量。而智能指针的作用就是调用new分配内存时,不必自己去调用delete,甚至不用调用new。 智能指针实际上就是对原始指针的包装。 unique_ptr 最简单的智能指针,是一种作用域指针,意思是当指针超出该作用域时,会自动调用delete。它名为unique的原因是这个

查看提交历史 —— Git 学习笔记 11

查看提交历史 查看提交历史 不带任何选项的git log-p选项--stat 选项--pretty=oneline选项--pretty=format选项git log常用选项列表参考资料 在提交了若干更新,又或者克隆了某个项目之后,你也许想回顾下提交历史。 完成这个任务最简单而又有效的 工具是 git log 命令。 接下来的例子会用一个用于演示的 simplegit

记录每次更新到仓库 —— Git 学习笔记 10

记录每次更新到仓库 文章目录 文件的状态三个区域检查当前文件状态跟踪新文件取消跟踪(un-tracking)文件重新跟踪(re-tracking)文件暂存已修改文件忽略某些文件查看已暂存和未暂存的修改提交更新跳过暂存区删除文件移动文件参考资料 咱们接着很多天以前的 取得Git仓库 这篇文章继续说。 文件的状态 不管是通过哪种方法,现在我们已经有了一个仓库,并从这个仓

忽略某些文件 —— Git 学习笔记 05

忽略某些文件 忽略某些文件 通过.gitignore文件其他规则源如何选择规则源参考资料 对于某些文件,我们不希望把它们纳入 Git 的管理,也不希望它们总出现在未跟踪文件列表。通常它们都是些自动生成的文件,比如日志文件、编译过程中创建的临时文件等。 通过.gitignore文件 假设我们要忽略 lib.a 文件,那我们可以在 lib.a 所在目录下创建一个名为 .gi

取得 Git 仓库 —— Git 学习笔记 04

取得 Git 仓库 —— Git 学习笔记 04 我认为, Git 的学习分为两大块:一是工作区、索引、本地版本库之间的交互;二是本地版本库和远程版本库之间的交互。第一块是基础,第二块是难点。 下面,我们就围绕着第一部分内容来学习,先不考虑远程仓库,只考虑本地仓库。 怎样取得项目的 Git 仓库? 有两种取得 Git 项目仓库的方法。第一种是在本地创建一个新的仓库,第二种是把其他地方的某个

Git 的特点—— Git 学习笔记 02

文章目录 Git 简史Git 的特点直接记录快照,而非差异比较近乎所有操作都是本地执行保证完整性一般只添加数据 参考资料 Git 简史 众所周知,Linux 内核开源项目有着为数众多的参与者。这么多人在世界各地为 Linux 编写代码,那Linux 的代码是如何管理的呢?事实是在 2002 年以前,世界各地的开发者把源代码通过 diff 的方式发给 Linus,然后由 Linus