Вариация Аллана

Часы легче всего проверить, сравнив их с гораздо более точными эталонными часами. В течение интервала времени τ , измеряемого эталонными часами, тестируемые часы продвигаются вперед на τy , где y - средняя (относительная) тактовая частота за этот интервал. Если мы измеряем два последовательных интервала, как показано, мы можем получить значение ( y - y ') ² - меньшее значение указывает на более стабильные и точные часы. Если мы повторим эту процедуру много раз, среднее значение ( y - y ') ² будет равно удвоенной дисперсии Аллана (или квадрату отклонения Аллана) для времени наблюдения.τ .

Дисперсия Аллана ( AVAR ), также известная как два образца дисперсия , является мерой стабильности частоты в часах , генераторах и усилителях . Он назван в честь Дэвида У. Аллана и математически выражается как . Отклонение Аллана ( ADEV ), также известное как сигма-тау , есть квадратный корень из дисперсии Аллана, .${\ Displaystyle \ sigma _ {у} ^ {2} (\ тау)}$${\ Displaystyle \ sigma _ {у} (\ тау)}$

М-выборочная дисперсия является мерой стабильности частоты с использованием М выборок, времени T между измерениями и времени наблюдения . M -выборочная дисперсия выражается как${\ Displaystyle \ тау}$

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

Дисперсия Аллана предназначена для оценки стабильности из-за шумовых процессов, а не из-за систематических ошибок или недостатков, таких как дрейф частоты или температурные эффекты. Дисперсия Аллана и девиация Аллана описывают стабильность частоты. См. Также раздел « Интерпретация стоимости» ниже.

Существуют также различные адаптации или изменения дисперсии Аллана, в частности, модифицированная дисперсия Аллана MAVAR или MVAR, общая дисперсия и дисперсия Адамара . Также существуют варианты временной стабильности, такие как временное отклонение (TDEV) или временное отклонение (TVAR). Дисперсия Аллана и ее варианты доказали свою полезность вне рамок хронометража и представляют собой набор улучшенных статистических инструментов, которые можно использовать всякий раз, когда шумовые процессы не являются безусловно стабильными, поэтому существует производная.

Общая дисперсия M- выборки остается важной, поскольку она допускает мертвое время в измерениях, а функции смещения позволяют преобразовывать в значения дисперсии Аллана. Тем не менее, для большинства приложений наибольший интерес представляет особый случай 2-выборки или «вариации Аллана» с .${\ Displaystyle Т = \ тау}$

Пример графика отклонения Аллана часов. При очень коротком времени наблюдения τ отклонение Аллана велико из-за шума. При больших τ он уменьшается, так как шум усредняется. При еще большем значении τ отклонение Аллана снова начинает увеличиваться, предполагая, что тактовая частота постепенно дрейфует из-за изменений температуры, старения компонентов или других подобных факторов. Планки погрешностей увеличиваются с увеличением τ просто потому, что получение большого количества точек данных для большого τ занимает много времени .

Фон [ править ]

При исследовании стабильности кварцевых генераторов и атомных часов было обнаружено, что они не имеют фазового шума, состоящего только из белого шума , но также из шума частоты мерцания . Эти формы шума становятся проблемой для традиционных статистических инструментов, таких как стандартное отклонение , поскольку оценка не сходится. Таким образом, шум считается расходящимся. Первые попытки анализа устойчивости включали как теоретический анализ, так и практические измерения. ^{[1] [2]}

Важным побочным следствием наличия этих типов шума было то, что, поскольку различные методы измерений не согласовывались друг с другом, ключевой аспект повторяемости измерения не мог быть достигнут. Это ограничивает возможность сравнения источников и составления значимых спецификаций, требуемых от поставщиков. Практически все формы научного и коммерческого использования тогда ограничивались специальными измерениями, которые, как мы надеемся, уловят потребность в этом приложении.

Чтобы решить эти проблемы, Дэвид Аллан ввел дисперсию M- выборки и (косвенно) дисперсию двух выборок. ^[3] Хотя двухвыборочная дисперсия не позволяла полностью различать все типы шума, она предоставила средства для значимого разделения многих форм шума для временных рядов измерений фазы или частоты между двумя или более генераторами. Аллан предоставил метод преобразования любой дисперсии M- выборки в любую дисперсию N- выборки с помощью общей дисперсии для двух выборок, что сделало все дисперсии M- выборки сопоставимыми. Механизм преобразования также доказал, что дисперсия M -выборки не сходится при больших M, что делает их менее полезными. Позже IEEE определил двухвыборочную дисперсию как предпочтительную меру. ^[4]

Первоначальное беспокойство было связано с приборами для измерения времени и частоты, которые имели мертвое время между измерениями. Такая серия измерений не обеспечивала непрерывного наблюдения за сигналом и, таким образом, вносила систематическую ошибку в измерение. Оценка этих предубеждений была проведена с большой осторожностью. Введение счетчиков с нулевым мертвым временем устранило необходимость, но инструменты анализа смещения оказались полезными.

Другой ранний аспект беспокойства был связан с тем, как полоса пропускания измерительного прибора будет влиять на измерение, так что это необходимо было отметить. Позже было обнаружено, что при алгоритмическом изменении наблюдения будут затронуты только низкие значения, в то время как более высокие значения не будут затронуты. Изменение осуществляется путем использования целого числа, кратного оси времени измерения :${\ Displaystyle \ тау}$${\ Displaystyle \ тау}$${\ Displaystyle \ тау}$${\ displaystyle n}$ ${\ displaystyle \ tau _ {0}}$

${\ Displaystyle \ тау = п \ тау _ {0}.}$

Физика кварцевых осцилляторов была проанализирована Д. Б. Лисоном ^[2], и теперь результат называется уравнением Лисона . Обратная связь в генераторе сделает белый шум и фликкер - шум усилителя обратной связи и кристалл становится степенными шумами от белого шума частоты и частоты мерцания шума соответственно. Эти формы шума приводят к тому, что стандартная дисперсия${\ displaystyle f ^ {- 2}}$${\ displaystyle f ^ {- 3}}$Оценщик не сходится при обработке выборок с временной ошибкой. Эта механика осцилляторов с обратной связью была неизвестна, когда началась работа над стабильностью осцилляторов, но была представлена ​​Лисоном в то же время, когда Дэвид У. Аллан предоставил набор статистических инструментов . Более подробное описание эффекта Лисона см. В современной литературе по фазовому шуму. ^[5]

Интерпретация значения [ править ]

Дисперсия Аллана определяется как половина среднего по времени квадратов разностей между последовательными показаниями отклонения частоты, выбранными за период выборки. Дисперсия Аллана зависит от периода времени, используемого между выборками, поэтому она является функцией периода выборки, обычно обозначаемой как τ , аналогично измеряемому распределению и отображается в виде графика, а не одного числа. Низкая дисперсия Аллана - характеристика часов с хорошей стабильностью в течение измеряемого периода.

Отклонение Аллана широко используется для построения графиков (обычно в формате log – log ) и представления чисел. Это предпочтительнее, так как дает относительную стабильность амплитуды, позволяя легко сравнивать с другими источниками ошибок.

Отклонение Аллана 1,3 × 10 ⁻⁹ при времени наблюдения 1 с (т. Е. Τ = 1 с) следует интерпретировать как нестабильность частоты между двумя наблюдениями с интервалом в 1 секунду с относительным среднеквадратичным значением (RMS) 1,3 × 10 ⁻⁹ . Для тактовой частоты 10 МГц это будет эквивалентно движению RMS на 13 МГц. Если необходима фазовая стабильность генератора, следует проконсультироваться и использовать варианты отклонения по времени .

Можно преобразовать дисперсию Аллана и другие дисперсии во временной области в измерения времени (фазы) и стабильности частоты в частотной области. ^[6]

Определения [ править ]

M- выборочная дисперсия [ править ]

-Sample дисперсия определяется ^[3] (здесь в модернизированном виде обозначений) , как${\ displaystyle M}$

${\displaystyle \sigma _{y}^{2}(M,T,\tau )={\frac {1}{M-1}}\left\{\sum _{i=0}^{M-1}\left[{\frac {x(iT+\tau )-x(iT)}{\tau }}\right]^{2}-{\frac {1}{M}}\left[\sum _{i=0}^{M-1}{\frac {x(iT+\tau )-x(iT)}{\tau }}\right]^{2}\right\},}$

где - показание часов (в секундах), измеренное во времени , или с временным рядом средней дробной частоты${\displaystyle x(t)}$${\displaystyle t}$

${\displaystyle \sigma _{y}^{2}(M,T,\tau )={\frac {1}{M-1}}\left\{\sum _{i=0}^{M-1}{\bar {y}}_{i}^{2}-{\frac {1}{M}}\left[\sum _{i=0}^{M-1}{\bar {y}}_{i}\right]^{2}\right\},}$

где - количество частотных выборок, используемых в дисперсии, - это время между каждой частотной выборкой, и - это продолжительность каждой оценки частоты.${\displaystyle M}$${\displaystyle T}$${\displaystyle \tau }$

Важным аспектом является то, что модель дисперсии выборки может включать мертвое время, позволяя времени отличаться от этого .${\displaystyle M}$${\displaystyle T}$${\displaystyle \tau }$

Вариант Аллана [ править ]

Дисперсия Аллана определяется как

${\displaystyle \sigma _{y}^{2}(\tau )=\langle \sigma _{y}^{2}(2,\tau ,\tau )\rangle ,}$

где обозначает оператор математического ожидания. Это удобно выразить как${\displaystyle \langle \dots \rangle }$

${\displaystyle \sigma _{y}^{2}(\tau )={\frac {1}{2}}\langle ({\bar {y}}_{n+1}-{\bar {y}}_{n})^{2}\rangle ={\frac {1}{2\tau ^{2}}}\langle (x_{n+2}-2x_{n+1}+x_{n})^{2}\rangle ,}$

где - период наблюдения, - n- е среднее значение дробной частоты за время наблюдения .${\displaystyle \tau }$${\displaystyle {\bar {y}}_{n}}$${\displaystyle \tau }$

Образцы отбираются без мертвого времени между ними, что достигается

${\displaystyle T=\tau .}$

Отклонение Аллана [ править ]

Так же, как стандартное отклонение и дисперсия , отклонение Аллана определяется как квадратный корень из дисперсии Аллана:

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

Дополнительные определения [ править ]

Модель осциллятора [ править ]

Предполагается, что анализируемый осциллятор следует базовой модели

${\displaystyle V(t)=V_{0}\sin(\Phi (t)).}$

Предполагается, что осциллятор имеет номинальную частоту , выраженную в циклах в секунду (единица СИ: герцы ). Номинальная угловая частота (в радианах в секунду) определяется выражением${\displaystyle \nu _{\text{n}}}$ ${\displaystyle \omega _{\text{n}}}$

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

Полная фаза может быть разделена на идеально циклическую составляющую вместе с флуктуирующей составляющей :${\displaystyle \omega _{\text{n}}t}$${\displaystyle \varphi (t)}$

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

Ошибка времени [ править ]

Функция временной ошибки x ( t ) - это разница между ожидаемым номинальным временем и фактическим нормальным временем:

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

Для измеренных значений ряд ошибок времени TE ( t ) определяется из функции эталонного времени T _REF ( t ) как

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

Функция частоты [ править ]

Функция частоты - это частота во времени, определяемая как${\displaystyle \nu (t)}$

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

Дробная частота [ править ]

Дробная частота y ( t ) - это нормализованная разница между частотой и номинальной частотой :${\displaystyle \nu (t)}$${\displaystyle \nu _{\text{n}}}$

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

Средняя дробная частота [ править ]

Средняя дробная частота определяется как

${\displaystyle {\bar {y}}(t,\tau )={\frac {1}{\tau }}\int \limits _{0}^{\tau }y(t+t_{v})\,dt_{v},}$

где среднее значение берется за время наблюдения τ , y ( t ) - это относительная частотная ошибка в момент времени t , а τ - время наблюдения.

Поскольку y ( t ) является производной от x ( t ), мы можем без ограничения общности переписать ее как

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

Оценщики [ править ]

Это определение основано на статистическом ожидаемом значении , интегрируемом за бесконечное время. Реальная ситуация не позволяет использовать такие временные ряды, и в этом случае необходимо использовать статистическую оценку . Будет представлен и обсужден ряд различных оценок.

Соглашения [ править ]

• Количество выборок частоты в серии дробных частот обозначаются M .
• Количество временных выборок ошибок в серии временных ошибок обозначаются N .

Связь между количеством отсчетов с дробной частотой и сериями временных ошибок фиксируется соотношением

${\displaystyle N=M+1.}$

• Для серии отсчетов временных ошибок x _i обозначает i-й отсчет функции непрерывного времени x ( t ), как задано

${\displaystyle x_{i}=x(iT),}$

где T - время между измерениями. Для дисперсии Аллана используемое время имеет значение T, равное времени наблюдения τ .

Последовательность отсчетов с временной ошибкой позволяет обозначать N количество отсчетов ( x ₀ ... x _{N -1} ) в серии. Традиционная конвенция использует индекс 1 через N .

• Для серии выборок средней дробной частоты , обозначает i- ю выборку средней непрерывной функции дробной частоты y ( t ), как задано${\displaystyle {\bar {y}}_{i}}$

${\displaystyle {\bar {y}}_{i}={\bar {y}}(Ti,\tau ),}$

который дает

${\displaystyle {\bar {y}}_{i}={\frac {1}{\tau }}\int \limits _{0}^{\tau }y(iT+t_{v})\,dt_{v}={\frac {x(iT+\tau )-x(iT)}{\tau }}.}$

Для предположения о вариации Аллана, что T равно τ, становится

${\displaystyle {\bar {y}}_{i}={\frac {x_{i+1}-x_{i}}{\tau }}.}$

Средняя дробно-частота выборки серии позволяет М обозначает число выборок ( ) в серии. Традиционная конвенция использует индекс 1 через M .${\displaystyle {\bar {y}}_{0}\ldots {\bar {y}}_{M-1}}$

Для краткости средняя дробная частота часто пишется без средней полосы над ней. Однако это формально неверно, поскольку дробная частота и средняя дробная частота - две разные функции. Измерительный прибор, способный производить оценки частоты без мертвого времени, фактически предоставляет временные ряды со средней частотой, которые нужно только преобразовать в среднюю дробную частоту и затем использовать напрямую.

• Кроме того, принято использовать τ для обозначения номинальной разницы во времени между соседними фазовыми или частотными отсчетами. Временной ряд, взятый для одной разницы во времени τ _0, может использоваться для генерации дисперсии Аллана для любого τ, являющегося целым кратным τ ₀ , и в этом случае используется τ = nτ ₀ , а n становится переменной для оценки.
• Время между измерениями обозначается буквой T , которая представляет собой сумму времени наблюдения τ и мертвого времени.

Фиксированные оценки τ [ править ]

Первым простым оценщиком было бы прямое преобразование определения в

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

или для временного ряда:

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

Эти формулы, однако, обеспечивают расчет только для случая τ = τ ₀ . Чтобы рассчитать другое значение τ , необходимо предоставить новый временной ряд.

Неперекрывающиеся оценки переменных τ [ править ]

Взяв временной ряд и пропустив n  - 1 выборку, возникнет новый (более короткий) временной ряд с τ ₀ как время между соседними выборками, для которого дисперсия Аллана может быть рассчитана с помощью простых оценок. Их можно было бы изменить, чтобы ввести новую переменную n , чтобы не нужно было создавать новые временные ряды, а вместо этого можно было бы повторно использовать исходные временные ряды для различных значений n . Оценщики становятся

${\displaystyle \sigma _{y}^{2}(n\tau _{0},M)={\text{AVAR}}(n\tau _{0},M)={\frac {1}{2{\frac {M-1}{n}}}}\sum _{i=0}^{{\frac {M-1}{n}}-1}({\bar {y}}_{ni+n}-{\bar {y}}_{ni})^{2}}$

с ,${\displaystyle n\leq M-1}$

и для временного ряда:

${\displaystyle \sigma _{y}^{2}(n\tau _{0},N)={\text{AVAR}}(n\tau _{0},N)={\frac {1}{2n^{2}\tau _{0}^{2}({\frac {N-1}{n}}-1)}}\sum _{i=0}^{{\frac {N-1}{n}}-2}(x_{ni+2n}-2x_{ni+n}+x_{ni})^{2}}$

с .${\displaystyle n\leq {\frac {N-1}{2}}}$

У этих оценщиков есть существенный недостаток, заключающийся в том, что они отбрасывают значительный объем выборочных данных, поскольку используется только 1 / n из доступных выборок.

Перекрывающиеся оценки переменных τ [ править ]

Методика, представленная Дж. Дж. Снайдером ^[7], предоставила улучшенный инструмент, поскольку измерения были перекрыты в n перекрывающихся сериях из исходной серии. Перекрывающаяся оценка дисперсии Аллана была введена Хоу, Алланом и Барнсом. ^[8] Можно показать, что это эквивалентно усреднению временных или нормированных частотных отсчетов в блоках по n отсчетов до обработки. Результирующий предиктор становится

${\displaystyle \sigma _{y}^{2}(n\tau _{0},M)={\text{AVAR}}(n\tau _{0},M)={\frac {1}{2n^{2}(M-2n+1)}}\sum _{j=0}^{M-2n}\left(\sum _{i=j}^{j+n-1}y_{i+n}-y_{i}\right)^{2}={\frac {1}{2(M-2n+1)}}\sum _{j=0}^{M-2n}\left({\bar {y}}_{j+n}-{\bar {y}}_{j}\right)^{2},}$

или для временного ряда:

${\displaystyle \sigma _{y}^{2}(n\tau _{0},N)={\text{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}.}$

Перекрывающиеся оценщики имеют гораздо лучшие характеристики по сравнению с непересекающимися оценщиками, поскольку n возрастает, а временной ряд имеет умеренную длину. Перекрывающиеся оценки были приняты в качестве предпочтительных оценок дисперсии Аллана в стандартах IEEE, ^[4] ITU-T ^[9] и ETSI ^[10] для сопоставимых измерений, например, необходимых для аттестации электросвязи.

Измененная дисперсия Аллана [ править ]

Чтобы решить проблему невозможности отделить модуляцию фазы белого от модуляции фазы мерцания с использованием традиционных средств оценки дисперсии Аллана, алгоритмическая фильтрация уменьшает полосу пропускания на n . Эта фильтрация обеспечивает модификацию определения и оценок и теперь идентифицируется как отдельный класс дисперсии, называемый измененной дисперсией Аллана . Модифицированная мера дисперсии Аллана является мерой стабильности частоты, так же как и дисперсия Аллана.

Оценщики временной стабильности [ править ]

Статистическая мера стабильности во времени (σ _x ), которую часто называют отклонением во времени (TDEV), может быть рассчитана на основе модифицированного отклонения Аллана (MDEV). TDEV основан на MDEV вместо исходного отклонения Аллана, потому что MDEV может различать фазовую модуляцию белого и мерцания (PM). Ниже приводится оценка временной дисперсии на основе модифицированной дисперсии Аллана:

${\displaystyle \sigma _{x}^{2}(\tau )={\frac {\tau ^{2}}{3}}{\bmod {\sigma }}_{y}^{2}(\tau ),}$

и аналогично для модифицированного отклонения Аллана к отклонению во времени :

${\displaystyle \sigma _{x}(\tau )={\frac {\tau }{\sqrt {3}}}{\bmod {\sigma }}_{y}(\tau ).}$

TDEV нормализован так, чтобы он был равен классическому отклонению для белого PM для постоянной времени τ =  τ ₀ . Чтобы понять масштабный коэффициент нормализации между статистическими показателями, используйте следующее статистическое правило: для независимых случайных величин X и Y дисперсия (σ _z ² ) суммы или разности ( z = x - y ) представляет собой квадрат суммы их дисперсий (σ _z ² = σ _x ² + σ _y ² ). Дисперсия суммы или разницы ( y = x _2τ - x _τ ) двух независимых выборок случайной величины вдвое превышает дисперсию случайной величины (σ _y ² = 2σ _x ² ). MDEV - это вторая разность независимых фазовых измерений ( x ), которые имеют дисперсию (σ _x ² ). Поскольку расчет представляет собой двойную разность, которая требует трех независимых фазовых измерений ( x _{2 τ} - 2 x _τ + x ), модифицированная дисперсия Аллана (MVAR) в три раза превышает дисперсию фазовых измерений.

Другие оценщики [ править ]

Дальнейшие разработки привели к созданию улучшенных методов оценки для той же меры стабильности, дисперсии / отклонения частоты, но они известны под разными названиями, такими как дисперсия Адамара , модифицированная дисперсия Адамара , общая дисперсия , модифицированная общая дисперсия и дисперсия Тео . Они отличаются лучшим использованием статистики для улучшения доверительных границ или способностью справляться с линейным дрейфом частоты.

Доверительные интервалы и эквивалентные степени свободы [ править ]

Статистические оценщики рассчитают оценочное значение для использованной выборки. Оценки могут отклоняться от истинного значения, и диапазон значений, который с некоторой вероятностью будет содержать истинное значение, называется доверительным интервалом . Доверительный интервал зависит от количества наблюдений в выборке, преобладающего типа шума и используемой оценки. Ширина также зависит от статистической достоверности, для которой значения доверительного интервала образуют ограниченный диапазон, таким образом, статистическая уверенность в том, что истинное значение находится в этом диапазоне значений. Для оценок переменной τ кратное τ ₀ n также является переменной.

Доверительный интервал [ править ]

Доверительный интервал может быть установлен с помощью критерия хи-квадрат распределение с использованием распределения образца дисперсии : ^{[4] [8]}

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

где s ² - выборочная дисперсия нашей оценки, σ ² - истинное значение дисперсии, df - степени свободы для оценки, а χ ² - степени свободы для определенной вероятности. Для вероятности 90%, охватывающей диапазон от 5% до 95% на кривой вероятности, верхний и нижний пределы могут быть найдены с помощью неравенства

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

который после перестановки для истинной дисперсии становится

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

Эффективные степени свободы [ править ]

В степени свободы представляет собой число свободных переменных , способных внести свой вклад в оценку. В зависимости от оценщика и типа шума эффективные степени свободы варьируются. Формулы оценки в зависимости от N и n были найдены эмпирически: ^[8]

Степени свободы по вариации Аллана

Тип шума	степени свободы
модуляция фазы белого (WPM)	${\displaystyle {\text{df}}\cong {\frac {(N+1)(N-2n)}{2(N-n)}}}$
фазовая модуляция мерцания (FPM)	${\displaystyle {\text{df}}\cong \exp \left[\left(\ln {\frac {N-1}{2n}}\ln {\frac {(2n+1)(N-1)}{4}}\right)^{-1/2}\right]}$
белая частотная модуляция (WFM)	${\displaystyle {\text{df}}\cong \left[{\frac {3(N-1)}{2n}}-{\frac {2(N-2)}{N}}\right]{\frac {4n^{2}}{4n^{2}+5}}}$
частотная модуляция мерцания (FFM)	${\displaystyle {\text{df}}\cong {\begin{cases}{\frac {2(N-2)}{2.3N-4.9}}&n=1\\{\frac {5N^{2}}{4n(N+3n)}}&n\geq 2\end{cases}}}$
частотная модуляция случайного блуждания (RWFM)	${\displaystyle {\text{df}}\cong {\frac {N-2}{n}}{\frac {(N-1)^{2}-3n(N-1)+4n^{2}}{(N-3)^{2}}}}$

Степенной шум [ править ]

Дисперсия Аллана будет обрабатывать различные типы степенного шума по- разному, что позволяет легко идентифицировать их и оценивать их интенсивность. По соглашению, ширина измерительной системы (высокая угловая частота) обозначается F _H .

Степенной закон дисперсии Аллана

Тип шума по степенному закону	Наклон фазового шума	Наклон частотного шума	Коэффициент мощности	Фазовый шум ${\displaystyle S_{x}(f)}$	Вариация Аллана ${\displaystyle \sigma _{y}^{2}(\tau )}$	Отклонение Аллана ${\displaystyle \sigma _{y}(\tau )}$
модуляция фазы белого (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}}}}$
фазовая модуляция мерцания (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}}}}$
белая частотная модуляция (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}}}}$
частотная модуляция мерцания (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}}}}$
частотная модуляция случайного блуждания (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}}}}$

Как найдено в ^{[11] [12]} и в современных формах. ^{[13] [14]}

Дисперсия Аллана не может различить WPM и FPM, но способна разрешить другие типы степенного шума. Чтобы различать WPM и FPM, необходимо использовать модифицированную дисперсию Аллана .

Приведенные выше формулы предполагают, что

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

и, таким образом, ширина полосы времени наблюдения намного меньше ширины полосы частот инструментов. Когда это условие не выполняется, все формы шума зависят от полосы пропускания прибора.

отображение α - μ [ править ]

Подробное отображение фазовой модуляции вида

${\displaystyle S_{x}(f)={\frac {1}{4\pi ^{2}}}h_{\alpha }f^{\alpha -2}={\frac {1}{4\pi ^{2}}}h_{\alpha }f^{\beta },}$

куда

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

или частотная модуляция формы

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

в вариацию Аллана формы

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

можно значительно упростить, предоставив отображение между α и μ . Для удобства также представлено отображение между α и K _α : ^[4]

Отображение вариации Аллана α - μ

α	β	μ	К α
−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}}}}$

Общее преобразование из фазового шума [ править ]

Сигнал со спектральным фазовым шумом в единицах рад ² / Гц может быть преобразован в дисперсию Аллана ^[14]${\displaystyle S_{\varphi }}$

${\displaystyle \sigma _{y}^{2}(\tau )={\frac {2}{\nu _{0}^{2}}}\int _{0}^{f_{b}}S_{\varphi }(f){\frac {\sin ^{4}(\pi \tau f)}{(\pi \tau )^{2}}}\,df.}$

Линейный ответ [ править ]

Хотя дисперсия Аллана предназначена для различения форм шума, она будет зависеть от некоторых, но не от всех линейных откликов на время. Они приведены в таблице:

Линейный отклик вариации Аллана

Линейный эффект	время ответа	частотный отклик	Вариация Аллана	Отклонение Аллана
фазовый сдвиг	${\displaystyle x_{0}}$	${\displaystyle 0}$	${\displaystyle 0}$	${\displaystyle 0}$
смещение частоты	${\displaystyle y_{0}t}$	${\displaystyle y_{0}}$	${\displaystyle 0}$	${\displaystyle 0}$
линейный дрейф	${\displaystyle {\frac {Dt^{2}}{2}}}$	${\displaystyle Dt}$	${\displaystyle {\frac {D^{2}\tau ^{2}}{2}}}$	${\displaystyle {\frac {D\tau }{\sqrt {2}}}}$

Таким образом, линейный дрейф будет способствовать выходному результату. При измерении реальной системы может потребоваться оценить линейный дрейф или другой механизм дрейфа и удалить его из временного ряда до расчета дисперсии Аллана. ^[13]

Свойства фильтра времени и частоты [ править ]

При анализе свойств дисперсии Аллана и др. Оказалось полезным рассмотреть свойства фильтра на нормированной частоте. Начиная с определения дисперсии Аллана для

${\displaystyle \sigma _{y}^{2}(\tau )={\frac {1}{2}}\langle ({\bar {y}}_{i+1}-{\bar {y}}_{i})^{2}\rangle ,}$

куда

${\displaystyle {\bar {y}}_{i}={\frac {1}{\tau }}\int \limits _{0}^{\tau }y(i\tau +t)\,dt.}$

Заменив временной ряд на вариант с преобразованием Фурье, дисперсию Аллана можно выразить в частотной области как${\displaystyle y_{i}}$${\displaystyle S_{y}(f)}$

${\displaystyle \sigma _{y}^{2}(\tau )=\int _{0}^{\infty }S_{y}(f){\frac {2\sin ^{4}\pi \tau f}{(\pi \tau f)^{2}}}\,df.}$

Таким образом, передаточная функция для дисперсии Аллана равна

${\displaystyle \left\vert H_{A}(f)\right\vert ^{2}={\frac {2\sin ^{4}\pi \tau f}{(\pi \tau f)^{2}}}.}$

Функции смещения [ править ]

М -sample дисперсия, и определяются частный случай Аллан, будут испытывать систематические смещения в зависимости от различного числа образцов М и различных отношений между Т и т . Для того , чтобы решить эти уклоны, смещения-функции B ₁ и B ₂ было определено ^[15] и позволяет преобразование между различными M и T значений.

Этих функций смещения недостаточно для обработки смещения, возникающего в результате конкатенации M отсчетов во время наблюдения Mτ _{0 в} течение MT ₀ с мертвым временем, распределенным между M блоками измерения, а не в конце измерения. Это привело к необходимости смещения B ₃ . ^[16]

Функции смещения оцениваются для конкретного значения µ, поэтому отображение α – µ должно быть выполнено для доминирующей формы шума, найденной с помощью идентификации шума . В качестве альтернативы ^{[3] [15]} значение µ преобладающей формы шума может быть выведено из измерений с использованием функций смещения.

B ₁ функция смещения [ править ]

Функция смещения B ₁ связывает дисперсию M -выборки с дисперсией 2-х выборок (дисперсия Аллана), сохраняя постоянным время между измерениями T и время для каждого измерения τ . Он определяется ^[15] как

${\displaystyle B_{1}(N,r,\mu )={\frac {\left\langle \sigma _{y}^{2}(N,T,\tau )\right\rangle }{\left\langle \sigma _{y}^{2}(2,T,\tau )\right\rangle }},}$

куда

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

Функция смещения становится после анализа

${\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]}}.}$

Функция смещения B ₂ [ править ]

Функция смещения B ₂ связывает двухвыборочную дисперсию для времени выборки T с двухвыборочной дисперсией (дисперсия Аллана), сохраняя количество выборок N = 2 и время наблюдения τ постоянным. Он определяется ^[15] как

${\displaystyle B_{2}(r,\mu )={\frac {\left\langle \sigma _{y}^{2}(2,T,\tau )\right\rangle }{\left\langle \sigma _{y}^{2}(2,\tau ,\tau )\right\rangle }},}$

куда

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

Функция смещения становится после анализа

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

Функция смещения B ₃ [ править ]

Функция смещения B ₃ связывает двухвыборочную дисперсию для времени выборки MT ₀ и времени наблюдения Mτ ₀ с двухвыборочной дисперсией (дисперсия Аллана) и определяется ^[16] как

${\displaystyle B_{3}(N,M,r,\mu )={\frac {\left\langle \sigma _{y}^{2}(N,M,T,\tau )\right\rangle }{\left\langle \sigma _{y}^{2}(N,T,\tau )\right\rangle }},}$

куда

${\displaystyle T=MT_{0},}$

${\displaystyle \tau =M\tau _{0}.}$

Функция смещения B ₃ полезна для корректировки неперекрывающихся и перекрывающихся значений оценки переменной τ на основе измерений мертвого времени времени τ ₀ наблюдения и времени между наблюдениями T ₀ до нормальных оценок мертвого времени.

Функция смещения становится после анализа (для  случая N = 2)

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

куда

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

функция смещения τ [ править ]

Хотя формально это не сформулировано, оно было косвенно выведено как следствие отображения α - µ . При сравнении двух показателей дисперсии Аллана для разных τ , предполагая один и тот же доминирующий шум в форме одного и того же коэффициента µ, смещение можно определить как

${\displaystyle B_{\tau }(\tau _{1},\tau _{2},\mu )={\frac {\left\langle \sigma _{y}^{2}(2,\tau _{2},\tau _{2})\right\rangle }{\left\langle \sigma _{y}^{2}(2,\tau _{1},\tau _{1})\right\rangle }}.}$

Функция смещения становится после анализа

${\displaystyle B_{\tau }(\tau _{1},\tau _{2},\mu )=\left({\frac {\tau _{2}}{\tau _{1}}}\right)^{\mu }.}$

Преобразование между значениями [ править ]

Для преобразования одного набора измерений в другой можно собрать функции смещения B ₁ , B ₂ и τ. Сначала функция B ₁ преобразует значение ( N ₁ ,  T ₁ ,  τ ₁ ) в (2,  T ₁ ,  τ ₁ ), из которого функция B ₂ преобразуется в значение (2,  τ ₁ ,  τ ₁ ), таким образом дисперсия Аллана при τ ₁ . Мера дисперсии Аллана может быть преобразована с использованием функции смещения τ из τ ₁до τ ₂ , из которого затем (2,  T ₂ ,  τ ₂ ) с использованием B ₂ и, наконец, с использованием B ₁ в дисперсию ( N ₂ ,  T ₂ ,  τ ₂ ). Полное преобразование становится

${\displaystyle \left\langle \sigma _{y}^{2}(N_{2},T_{2},\tau _{2})\right\rangle =\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]\left\langle \sigma _{y}^{2}(N_{1},T_{1},\tau _{1})\right\rangle ,}$

куда

${\displaystyle r_{1}={\frac {T_{1}}{r_{1}}},}$

${\displaystyle r_{2}={\frac {T_{2}}{r_{2}}}.}$

Аналогично, для объединенных измерений с использованием M секций логическое расширение становится

${\displaystyle \left\langle \sigma _{y}^{2}(N_{2},M_{2},T_{2},\tau _{2})\right\rangle =\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]\left\langle \sigma _{y}^{2}(N_{1},M_{1},T_{1},\tau _{1})\right\rangle .}$

Проблемы с измерениями [ править ]

При выполнении измерений для расчета дисперсии Аллана или отклонения Аллана ряд проблем может привести к ухудшению результатов измерений. Здесь рассматриваются эффекты, характерные для дисперсии Аллана, результаты которых могут быть смещенными.

Пределы полосы пропускания измерения [ править ]

Ожидается, что измерительная система будет иметь полосу пропускания, равную или ниже скорости Найквиста , как описано в теореме Шеннона – Хартли . Как видно из формул степенного шума, модуляция белого и фликкер-шума зависит от верхней угловой частоты (предполагается, что эти системы фильтруются только нижними частотами ). Учитывая свойство частотного фильтра, ясно видно, что низкочастотный шум оказывает большее влияние на результат. Для типов шума с относительно плоской фазовой модуляцией (например, WPM и FPM) фильтрация имеет значение, тогда как для типов шума с большим наклоном верхний предел частоты становится менее важным, если предположить, что полоса пропускания системы измерения является широкой по сравнению с заданным соотношением${\displaystyle f_{H}}$${\displaystyle \tau }$

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

Когда это предположение не выполняется, эффективную полосу пропускания необходимо записывать вместе с измерением. Заинтересованным следует проконсультироваться по NBS TN394. ^[11]${\displaystyle f_{H}}$

Если, однако, отрегулировать полосу пропускания оценщика, используя целые числа, кратные времени выборки , то влияние полосы пропускания системы может быть уменьшено до незначительного уровня. Для нужд электросвязи такие методы требовались, чтобы гарантировать сопоставимость измерений и дать производителям некоторую свободу при различных реализациях. Рекомендация МСЭ-Т Rec. G.813 ^[17] для измерения TDEV.${\displaystyle n\tau _{0}}$

Можно рекомендовать игнорировать первые кратные, чтобы большая часть обнаруженного шума находилась в пределах полосы пропускания полосы пропускания измерительной системы.${\displaystyle \tau _{0}}$

Дальнейшие разработки по дисперсии Аллана были выполнены, чтобы позволить аппаратной полосе пропускания быть уменьшенной программными средствами. Эта разработка программного обеспечения полосы пропускания позволила устранить остающийся шум, и теперь метод называется модифицированной дисперсией Аллана . Этот метод уменьшения полосы пропускания не следует путать с расширенным вариантом модифицированной дисперсии Аллана , который также изменяет полосу пропускания сглаживающего фильтра.

Мертвое время в измерениях [ править ]

Многие приборы для измерения времени и частоты имеют стадии постановки на охрану, временной развертки, времени обработки и могут затем повторно запускать постановку на охрану. Время постановки на охрану - от момента срабатывания постановки на охрану до момента возникновения стартового события на стартовом канале. Таким образом, временная база гарантирует, что до принятия события на канале остановки в качестве события остановки проходит минимальное количество времени. Количество событий и время, прошедшее между событием запуска и событием остановки, записываются и представляются во время обработки. Когда происходит обработка (также известная как время задержки), прибор обычно не может выполнить другое измерение. После завершения обработки прибор в непрерывном режиме снова запускает схему плеча. Время между событием остановки и следующим событием запуска становится мертвым временем., во время которого сигнал не наблюдается. Такое мертвое время вносит систематические погрешности измерения, которые необходимо компенсировать, чтобы получить правильные результаты. Для таких систем измерения время T будет обозначать время между соседними стартовыми событиями (и, следовательно, измерениями), а обозначать длину временной базы, то есть номинальную продолжительность между стартовым и конечным событиями любого измерения.${\displaystyle \tau }$

Влияние мертвого времени на измерения оказывает такое влияние на получаемый результат, что было проведено много исследований поля для правильной количественной оценки его свойств. Введение счетчиков с нулевым мертвым временем устранило необходимость в этом анализе. Счетчик с нулевым мертвым временем обладает тем свойством, что событие остановки одного измерения также используется в качестве стартового события следующего события. Такие счетчики создают серию пар событий и временных меток, по одной для каждого канала, разделенных временной базой. Такие измерения также оказались полезными при упорядоченном анализе временных рядов.

Измерения, выполняемые с мертвым временем, можно скорректировать с помощью функции смещения B ₁ , B ₂ и B ₃ . Таким образом, мертвое время как таковое не запрещает доступ к дисперсии Аллана, но делает его более проблематичным. Необходимо знать мертвое время, чтобы можно было установить время между выборками T.

Длина измерения и эффективное использование образцов [ править ]

Изучая влияние на доверительные интервалы, которые имеет длина N серии выборок, и влияние переменного параметра τ n, доверительные интервалы могут стать очень большими, поскольку эффективная степень свободы может стать небольшой для некоторой комбинации N и n для доминирующей формы шума (для этого τ ).

Эффект может заключаться в том, что оценочное значение может быть намного меньше или намного больше, чем реальное значение, что может привести к ложным выводам о результате.

Рекомендуется нанести доверительный интервал на график вместе с данными, чтобы читатель графика мог быть осведомлен о статистической неопределенности значений.

Рекомендуется, чтобы длина последовательности выборок, т. Е. Количество выборок N , поддерживалась на высоком уровне, чтобы гарантировать, что доверительный интервал невелик в интересующем диапазоне τ .

Рекомендуется, чтобы τ диапазон , как прокатилась по τ ₀ умножителя п ограничена в верхнем конце относительно N , таким образом, что чтение сюжета не спутать высоко нестабильными значениями оценщик.

Рекомендуется использовать оценщики, обеспечивающие лучшие значения степеней свободы, вместо оценок дисперсии Аллана или в качестве дополнения к ним там, где они превосходят оценщики дисперсии Аллана. Среди них следует учитывать оценки общей дисперсии и дисперсии Тео .

Доминирующий тип шума [ править ]

Большое количество констант преобразования, поправок смещения и доверительных интервалов зависит от доминирующего типа шума. Для правильной интерпретации доминирующий тип шума для конкретного интересующего τ должен быть идентифицирован посредством идентификации шума. Неспособность определить доминирующий тип шума приведет к смещению значений. Некоторые из этих смещений могут иметь величину нескольких порядков, поэтому они могут иметь большое значение.

Линейный дрейф [ править ]

Систематическое воздействие на сигнал отменяется лишь частично. Сдвиг фазы и частоты отменяется, но линейный дрейф или другие формы полиномиальных фазовых кривых высокой степени не отменяются и, таким образом, формируют ограничение измерения. Можно использовать подгонку кривой и удаление систематического смещения. Часто бывает достаточно устранения линейного дрейфа. Также можно использовать оценки линейного дрейфа, такие как дисперсия Адамара . Удаление линейного дрейфа может быть использовано с использованием оценки на основе момента.

Смещение оценщика измерительного прибора [ править ]

Традиционные инструменты обеспечивали измерение только отдельных событий или пар событий. Внедрение улучшенного статистического инструмента перекрывающихся измерений Дж. Дж. Снайдером ^[7] позволило значительно улучшить разрешение при считывании частоты, нарушив традиционный баланс цифр / временной развертки. Хотя такие методы полезны по своему прямому назначению, использование таких сглаженных измерений для расчетов дисперсии Аллана могло бы создать ложное впечатление о высоком разрешении ^{[18] [19] [20],} но при более длительном значении τ эффект постепенно устраняется, а более низкие значения τобласть измерения имеет смещенные значения. Это смещение дает более низкие значения, чем следовало бы, поэтому это чрезмерно оптимистичное (предполагающее, что низкие числа - это то, что нужно) смещение, которое снижает удобство использования измерения, а не улучшает его. Такие интеллектуальные алгоритмы обычно можно отключить или иным образом обойти с помощью режима отметки времени, который гораздо предпочтительнее, если он доступен.

Практические измерения [ править ]

Хотя можно разработать несколько подходов к измерению дисперсии Аллана, простой пример может проиллюстрировать, как можно проводить измерения.

Измерение [ править ]

Все измерения дисперсии Аллана будут фактически сравнением двух разных часов. Рассмотрим опорные часы и тестируемое устройство (DUT), и оба имеют общую номинальную частоту 10 МГц. Временной интервал счетчика используется для измерения времени между передним фронтом опорного канала (А) и по переднему фронту тестируемого устройства.

Чтобы обеспечить равномерно распределенные измерения, опорные часы будут разделены на частоту измерения, запустив счетчик временных интервалов (вход ARM). Эта скорость может быть 1 Гц ( с использованием 1 PPS выхода опорного тактового сигнала), но другие скорости , как 10 Гц и 100 Гц , также могут быть использованы. Скорость, с которой счетчик временных интервалов может завершить измерение, вывести результат и подготовиться к следующей руке, ограничивает частоту запуска.

Затем компьютер полезен для записи серии наблюдаемых разниц во времени.

Постобработка [ править ]

Записанный временной ряд требует постобработки, чтобы развернуть свернутую фазу, чтобы обеспечить непрерывную фазовую ошибку. При необходимости также следует исправить ошибки записи и измерения. Должны быть выполнены оценка сноса и удаление сноса, механизм сноса должен быть идентифицирован и понят для источников. Ограничения по дрейфу при измерениях могут быть серьезными, поэтому необходимо дать генераторам возможность стабилизироваться за счет достаточного времени включения питания.

Затем можно рассчитать дисперсию Аллана с использованием данных оценщиков, и для практических целей следует использовать оценщик с перекрытием, поскольку он использует данные лучше, чем оценщик без перекрытия. Другие оценщики, такие как оценщики общей или Тео-дисперсии, также могут быть использованы, если применяются поправки на смещение, так что они дают результаты, совместимые с дисперсией Аллана.

Для построения классических графиков отклонение Аллана (квадратный корень из дисперсии Аллана) строится в логарифмическом формате в зависимости от интервала наблюдения  τ .

Оборудование и программное обеспечение [ править ]

Счетчик интервалов времени обычно представляет собой серийный счетчик, имеющийся в продаже. Ограничивающие факторы включают однократное разрешение, джиттер триггера, скорость измерений и стабильность эталонных часов. Компьютерный сбор и последующая обработка могут быть выполнены с использованием существующего коммерческого или общедоступного программного обеспечения. Существуют высокоразвитые решения, позволяющие выполнять измерения и вычисления в одном устройстве.

История исследований [ править ]

Область частотной стабильности изучается давно. Однако в 1960-х годах было обнаружено, что последовательных определений не было. Симпозиум NASA-IEEE по краткосрочной стабильности в ноябре 1964 года ^[21] привел к выпуску специального выпуска журнала IEEE Proceedings on Frequency Stability в феврале 1966 года.

Симпозиум NASA-IEEE объединил множество областей и применений краткосрочной и долгосрочной стабильности с докладами многих разных участников. Статьи и групповые обсуждения сходятся во мнении о существовании частотного фликкер-шума и стремлении достичь общего определения как краткосрочной, так и долгосрочной стабильности.

Важные статьи, в том числе статьи Дэвида Аллана, ^[3] Джеймса А. Барнса, ^[22] Л.С. Катлера и К.Л. Сирла ^[1] и Д.Б. Лисона, ^[2], появились в Трудах IEEE по стабильности частоты и помогли сформировать эту область.

В статье Дэвида Аллана анализируется классическая дисперсия частоты M- выборки, решается проблема мертвого времени между измерениями вместе с начальной функцией смещения. ^[3] Хотя исходная функция смещения Аллана не предполагает мертвого времени, его формулы включают вычисления мертвого времени. В его статье анализируется случай M выборок частоты (называемых в статье N) и оценок дисперсии. Он обеспечивает теперь стандартное отображение α – µ, явно основанное на работе Джеймса Барнса ^[22] в том же выпуске.

Случай двухвыборочной дисперсии является частным случаем дисперсии M- выборки, которая дает среднее значение производной частоты. Аллан неявно использует двухвыборочную дисперсию в качестве базового случая, поскольку для произвольно выбранного M значения могут быть перенесены через двухвыборочную дисперсию в M -выборочную дисперсию. Не было четко заявлено о предпочтении двухвыборочной дисперсии, даже если бы инструменты были предоставлены. Однако эта статья заложила основу для использования двухвыборочной дисперсии как способа сравнения других M- выборочных дисперсий.

Джеймс Барнс значительно расширил работу по функциям смещения ^[15], введя современные функции смещения B ₁ и B ₂ . Достаточно любопытно, что это относится к дисперсии M -выборки как «дисперсия Аллана», при этом ссылаясь на статью Аллана «Статистика атомных стандартов частоты». ^[3] С помощью этих современных функций смещения может быть выполнено полное преобразование между мерами дисперсии M- выборки различных значений M , T и τ путем преобразования дисперсии 2-выборки.

Джеймс Барнс и Дэвид Алан продлена функция смещения с B ₃ функцией ^[16] , чтобы обрабатывать конкатенированную образцы оценщика смещение. Это было необходимо для того, чтобы справиться с новым использованием объединенных наблюдений выборки с мертвым временем между ними.

В 1970 году Технический комитет IEEE по частоте и времени в составе группы IEEE по приборам и измерениям представил краткое изложение этой области, опубликованное как NBS Technical Notice 394. ^[11] Эта статья была первой из серии более образовательных и практических материалов. документы, помогающие коллегам-инженерам понять эту область. В этой статье рекомендуется двухвыборочная дисперсия с T = τ , называемая дисперсией Аллана (теперь без кавычек). Выбор такой параметризации позволяет хорошо обрабатывать некоторые формы шума и получать сопоставимые измерения; по существу, это наименьший общий знаменатель с помощью функций смещения B ₁ и B ₂ .

Дж. Дж. Снайдер предложил улучшенный метод оценки частоты или дисперсии с использованием выборочной статистики для частотомеров. ^[7] Чтобы получить более эффективные степени свободы из доступного набора данных, хитрость заключается в использовании перекрывающихся периодов наблюдения. Это обеспечивает улучшение √ n и было включено в перекрывающуюся оценку дисперсии Аллана . ^[8] Также была включена программная обработка переменных τ. ^[8] Эта разработка улучшила классические оценки дисперсии Аллана, что также послужило прямым вдохновением для работы над модифицированной дисперсией Аллана .

Хау, Аллан и Барнс представили анализ доверительных интервалов, степеней свободы и установленных оценок. ^[8]

Образовательные и практические ресурсы [ править ]

Поле времени и частоты, а также использование вариации Аллана, девиации Аллана и друзей - это область, включающая множество аспектов, для которых как понимание концепций, так и практические измерения и постобработка требуют внимательности и понимания. Таким образом, существует область учебных материалов, охватывающая около 40 лет. Поскольку они отражают развитие исследований своего времени, они сосредоточены на обучении различным аспектам с течением времени, и в этом случае обзор доступных ресурсов может быть подходящим способом поиска нужного ресурса.

Первым значимым резюме является Техническая записка NBS № 394 «Характеристики стабильности частоты». ^[11] Это продукт Технического комитета по частоте и времени Группы IEEE по контрольно-измерительным приборам и измерениям. Он дает первый обзор области, формулируя проблемы, определяя основные вспомогательные определения и углубляясь в дисперсию Аллана, функции смещения B ₁ и B ₂ , преобразование измерений во временной области. Это полезно, так как это одна из первых ссылок для табулирования дисперсии Аллана для пяти основных типов шума.

Классической ссылкой является Монография NBS 140 ^[23] 1974 года, в которой в главе 8 есть «Статистика анализа временных и частотных данных». ^[24] Это расширенный вариант Технической записки NBS 394, существенно добавляющий методы измерения и практическую обработку значений.

Важным дополнением станут Свойства источников сигналов и методов измерения . ^[8] Он охватывает эффективное использование данных, доверительные интервалы, эффективную степень свободы, а также вводит перекрывающуюся оценку дисперсии Аллана. Это настоятельно рекомендуется для чтения по этим темам.

Стандарт IEEE 1139 Стандартные определения физических величин для фундаментальной метрологии частоты и времени ^[4] выходит за рамки стандартного всеобъемлющего справочного и образовательного ресурса.

Современная книга, посвященная телекоммуникациям, - это Стефано Брегни «Синхронизация цифровых телекоммуникационных сетей». ^[13] Это резюмирует не только область, но и большую часть его исследований в этой области до того момента. Он призван включать как классические меры, так и меры, специфичные для электросвязи, такие как MTIE. Это удобный помощник при проведении измерений, связанных со стандартами электросвязи.

Специальная публикация NIST 1065 «Справочник по анализу стабильности частоты» WJ Riley ^[14] рекомендуется к прочтению всем, кто хочет продолжить изучение этой области. Он богат ссылками, а также охватывает широкий спектр показателей, предубеждений и связанных функций, которые должны быть доступны современному аналитику. Далее в нем описывается общая обработка, необходимая для современного инструмента.

Использует [ редактировать ]

Дисперсия Аллана используется в качестве меры стабильности частоты в различных прецизионных генераторах, таких как кварцевые генераторы , атомные часы и лазеры со стабилизированной частотой в течение секунды или более. Кратковременная стабильность (менее секунды) обычно выражается как фазовый шум . Дисперсии Аллана также используются для характеристики стабильности смещения гироскопов , в том числе волоконно - оптических гироскопов , полусферический резонаторы гироскопов и MEMS гироскопов и акселерометров. ^{[25] [26]}

50 лет [ править ]

В 2016 году IEEE-UFFC собирается опубликовать «Специальный выпуск, посвященный 50-летию Allan Variance (1966–2016)». ^[27] Приглашенным редактором этого выпуска будет бывший коллега Дэвида по NIST , Джуда Левайн, который недавно стал лауреатом II премии Раби .

См. Также [ править ]

Ссылки [ править ]

Внешние ссылки [ править ]

• Учебные ресурсы UFFC по контролю частоты
• Инструмент поиска публикаций NIST
• Обзор дисперсии Аллана Дэвида В. Аллана
• Официальный веб-сайт Дэвида В. Аллана
• Публикации JPL - Анализ и статистика шума
• Публикации Уильяма Райли
• Stable32 , Программное обеспечение для анализа стабильности частоты, Уильям Райли
• Публикации Стефано Брегни
• Публикации Энрико Рубиолы
• Allanvar: пакет R для определения характеристик ошибок датчика с помощью дисперсии Аллана
• Программное обеспечение Alavar для Windows с инструментами отчетности; Бесплатное ПО
• Библиотека Python с открытым исходным кодом AllanTools для вариации Аллана
• MATLAB AVAR приложение MATLAB с открытым исходным кодом