Remove DC component from EEG signals

15 visualizzazioni (ultimi 30 giorni)
altaf adil
altaf adil il 23 Feb 2011
Commentato: Jinrong Tang il 11 Set 2019
I have tried two methods to remove the DC component from the EEG signal.
Method 1: Data_without_dc = Data_with_dc - mean (Data_with_dc)
Method 2: Data_fft = fft(Data_with_dc); Data_fft(1) = 0; Data_without_dc = ifft (Data_fft);
I think the first component of signal's FFT is actually the mean of whole signal and it represents signal amplitude at zero frequency, so by making it to zero, we can remove DC component.
But from both the methods, I am getting different results. Please tell me which method is preferable if I have to remove the DC component whose amplitude is unknown in case of EEG data.
  1 Commento
Edith
Edith il 28 Mar 2015
Please why do you want to remove the DC component from the EEG signal? I also trying to show in a plot my signals from 16 electrodes.

Accedi per commentare.

Risposta accettata

Jan
Jan il 23 Feb 2011
You can remove the DC component by removing the MEAN. This is fast, accurate and cheap.
I cannot see a difference between the methods, if you remove the imaginary roundoff errors:
x = rand(100, 1);
f = fft(x);
f(1) = 0;
x_ac = real(ifft(f));
[x_ac, x - mean(x)]
>> -0.4723 -0.4723
0.4304 0.4304
-0.4166 -0.4166
etc.
The difference is as expected in the magnitude of EPS(1).
BTW: I prefer "SUM(x)./SIZE(x,1)", because it is faster than MEAN.
  1 Commento
Jinrong Tang
Jinrong Tang il 11 Set 2019
Hi, I don't think remove the Mean is suitable for the ECG signal since it's asymmetrical and not like the sine wave. So removing the Mean will make the signal's baseline a little below the x-axis.

Accedi per commentare.

Più risposte (2)

Narayan Puthanmadam Subramaniyam
Hey altaf, if you data is in form of a matrix, then make sure to remove the entire first row and not just the first element. For example, make Data_fft(1,:) = 0. Then your answers should be identical!

VENKATA PHANIKRISHNA B
VENKATA PHANIKRISHNA B il 7 Giu 2017
Removal of DC component means Mean removal from the signal. By doing this we can eliminate the dc, i.e., 0Hz noise. Here DC component means, the signal positive half cycles average and the negative half cycle is not zero. if x is the signal then xm=x-mean(x).

Community

Più risposte nel  Power Electronics Control

Categorie

Scopri di più su EEG/MEG/ECoG in Help Center e File Exchange

Community Treasure Hunt

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

Start Hunting!

Translated by