Использование Matlab FFT для извлечения частот из сигнала ЭЭГ

Я новичок в BCI. У меня есть ЭЭГ-устройство Mindset от Neurosky, и я записываю необработанные значения данных, поступающие с устройства, в CSV-файл. Я могу читать и извлекать данные из csv в Matlab и применять БПФ. Теперь мне нужно извлечь определенные частоты (альфа, бета, тета, гамма) из БПФ.

Где Дельта = 1–3 Гц

Тета = 4–7 Гц

Альфа = 8–12 Гц

Бета = 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');

Кто-нибудь может сказать мне, как извлечь эти конкретные частотные диапазоны из сигнала?? Большое тебе спасибо!