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