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' 返回右侧的数据点值
axis y 数组中对应于 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,'-')

相关推荐

  1. 数据科学:Scipy、Scikit-Learn笔记

    2023-12-06 14:34:02       56 阅读
  2. pytorch笔记:ReplicationPad1d

    2023-12-06 14:34:02       34 阅读
  3. scipy.sparse.lil_matrix

    2023-12-06 14:34:02       31 阅读

最近更新

  1. docker php8.1+nginx base 镜像 dockerfile 配置

    2023-12-06 14:34:02       98 阅读
  2. Could not load dynamic library ‘cudart64_100.dll‘

    2023-12-06 14:34:02       106 阅读
  3. 在Django里面运行非项目文件

    2023-12-06 14:34:02       87 阅读
  4. Python语言-面向对象

    2023-12-06 14:34:02       96 阅读

热门阅读

  1. 【C/C++】可变参数va_list与格式化输出

    2023-12-06 14:34:02       58 阅读
  2. Windows如何启动MySQL

    2023-12-06 14:34:02       52 阅读
  3. Pytorch:torch.utils.data.random_split()

    2023-12-06 14:34:02       57 阅读
  4. VB.NET二维数组的组合

    2023-12-06 14:34:02       73 阅读
  5. 通俗讲解分布式锁:场景和使用方法

    2023-12-06 14:34:02       46 阅读
  6. 2312skia,12画布包与路径包

    2023-12-06 14:34:02       48 阅读
  7. 大数据(十一):概率统计基础

    2023-12-06 14:34:02       60 阅读
  8. 关于SQL注入问题及解决--小记

    2023-12-06 14:34:02       50 阅读
  9. html复习

    2023-12-06 14:34:02       51 阅读
  10. Python中的global关键字

    2023-12-06 14:34:02       52 阅读
  11. MIT 6.S081学习笔记(第五章)

    2023-12-06 14:34:02       39 阅读
  12. js实现tab栏切换

    2023-12-06 14:34:02       56 阅读