How to correctly apply bandpass filter to EEG data?

66 visualizaciones (últimos 30 días)
Ali Al-Saegh
Ali Al-Saegh el 29 de Jun. de 2020
Editada: Star Strider el 29 de Oct. de 2023
I have 25-electrodes EEG data, I used bandpass(myData, [0.5 40], 250) to apply a simple bandpass filtering on the data I have. I do not know why some time delay is exist after the application of the filter.

Respuestas (1)

Star Strider
Star Strider el 29 de Jun. de 2020
Editada: Star Strider el 29 de Oct. de 2023
The bandpass function (and its friends) in my experience have all used filtfilt to do the actual filtering, so there should not be any delay or phase distortion.
You can check it with:
[myFilteredData, bpdf] = bandpass(myData, [0.5 40], 250);
then:
myFilteredData2 = filtfilt(bpdf, myData);
and see if there is a difference when you plot them.
EDIT — (29 Oct 2023 at 21:58)
The bandpass function (and its friends) will always return an efficient elliptic filter if the 'ImpulseResponse','iir' name-value pair is added to the argument list:
[myFilteredData, bpdf] = bandpass(myData, [0.5 40], 250, 'ImpulseResponse','iir');
.

Categorías

Más información sobre EEG/MEG/ECoG en Help Center y File Exchange.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by