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)')