TheAllan variance (AVAR), also known astwo-sample variance, is a measure offrequency stability inclocks,oscillators andamplifiers. It is named afterDavid W. Allan and expressed mathematically as${\displaystyle {\sigma }_y^2(\tau )}$.TheAllan deviation (ADEV), also known assigma-tau, is the square root of the Allan variance,${\displaystyle {\sigma }_y(\tau )}$.

TheM-sample variance is a measure of frequency stability usingM samples, timeT between measurements and observation time${\displaystyle \tau }$.M-samplevariance is expressed as

${\displaystyle {\sigma }_y^2(M,T,\tau ).}$

The Allan variance is intended to estimate stability due to noise processes and not that of systematic errors or imperfections such as frequency drift or temperature effects. The Allan variance and Allan deviation describe frequency stability. See also the sectionInterpretation of value below.

There are also different adaptations or alterations of Allan variance, notably themodified Allan variance MAVAR or MVAR, thetotal variance, and theHadamard variance. There also exist time-stability variants such astime deviation (TDEV) or time variance (TVAR). Allan variance and its variants have proven useful outside the scope oftimekeeping and are a set of improved statistical tools to use whenever the noise processes are not unconditionally stable, thus a derivative exists.

The generalM-sample variance remains important, since it allowsdead time in measurements, and bias functions allow conversion into Allan variance values. Nevertheless, for most applications the special case of 2-sample, or "Allan variance" with${\displaystyle T=\tau }$ is of greatest interest.

When investigating the stability ofcrystal oscillators andatomic clocks, it was found that they did not have aphase noise consisting only ofwhite noise, but also offlicker frequency noise. These noise forms become a challenge for traditional statistical tools such asstandard deviation, as the estimator will not converge. The noise is thus said to be divergent. Early efforts in analyzing the stability included both theoretical analysis and practical measurements.^[2][3]

An important side consequence of having these types of noise was that, since the various methods of measurements did not agree with each other, the key aspect of repeatability of a measurement could not be achieved. This limits the possibility to compare sources and make meaningful specifications to require from suppliers. Essentially all forms of scientific and commercial uses were then limited to dedicated measurements, which hopefully would capture the need for that application.

To address these problems, David Allan introduced theM-sample variance and (indirectly) the two-sample variance.^[4] While the two-sample variance did not completely allow all types of noise to be distinguished, it provided a means to meaningfully separate many noise-forms for time-series of phase or frequency measurements between two or more oscillators. Allan provided a method to convert between anyM-sample variance to anyN-sample variance via the common 2-sample variance, thus making allM-sample variances comparable. The conversion mechanism also proved thatM-sample variance does not converge for largeM, thus making them less useful.IEEE later identified the 2-sample variance as the preferred measure.^[5]

An early concern was related to time- and frequency-measurement instruments that had adead time between measurements. Such a series of measurements did not form a continuous observation of the signal and thus introduced asystematic bias into the measurement. Great care was spent in estimating these biases. The introduction of zero-dead-time counters removed the need, but the bias-analysis tools have proved useful.

Another early aspect of concern was related to how thebandwidth of the measurement instrument would influence the measurement, such that it needed to be noted. It was later found that by algorithmically changing the observation${\displaystyle \tau }$, only low${\displaystyle \tau }$ values would be affected, while higher values would be unaffected. The change of${\displaystyle \tau }$ is done by letting it be an integer multiple${\displaystyle n}$ of the measurementtimebase${\displaystyle {\tau }_0}$:

${\displaystyle \tau =n{\tau }_0.}$

The physics ofcrystal oscillators were analyzed by D. B. Leeson,^[3] and the result is now referred to asLeeson's equation. The feedback in theoscillator will make thewhite noise andflicker noise of the feedback amplifier and crystal become thepower-law noises of${\displaystyle f^{-2}}$ white frequency noise and${\displaystyle f^{-3}}$ flicker frequency noise respectively. These noise forms have the effect that thestandard variance estimator does not converge when processing time-error samples. The mechanics of the feedback oscillators was unknown when the work on oscillator stability started, but was presented by Leeson at the same time as the set of statistical tools was made available byDavid W. Allan. For a more thorough presentation on theLeeson effect, see modern phase-noise literature.^[6]

Allan variance is defined as one half of the time average of the squares of the differences between successive readings of thefrequency deviation sampled over the sampling period. The Allan variance depends on the time period used between samples, therefore, it is a function of the sample period, commonly denoted asτ, likewise the distribution being measured, and is displayed as a graph rather than a single number. A low Allan variance is a characteristic of a clock with good stability over the measured period.

Allan deviation is widely used for plots (conventionally inlog–log format) and presentation of numbers. It is preferred, as it gives the relative amplitude stability, allowing ease of comparison with other sources of errors.

An Allan deviation of 1.3×10⁻⁹ at observation time 1 s (i.e.τ = 1 s) should be interpreted as there being an instability in frequency between two observations 1 second apart with a relativeroot mean square (RMS) value of 1.3×10⁻⁹. For a 10 MHz clock, this would be equivalent to 13 mHz RMS movement. If the phase stability of an oscillator is needed, then thetime deviation variants should be consulted and used.

One may convert the Allan variance and other time-domain variances into frequency-domain measures of time (phase) and frequency stability.^[7]

Given a time-series${\displaystyle x(t)}$, for any positive real numbers${\displaystyle T,\tau }$, define the real number sequence$${\displaystyle {\overline{y}}_i=\frac{x(iT+\tau )-x(iT)}{\tau }i=0,1,2,...}$$Then the${\displaystyle M}$-sample variance is defined^[4] (here in a modernized notation form) as theBessel-corrected variance of the sequence${\displaystyle {\overline{y}}_0,...,{\overline{y}}_{M-1}}$:$${\displaystyle {\sigma }_y^2(M,T,\tau )=\frac{M}{M-1}\left(\frac{1}{M}\sum _{i=0}^{M-1}{\overline{y}}_i^2-{\left[\frac{1}{M}\sum _{i=0}^{M-1}{\overline{y}}_i\right]}^2\right),}$$The interpretation of the symbols is as follows:

• ${\displaystyle t}$ is the reading on a reference clock (in arbitrary units).
• ${\displaystyle x(t)}$ is the reading of a clock we are testing (in arbitrary units), as a function of the reference clock's reading. It can also be interpreted as theaverage fractional frequencytime series.
• ${\displaystyle {\overline{y}}_n}$ is then^thfractional frequency average over the observation time${\displaystyle \tau }$.
• ${\displaystyle M}$ is the number of clock reading intervals used in computing the${\displaystyle M}$-sample variance,
• ${\displaystyle T}$ is the time between each frequency sample,
• ${\displaystyle \tau }$ is the time length of each frequency estimate, or the observation period.

Dead-time can be accounted for by letting the time${\displaystyle T}$ be different from that of${\displaystyle \tau }$.

The Allan variance is defined as

${\displaystyle {\sigma }_y^2(\tau )=⟨{\sigma }_y^2(2,\tau ,\tau )⟩=\frac{1}{2}⟨{\left({\overline{y}}_{n+1}-{\overline{y}}_n\right)}^2⟩=\frac{1}{2{\tau }^2}⟨{\left(x_{n+2}-2x_{n+1}+x_n\right)}^2⟩}$

where${\displaystyle x_n:=x(n\tau )}$ and${\displaystyle ⟨\cdots ⟩}$ denotes the expectation operator.

The condition$T=\tau$ means the samples are taken with no dead-time between them.

Just as withstandard deviation andvariance, the Allan deviation is defined as the square root of the Allan variance:

${\displaystyle {\sigma }_y(\tau )=\sqrt{{\sigma }_y^2(\tau )}.}$

The oscillator being analysed is assumed to follow the basic model of

${\displaystyle V(t)=V_0\sin (\mathrm{\Phi }(t)).}$

The oscillator is assumed to have a nominal frequency of${\displaystyle {\nu }_{\text{n}}}$, given in cycles per second (SI unit:hertz). The nominalangular frequency${\displaystyle {\omega }_{\text{n}}}$ (in radians per second) is given by

${\displaystyle {\omega }_{\text{n}}=2\pi {\nu }_{\text{n}}.}$

The total phase can be separated into a perfectly cyclic component${\displaystyle {\omega }_{\text{n}}t}$, along with a fluctuating component${\displaystyle \phi (t)}$:

${\displaystyle \mathrm{\Phi }(t)={\omega }_{\text{n}}t+\phi (t)=2\pi {\nu }_{\text{n}}t+\phi (t).}$

The time-error functionx(t) is the difference between expected nominal time and actual normal time:

${\displaystyle x(t)=\frac{\phi (t)}{2\pi {\nu }_{\text{n}}}=\frac{\mathrm{\Phi }(t)}{2\pi {\nu }_{\text{n}}}-t=T(t)-t.}$

For measured values a time-error series TE(t) is defined from the reference time functionT_ref(t) as

${\displaystyle TE(t)=T(t)-T_{\text{ref}}(t).}$

The frequency function${\displaystyle \nu (t)}$ is the frequency over time, defined as

${\displaystyle \nu (t)=\frac{1}{2\pi }\frac{d\mathrm{\Phi }(t)}{dt}.}$

The fractional frequencyy(t) is the normalized difference between the frequency${\displaystyle \nu (t)}$ and the nominal frequency${\displaystyle {\nu }_{\text{n}}}$:

${\displaystyle y(t)=\frac{\nu (t)-{\nu }_{\text{n}}}{{\nu }_{\text{n}}}=\frac{\nu (t)}{{\nu }_{\text{n}}}-1.}$

The average fractional frequency is defined as

${\displaystyle \overline{y}(t,\tau )=\frac{1}{\tau }{\int }_0^{\tau }y(t+t_v)d{t}_v,}$

where the average is taken over observation timeτ, they(t) is the fractional-frequency error at timet, andτ is the observation time.

Sincey(t) is the derivative ofx(t), we canwithout loss of generality rewrite it as

${\displaystyle \overline{y}(t,\tau )=\frac{x(t+\tau )-x(t)}{\tau }.}$

This definition is based on the statisticalexpected value, integrating over infinite time. The real-world situation does not allow for such time-series, in which case a statisticalestimator needs to be used in its place. A number of different estimators will be presented and discussed.

A first simple estimator would be to directly translate the definition into

${\displaystyle {\sigma }_y^2(\tau ,M)=\mathrm{AVAR}(\tau ,M)=\frac{1}{2(M-1)}\sum _{i=0}^{M-2}({\overline{y}}_{i+1}-{\overline{y}}_i{)}^2,}$

or for the time series:

${\displaystyle {\sigma }_y^2(\tau ,N)=\mathrm{AVAR}(\tau ,N)=\frac{1}{2{\tau }^2(N-2)}\sum _{i=0}^{N-3}(x_{i+2}-2x_{i+1}+x_i{)}^2.}$

These formulas, however, only provide the calculation for theτ =τ₀ case. To calculate for a different value ofτ, a new time-series needs to be provided.

Taking the time-series and skipping pastn − 1 samples, a new (shorter) time-series would occur withτ₀ as the time between the adjacent samples, for which the Allan variance could be calculated with the simple estimators. These could be modified to introduce the new variablen such that no new time-series would have to be generated, but rather the original time series could be reused for various values ofn. The estimators become

${\displaystyle {\sigma }_y^2(n{\tau }_0,M)=\mathrm{AVAR}(n{\tau }_0,M)=\frac{1}{2\frac{M-1}{n}}\sum _{i=0}^{\frac{M-1}{n}-1}{\left({\overline{y}}_{ni+n}-{\overline{y}}_{ni}\right)}^2}$

with${\displaystyle n\le \frac{M-1}{2}}$,

and for the time series:

${\displaystyle {\sigma }_y^2(n{\tau }_0,N)=\mathrm{AVAR}(n{\tau }_0,N)=\frac{1}{2n^2{\tau }_0^2\left(\frac{N-1}{n}-1\right)}\sum _{i=0}^{\frac{N-1}{n}-2}{\left(x_{ni+2n}-2x_{ni+n}+x_{ni}\right)}^2}$

with${\displaystyle n\le \frac{N-1}{2}}$.

These estimators have a significant drawback in that they will drop a significant amount of sample data, as only 1/n of the available samples is being used.

A technique presented by J. J. Snyder^[8] provided an improved tool, as measurements were overlapped inn overlapped series out of the original series. The overlapping Allan variance estimator was introduced by Howe, Allan and Barnes.^[9] This can be shown to be equivalent to averaging the time or normalized frequency samples in blocks ofn samples prior to processing. The resulting predictor becomes

or for the time series:

${\displaystyle {\sigma }_y^2(n{\tau }_0,N)=\mathrm{AVAR}(n{\tau }_0,N)=\frac{1}{2n^2{\tau }_0^2(N-2n)}\sum _{i=0}^{N-2n-1}(x_{i+2n}-2x_{i+n}+x_i{)}^2.}$

The overlapping estimators have far superior performance over the non-overlapping estimators, asn rises and the time-series is of moderate length. The overlapped estimators have been accepted as the preferred Allan variance estimators in IEEE,^[5] ITU-T^[10] and ETSI^[11] standards for comparable measurements such as needed for telecommunication qualification.

In order to address the inability to separate white phase modulation from flicker phase modulation using traditional Allan variance estimators, an algorithmic filtering reduces the bandwidth byn. This filtering provides a modification to the definition and estimators and it now identifies as a separate class of variance calledmodified Allan variance. The modified Allan variance measure is a frequency stability measure, just as is the Allan variance.

A time stability (σ_x) statistical measure, which is often called the time deviation (TDEV), can be calculated from the modified Allan deviation (MDEV). The TDEV is based on the MDEV instead of the original Allan deviation, because the MDEV can discriminate between white and flicker phase modulation (PM). The following is the time variance estimation based on the modified Allan variance:

${\displaystyle {\sigma }_x^2(\tau )=\frac{{\tau }^2}{3}{mod\sigma }_y^2(\tau ),}$

and similarly for modified Allan deviation totime deviation:

${\displaystyle {\sigma }_x(\tau )=\frac{\tau }{\sqrt{3}}{mod\sigma }_y(\tau ).}$

The TDEV is normalized so that it is equal to the classical deviation for white PM for time constantτ = τ₀. To understand the normalization scale factor between the statistical measures, the following is the relevant statistical rule: For independent random variablesX andY, the variance (σ_z²) of a sum or difference (z =x −y) is the sum square of their variances (σ_z² = σ_x² + σ_y²). The variance of the sum or difference (y =x_2τ −x_τ) of two independent samples of a random variable is twice the variance of the random variable (σ_y² = 2σ_x²). The MDEV is the second difference of independent phase measurements (x) that have a variance (σ_x²). Since the calculation is the double difference, which requires three independent phase measurements (x_2τ − 2x_τ +x), the modified Allan variance (MVAR) is three times the variances of the phase measurements.

Further developments have produced improved estimation methods for the same stability measure, the variance/deviation of frequency, but these are known by separate names such as theHadamard variance,modified Hadamard variance, thetotal variance,modified total variance and theTheo variance. These distinguish themselves in better use of statistics for improved confidence bounds or ability to handle linear frequency drift.

Statistical estimators will calculate an estimated value on the sample series used. The estimates may deviate from the true value and the range of values which for some probability will contain the true value is referred to as theconfidence interval. The confidence interval depends on the number of observations in the sample series, the dominant noise type, and the estimator being used. The width is also dependent on the statistical certainty for which the confidence interval values forms a bounded range, thus the statistical certainty that the true value is within that range of values. For variable-τ estimators, theτ₀ multiplen is also a variable.

Theconfidence interval can be established usingchi-squared distribution with df degrees of freedom by using thedistribution of the sample variance:^[5][9]

${\displaystyle {\chi }^2=\frac{\text{df}{s}^2}{{\sigma }^2},}$

wheres² is the sample variance of our estimate,σ² is the true variance value, df is the degrees of freedom for the estimator, andχ² is calculated based on the inverse cumulative density distribution of aχ² with df degrees of freedom. For a 90% probability, covering the range from the 5% to the 95% range on the probability curve, the upper and lower limits can be found using the inequality

${\displaystyle {\chi }^2(0.05)\le \frac{\text{df}{s}^2}{{\sigma }^2}\le {\chi }^2(0.95),}$

which after rearrangement for the true variance becomes

${\displaystyle \frac{\text{df}{s}^2}{{\chi }^2(0.95)}\le {\sigma }^2\le \frac{\text{df}{s}^2}{{\chi }^2(0.05)}.}$

Thedegrees of freedom represents the number of free variables capable of contributing to the estimate. Depending on the estimator and noise type, the effective degrees of freedom varies. Estimator formulas depending onN (number of total sample points) andn (integer multiple ofτ₀) has been found empirically:^[9]

Allan variance degrees of freedom

Noise type	degrees of freedom
white phase modulation (WPM)	${\displaystyle \text{df}\cong \frac{(N+1)(N-2n)}{2(N-n)}}$
flicker phase modulation (FPM)	${\displaystyle \text{df}\cong \exp \left[{\left(\ln \frac{N-1}{2n}\ln \frac{(2n+1)(N-1)}{4}\right)}^{-1/2}\right]}$
white frequency modulation (WFM)	${\displaystyle \text{df}\cong \left[\frac{3(N-1)}{2n}-\frac{2(N-2)}{N}\right]\frac{4n^2}{4n^2+5}}$
flicker frequency modulation (FFM)
random-walk frequency modulation (RWFM)	${\displaystyle \text{df}\cong \frac{N-2}{n}\frac{(N-1{)}^2-3n(N-1)+4n^2}{(N-3{)}^2}}$

The Allan variance will treat variouspower-law noise types differently, conveniently allowing them to be identified and their strength estimated. As a convention, the measurement system width (high corner frequency) is denotedf_H.

Allan variance power-law response

Power-law noise type	Phase noise slope	Frequency noise slope	Power coefficient	Phase noise ${\displaystyle S_x(f)}$	Allan variance ${\displaystyle {\sigma }_y^2(\tau )}$	Allan deviation ${\displaystyle {\sigma }_y(\tau )}$
white phase modulation (WPM)	${\displaystyle f^0=1}$	${\displaystyle f^2}$	${\displaystyle h_2}$	${\displaystyle \frac{1}{(2\pi {)}^2}{h}_2}$	${\displaystyle \frac{3f_H}{4{\pi }^2{\tau }^2}{h}_2}$	${\displaystyle \frac{\sqrt{3f_H}}{2\pi \tau }\sqrt{h_2}}$
flicker phase modulation (FPM)	${\displaystyle f^{-1}}$	${\displaystyle f^1=f}$	${\displaystyle h_1}$	${\displaystyle \frac{1}{(2\pi {)}^{2}f}{h}_1}$	${\displaystyle \frac{3[\gamma +\ln (2\pi f_H\tau )]-\ln 2}{4{\pi }^2{\tau }^2}{h}_1}$	${\displaystyle \frac{\sqrt{3[\gamma +\ln (2\pi f_H\tau )]-\ln 2}}{2\pi \tau }\sqrt{h_1}}$
white frequency modulation (WFM)	${\displaystyle f^{-2}}$	${\displaystyle f^0=1}$	${\displaystyle h_0}$	${\displaystyle \frac{1}{(2\pi {)}^2{f}^2}{h}_0}$	${\displaystyle \frac{1}{2\tau }{h}_0}$	${\displaystyle \frac{1}{\sqrt{2\tau }}\sqrt{h_0}}$
flicker frequency modulation (FFM)	${\displaystyle f^{-3}}$	${\displaystyle f^{-1}}$	${\displaystyle h_{-1}}$	${\displaystyle \frac{1}{(2\pi {)}^2{f}^3}{h}_{-1}}$	${\displaystyle 2\ln (2)h_{-1}}$	${\displaystyle \sqrt{2\ln (2)}\sqrt{h_{-1}}}$
random walk frequency modulation (RWFM)	${\displaystyle f^{-4}}$	${\displaystyle f^{-2}}$	${\displaystyle h_{-2}}$	${\displaystyle \frac{1}{(2\pi {)}^2{f}^4}{h}_{-2}}$	${\displaystyle \frac{2{\pi }^2\tau }{3}{h}_{-2}}$	${\displaystyle \frac{\pi \sqrt{2\tau }}{\sqrt{3}}\sqrt{h_{-2}}}$

As found in^[12][13] and in modern forms.^[14][15]

The Allan variance is unable to distinguish between WPM and FPM, but is able to resolve the other power-law noise types. In order to distinguish WPM and FPM, themodified Allan variance needs to be employed.

The above formulas assume that

${\displaystyle \tau \gg \frac{1}{2\pi f_H},}$

and thus that thebandwidth of the observation time is much lower than the instruments bandwidth. When this condition is not met, all noise forms depend on the instrument's bandwidth.

The detailed mapping of a phase modulation of the form

${\displaystyle S_x(f)=\frac{1}{4{\pi }^2}{h}_{\alpha }{f}^{\alpha -2}=\frac{1}{4{\pi }^2}{h}_{\alpha }{f}^{\beta },}$

where

${\displaystyle \beta \equiv \alpha -2,}$

or frequency modulation of the form

${\displaystyle S_y(f)=h_{\alpha }{f}^{\alpha }}$

into the Allan variance of the form

${\displaystyle {\sigma }_y^2(\tau )=K_{\alpha }{h}_{\alpha }{\tau }^{\mu }}$

can be significantly simplified by providing a mapping betweenα andμ. A mapping betweenα andK_α is also presented for convenience:^[5]

Allan varianceα–μ mapping

α	β	μ	Kα
−2	−4	1	${\displaystyle \frac{2{\pi }^2}{3}}$
−1	−3	0	${\displaystyle 2\ln 2}$
0	−2	−1	${\displaystyle \frac{1}{2}}$
1	−1	−2	${\displaystyle \frac{3[\gamma +\ln (2\pi f_H\tau )]-\ln 2}{4{\pi }^2}}$
2	0	−2	${\displaystyle \frac{3f_H}{4{\pi }^2}}$

A signal with spectral phase noise${\displaystyle S_{\phi }}$ with units rad²/Hz can be converted to Allan Variance by^[15]

${\displaystyle {\sigma }_y^2(\tau )=\frac{2}{{\nu }_0^2}{\int }_0^{f_b}{S}_{\phi }(f)\frac{{\sin }^4(\pi \tau f)}{(\pi \tau {)}^2}df.}$

While Allan variance is intended to be used to distinguish noise forms, it will depend on some but not all linear responses to time. They are given in the table:

Allan variance linear response

Linear effect	time response	frequency response	Allan variance	Allan deviation
phase offset	${\displaystyle x_0}$	${\displaystyle 0}$	${\displaystyle 0}$	${\displaystyle 0}$
frequency offset	${\displaystyle y_{0}t}$	${\displaystyle y_0}$	${\displaystyle 0}$	${\displaystyle 0}$
linear drift	${\displaystyle \frac{D{t}^2}{2}}$	${\displaystyle Dt}$	${\displaystyle \frac{D^2{\tau }^2}{2}}$	${\displaystyle \frac{D\tau }{\sqrt{2}}}$

Thus, linear drift will contribute to output result. When measuring a real system, the linear drift or other drift mechanism may need to be estimated and removed from the time-series prior to calculating the Allan variance.^[14]

In analyzing the properties of Allan variance and friends, it has proven useful to consider the filter properties on the normalize frequency. Starting with the definition for Allan variance for

${\displaystyle {\sigma }_y^2(\tau )=\frac{1}{2}⟨{\left({\overline{y}}_{i+1}-{\overline{y}}_i\right)}^2⟩,}$

where

${\displaystyle {\overline{y}}_i=\frac{1}{\tau }{\int }_0^{\tau }y(i\tau +t)dt.}$

Replacing the time series of${\displaystyle y_i}$ with the Fourier-transformed variant${\displaystyle S_y(f)}$ the Allan variance can be expressed in thefrequency domain as

${\displaystyle {\sigma }_y^2(\tau )={\int }_0^{\mathrm{\infty }}{S}_y(f)\frac{2{\sin }^4\pi \tau f}{(\pi \tau f{)}^2}df.}$

Thus the transfer function for Allan variance is

${\displaystyle {\left|H_A(f)\right|}^2=\frac{2{\sin }^4\pi \tau f}{(\pi \tau f{)}^2}.}$

TheM-sample variance, and the defined special case Allan variance, will experiencesystematic bias depending on different number of samplesM and different relationship betweenT andτ. In order to address these biases the bias-functionsB₁ andB₂ has been defined^[16] and allows conversion between differentM andT values.

These bias functions are not sufficient for handling the bias resulting from concatenatingM samples to theMτ₀ observation time over theMT₀ with the dead-time distributed among theM measurement blocks rather than at the end of the measurement. This rendered the need for theB₃ bias.^[17]

The bias functions are evaluated for a particular μ value, so the α–μ mapping needs to be done for the dominant noise form as found usingnoise identification. Alternatively,^[4][16] the μ value of the dominant noise form may be inferred from the measurements using the bias functions.

TheB₁ bias function relates theM-sample variance with the 2-sample variance, keeping the time between measurementsT and time for each measurementsτ constant. It is defined^[16] as

${\displaystyle B_1(N,r,\mu )=\frac{⟨{\sigma }_y^2(N,T,\tau )⟩}{⟨{\sigma }_y^2(2,T,\tau )⟩},}$

where

${\displaystyle r=\frac{T}{\tau }.}$

The bias function becomes after analysis

${\displaystyle B_1(N,r,\mu )=\frac{1+\sum _{n=1}^{N-1}\frac{N-n}{N(N-1)}\left[2(rn{)}^{\mu +2}-(rn+1{)}^{\mu +2}-|rn-1{|}^{\mu +2}\right]}{1+\frac{1}{2}\left[2r^{\mu +2}-(r+1{)}^{\mu +2}-|r-1{|}^{\mu +2}\right]}.}$

TheB₂ bias function relates the 2-sample variance for sample timeT with the 2-sample variance (Allan variance), keeping the number of samplesN = 2 and the observation timeτ constant. It is defined^[16] as

${\displaystyle B_2(r,\mu )=\frac{⟨{\sigma }_y^2(2,T,\tau )⟩}{⟨{\sigma }_y^2(2,\tau ,\tau )⟩},}$

where

${\displaystyle r=\frac{T}{\tau }.}$

The bias function becomes after analysis

${\displaystyle B_2(r,\mu )=\frac{1+\frac{1}{2}\left[2r^{\mu +2}-(r+1{)}^{\mu +2}-|r-1{|}^{\mu +2}\right]}{2\left(1-2^{\mu }\right)}.}$

TheB₃ bias function relates the 2-sample variance for sample timeMT₀ and observation timeMτ₀ with the 2-sample variance (Allan variance) and is defined^[17] as

${\displaystyle B_3(N,M,r,\mu )=\frac{⟨{\sigma }_y^2(N,M,T,\tau )⟩}{⟨{\sigma }_y^2(N,T,\tau )⟩},}$

where

${\displaystyle T=M{T}_0,}$

${\displaystyle \tau =M{\tau }_0.}$

TheB₃ bias function is useful to adjust non-overlapping and overlapping variableτ estimator values based on dead-time measurements of observation timeτ₀ and time between observationsT₀ to normal dead-time estimates.

The bias function becomes after analysis (for theN = 2 case)

${\displaystyle B_3(2,M,r,\mu )=\frac{2M+MF(Mr)-\sum _{n=1}^{M-1}(M-n)\left[2F(nr)-F((M+n)r)+F((M-n)r)\right]}{M^{\mu +2}[F(r)+2]},}$

where

${\displaystyle F(A)=2A^{\mu +2}-(A+1{)}^{\mu +2}-|A-1{|}^{\mu +2}.}$

While formally not formulated, it has been indirectly inferred as a consequence of theα–μ mapping. When comparing two Allan variance measure for differentτ, assuming same dominant noise in the form of same μ coefficient, a bias can be defined as

${\displaystyle B_{\tau }({\tau }_1,{\tau }_2,\mu )=\frac{⟨{\sigma }_y^2(2,{\tau }_2,{\tau }_2)⟩}{⟨{\sigma }_y^2(2,{\tau }_1,{\tau }_1)⟩}.}$

The bias function becomes after analysis

${\displaystyle B_{\tau }({\tau }_1,{\tau }_2,\mu )={\left(\frac{{\tau }_2}{{\tau }_1}\right)}^{\mu }.}$

In order to convert from one set of measurements to another theB₁,B₂ and τ bias functions can be assembled. First theB₁ function converts the (N₁, T₁, τ₁) value into (2, T₁, τ₁), from which theB₂ function converts into a (2, τ₁, τ₁) value, thus the Allan variance atτ₁. The Allan variance measure can be converted using the τ bias function fromτ₁ toτ₂, from which then the (2, T₂, τ₂) usingB₂ and then finally usingB₁ into the (N₂, T₂, τ₂) variance. The complete conversion becomes

${\displaystyle ⟨{\sigma }_y^2(N_2,T_2,{\tau }_2)⟩={\left(\frac{{\tau }_2}{{\tau }_1}\right)}^{\mu }\left[\frac{B_1(N_2,r_2,\mu )B_2(r_2,\mu )}{B_1(N_1,r_1,\mu )B_2(r_1,\mu )}\right]⟨{\sigma }_y^2(N_1,T_1,{\tau }_1)⟩,}$

where

${\displaystyle r_1=\frac{T_1}{r_1},}$

${\displaystyle r_2=\frac{T_2}{r_2}.}$

Similarly, for concatenated measurements usingM sections, the logical extension becomes

${\displaystyle ⟨{\sigma }_y^2(N_2,M_2,T_2,{\tau }_2)⟩={\left(\frac{{\tau }_2}{{\tau }_1}\right)}^{\mu }\left[\frac{B_3(N_2,M_2,r_2,\mu )B_1(N_2,r_2,\mu )B_2(r_2,\mu )}{B_3(N_1,M_1,r_1,\mu )B_1(N_1,r_1,\mu )B_2(r_1,\mu )}\right]⟨{\sigma }_y^2(N_1,M_1,T_1,{\tau }_1)⟩.}$

When making measurements to calculate Allan variance or Allan deviation, a number of issues may cause the measurements to degenerate. Covered here are the effects specific to Allan variance, where results would be biased.

A measurement system is expected to have a bandwidth at or below that of theNyquist rate, as described within theShannon–Hartley theorem. As can be seen in the power-law noise formulas, the white and flicker noise modulations both depend on the upper corner frequency${\displaystyle f_H}$ (these systems are assumed to be low-pass filtered only). Considering the frequency filter property, it can be clearly seen that low-frequency noise has greater impact on the result. For relatively flat phase-modulation noise types (e.g. WPM and FPM), the filtering has relevance, whereas for noise types with greater slope the upper frequency limit becomes of less importance, assuming that the measurement system bandwidth is wide relative the${\displaystyle \tau }$ as given by

${\displaystyle \tau \gg \frac{1}{2\pi f_H}.}$

When this assumption is not met, the effective bandwidth${\displaystyle f_H}$ needs to be notated alongside the measurement. The interested should consult NBS TN394.^[12]

If, however, one adjust the bandwidth of the estimator by using integer multiples of the sample time${\displaystyle n{\tau }_0}$, then the system bandwidth impact can be reduced to insignificant levels. For telecommunication needs, such methods have been required in order to ensure comparability of measurements and allow some freedom for vendors to do different implementations. The ITU-T Rec. G.813^[18] for the TDEV measurement.

It can be recommended that the first${\displaystyle {\tau }_0}$ multiples be ignored, such that the majority of the detected noise is well within the passband of the measurement systems bandwidth.

Further developments on the Allan variance was performed to let the hardware bandwidth be reduced by software means. This development of a software bandwidth allowed addressing the remaining noise, and the method is now referred tomodified Allan variance. This bandwidth reduction technique should not be confused with the enhanced variant ofmodified Allan variance, which also changes a smoothing filter bandwidth.

Many measurement instruments of time and frequency have the stages of arming time, time-base time, processing time and may then re-trigger the arming. The arming time is from the time the arming is triggered to when the start event occurs on the start channel. The time-base then ensures that minimal amount of time goes prior to accepting an event on the stop channel as the stop event. The number of events and time elapsed between the start event and stop event is recorded and presented during the processing time. When the processing occurs (also known as thedwell time), the instrument is usually unable to do another measurement. After the processing has occurred, an instrument in continuous mode triggers the arm circuit again. The time between the stop event and the following start event becomesdead time, during which the signal is not being observed. Such dead time introduces systematic measurement biases, which needs to be compensated for in order to get proper results. For such measurement systems will the timeT denote the time between the adjacent start events (and thus measurements), while${\displaystyle \tau }$ denote the time-base length, i.e. the nominal length between the start and stop event of any measurement.

Dead-time effects on measurements have such an impact on the produced result that much study of the field have been done in order to quantify its properties properly. The introduction of zero-dead-time counters removed the need for this analysis. A zero-dead-time counter has the property that the stop event of one measurement is also being used as the start event of the following event. Such counters create a series of event and time timestamp pairs, one for each channel spaced by the time-base. Such measurements have also proved useful in order forms of time-series analysis.

Measurements being performed with dead time can be corrected using the bias functionB₁,B₂ andB₃. Thus, dead time as such is not prohibiting the access to the Allan variance, but it makes it more problematic. The dead time must be known, such that the time between samplesT can be established.

Studying the effect on theconfidence intervals that the lengthN of the sample series have and the effect of the variableτ parametern, the confidence intervals may become very large since theeffective degree of freedom may become small for some combination ofN andn for the dominant noise form (for thatτ).

The effect may be that the estimated value may be much smaller or much greater than the real value, which may lead to false conclusions of the result.

It is recommended that:

• The confidence interval be plotted along with the data, such that the reader of the plot knows of the statistical uncertainty of the values.
• The length of the sample sequence (i.e. the number of samplesN) must be kept as high as possible to ensure that confidence interval is small over theτ range of interest.
• Estimators providing better degrees of freedom values be used in replacement of the Allan variance estimators or as complementing them where they outperform the Allan variance estimators. Among those thetotal variance andTheo variance estimators should be considered.
• Theτ range as swept by theτ₀ multipliern is limited in the upper end relativeN, such that the reader of the plot may not be confused by highly unstable estimator values.

A large number of conversion constants, bias corrections and confidence intervals depends on the dominant noise type. For proper interpretation shall the dominant noise type for the particularτ of interest be identified through noise identification. Failing to identify the dominant noise type will produce biased values. Some of these biases may be of several order of magnitude, so it may be of large significance.

Systematic effects on the signal is only partly cancelled. Phase and frequency offset is cancelled, but linear drift or other high-degree forms of polynomial phase curves will not be cancelled and thus form a measurement limitation.Curve fitting and removal of systematic offset could be employed. Often removal of linear drift can be sufficient. Use of linear-drift estimators such as theHadamard variance could also be employed. A linear drift removal could be employed using a moment-based estimator.

Traditional instruments provided only the measurement of single events or event pairs. The introduction of the improved statistical tool of overlapping measurements by J. J. Snyder^[8] allowed much improved resolution in frequency readouts, breaking the traditional digits/time-base balance. While such methods is useful for their intended purpose, using such smoothed measurements for Allan variance calculations would give a false impression of high resolution,^[19][20][21] but for longerτ the effect is gradually removed, and the lower-τ region of the measurement has biased values. This bias is providing lower values than it should, so it is an overoptimistic (assuming that low numbers is what one wishes) bias, reducing the usability of the measurement rather than improving it. Such smart algorithms can usually be disabled or otherwise circumvented by using time-stamp mode, which is much preferred if available.

This sectiondoes notcite anysources. Please helpimprove this section byadding citations to reliable sources. Unsourced material may be challenged andremoved.(January 2018) (Learn how and when to remove this message)

While several approaches to measurement of Allan variance can be devised, a simple example may illustrate how measurements can be performed.

All measurements of Allan variance will in effect be the comparison of two different clocks. Consider a reference clock and a device under test (DUT), and both having a common nominal frequency of 10 MHz. A time-interval counter is being used to measure the time between the rising edge of the reference (channel A) and the rising edge of the device under test.

In order to provide evenly spaced measurements, the reference clock will be divided down to form the measurement rate, triggering the time-interval counter (ARM input). This rate can be 1 Hz (using the1 PPS output of a reference clock), but other rates like 10 Hz and 100 Hz can also be used. The speed of which the time-interval counter can complete the measurement, output the result and prepare itself for the next arm will limit the trigger frequency.

A computer is then useful to record the series of time differences being observed.

The recorded time-series require post-processing to unwrap the wrapped phase, such that a continuous phase error is being provided. If necessary, logging and measurement mistakes should also be fixed. Drift estimation and drift removal should be performed, the drift mechanism needs to be identified and understood for the sources. Drift limitations in measurements can be severe, so letting the oscillators become stabilized, by long enough time being powered on, is necessary.

The Allan variance can then be calculated using the estimators given, and for practical purposes the overlapping estimator should be used due to its superior use of data over the non-overlapping estimator. Other estimators such as total or Theo variance estimators could also be used if bias corrections is applied such that they provide Allan variance-compatible results.

To form the classical plots, the Allan deviation (square root of Allan variance) is plotted in log–log format against the observation interval τ.

The time-interval counter is typically an off-the-shelf counter commercially available. Limiting factors involve single-shot resolution, trigger jitter, speed of measurements and stability of reference clock. The computer collection and post-processing can be done using existing commercial or public-domain software. Highly advanced solutions exists, which will provide measurement and computation in one box.

The field of frequency stability has been studied for a long time. However, during the 1960s it was found that coherent definitions were lacking. ANASA-IEEE Symposium on Short-Term Stability in November 1964^[22] resulted in the special February 1966 issue of the IEEE Proceedings on Frequency Stability.

The NASA-IEEE Symposium brought together many fields and uses of short- and long-term stability, with papers from many different contributors. The articles and panel discussions concur on the existence of the frequency flicker noise and the wish to achieve a common definition for both short-term and long-term stability.

Important papers, including those of David Allan,^[4] James A. Barnes,^[23] L. S. Cutler and C. L. Searle^[2] and D. B. Leeson,^[3] appeared in the IEEE Proceedings on Frequency Stability and helped shape the field.

David Allan's article analyses the classicalM-sample variance of frequency, tackling the issue of dead-time between measurements along with an initial bias function.^[4] Although Allan's initial bias function assumes no dead-time, his formulas do include dead-time calculations. His article analyses the case of M frequency samples (called N in the article) and variance estimators. It provides the now standard α–μ mapping, clearly building on James Barnes' work^[23] in the same issue.

The 2-sample variance case is a special case of theM-sample variance, which produces an average of the frequency derivative. Allan implicitly uses the 2-sample variance as a base case, since for arbitrary chosenM, values may be transferred via the 2-sample variance to theM-sample variance. No preference was clearly stated for the 2-sample variance, even if the tools were provided. However, this article laid the foundation for using the 2-sample variance as a way of comparing otherM-sample variances.

James Barnes significantly extended the work on bias functions,^[16] introducing the modernB₁ andB₂ bias functions. It refers to theM-sample variance as "Allan variance", while referring to Allan's article "Statistics of Atomic Frequency Standards".^[4] With these modern bias functions, full conversion amongM-sample variance measures of variousM,T andτ values could be performed, by conversion through the 2-sample variance.

James Barnes and David Allan further extended the bias functions with theB₃ function^[17] to handle the concatenated samples estimator bias. This was necessary to handle the new use of concatenated sample observations with dead-time in between.

In 1970, the IEEE Technical Committee on Frequency and Time, within the IEEE Group on Instrumentation & Measurements, provided a summary of the field, published as NBS Technical Notice 394.^[12] This paper was first in a line of more educational and practical papers helping fellow engineers grasp the field. This paper recommended the 2-sample variance withT =τ, referring to it asAllan variance (now without the quotes). The choice of such parametrisation allows good handling of some noise forms and getting comparable measurements; it is essentially the least common denominator with the aid of the bias functionsB₁ andB₂.

J. J. Snyder proposed an improved method for frequency or variance estimation, using sample statistics for frequency counters.^[8] To get more effective degrees of freedom out of the available dataset, the trick is to use overlapping observation periods. This provides a√n improvement, and was incorporated in theoverlapping Allan variance estimator.^[9] Variable-τ software processing was also incorporated.^[9] This development improved the classical Allan variance estimators, likewise providing a direct inspiration for the work onmodified Allan variance.

Howe, Allan and Barnes presented the analysis of confidence intervals, degrees of freedom, and the established estimators.^[9]

The field of time and frequency and its use of Allan variance,Allan deviation and friends is a field involving many aspects, for which both understanding of concepts and practical measurements and post-processing requires care and understanding. Thus, there is a realm of educational material stretching about 40 years available. Since these reflect the developments in the research of their time, they focus on teaching different aspect over time, in which case a survey of available resources may be a suitable way of finding the right resource.

The first meaningful summary is the NBS Technical Note 394 "Characterization of Frequency Stability".^[12] This is the product of the Technical Committee on Frequency and Time of the IEEE Group on Instrumentation & Measurement. It gives the first overview of the field, stating the problems, defining the basic supporting definitions and getting into Allan variance, the bias functionsB₁ andB₂, the conversion of time-domain measures. This is useful, as it is among the first references to tabulate the Allan variance for the five basic noise types.

A classical reference is the NBS Monograph 140^[24] from 1974, which in chapter 8 has "Statistics of Time and Frequency Data Analysis".^[25] This is the extended variant of NBS Technical Note 394 and adds essentially in measurement techniques and practical processing of values.

An important addition will be theProperties of signal sources and measurement methods.^[9] It covers the effective use of data, confidence intervals, effective degree of freedom, likewise introducing the overlapping Allan variance estimator. It is a highly recommended reading for those topics.

The IEEE standard 1139Standard definitions of Physical Quantities for Fundamental Frequency and Time Metrology^[5] is beyond that of a standard a comprehensive reference and educational resource.

A modern book aimed towards telecommunication is Stefano Bregni "Synchronisation of Digital Telecommunication Networks".^[14] This summarises not only the field, but also much of his research in the field up to that point. It aims to include both classical measures and telecommunication-specific measures such as MTIE. It is a handy companion when looking at measurements related to telecommunication standards.

The NIST Special Publication 1065 "Handbook of Frequency Stability Analysis" of W. J. Riley^[15] is a recommended reading for anyone wanting to pursue the field. It is rich of references and also covers a wide range of measures, biases and related functions that a modern analyst should have available. Further it describes the overall processing needed for a modern tool.

Allan variance is used as a measure of frequency stability in a variety of precision oscillators, such ascrystal oscillators,atomic clocks and frequency-stabilizedlasers over a period of a second or more. Short-term stability (under a second) is typically expressed asphase noise. The Allan variance is also used to characterize the bias stability ofgyroscopes, includingfiber optic gyroscopes,hemispherical resonator gyroscopes andMEMS gyroscopes andaccelerometers.^[26][27]

In 2016,IEEE-UFFC is going to be publishing a "Special Issue to celebrate the 50th anniversary of the Allan Variance (1966–2016)".^[28] A guest editor for that issue will be David's former colleague atNIST, Judah Levine, who is the most recent recipient of theI. I. Rabi Award.

• UFFC Frequency Control Teaching Resources
• NIST Publication search tool
• David W. Allan's Allan Variance Overview
• David W. Allan's official web site
• JPL Publications – Noise Analysis and Statistics
• William Riley publications
• Stable32, Software for Frequency Stability Analysis, by William Riley
• Stefano Bregni publications
• Enrico Rubiola publications
• Allanvar: R package for sensor error characterization using the Allan Variance
• Alavar windows software with reporting tools; Freeware
• AllanTools open-source python library for Allan variance
• MATLAB AVAR open-source MATLAB application

• Clocks
• Signal processing metrics
• Measurement

• Webarchive template wayback links
• Articles with short description
• Short description matches Wikidata
• Use dmy dates from June 2023
• Articles needing additional references from January 2018
• All articles needing additional references