scipy.signal.

cspline1d_eval#

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

在新的点集上计算三次样条。

dx 是旧的采样间距,而 x0 是旧的原点。换句话说,cj 代表样条系数的旧采样点(节点点)位于等间距的点上:

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

边缘使用镜像对称边界条件处理。

参数:
cjndarray

三次样条系数

newxndarray

新的点集。

dxfloat, 可选

旧的采样间距,默认值为 1.0。

x0int, 可选

旧的原点,默认值为 0。

返回:
resndarray

计算出的三次样条点。

另请参阅

cspline1d

计算秩为 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