<div dir="ltr">Dear experts,<div><br></div><div>I'm basically a beginner in coding with Fieldtrip, and I'm having a problem with line noise filtering.</div><div>I would like to remove it from the heart rate signal, and I tried just typing 'yes' at the option cfg.dftfilter.</div><div>But is doesn't seem to work, since the power spectrum contains the 50Hz and her harmonics even after the filtering.</div><div><br></div><div>You can find the code in the following lines:</div><div><br></div><div><div><font face="monospace, monospace" size="1" color="#0b5394">%% ===FREQUENCY Pre Filtering VERIFICATION===</font></div><div><font face="monospace, monospace" size="1" color="#0b5394">fs=data.fsample;</font></div><div><font face="monospace, monospace" size="1" color="#0b5394"><br></font></div><div><font face="monospace, monospace" size="1" color="#0b5394">cfg = [];</font></div><div><font face="monospace, monospace" size="1" color="#0b5394">cfg.output        = 'pow';</font></div><div><font face="monospace, monospace" size="1" color="#0b5394">cfg.channel       = 'ECG 2';</font></div><div><font face="monospace, monospace" size="1" color="#0b5394">cfg.method        = 'mtmfft';</font></div><div><font face="monospace, monospace" size="1" color="#0b5394">cfg.taper         = 'hanning';</font></div><div><font face="monospace, monospace" size="1" color="#0b5394">cfg.pad           = 'nextpow2';</font></div><div><font face="monospace, monospace" size="1" color="#0b5394">cfg.foilim        = [0.3 fs/2]; % teorema di Nyquist</font></div><div><font face="monospace, monospace" size="1" color="#0b5394">freq_epoched_Cz   = ft_freqanalysis(cfg, data_epoched_Cz);</font></div><div><font face="monospace, monospace" size="1" color="#0b5394"><br></font></div><div><font face="monospace, monospace" size="1" color="#0b5394">%% ===Plot Spectra===</font></div><div><font face="monospace, monospace" size="1" color="#0b5394">figure;</font></div><div><font face="monospace, monospace" size="1" color="#0b5394">plot(freq_epoched_Cz.freq,mag2db(freq_epoched_Cz.powspctrm));</font></div><div><font face="monospace, monospace" size="1" color="#0b5394"><br></font></div><div><font face="monospace, monospace" size="1" color="#0b5394">%% ===Filtering===</font></div><div><font face="monospace, monospace" size="1" color="#0b5394">cfg            = [];</font></div><div><font face="monospace, monospace" size="1" color="#0b5394">cfg.channel    = {'ECG 2'};</font></div><div><font face="monospace, monospace" size="1" color="#0b5394">cfg.dftfilter  = 'yes';</font></div><div><font face="monospace, monospace" size="1" color="#0b5394">ECG_Filtered   = ft_preprocessing(cfg,data_Cz);</font></div><div><font face="monospace, monospace" size="1" color="#0b5394"><br></font></div><div><font face="monospace, monospace" size="1" color="#0b5394">cfg          = [];</font></div><div><font face="monospace, monospace" size="1" color="#0b5394">cfg.length   = 30; % in seconds;</font></div><div><font face="monospace, monospace" size="1" color="#0b5394">cfg.overlap  = 0;</font></div><div><font face="monospace, monospace" size="1" color="#0b5394">data_epoched_ECG_Filtered = ft_redefinetrial(cfg, ECG_Filtered);</font></div><div><font face="monospace, monospace" size="1" color="#0b5394"><br></font></div><div><font face="monospace, monospace" size="1" color="#0b5394">cfg = [];</font></div><div><font face="monospace, monospace" size="1" color="#0b5394">cfg.output        = 'pow';</font></div><div><font face="monospace, monospace" size="1" color="#0b5394">cfg.method        = 'mtmfft';</font></div><div><font face="monospace, monospace" size="1" color="#0b5394">cfg.taper         = 'hanning';</font></div><div><font face="monospace, monospace" size="1" color="#0b5394">cfg.pad           = 'nextpow2';</font></div><div><font face="monospace, monospace" size="1" color="#0b5394">cfg.foilim        = [0.3 fs/2]; % teorema di Nyquist</font></div><div><font face="monospace, monospace" size="1" color="#0b5394">freq_epoched_Cz_Filtered = ft_freqanalysis(cfg, data_epoched_ECG_Filtered);</font></div><div><font face="monospace, monospace" size="1" color="#0b5394"><br></font></div><div><font face="monospace, monospace" size="1" color="#0b5394">%% ===Plot Spectra===</font></div><div><font face="monospace, monospace" size="1" color="#0b5394">plot(freq_epoched_Cz_Filtered.freq,mag2db(freq_epoched_Cz_Filtered.powspctrm));</font></div></div><div><br></div><div>Thanks a lot.</div><div><br></div><div>Kind regards,</div><div>Alessandro Orticoni</div></div>