EEG preprocessing using bandpass filter
14 views (last 30 days)
Show older comments
I have EEG data obtained from BCI Competition III dataset IIIA. I want to band pass this data in 8-13 HZ. Currently, the data is in the frequency range of 0-50HZ.How do I apply Butterworth filter on this data? What should be the ideal filter order?can somebody provide the code?
0 Comments
Answers (1)
Ashish Tiwari
on 28 Dec 2019
First of all let me know the sampling frequency you are working on. if your sampling frequency is 128HZ then normalise your lower cut off frequency and upper cut off frequency like this:
Fs=128;
for signals between 8-13Hz:--
[8 13]/Fs/2
then put this in butter function
[b,a]=butter(10,[8 13]/Fs/2);%where 10 is the order of the filter.The resulting bandpass and bandstop designs are of order 2n.i.e. 20th order
2 Comments
vo
on 12 Dec 2022
@Ashish Tiwari Why the bandpass [flo fhi]/128/2 ? I don't know why? why must divide 128/2? Can you explain for me? thank you
Star Strider
on 12 Dec 2022
Edited: Star Strider
on 12 Dec 2022
@vo — It is necessary to divide the frequencies in Hz by the Nyquist frequency (the highest uniquely identifiable frequency in a sampled signal, one-half the sampling frequency) in order to normalise the signal on the interval
radians/independent variable unit (seconds, centimetres, etc.), although actually to the interval
as arguments to the MATLAB filter design functions.
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1228232/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1228237/image.png)
See Also
Categories
Find more on EEG/MEG/ECoG in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!