python - 低频带通滤波

Question

我想知道是否有任何最佳实践来解决我认为可能引起普遍兴趣的以下问题。

我有一个以 40 kHz 采样约 60 秒的信号。该信号包含高达 5kHz 的频率，但我只对它的 40Hz 分量感兴趣。

如何获得稳定的、窄带通滤波（大约 40 Hz）信号，它也与原始信号具有相同的时基（零滞后）。滤波器还应具有线性相位。

计算工作不是问题。

标签： pythonpython-2.7matlabsignal-processingbandpass-filter

Answer 1

使用： https ://docs.scipy.org/doc/scipy/reference/generated/scipy.signal.iirpeak.html设计一个窄带通滤波器 ，然后使用前向后向滤波器 filtfilt 在信号上实现它： https://docs.scipy.org/doc/scipy-0.14.0/reference/generated/scipy.signal.filtfilt.html。

前向后向滤波需要对信号进行分段，但您始终可以生成具有较大重叠的部分。