Python SciPy signal.cspline1d用法及代码示例

本文简要介绍 python 语言中 scipy.signal.cspline1d 的用法。

用法:

scipy.signal.cspline1d(signal, lamb=0.0)#

计算 rank-1 数组的三次样条系数。

假设 mirror-symmetric 边界条件，求一维信号的三次样条系数。为了从样条表示 mirror-symmetric-convolve 中获得信号，这些系数的长度为 3 FIR 窗口 [1.0, 4.0, 1.0]/6.0 。

参数 ：：

signal： ndarray

表示信号样本的 rank-1 数组。

lamb： 浮点数，可选

平滑系数，默认为 0.0。

返回 ：：

c： ndarray

三次样条系数。

例子：

我们可以使用三次样条过滤信号以减少和平滑high-frequency 噪声：

>>> import numpy as np
>>> import matplotlib.pyplot as plt
>>> from scipy.signal import cspline1d, cspline1d_eval
>>> rng = np.random.default_rng()
>>> sig = np.repeat([0., 1., 0.], 100)
>>> sig += rng.standard_normal(len(sig))*0.05  # add noise
>>> time = np.linspace(0, len(sig))
>>> filtered = cspline1d_eval(cspline1d(sig), time)
>>> plt.plot(sig, label="signal")
>>> plt.plot(time, filtered, label="filtered")
>>> plt.legend()
>>> plt.show()