clc
  close all
  clear all
  
  Ts = 0.2;
  N = 1e6;
  D = 5;
  
  data_raw.rate = 1/Ts;
  data_raw.freq = 1/2.*noise(N, 1, 'white') + 1.*noise(N, 1, 'pink') + 1.*noise(N, 1, 'brown') + 0e-8.*[0:N-1]';
  t_raw = [0:length(data_raw.freq)-1];
  [ad_raw, S, err_raw, tau_raw] = allan(data_raw, 2.^[0:20].*Ts, 'raw',0);
  
  data_avg.rate = 1/(Ts*D);
  data_avg.freq = mean(reshape(data_raw.freq', [D, N/D]))';
  t_avg = D.*[1:length(data_avg.freq)];
  [ad_avg, S, err_avg, tau_avg] = allan(data_avg, 2.^[0:20].*Ts*D, 'avg',0);
  
  data_decim.rate = 1/(Ts*D);
  data_decim.freq = reshape(data_raw.freq', [D, N/D])(1,:)';
  t_decim = D.*[0:length(data_decim.freq)-1];
  [ad_decim, S, err_decim, tau_decim] = allan(data_decim, 2.^[0:20].*Ts*D, 'decim',0);
  
  figure
  hold on
  loglogerr(tau_raw, ad_raw, err_raw, 'b')
  loglogerr(tau_decim, ad_decim, err_decim, 'm')
  loglogerr(tau_avg, ad_avg, err_avg, 'r')
  grid on
  
  figure
  hold on
  plot(t_raw, data_raw.freq, 'b')
  plot(t_decim, data_decim.freq, 'm')
  plot(t_avg, data_avg.freq, 'r')
  grid on