ALLAN方差的matlab代码

function [ AB ] = Allan（ c ）
%UNtitleD 此处显示有关此函数的摘要
%   此处显示详细说明
%Allan方差分析
%clear;
%clc
%load c.dat;
%c=a;   %转换成角速率，标度因数11300
tt=0.001;       %数据输出周期tt=0.01
% c=kalman（33000.1c）;     %kalman（RQmeasuredValue ）
c=RemoveNaN（c）;  %去除Nan
   c=kalman（330.1c）;  %进行卡尔曼滤波

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;%sig
end

figure（2）;
loglog（（1:NN）*ttsig（1:NN））;
%grid on;%打开分格线

tau=（1:NN）*tt;
M=3;%阶数此处是3阶拟合

A（: