由于DAVAR源代码是意大利Lorenzo Galleani教授写的，所以不方面上传网络，如果您有任何需求，请联系：何成文（cwhe_10@163.com），下面上传一个自己写的Allan方差源代码：

%--------------------------------------------------------------------------

%功能：                Allan方差（时域）-何成文(2016-8-24)

%参考论文：‘Analysis and Modeling of Inertial Sensors Using Allan Variance'

%example:

%   y=randn(100000,1)+0.00001*[1:100000]';

%   [sigma,tau,err]=avar(y,0.1);

%--------------------------------------------------------------------------

function [sigma,tau]=cwhe_AVAR(x0,tau0)

N=length(x0);

for i=1:inf

   m=2^(i-1);  

   if  (N-2*m) < 0

       break;

   end

   tau(i,1)=m*tau0;

   k=[1:N-2*m];

   sigma(i,1)=sqrt(1/(2*(N-2*m)*(m*tau0)^2) * sum([x0(k+2*m)-2*x0(k+m)+x0(k)].^2));

end

figure,plot(tau0*[1:N],x0);grid on;

xlabel('ittrm/s');ylabel('ity');

figure;loglog(tau,sigma,'b');grid on;hold on;loglog(tau,sigma,'ro');

xlabel('ittrm/s');ylabel('itsigma_Arm(tau)');

%--------------------------------------------------------------------------测试成功！

下面是动态阿兰方差的测试图：

尽管动态Allan方差的源代码没有上传，但是具体公式可参考：        

         《光学学报》论文‘基于动态Allan方差的光纤陀螺随机误差分析及算法改进’中的公式（6）

致谢：            

         对意大利Lorenzo Galleani教授的帮助表示感谢！