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,'-')

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/198747.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Linux:使用pv实现执行进度监控

pv全称&#xff1a;Pipe Viewer&#xff0c;通过管道显示数据处理进度的信息 安装 yum install pv -y示例 复制文件 # 显示进度 pv data.sql > ./data-new.sql330MiB 0:00:00 [1.32GiB/s] [>] 100%限制mysql数据导出速率 mysqldump | pv -L10m > data.sql# -L, -…

gitlab注册无中国区电话验证问题

众所周知gitlab对中国区不友好&#xff0c;无法直接注册&#xff0c;页面无法选择86的手机号进行验证码发送。 Google上众多的方案是修改dom&#xff0c;而且时间大约是21年以前。 修改dom&#xff0c;对于现在的VUE、React框架来说是没有用的&#xff0c;所以不用尝试。 直接看…

Docker 安装 Nacos

Docker 安装 Nacos tags: docker Nacos 文章目录 Docker 安装 Nacostags: docker Nacos 下载镜像docker 运行命令说明 Nacos 端口说明 下载镜像 下载镜像 # 搜索镜像 docker search nacos # 下载镜像 dockers pull nacos/nacos-server创建挂载文件夹 # 存放日志 mkdir naco…

Linux结束程序运行的命令

kill 通过进程 ID&#xff08;PID&#xff09;结束一个程序的运行。例如&#xff0c;要结束进程 ID 为 1234 的进程&#xff1a; kill 1234 pkill 通过进程名称结束一个程序的运行。例如&#xff0c;要结束名称为example_process的进程&#xff1a; pkill example_process …

postman参数为D:\\audio\\test.mp3请求报错

报错信息 报错 java.lang.IllegalArgumentException: Invalid character found in the request target [/v1/audio/transcriptions?audioPathD:\\audio\\test.mp3 ]. The valid characters are defined in RFC 7230 and RFC 3986 解决方式 yml文件上放行指定字符 relaxed-pa…

Mac电脑每次修改完java的版本后,没有成功

问题&#xff0c;本地有多个java版本8,11,15,17但是每次执行代码后&#xff0c;版本没有变化。在环境变量文件.bash_profile中设置无效 export JAVA_HOME/Library/Java/JavaVirtualMachines/<Java版本目录>/Contents/Home 这个问题通常是由于系统默认使用的Shell不同导…

安装获取mongodb

目录 本地安装 获取云上资源 获取Atlas免费数据库 本地连接数据库 在Atlas中连接数据库 本文适合初学者或mongodb感兴趣的同学来准备学习测试环境&#xff0c;或本地临时开发环境。mongodb是一个对用户非常友好的数据库。这种友好&#xff0c;不仅仅体现在灵活的数据结构和…

力扣:191. 位1的个数(Python3)

题目&#xff1a; 编写一个函数&#xff0c;输入是一个无符号整数&#xff08;以二进制串的形式&#xff09;&#xff0c;返回其二进制表达式中数字位数为 1 的个数&#xff08;也被称为汉明重量&#xff09;。 提示&#xff1a; 请注意&#xff0c;在某些语言&#xff08;如 J…

评论功能实现方案

构建高效且安全的评论功能&#xff1a;实现方案探讨。 1、分析 我们以b站的评论为例&#xff0c;用下图来解释我们评论的分级。 我们可以抽出存储评论的数据表属性 评论id父级id评论作者id被回复用户ID评论帖子ID评论内容创建时间 可以设计如下的数据表 其中pid表示父id。 …

考研失利后,我是如何零基础转行测试开发 ,成功拿下独角兽公司offer?

想当年&#xff0c;从一个什么都不懂的非科班测试小白&#xff0c;考研失利后&#xff0c;转行到K12教育知名互联网公司做测试开发工程师&#xff0c;我用了大概半年的时间。 这个过程中我自己也摸索出了一条学习路线&#xff0c;在这里想给大家分享一下我的学习路线&#xff…

EasyExcel list<Map>批量导出多个sheet

1 列表List<Map<String,Object>> list 按类型分模块分别导出到各自sheet 2 首先获取列表 List<Map<String,Object>> list tzBusiTaskUserListService.getTaskUserList(tzBusiTaskUserList); 3 分组 Map<String,List<Map<String,Object&g…

Hadoop学习笔记(HDP)-Part.16 安装HBase

目录 Part.01 关于HDP Part.02 核心组件原理 Part.03 资源规划 Part.04 基础环境配置 Part.05 Yum源配置 Part.06 安装OracleJDK Part.07 安装MySQL Part.08 部署Ambari集群 Part.09 安装OpenLDAP Part.10 创建集群 Part.11 安装Kerberos Part.12 安装HDFS Part.13 安装Ranger …

vue3 vue-router过渡动效 滚动行为 (四)

文章目录 一、过渡动效1.1安装animate.css1.2 利用元信息存储过渡名称1.3 在组件中使用 二、滚动行为2.1 始终滚动到顶部2.2 相对于某个元素的偏移量2.3 保持之前的滚动位置 一、过渡动效 1.1安装animate.css npm install animate.css --save1.2 利用元信息存储过渡名称 {pa…

ROS opencv PCL Ceres-solver之间版本对应关系

ROS1 : neotic Opencv : 4.6.0 Ceres-solver : 2.0.0

ABAP 报表工具栏缺少小计按钮

解决方案&#xff1a; 在sap标准程序 SAPLKKBL 中有多个标准的的状态栏 都有小计按钮 复制过来之后却不显示&#xff0c;调试发现&#xff0c; 在 pf_status_alv里面做了excluding &#xff0c;需要把小计排除 调试RT_EXTAB. 说明程序默认给隐藏了 不显示&#xff0c;删除调…

教你用Python+selenium搭建自动化测试环境

一、环境搭建 1、安装pythonpycharm软件 。python安装网址官网&#xff1a;About Python™ | Python.org 根据自己的电脑系统选择最新版本 下载到本地&#xff0c;选择安装路径并配置好环境变量 验证安装是否成功 搜索中录入cmd 打开命令窗口 录入python显示一下版本号表示…

微服务实战系列之Cache(技巧篇)

前言 凡工具必带使用说明书&#xff0c;如不合理的使用&#xff0c;可能得到“意外收获”。这就好比每个人擅长的领域有所差异&#xff0c;如果放错了位置或用错了人&#xff0c;也一定会让 Leader 们陷入两难之地&#xff1a;“上无法肩负领导之重托&#xff0c;下难免失去伙伴…

软件测试基础知识总结

之前有将基础的软件测试知识做了一个总结&#xff0c;但比较潦草&#xff0c;很多内容只是一笔带过&#xff0c;快到年底了&#xff0c;自己也有个写年终知识总结文档的计划&#xff0c;就将基础的理论知识重新整理一番。 有人问我&#xff0c;这些都是能搜索到的知识&#xf…

用IDEA创建Java类时自动生成作者注释等信息

File->settings-->Editor->File and Code Templates->Includes->File Header(双击)&#xff0c;然后在右边输入框内输入代码即可 /*** Author 作者名* Date ${DATE} ${TIME}* version 1.0* 注释*/

C++日常遇到的一些坑的总结

一、const 相关 C中const的不同位置的用法 const 修饰符用法总结 二、函数形参没有变量名 三、指针偏移问题 笔记&#xff1a; 包含来自C标准库的头文件&#xff0c;用#inlcude<xxx>&#xff0c;包含不来自C标准库的头文件&#xff0c;用#include"xxx"最…