使用 Matlab FFT 从脑电信号中提取频率答案

【问题标题】：Using Matlab FFT to extract frequencies from EEG signal使用 Matlab FFT 从脑电信号中提取频率
【发布时间】：2012-05-18 00:30:06
【问题描述】：

我是 BCI 的新手。我有一个来自Neurosky 的 Mindset EEG 设备，我将来自设备的原始数据值记录在一个 csv 文件中。我可以从 csv 读取数据并将其提取到 Matlab 中，然后应用 FFT。我现在需要从 FFT 中提取某些频率（Alpha、Beta、Theta、Gamma）。

其中 Delta = 1-3 Hz

Theta= 4-7 赫兹

阿尔法 = 8-12 赫兹

Beta = 13-30 赫兹

伽马 = 31-40 赫兹

这是我到目前为止所做的：

f = (0:N-1)*(Fs/N);
plot(rawDouble);
title ('Raw Signal');
p = abs(fft(rawDouble));
figure,plot (f,p);
title('Magnitude of FFT of Raw Signal');

谁能告诉我如何从信号中提取那些特定的频率范围？非常感谢！

【问题讨论】：

查看 MATLAB 文档以了解带宽滤波器。这应该比先通过 fft 更容易。
请大家在这里回答我的问题stackoverflow.com/questions/23167374/…

标签： matlab fft frequency

【解决方案1】：

我认为最简单的方法是在加载数据后过滤这些范围内的信号。例如

band=[30 100] eeglocal.lowpass(band(2)).highpass(band(1));

然后您可以使用选择您想要处理的时间。

这应该就是你所需要的。

【讨论】：

【解决方案2】：

您可能会发现从 MATLAB 的 periodogram 函数开始比尝试直接使用 FFT 更容易。这负责为您和各种其他实现细节窗口化数据。

【讨论】：

【解决方案3】：

为了方便使用 MatLab 分析 EEG 数据，您可以考虑使用 EEGLAB 工具箱 (http://sccn.ucsd.edu/eeglab/) 或 fieldtrip 工具箱 (http://fieldtrip.fcdonders.nl/start)。

两个工具箱都附有很好的教程：

http://sccn.ucsd.edu/eeglab/eeglabtut.html

http://fieldtrip.fcdonders.nl/tutorial

【讨论】：