scipy.signal.

qspline1d_eval

scipy.signal.qspline1d_eval(cj, newx, dx=1.0, x0=0)

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

参数:

cjndarray

二次样条系数

newxndarray

新的点集。

dxfloat，可选

旧的采样间隔，默认值为 1.0。

x0int，可选

旧的起点，默认值为 0。

返回:

resndarray

计算的二次样条点。

另请参阅

qspline1d

计算 1 阶阵列的二次样条系数。

注释

dx 是旧的采样间隔，而 x0 是旧的起点。换句话说，其中的 cj 表示样条系数的旧样例点（结点）在等间距点

oldx = x0 + j*dx  j=0...N-1, with N=len(cj)

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

示例

在浏览器中尝试！

我们可以用二次样条滤波一个信号，以减少和消除高频噪声

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

返回在新标签页中打开