scipy.signal.
qspline1d#
- scipy.signal.qspline1d(signal, lamb=0.0)[原始碼]#
計算 rank-1 陣列的二次樣條係數。
- 參數:
- signalndarray
代表訊號樣本的 rank-1 陣列。
- lambfloat, optional
平滑係數 (目前必須為零)。
- 返回:
- cndarray
二次樣條係數。
另請參閱
qspline1d_eval
在新點集評估二次樣條。
筆記
尋找 1-D 訊號的二次樣條係數,假設鏡像對稱邊界條件。為了從樣條表示法中取回訊號,請將這些係數與長度為 3 的 FIR 視窗 [1.0, 6.0, 1.0]/ 8.0 進行鏡像對稱卷積。
範例
我們可以使用二次樣條過濾訊號,以減少並平滑高頻雜訊
>>> import numpy as np >>> import matplotlib.pyplot as plt >>> from scipy.signal import qspline1d, qspline1d_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 = qspline1d_eval(qspline1d(sig), time) >>> plt.plot(sig, label="signal") >>> plt.plot(time, filtered, label="filtered") >>> plt.legend() >>> plt.show()