cwhe10的个人博客分享 http://blog.sciencenet.cn/u/cwhe10

博文

0016:动态Allan方差源代码及效果图展示

已有 6789 次阅读 2016-9-1 23:06 |个人分类:科学研究|系统分类:论文交流| DAVAR, 动态Allan方差

     由于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教授的帮助表示感谢!







https://blog.sciencenet.cn/blog-2777453-1000248.html

上一篇:0016:现代信号谱分析(下篇:源代码与图)
下一篇:0017:全球六大卫星导航系统ICD文件汇总
收藏 IP: 202.193.56.*| 热度|

1 黄仁勇

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-3-29 21:17

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部