scipy.signal.

cspline1d_eval#

scipy.signal.cspline1d_eval(cj, newx, dx=1.0, x0=0)[source]#

在新點集上評估三次樣條。

dx 是舊的樣本間距,而 x0 是舊的原點。換句話說,cj 代表樣條係數的舊樣本點(節點)位於等間隔的點上:

oldx = x0 + j*dx j=0…N-1, 其中 N=len(cj)

邊緣使用鏡像對稱邊界條件處理。

參數:
cjndarray

三次樣條係數

newxndarray

新的點集。

dxfloat, optional

舊的樣本間距,預設值為 1.0。

x0int, optional

舊的原點,預設值為 0。

回傳值:
resndarray

評估後的三次樣條點。

另請參閱

cspline1d

計算 rank-1 陣列的三次樣條係數。

範例

我們可以過濾訊號,以三次樣條減少並平滑高頻雜訊

>>> 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()
../../_images/scipy-signal-cspline1d_eval-1.png