allan.m

光纤陀螺输出误差的allan方差分析

%Allan方差分析
clear;
clc
load TG400_5.dat;
c=TG400_5/11399*3600;   %转换成角速率，标度因数11300
tt=0.3;       %数据输出周期
NA=length(c);
NN=NA/7;
for M=1:NN;    %
    K=NA/M;
    for i=1:K
       b_aver(i)=mean(c(M*(i-1)+1:M*i)); 
    end
    si=0;
    for j=1:K-1
       si=si+(b_aver(j+1)-b_aver(j))^2;
    end
    sig(M)=si/(K-1)/2;
end
figure(2);
loglog((1:NN)*tt,sig(1:NN));
%grid on;%打开分格线
tau=(1:NN)*tt;
M=2;%阶数
C=nihe(tau,sig,M);
C
Q=sqrt(abs(C(M-1)/3))*(1000000/180/3600);
N=sqrt(abs(C(M)/1))/60; 
Bs=sqrt(abs(C(M+1))/(log(2)*2/pi));
K=sqrt(abs(C(M+2)*3))*60; 
R=sqrt(abs(C(M+3)*2))*3600; 
Q
N
Bs
K
R
X=tau';Y=sig';
YY=zeros(length(X),1);
for i=-M:M
YY=C(M+i+1).*X.^(i)+YY;
end
YY=sqrt(YY);sig=sqrt(sig);
loglog(tau,sig);%grid on;
hold on;loglog(X,YY,'r');xlabel('相关时间对数\tau');ylabel('Allan方差对数\sigma');title('Allan Analysis');
figure(1);
plot(1:NA,c);
xlabel('采样点');ylabel('输出信号');title('输出信号波形');