if (exist('x1')==0)
    x1=load('220121_18h48.hp53132a');
    x1(1:29641)=x1(1:29641)+x1(29641+3)-x1(29641-3);
    x1(29639:29642)=x1(29643);
    x1=x1(1:67000);
    x2=load('220125_18h20_integer_driver.hp53132a');
    [sig1,sig12,osig1,msig1,tsig1,tau1]=avar(x1,1);
  %  [sig1,sig22,osig1,msig1,tsig1,tau2]=avar(x1(69000:end),1);
    [sig1,sig22,osig1,msig1,tsig1,tau2]=avar(x2,1);
  end
  
  loglog(tau2,sig22);hold on
  loglog(tau1,sig12);
  legend('controlled','free','location','southeast')
  xlabel('integration time (s)')
  ylabel('Allan deviation (no unit)')