MAMEM / eeg-processing-toolbox

Matlab code for proccesing EEG signals.
http://www.mamem.eu/
Apache License 2.0
74 stars 36 forks source link

"fft method" for dataset II #63

Open ShivaNaghsh opened 5 years ago

ShivaNaghsh commented 5 years ago

Hi, I am working with Dataset II, and I want to see "the effct of Stimulus frequency" on "each Trial" using "fft method" in MATLAB. I mean, I know the label of each Trial, so I know the stimulus frequency used, and I expect to see a significant Peak regarding "the effect of Stimulus frequency on the EEG signal of each Trial " on the "frequency plot" (for example: a peak on 10 Hz or 6.66 Hzz or ...); however, I can't see such a peak at all. I am wondering why I can't see such a peak!!! The code I used to apply fft to signals of each trial is this:

close all; lbl=sess.trials{1,55}.label; %oz num is 126 sig=sess.trials{1,55}.signal(126,:); %oz num is 126 fs = 250; x = sig; nfft = 2^( nextpow2(length(x)) ); df = fs/nfft; f = 0:df:fs/2; X = fft(x,nfft); X = X(1:nfft/2+1); plot(f,abs(X)) xlim([0 20]);

Thank you in advance for any help. Best regards

ShivaNaghsh commented 5 years ago

Dear @liarosge and @vangelis2015, I have to work with fft methods to extract the frequency components of SSVEP for this stage of my project. I would really appreciate it if you guide me to solve my problem.

vangelis2015 commented 5 years ago

Hello,

In general it is difficult to see, by naked eye, the effect of stimulus frequency by adopting this simple approach. Perhaps in some trials and subjects, you can observe it. You must not forgot that the EEG trials are noisy. In the case that you need some trials where the effect can be observed by naked eye then i suggest to apply your fft procedure in all trials and choose the desired trials for further processing. Also, since the SSVEP responses includes, also, multiples of frequency of the visual stimulus, i suggest to remove the last line of your program (ie. xlim). Finally, instead of using the fft function to find the spectrum you can use as an alternative the pwelch function.

ShivaNaghsh commented 5 years ago

@vangelis2015, Thank you for your response. I will consider the points you mentioned.