使用Scipy选择与特定频率相关的样本

Question

我正在研究以1kHz采样的EEG数据集。我已经应用了4-40 Hz的带通滤波器，并使用2秒的滑动窗口提取段进行FFT。现在我想选择对应于5-35 Hz的200个样本。有人请指导我如何使用Scipy或Numpy做到这一点。

提前谢谢你。

Answer 1

对您提取的数据应用5-35 Hz的带通滤波器。您可以使用scipy.signal.filtfilt函数（它可以避免为信号添加相位）：

https://docs.scipy.org/doc/scipy-0.14.0/reference/generated/scipy.signal.filtfilt.html

您可以使用Butterworth滤波器系数（使用原始采样率）提供此功能。您可以使用scipy.signal.butter函数。