How to correctly apply bandpass filter to EEG data?

56 vues (au cours des 30 derniers jours)
Ali Al-Saegh
Ali Al-Saegh le 29 Juin 2020
Modifié(e) : Star Strider le 29 Oct 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.

Réponses (1)

Star Strider
Star Strider le 29 Juin 2020
Modifié(e) : Star Strider le 29 Oct 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');
.

Catégories

En savoir plus sur EEG/MEG/ECoG dans Help Center et File Exchange

Community Treasure Hunt

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

Start Hunting!

Translated by