scipy.signal.
qspline1d#
- scipy.signal.qspline1d(signal, lamb=0.0)[源代码][源代码]#
计算秩为1数组的二次样条系数。
- 参数:
- 信号ndarray
表示信号样本的秩为1的数组。
- lambfloat, 可选
平滑系数(目前必须为零)。
- 返回:
- cndarray
二次样条系数。
参见
qspline1d_eval
在新的一组点上评估二次样条。
注释
假设镜像对称边界条件,找到一维信号的二次样条系数。要从样条表示中恢复信号,请使用长度为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()