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()