Contents

Signal params

fs=2048%samples/s
t=linspace(0,8,8*2048);%milliseconds
f=6%Hz

fs =

        2048


f =

     6

Signal Generation Two tones

x=.05*sin(2*pi*f*t);
x1=.05*sin(2*pi*f*t*2);
%X=cell([1, 2048*8,22]);
%X(1,:,:)={repmat(x,[22,1])'};

% X=zeros(1, 2048*8,22);
% X(1,:,:)=repmat(x,[22,1])';

% half 6Hz half 12Hz
x=[x(1:8*2048/2) x1(1:8*2048/2)];

compare two conditions

t_pretrig=EEG.xmin; t_postrig=EEG.xmax; newtimef({EEG.data(1,:,:) EEG.data(2,:,:)},EEG.pnts, [ t_pretrig t_postrig]*1000, EEG.srate,[3 .5],'baseline',[3000 4000 ],'commonbase','off' ,'basenorm','on', 'maxfreq' ,[],'nfreqs',50, 'freqscale' ,'log');

%figure; pop_newtimef( EEG, 1, 1, [-1000  5000], [3         0.5] , 'topovec', 1, 'elocs', EEG.chanlocs, 'chaninfo', EEG.chaninfo, 'baseline',[0], 'plotphase', 'off', 'padratio', 1);

Baseline from [-500 500] hence only on the 6Hz part,wavelet 3cycles increase with each freq by .5cycle

figure
[ersp,itc,powbase,times,freqs,erspboot,itcboot] = newtimef(x, 8*2048, [-500 6000], 2048, [3 .5],'baseline',[-500 500 ],'commonbase','off' ,'basenorm','on', 'maxfreq' ,20,'nfreqs',20,'padratio', 32);

Computing Event-Related Spectral Perturbation (ERSP) and
  Inter-Trial Phase Coherence (ITC) images based on 1 trials
  of 16384 frames sampled at 2048 Hz.
Each trial contains samples from -500 ms before to
  6000 ms after the timelocking event.
  Image frequency direction: normal
Using 3 cycles at lowest frequency to 10 at highest.
Generating 200 time points (-47.5 to 5547.5 ms)
Finding closest points for time variable
Time values for time/freq decomposition is not perfectly uniformly distributed
The window size used is 2281 samples (1113.77 ms) wide.
Estimating 20 linear-spaced frequencies from 3.0 Hz to 20.0 Hz.
Processing time point (of 200): 10 20 30 40 50 60 70 80 90 100 110 120
 130 140 150 160 170 180 190 200
Computing the mean baseline spectrum
Note: Add output variables to command line call in history to
      retrieve results and use the tftopo function to replot them

Baseline on [5000 6000] hence only on the 12Hz part,wavelet 3cycles increase with each freq by .5cycle

figure
[ersp,itc,powbase,times,freqs,erspboot,itcboot] = newtimef(x, 8*2048, [-500 6000], 2048, [3 .5],'baseline',[5000 6000],'commonbase','off' ,'basenorm','on', 'maxfreq' ,20,'nfreqs',50,'padratio',32);

Computing Event-Related Spectral Perturbation (ERSP) and
  Inter-Trial Phase Coherence (ITC) images based on 1 trials
  of 16384 frames sampled at 2048 Hz.
Each trial contains samples from -500 ms before to
  6000 ms after the timelocking event.
  Image frequency direction: normal
Using 3 cycles at lowest frequency to 10 at highest.
Generating 200 time points (-47.5 to 5547.5 ms)
Finding closest points for time variable
Time values for time/freq decomposition is not perfectly uniformly distributed
The window size used is 2281 samples (1113.77 ms) wide.
Estimating 50 linear-spaced frequencies from 3.0 Hz to 20.0 Hz.
Processing time point (of 200): 10 20 30 40 50 60 70 80 90 100 110 120
 130 140 150 160 170 180 190 200
Computing the mean baseline spectrum
Note: Add output variables to command line call in history to
      retrieve results and use the tftopo function to replot them

Baseline on NONE,wavelet 3cycles increase with freq by .5cycle

figure
[ersp,itc,powbase,times,freqs,erspboot,itcboot] = newtimef(x, 8*2048, [-500 6000], 2048, [3 .5],'baseline',[0],'commonbase','off' ,'basenorm','on', 'maxfreq' ,20,'nfreqs',50,'padratio', 32);

Computing Event-Related Spectral Perturbation (ERSP) and
  Inter-Trial Phase Coherence (ITC) images based on 1 trials
  of 16384 frames sampled at 2048 Hz.
Each trial contains samples from -500 ms before to
  6000 ms after the timelocking event.
  Image frequency direction: normal
Using 3 cycles at lowest frequency to 10 at highest.
Generating 200 time points (-47.5 to 5547.5 ms)
Finding closest points for time variable
Time values for time/freq decomposition is not perfectly uniformly distributed
The window size used is 2281 samples (1113.77 ms) wide.
Estimating 50 linear-spaced frequencies from 3.0 Hz to 20.0 Hz.
Processing time point (of 200): 10 20 30 40 50 60 70 80 90 100 110 120
 130 140 150 160 170 180 190 200
Computing the mean baseline spectrum
Note: Add output variables to command line call in history to
      retrieve results and use the tftopo function to replot them

Baseline from [-500 500] hence only on the 6Hz part,FFT

figure
[ersp,itc,powbase,times,freqs,erspboot,itcboot] = newtimef(x, 8*2048, [-500 6000], 2048, 0,'baseline',[-500 500 ],'commonbase','off' ,'basenorm','on', 'maxfreq' ,20,'nfreqs',50,'padratio', 32);

Computing Event-Related Spectral Perturbation (ERSP) and
  Inter-Trial Phase Coherence (ITC) images based on 1 trials
  of 16384 frames sampled at 2048 Hz.
Each trial contains samples from -500 ms before to
  6000 ms after the timelocking event.
  Image frequency direction: normal
Generating 200 time points (-93.7 to 5593.7 ms)
Finding closest points for time variable
Time values for time/freq decomposition is not perfectly uniformly distributed
The window size used is 2048 samples (1000 ms) wide.
Estimating 50 linear-spaced frequencies from 0.0 Hz to 20.0 Hz.
Processing trial (of 1):
finding closest frequencies: 32718 freqs removed
Computing the mean baseline spectrum
Note: Add output variables to command line call in history to
      retrieve results and use the tftopo function to replot them

Baseline on [5000 6000] hence only on the 12Hz part,FFT

figure
[ersp,itc,powbase,times,freqs,erspboot,itcboot] = newtimef(x, 8*2048, [-500 6000], 2048, 0,'baseline',[5000 6000 ],'commonbase','off' ,'basenorm','on', 'maxfreq' ,20,'nfreqs',50,'padratio',32);

Computing Event-Related Spectral Perturbation (ERSP) and
  Inter-Trial Phase Coherence (ITC) images based on 1 trials
  of 16384 frames sampled at 2048 Hz.
Each trial contains samples from -500 ms before to
  6000 ms after the timelocking event.
  Image frequency direction: normal
Generating 200 time points (-93.7 to 5593.7 ms)
Finding closest points for time variable
Time values for time/freq decomposition is not perfectly uniformly distributed
The window size used is 2048 samples (1000 ms) wide.
Estimating 50 linear-spaced frequencies from 0.0 Hz to 20.0 Hz.
Processing trial (of 1):
finding closest frequencies: 32718 freqs removed
Computing the mean baseline spectrum
Note: Add output variables to command line call in history to
      retrieve results and use the tftopo function to replot them

Baseline on NONE,FFT

figure
[ersp,itc,powbase,times,freqs,erspboot,itcboot] = newtimef(x, 8*2048, [-500 6000], 2048, 0,'baseline',[0 ],'commonbase','off' ,'basenorm','on', 'maxfreq' ,20,'nfreqs',50,'padratio',32);

Computing Event-Related Spectral Perturbation (ERSP) and
  Inter-Trial Phase Coherence (ITC) images based on 1 trials
  of 16384 frames sampled at 2048 Hz.
Each trial contains samples from -500 ms before to
  6000 ms after the timelocking event.
  Image frequency direction: normal
Generating 200 time points (-93.7 to 5593.7 ms)
Finding closest points for time variable
Time values for time/freq decomposition is not perfectly uniformly distributed
The window size used is 2048 samples (1000 ms) wide.
Estimating 50 linear-spaced frequencies from 0.0 Hz to 20.0 Hz.
Processing trial (of 1):
finding closest frequencies: 32718 freqs removed
Computing the mean baseline spectrum
Note: Add output variables to command line call in history to
      retrieve results and use the tftopo function to replot them

only 6Hz signal over the whole 8sec period, Baseline NONE,wavelet

3cycles increase with each freq by .5cycle

figure
x=.05*sin(2*pi*f*t);
[ersp,itc,powbase,times,freqs,erspboot,itcboot] = newtimef(x, 8*2048, [-500 6000], 2048, [3 .5],'baseline',[0 ],'commonbase','off' ,'basenorm','on', 'maxfreq' ,20,'nfreqs',50,'padratio', 32);

Computing Event-Related Spectral Perturbation (ERSP) and
  Inter-Trial Phase Coherence (ITC) images based on 1 trials
  of 16384 frames sampled at 2048 Hz.
Each trial contains samples from -500 ms before to
  6000 ms after the timelocking event.
  Image frequency direction: normal
Using 3 cycles at lowest frequency to 10 at highest.
Generating 200 time points (-47.5 to 5547.5 ms)
Finding closest points for time variable
Time values for time/freq decomposition is not perfectly uniformly distributed
The window size used is 2281 samples (1113.77 ms) wide.
Estimating 50 linear-spaced frequencies from 3.0 Hz to 20.0 Hz.
Processing time point (of 200): 10 20 30 40 50 60 70 80 90 100 110 120
 130 140 150 160 170 180 190 200
Computing the mean baseline spectrum
Note: Add output variables to command line call in history to
      retrieve results and use the tftopo function to replot them

only 6Hz signal over the whole 8sec period, Baseline [-500 500],wavelet 3cycles increase with each freq by .5cycle

figure
%x=.05*sin(2*pi*f*t);
[ersp,itc,powbase,times,freqs,erspboot,itcboot] = newtimef(x, 8*2048, [-500 6000], 2048, [3 .5],'baseline',[-500 500 ],'commonbase','off' ,'basenorm','on', 'maxfreq' ,20,'nfreqs',50,'padratio', 32);

Computing Event-Related Spectral Perturbation (ERSP) and
  Inter-Trial Phase Coherence (ITC) images based on 1 trials
  of 16384 frames sampled at 2048 Hz.
Each trial contains samples from -500 ms before to
  6000 ms after the timelocking event.
  Image frequency direction: normal
Using 3 cycles at lowest frequency to 10 at highest.
Generating 200 time points (-47.5 to 5547.5 ms)
Finding closest points for time variable
Time values for time/freq decomposition is not perfectly uniformly distributed
The window size used is 2281 samples (1113.77 ms) wide.
Estimating 50 linear-spaced frequencies from 3.0 Hz to 20.0 Hz.
Processing time point (of 200): 10 20 30 40 50 60 70 80 90 100 110 120
 130 140 150 160 170 180 190 200
Computing the mean baseline spectrum
Note: Add output variables to command line call in history to
      retrieve results and use the tftopo function to replot them

only 6Hz signal over the whole 8sec period, Baseline NONE, FFT

figure
x=.05*sin(2*pi*f*t);
[ersp,itc,powbase,times,freqs,erspboot,itcboot] = newtimef(x, 8*2048, [-500 6000], 2048, 0,'baseline',[0 ],'commonbase','off' ,'basenorm','on', 'maxfreq' ,20,'nfreqs',50,'padratio', 32);

Computing Event-Related Spectral Perturbation (ERSP) and
  Inter-Trial Phase Coherence (ITC) images based on 1 trials
  of 16384 frames sampled at 2048 Hz.
Each trial contains samples from -500 ms before to
  6000 ms after the timelocking event.
  Image frequency direction: normal
Generating 200 time points (-93.7 to 5593.7 ms)
Finding closest points for time variable
Time values for time/freq decomposition is not perfectly uniformly distributed
The window size used is 2048 samples (1000 ms) wide.
Estimating 50 linear-spaced frequencies from 0.0 Hz to 20.0 Hz.
Processing trial (of 1):
finding closest frequencies: 32718 freqs removed
Computing the mean baseline spectrum
Note: Add output variables to command line call in history to
      retrieve results and use the tftopo function to replot them

only 6Hz signal over the whole 8sec period, Baseline [-500 500],FFT

figure
%x=.05*sin(2*pi*f*t);
[ersp,itc,powbase,times,freqs,erspboot,itcboot] = newtimef(x, 8*2048, [-500 6000], 2048, 0,'baseline',[-500 500 ],'commonbase','off' ,'basenorm','on', 'maxfreq' ,20,'nfreqs',50,'padratio', 32);

Computing Event-Related Spectral Perturbation (ERSP) and
  Inter-Trial Phase Coherence (ITC) images based on 1 trials
  of 16384 frames sampled at 2048 Hz.
Each trial contains samples from -500 ms before to
  6000 ms after the timelocking event.
  Image frequency direction: normal
Generating 200 time points (-93.7 to 5593.7 ms)
Finding closest points for time variable
Time values for time/freq decomposition is not perfectly uniformly distributed
The window size used is 2048 samples (1000 ms) wide.
Estimating 50 linear-spaced frequencies from 0.0 Hz to 20.0 Hz.
Processing trial (of 1):
finding closest frequencies: 32718 freqs removed
Computing the mean baseline spectrum
Note: Add output variables to command line call in history to
      retrieve results and use the tftopo function to replot them


Published with MATLAB® 7.10