matlab Allan方差分析文件

Allan方差分析的m代码文件，里面有可以测试的data.mat文件，主要步骤在allan.m和nihe.m文件中

function [retval s errorb] = allan（datatauname）

% Compute the Allan deviation for a set of time-domain frequency data
% [RETVAL S ERRORB] = ALLAN（DATATAUNAME）
% DATA should be a struct and have the following fields:
%  DATA.freq    The frequency measurements in Hz
%  DATA.rate or DATA.time
%               The sampling rate in Hertz （DATA.rate） or a timestamp for
%               each measurement in seconds （DATA.time）. Computation is
%               shorter when the rate is known but if the rate is
%               inaccurate then the Allan plot will be skewed.
%               DATA.rate is used if both fields are present.
%               If DATA.rate == 0 then the timestamps are used.
% TAU is an array of tau values for computing Allan deviation.
% NAME is a label that is added to the plot titles.
%
% RETVAL is the array of Allan deviation values at each TAU.
% S is an optional output of other statistical measures of the data （mean std etc）.
% ERRORB is an optional output containing the error estimates for a 1-sigma
%   confidence interval. Error bars are plotted as vertical lines at each point
%   （MATLAB-style error bars cannot be used due to a known bug in MATLAB R14SP3）.
%
% Example:
%
% To compute the Allan deviation for the data in the variable “lt“:
% >> lt
% lt = 
%     freq: [1x86400 double]
%     rate: 0.50
%
% Use:
%
% >> ad = allan（lt[1 10 100]‘lt data‘）;
%
% The Allan deviation will be computed and plotted at tau = 110100 seconds.
%  1-sigma confidence intervals will be indicated by vertical lines.
%
% Notes:
%  No pre-processing of the data is performed. 
%  For rate-based data AD is computed only for tau values greater than the
%   minimum time between samples and less than the half the total time. For
%   time-stamped data only tau values greater than the maximum gap between
%   samples and less than half the total time are used.
%  The calculation for rate-based data is *much* faster than for time-stamp
%   data. You may wish to run the rate-based calculation first then
%   compare with time-stamp-based. Often the differences are insignificant.
%  To plot the “tau bins“ uncomment the code at the beginning of the
%   “plot“ section （last section of code search for “TAUBIN“）.
%  This function has been validated using the test data from NBS Monograph
%   140 and the 1000-point test data set given by Riley [1].
%   If you have other validation results please let me know!
%
% For more information see:
% [1] W. J. Riley “Addendum to a test suite for the calculation of time domain
%  frequency stability“ presented at IEEE Frequency Control Symposium
%  1996.
% Available on the web:
%  http://www.ieee-uffc.org/freqcontrol/paper1ht.html
%
%
% M.A. Hopcroft
%      hopcroft at mems stanford edu
%
%
% MH APR2008
% v1.61  improve error handling plotting
%        fix bug in regular data calc for high-rate data
%        fix bug in timestamp data calc for large starting gap
%         （thanks to C. B. Ruiz for identifyi

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1788  2018-05-14 19:52  allen1\allan_sim_3.m
     文件         846  2018-05-14 19:52  allen1\allantrend.m
     文件       14112  2018-05-14 19:52  allen1\allan-w.m
     文件        1282  2018-05-14 19:52  allen1\allen.m
     文件        1459  2018-05-14 19:52  allen1\allen2.m
     文件         443  2018-05-14 19:52  allen1\D.m
     文件      644654  2018-05-14 19:52  allen1\data.mat
     文件     1440162  2018-05-14 19:52  allen1\data.txt
     文件         130  2018-05-14 19:52  allen1\dataavg.m
     文件          85  2018-05-14 19:52  allen1\data-avg.m
     文件         616  2018-05-14 19:52  allen1\datacraet.m
     文件         247  2018-05-14 19:52  allen1\k.mat
     文件         350  2018-05-14 19:52  allen1\k.txt
     文件         145  2018-05-14 19:52  allen1\kCreat.m
     文件         222  2018-05-14 19:52  allen1\nihe.asv
     文件         224  2018-05-14 19:52  allen1\nihe.m
     文件          34  2018-05-14 19:52  allen1\t.m
     文件       90000  2018-05-14 19:52  allen1\xacc.txt
     文件       90000  2018-05-14 19:52  allen1\xraw.txt
     文件       90132  2018-05-14 19:52  allen1\yacc.txt
     文件       90000  2018-05-14 19:52  allen1\yraw.txt
     文件       90000  2018-05-14 19:52  allen1\zacc.txt
     文件       90000  2018-05-14 19:52  allen1\zraw.txt