[FieldTrip] saw-tooth pattern of frequency spectrum

hui zhang hzhang.lib at gmail.com
Thu Sep 24 15:44:22 CEST 2015

Dear fieldtrip community,

I have encountered a problem of time-frequency transformation on continuous data. Basically, the frequency spectrum established a rhythmic pattern as indicated in figures below (the first figure and the blue line in the second figure).
The code I use are as below:

cfg = [];
    cfg.reref = 'yes';
    cfg.refchannel = 'all';
    cfg.datafile = 'subj1.dat';
    cfg.headerfile = 'subj1.vhdr';
    cfg.bsfilter = 'yes';
    cfg.bsfreq = [49 51; 99 101; 149 151; 199 201];
    data_all = ft_preprocessing(cfg);

    cfg = [];
    cfg.method = 'wavelet';
    cfg.output = 'pow';
    cfg.keeptrials = 'yes';
    cfg.foi = 1:1:200;
    cfg.toi = 0:0.01:size(data_all.time{1},2)/data_all.fsample;
    cfg.trials = 'all';
    freq = ft_freqanalysis(cfg,data_all);

There are different methods I tried to improve the result
1). highpass filter the data (0.1 Hz)
2). bandpass filter the data ([0.1 300Hz])
3). detrend the data

Among these methods, 1) and 2) works better and generate the same result, but can still see the rhythmic pattern of the frequency spectrum with lower amplitude and slower frequencies of the saw-tooth pattern (The green line in the second figure).
3) does not work at all.

Any idea of why the saw-tooth pattern appears in the frequency spectrum? How to efficiently improve the result?

[Inline image 3]


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20150924/2c3b64b6/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 191619 bytes
Desc: image.png
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20150924/2c3b64b6/attachment-0001.png>

More information about the fieldtrip mailing list