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

редактировать
Измерение стабильности частоты в часах и генераторах

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

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

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

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

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

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

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

Пример графика отклонения Аллана часов. При очень коротком времени наблюдения τ отклонение Аллана велико из-за шума. При больших τ он уменьшается, так как шум усредняется. При еще большем значении τ отклонение Аллана снова начинает увеличиваться, предполагая, что тактовая частота постепенно дрейфует из-за изменений температуры, старения компонентов или других подобных факторов. Планки погрешностей увеличиваются с увеличением τ просто потому, что получение большого количества точек данных для больших τ занимает много времени.
Содержание
  • 1 Предпосылки
  • 2 Интерпретация значения
  • 3 Определения
    • 3,1 млн -выборочная дисперсия
    • 3.2 Вариация Аллана
    • 3.3 Отклонение Аллана
  • 4 Дополнительные определения
    • 4.1 Модель осциллятора
    • 4.2 Временная ошибка
    • 4.3 Функция частоты
    • 4.4 Дробная частота
    • 4.5 Среднее относительная частота
  • 5 Оценщиков
    • 5.1 Условные обозначения
    • 5.2 Фиксированные оценки τ
    • 5.3 Неперекрывающиеся оценки τ переменных
    • 5.4 Перекрывающиеся переменные оценки τ
    • 5.5 Модифицированная дисперсия Аллана
    • 5.6 Временная стабильность Оценщики
    • 5.7 Другие оценщики
  • 6 Доверительные интервалы и эквивалентные степени свободы
    • 6.1 Доверительный интервал
    • 6.2 Эффективные степени свободы
  • 7 Степенный шум
    • 7.1 Отображение α – μ
    • 7.2 Общее преобразование из фазового шума
  • 8 Линейная характеристика
  • 9 Свойства частотно-временного фильтра
  • 10 Функции смещения
    • 10.1 B 1 Функция смещения
    • 10,2 B 2 функция смещения
    • 10,3 B 3 функция смещения
    • 10,4 функция смещения τ
    • 10,5 Преобразование между значениями
  • 11 Проблемы измерения
    • 11.1 Пределы ширины полосы измерения
    • 11.2 Время простоя при измерениях
    • 11.3 Длина измерения и эффективное использование образцов
    • 11.4 Доминирующий тип шума
    • 11.5 Линейный дрейф
    • 11.6 Смещение оценки измерительного прибора
  • 12 Практические измерения
    • 12.1 Измерение
    • 12.2 Постобработка
    • 12.3 Оборудование и программное обеспечение
  • 13 История исследований
  • 14 Образовательные и практические ресурсы
  • 15 Использование
  • 16 50-летие
  • 17 См. Также
  • 18 Источники
  • 19 Внешние ссылки
Предпосылки

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

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

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

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

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

τ = n τ 0. {\ displaystyle \ tau = n \ tau _ {0}.}{\displaystyle \tau =n\tau _{0}.}

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

Интерпретация значения

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

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

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

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

Определения

Дисперсия M-выборки

M {\ displaystyle M}M -выборочная дисперсия определяется (здесь в модернизированной форме обозначений) как

σ y 2 (M, T, τ) = 1 M - 1 {∑ i = 0 M - 1 [x (i T + τ) - x (i T) τ] 2 - 1 M [∑ i = 0 M - 1 x (i T + τ) - x (i T) τ] 2}, {\ 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 \ 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 (i T + \ tau) -x (iT)} {\ tau}} \ right] ^ {2} \ right \},}

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

σ y 2 (M, T, τ) = 1 M - 1 {∑ i = 0 M - 1 y ¯ i 2 - 1 M [∑ я знак равно 0 M - 1 y ¯ я] 2}, {\ displaystyle \ sig ma _ {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 \ 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 \},}

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

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

дисперсия Аллана

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

σ y 2 (τ) = ⟨σ y 2 (2, τ, τ)⟩, {\ Displaystyle \ sigma _ {y} ^ {2} (\ tau) = \ langle \ sigma _ {y} ^ {2} (2, \ tau, \ tau) \ rangle,}{\ displaystyle \ sigma _ {y} ^ {2} (\ tau) = \ langle \ sigma _ {y} ^ {2} (2, \ tau, \ tau) \ rangle,}

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

σ y 2 (τ) = 1 2 ⟨(y ¯ n + 1 - y ¯ n) 2⟩ = 1 2 τ 2 ⟨(xn + 2 - 2 xn + 1 + xn) 2⟩, {\ 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,}{\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,}

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

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

T = τ. {\ displaystyle T = \ tau.}{\displaystyle T=\tau.}

Отклонение Аллана

Как и в случае с стандартным отклонением и дисперсией, отклонение Аллана определяется как квадратный корень из дисперсия Аллана:

σ y (τ) = σ y 2 (τ). {\ displaystyle \ sigma _ {y} (\ tau) = {\ sqrt {\ sigma _ {y} ^ {2} (\ tau)}}.}{\ displaystyle \ sigma _ {y} (\ tau) = {\ sqrt {\ sigma _ {y} ^ {2} (\ tau)}}.}
Дополнительные определения

Модель осциллятора

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

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

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

ω n = 2 π ν n. {\ displaystyle \ omega _ {\ text {n}} = 2 \ pi \ nu _ {\ text {n}}.}{\displaystyle \omega _{\text{n}}=2\pi \nu _{\text{n}}.}

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

Φ (t) = ω nt + ϕ (t) = 2 π ν nt + ϕ (t). {\ displaystyle \ Phi (t) = \ omega _ {\ text {n}} t + \ phi (t) = 2 \ pi \ nu _ {\ text {n}} t + \ phi (t).}{\ displaystyle \ Phi (t) = \ omega _ {\ text {n}} t + \ phi (t) = 2 \ pi \ nu _ {\ text {n}} t + \ phi (t).}

Ошибка времени

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

x (t) = ϕ (t) 2 π ν n = Φ (t) 2 π ν n - t = T (t) - t. {\ displaystyle x (t) = {\ frac {\ phi (t)} {2 \ pi \ nu _ {\ text {n}}}} = {\ frac {\ Phi (t)} {2 \ pi \ nu _ {\ text {n}}}} - t = T (t) -t.}{\ displaystyle x (t) = {\ frac {\ phi (t)} {2 \ pi \ nu _ {\ text {n}}}} = {\ frac {\ Phi (t)} {2 \ pi \ nu _ {\ text {n}}}} -t = T (t) -t.}

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

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

Функция частоты

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

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

Дробная частота

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

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

Средняя относительная частота

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

y ¯ (t, τ) = 1 τ ∫ 0 τ y (t + tv) dtv, {\ displaystyle {\ bar {y}} (t, \ tau) = {\ frac {1} {\ tau}} \ int \ limits _ {0} ^ {\ tau} y (t + t_ {v}) \, dt_ {v},}{\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), мы можем без ограничения общности переписать ее как

y ¯ (t, τ) = x (t + τ) - x (t) τ. {\ displaystyle {\ bar {y}} (t, \ tau) = {\ frac {x (t + \ tau) -x (t)} {\ tau}}.}{\ displaystyle {\ bar { y}} (t, \ tau) = {\ frac {x (t + \ tau) -x (t)} { \ tau}}.}
Оценщики

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

Условные обозначения

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

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

N = M + 1. {\ Displaystyle N = M + 1.}{\displaystyle N=M+1.}
  • Для time- error sample series, x i обозначает i-ую выборку функции непрерывного времени x (t), как задано
xi = x (i T), {\ displaystyle x_ {i } = x (iT),}{\displaystyle x_{i}=x(iT),}

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

В серии выборок ошибка времени пусть N обозначает количество выборок (x 0... x N-1) в серии. В традиционном соглашении используются индексы от 1 до N.

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

, что дает

y ¯ i = 1 τ ∫ 0 τ y (i T + tv) dtv = x (i T + τ) - x (i T) τ. {\ displaystyle {\ bar {y}} _ {i} = {\ frac {1} {\ tau}} \ int \ limits _ {0} ^ {\ tau} y (iT + t_ {v}) \, dt_ {v} = {\ frac {x (iT + \ tau) -x (iT)} {\ 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 равно τ, это становится

y ¯ i = xi + 1 - xi τ. {\ displaystyle {\ bar {y}} _ {i} = {\ frac {x_ {i + 1} -x_ {i}} {\ tau}}.}{\ displaystyle {\ bar {y}} _ {i} = {\ гидроразрыва {x_ {i + 1} -x_ {i}} {\ tau}}.}

средняя дробная частота выборочная серия позволяет M обозначать количество выборок (y ¯ 0… y ¯ M - 1 {\ displaystyle {\ bar {y}} _ {0} \ ldots {\ bar {y}} _ {M -1}}{\ bar {y}} _ {0} \ ldots {\ bar {y}} _ {{M-1}} ) в серии. В традиционном соглашении используются индексы от 1 до M.

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

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

Фиксированное время τ оценки

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

σ y 2 (τ, M) = AVAR (τ, M) = 1 2 (M - 1) ∑ i = 0 M - 2 (Y ¯ я + 1 - Y ¯ я) 2, {\ 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, M) ={\text{AVAR}}(\tau,M)={\frac {1}{2(M-1)}}\sum _{i=0}^{M-2}({\bar {y }}_{i+1}-{\bar {y}}_{i})^{2},}

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

σ y 2 (τ, N) = AVAR (τ, N) = 1 2 τ 2 (N - 2) ∑ i Знак равно 0 N - 3 (xi + 2-2 xi + 1 + xi) 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}.}{\ 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}.}

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

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

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

σ y 2 (n τ 0, M) = AVAR (n τ 0, M) = 1 2 M - 1 n ∑ i = 0 M - 1 n - 1 (y ¯ ni + n - y ¯ ni) 2 {\ 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 \ 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}}

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

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

σ y 2 (n τ 0, N) = AVAR (n τ 0, N) = 1 2 n 2 τ 0 2 (N - 1 n - 1) ∑ i = 0 N - 1 n - 2 (xni + 2 n - 2 xni + n + xni) 2 {\ 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 \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}}

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

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

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

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

σ y 2 (n τ 0, M) = AVAR (n τ 0, M) = 1 2 n 2 (M - 2 n + 1) ∑ j = 0 M - 2 n (∑ i знак равно jj + n - 1 yi + n - yi) 2 знак равно 1 2 (M - 2 n + 1) ∑ j = 0 M - 2 n (y ¯ j + n - y ¯ j) 2, {\ 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},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},}

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

σ y 2 (n τ 0, N) = AVAR (n τ 0, N) = 1 2 n 2 τ 0 2 (N - 2 n) ∑ i = 0 N - 2 n - 1 (xi + 2 n - 2 xi + n + xi) 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}.}{\ 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}.}

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

Модифицированная дисперсия Аллана

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

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

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

σ x 2 (τ) = τ 2 3 Mod σ y 2 (τ), {\ displaystyle \ sigma _ {x} ^ {2} (\ tau) = {\ frac {\ tau ^ {2}} {3}} {\ text {Mod}} \ sigma _ {y} ^ {2} (\ tau),}{\ displaystyle \ sigma _ {x} ^ {2} (\ tau) = {\ frac {\ tau ^ {2}} {3}} {\ text {Mod}} \ sigm a _ {y} ^ {2} (\ tau),}

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

σ x (τ) = τ 3 Mod σ y (τ). {\ displaystyle \ sigma _ {x} (\ tau) = {\ frac {\ tau} {\ sqrt {3}}} {\ text {Mod}} \ sigma _ {y} (\ tau).}{\displaystyle \sigma _{x}(\tau)={\frac {\tau }{\sqrt {3}}}{\text{Mod}}\sigma _{y}(\tau).}

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

Другие средства оценки

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

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

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

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

Доверительный интервал может быть установлен с использованием распределения хи-квадрат с использованием распределения дисперсии выборки :

χ 2 знак равно df s 2 σ 2, {\ displaystyle \ chi ^ {2} = {\ frac {{\ text {df}} \, s ^ ​​{2}} {\ sigma ^ {2}}},}{\ displaystyle \ chi ^ {2} = {\ frac {{\ text {df}} \, s ^ {2}} {\ sigma ^ {2}}},}

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

χ 2 (0,05) ≤ df s 2 σ 2 ≤ χ 2 (0,95), {\ Displaystyle \ чи ^ {2} (0,05) \ leq {\ frac {{\ text {df}} \, s ^ ​​{2}} {\ sigma ^ {2}}} \ leq \ chi ^ {2} (0,95),}{\displaystyle \chi ^{2}(0.05)\leq {\frac {{\text{df}}\,s^{2}}{\sigma ^{2}}}\leq \chi ^{2}(0.95),}

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

df s 2 χ 2 (0,95) ≤ σ 2 ≤ df s 2 χ 2 (0,05). {\ displaystyle {\ frac {{\ text {df}} \, s ^ ​​{2}} {\ chi ^ {2} (0.95)}} \ leq \ sigma ^ {2} \ leq {\ frac { {\ text {df}} \, s ^ ​​{2}} {\ chi ^ {2} (0.05)}}.}{\ displaystyle {\ frac {{\ text {df}} \, s ^ ​​{2}} {\ chi ^ {2} (0.95)}} \ leq \ sigma ^ {2} \ leq {\ frac {{\ text {df}} \, s ^ ​​{2}} {\ чи ^ {2} (0,05)}}.}

Эффективные степени свободы

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

Степени свободы по дисперсии Аллана
Тип шумаСтепени свободы
Фазовая модуляция белого (WPM)df ≅ (N + 1) (N - 2 n) 2 (N - n) {\ displaystyle {\ text {df}} \ cong {\ frac {(N + 1) (N-2n)} {2 (Nn)}}}{\displaystyle {\text{df}}\cong {\frac {(N+1)(N-2n)}{2(N-n)}}}
фазовая модуляция мерцания (FPM)df ≅ exp ⁡ [(ln ⁡ N - 1 2 n ln ⁡ (2 n + 1) (N - 1) 4) - 1/2] {\ displaystyle {\ текст {df}} \ cong \ exp \ left [\ left (\ ln {\ frac {N-1} {2n}} \ ln {\ frac {(2n + 1) (N-1)} {4}} \ right) ^ {- 1/2} \ right]}{\ displaystyle {\ text {df}} \ cong \ exp \ left [\ left (\ ln {\ frac {N-1} {2n}) } \ ln {\ frac {(2n + 1) (N-1)} {4}} \ right) ^ {- 1/2} \ right]}
частотная модуляция белого (WFM)df ≅ [3 (N - 1) 2 n - 2 (N - 2) N] 4 n 2 4 п 2 + 5 {\ displaystyle {\ text {df}} \ cong \ left [{\ frac {3 (N-1)} {2n}} - {\ frac {2 (N-2)} {N}} \ right] {\ frac {4n ^ {2}} {4n ^ {2} +5}}}{\displaystyle {\text{df}}\cong \left[{\frac {3(N-1)}{2n}}-{\frac {2(N-2)}{N}}\right]{\frac {4n^{2}}{4n^{2}+5}}}
частотная модуляция мерцания (FFM)df ≅ {2 (N - 2) 2.3 N - 4.9 n Знак равно 1 5 N 2 4 n (N + 3 n) n ≥ 2 {\ 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}}}{\ 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)df ≅ N - 2 n (N - 1) 2 - 3 N (N - 1) + 4 N 2 (N - 3) 2 {\ displaystyle {\ text {df}} \ cong {\ frac {N-2} {n}} {\ frac {(N-1) ^ {2} -3n (N-1) + 4n ^ {2}} {(N-3) ^ {2}}}}{\ displaystyle {\ text {df}} \ cong {\ frac {N-2} {n}} {\ frac {(N-1) ^ {2} - 3n (N-1) + 4n ^ {2}} {(N-3) ^ {2}}}}
Степенной шум

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

характеристика степенного шума дисперсии Аллана
Тип степенного шумаНаклон фазового шумаЧастотный шум наклонКоэффициент мощностиФазовый шум. S x (f) {\ displaystyle S_ {x} (f)}{\ displaystyle S_ {x} (f)} Дисперсия Аллана. σ Y 2 (τ) {\ displaystyle \ sigma _ {y} ^ {2} (\ tau)}{\displaystyle \sigma _{y}^{2}(\tau)}отклонение Аллана. σ y (τ) {\ displaystyle \ sigma _ {y} (\ tau)}{\ displaystyle \ sigma _ {y} (\ tau)}
модуляция фазы белого (WPM)f 0 = 1 {\ displaystyle f ^ {0} = 1}f^{0}=1f 2 {\ displaystyle f ^ {2}}f^{2}h 2 {\ displaystyle h_ {2}}h_{2}1 (2 π) 2 час 2 {\ displaystyle {\ frac {1} {(2 \ pi) ^ {2}}} h_ {2}}{\displaystyle {\frac {1}{(2\pi)^{2}}}h_{2}}3 f H 4 π 2 τ 2 час 2 {\ displaystyle {\ frac {3f_ {H}} {4 \ pi ^ {2} \ tau ^ {2}}} h_ {2}}{\displaystyle {\frac {3f_{H}}{4\pi ^{2}\tau ^{2}}}h_{2}}3 f H 2 π τ h 2 {\ displaystyle {\ frac {\ sqrt {3f_ {H}}} {2 \ pi \ tau}} {\ sqrt {h_ {2}}}}{\ displaystyle {\ frac {\ sqrt {3f_ {H}}} {2 \ pi \ tau}} {\ sqrt {h_ {2}}}}
фазовая модуляция мерцания (FPM)f - 1 {\ displaystyle f ^ { - 1}}е ^ {- 1} f 1 = f {\ displaystyl ef ^ {1} = f}f ^ {1} = f час 1 {\ displaystyle h_ {1}}h_ {1} 1 (2 π) 2 fh 1 {\ displayst yle {\ frac {1} {(2 \ pi) ^ { 2} f}} h_ {1}}{\ displaystyle {\ frac {1} {(2 \ pi) ^ {2} f}} h_ {1}} 3 [γ + ln ⁡ (2 π f H τ)] - ln ⁡ 2 4 π 2 τ 2 час 1 {\ displaystyle {\ frac {3 [\ gamma + \ ln (2 \ pi f_ {H} \ tau)] - \ ln 2} {4 \ pi ^ {2} \ tau ^ {2}}} час_ {1}}{\displaystyle {\frac {3[\gamma +\ln(2\pi f_{H}\tau)]-\ln 2}{4\pi ^{2}\tau ^{2}}}h_{1}}3 [γ + пер ⁡ ( 2 π е час τ)] - пер ⁡ 2 2 π τ час 1 {\ displaystyle {\ frac {\ sqrt {3 [\ gamma + \ ln (2 \ pi f_ {H} \ tau)] - \ ln 2} } {2 \ pi \ tau}} {\ sqrt {h_ {1}}}}{\displaystyle {\frac {\sqrt {3[\gamma +\ln(2\pi f_{H}\tau)]-\ln 2}}{2\pi \tau }}{\sqrt {h_{1}}}}
частотная модуляция белого (WFM)f - 2 {\ displaystyle f ^ {- 2}}f^{{-2}}f 0 = 1 {\ displaystyle f ^ {0} = 1}f^{0}=1час 0 {\ displaystyle h_ {0}}h_{0}1 (2 π) 2 е 2 час 0 {\ displaystyle {\ frac {1} {( 2 \ pi) ^ {2} f ^ {2}}} h_ {0}}{\displaystyle {\frac {1}{(2\pi)^{2}f^{2}}}h_{0}}1 2 τ час 0 {\ displaystyle {\ frac {1} {2 \ tau}} h_ {0}}{\ displaystyle {\ frac {1} {2 \ tau}} h_ {0}} 1 2 τ h 0 {\ displaystyle {\ frac {1} {\ sqrt {2 \ tau}}} {\ sqrt {h_ {0}}}}{\displaystyle {\frac {1}{\sqrt {2\tau }}}{\sqrt {h_{0}}}}
частотная модуляция мерцания (FFM)f - 3 {\ displaystyle f ^ {- 3}}е ^ {{- 3}} f - 1 {\ displaystyle f ^ {- 1}}е ^ {- 1} h - 1 {\ displaystyle h _ {- 1}}ч _ {{- 1}} 1 (2 π) 2 е 3 час - 1 {\ displaystyle {\ frac {1} {(2 \ pi) ^ {2} f ^ {3}}} h _ {- 1 }}{\displaystyle {\frac {1}{(2\pi)^{2}f^{3}}}h_{-1}}2 пер (2) час - 1 {\ displaystyle 2 \ ln (2) час _ {- 1}}{\displaystyle 2\ln(2)h_{-1}}2 пер ⁡ (2) час - 1 {\ displaystyle {\ sqrt {2 \ ln (2)}} {\ sqrt {h _ {- 1}}}}{\ displaystyle {\ sqrt {2 \ ln (2)}} {\ sqrt {h _ {- 1}}}}
частотная модуляция случайного блуждания (RWFM)f - 4 {\ displaystyle f ^ {- 4}}f^{{-4}}е - 2 {\ displaystyle f ^ {- 2}}f^{{-2}}час - 2 {\ displaystyle h _ {- 2}}h_{{-2}}1 (2 π) 2 f 4 час - 2 {\ displaystyle {\ frac {1} {(2 \ pi) ^ {2} f ^ {4}}} h _ {- 2}}{\ displaystyle {\ frac {1} {(2 \ pi) ^ {2} f ^ {4}}} h_ {-2}} 2 π 2 τ 3 h - 2 {\ displaystyle {\ frac {2 \ pi ^ {2} \ tau } {3}} час _ {- 2}}{\displaystyle {\frac {2\pi ^{2}\tau }{3}}h_{-2}}π 2 τ 3 час - 2 {\ displaystyle {\ frac {\ pi {\ sqrt {2 \ tau}}} {\ sqrt {3}}} {\ sqrt {h _ {- 2}}}}{\displaystyle {\frac {\pi {\sqrt {2\tau }}}{\sqrt {3}}}{\sqrt {h_{-2}}}}

Как в современных формах.

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

В приведенных выше формулах предполагается, что

τ ≫ 1 2 π f H, {\ displaystyle \ tau \ gg {\ frac {1} {2 \ pi f_ {H}}},}{\displa ystyl e \tau \gg {\frac {1}{2\pi f_{H}}},}

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

отображение α – μ

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

S x (f) = 1 4 π 2 h α f α - 2 = 1 4 π 2 час α е β, {\ Displaystyle S_ {x} (f) = {\ гидроразрыва {1} {4 \ pi ^ {2}}} h _ {\ alpha} f ^ {\ alpha -2} = {\ frac { 1} {4 \ pi ^ {2}}} h _ {\ alpha} f ^ {\ beta},}{\displaystyle S_{x}(f)={\frac {1}{4\pi ^{2}}}h_{\alpha }f^{\alpha -2}={\frac {1}{4\pi ^{2}}}h_{\alpha }f^{\beta },}

где

β ≡ α - 2, {\ displaystyle \ beta \ Equiv \ alpha -2, }{\displaystyle \beta \equiv \alpha -2,}

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

S y (f) = h α f α {\ displaystyle S_ {y} (f) = h _ {\ alpha} f ^ {\ alpha}}{\ displaystyle S_ {y} (f) = h _ {\ alpha} f ^ { \ альфа}}

в дисперсия Аллана в форме

σ y 2 (τ) = К α час α τ μ {\ displaystyle \ sigma _ {y} ^ {2} (\ tau) = K _ {\ alpha} h _ {\ alpha} \ tau ^ {\ mu}}{\displaystyle \sigma _{y}^{2}(\tau)=K_{\alpha }h_{\alpha }\tau ^{\mu }}

можно значительно упростить, предоставив отображение между α и μ. Отображение между α и K α также представлено для удобства:

отображение дисперсии Аллана α – μ
αβμ
−2−412 π 2 3 {\ displaystyle {\ frac {2 \ pi ^ {2}} {3}}}{\frac {2\pi ^{2}}{3}}
−1−302 ln ⁡ 2 {\ displaystyle 2 \ ln {2}}2\ln {2}
0−2−11 2 {\ displaystyle {\ frac {1} {2}}}{\ frac {1} {2 }}
1−1−23 [γ + пер ⁡ (2 π е час τ)] - пер ⁡ 2 4 π 2 {\ displaystyle {\ frac {3 [\ gamma + \ ln (2 \ pi f_ {H} \ тау)] - \ ln 2} {4 \ pi ^ {2}}}}{\frac {3[\gamma +\ln(2\pi f_{H}\tau)]-\ln 2}{4\pi ^{2}}}
20−23 f H 4 π 2 {\ displaystyle {\ frac {3f_ {H}} {4 \ pi ^ {2}}}}{\ frac {3f_ {H}} {4 \ pi ^ {2}}}

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

Сигнал со спектральным фазовым шумом S ϕ {\ displaystyle S _ {\ phi}}S_{\phi }с единицами измерения рад / Гц может быть преобразован в дисперсию Аллана с точностью

σ y 2 (τ) = 2 ν 0 2 ∫ 0 fb S ϕ (f) sin 4 ⁡ (π τ f) (π τ) 2 df. {\ 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.}{\ displaystyle \ sigma _ {y} ^ {2} (\ tau) = {\ frac {2} {\ nu _ {0} ^ {2}}} \ int _ {0} ^ {f_ {b}} S _ {\ phi} (f) {\ гидроразрыва {\ грех ^ {4} (\ пи \ тау f)} {(\ пи \ тау) ^ {2}}} \, df.}
Линейный отклик

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

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

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

Свойства частотно-временного фильтра

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

σ y 2 (τ) = 1 2 ⟨(y ¯ i + 1 - y ¯ i) 2⟩, {\ displaystyle \ sigma _ {y} ^ {2} (\ tau) = {\ frac {1} {2}} \ langle ({\ bar {y}} _ {i + 1} - {\ bar {y}} _ {i}) ^ {2} \ rangle,}{\displaystyle \sigma _{y}^{2}(\tau)={\frac {1}{2}}\langle ({\bar {y}}_{i+1}-{\bar {y}}_{i})^{2}\rangle,}

где

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

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

σ y 2 (τ) = ∫ 0 ∞ S y (f) 2 sin 4 ⁡ π τ f (π τ f) 2 дф. {\ displaystyle \ sigma _ {y} ^ {2} (\ tau) = \ int _ {0} ^ {\ infty} S_ {y} (f) {\ frac {2 \ sin ^ {4} \ pi \ tau f} {(\ pi \ tau f) ^ {2}}} \, df.}{\displaystyle \sigma _{y}^{2}(\tau)=\int _{0}^{\infty }S_{y}(f){\frac {2\sin ^{4}\pi \tau f}{(\pi \tau f)^{2}}}\,df.}

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

| H A (f) | 2 = 2 sin 4 ⁡ π τ f (π τ f) 2. {\ displaystyle \ left \ vert H_ {A} (f) \ right \ vert ^ {2} = {\ frac {2 \ sin ^ {4} \ pi \ tau f} {(\ pi \ tau f) ^ { 2}}}.}{\displaystyle \left\vert H_{A}(f)\right\vert ^{2}={\frac {2\sin ^{4}\pi \tau f}{(\pi \tau f)^{2}}}.}
Функции смещения

Дисперсия M-выборки и определенная дисперсия Аллана для особого случая будут испытывать систематическое смещение в зависимости от разного количества выборок M и различного отношения между T и τ. Чтобы устранить эти функции, были реализованы функции с ущербом B 1 B 2, позволяют преобразовывать между различными значениями M и T.

Этих функции ущерба недостаточно для обработки, развивающего в результате объединения M со временем наблюдения Mτ 0 в MT 0 с мертвым временем, распределенным блоком измерения M, а не в конце измерения. Это привело к необходимости причинения ущерба B 3.

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

B1функция с ущербом

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

В 1 (N, r, μ) = ⟨σ y 2 (N, T, τ)⟩ ⟨σ y 2 (2, T, τ)⟩, {\ 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 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 }},}

где

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

После анализа с ущербом

B 1 (N, r, μ) = 1 + ∑ n = 1 N - 1 N - n N (N - 1) [2 (rn) μ + 2 - (rn + 1) μ + 2 - | r n - 1 | μ + 2] 1 + 1 2 [2 r μ + 2 - (r + 1) μ + 2 - | г - 1 | μ + 2]. {\ displaystyle B_ {1} (N, r, \ mu) = {\ frac {1+ \ sum _ {n = 1} ^ {N-1} {\ frac {Nn} {N (N-1)} } \ left [2 (rn) ^ {\ mu +2} - (rn + 1) ^ {\ mu +2} - | рН-1 | ^ {\ mu +2} \ right]} {1 + {\ гидроразрыв {1} {2}} \ left [2r ^ {\ mu +2} - (r + 1) ^ {\ mu +2} - | р-1 | ^ {\ mu +2} \ right]}}.}{\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]}}.}

B2функция с ущербом

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

B 2 (r, μ) = ⟨σ y 2 (2, T, τ)⟩ ⟨σ y 2 (2, τ, τ)⟩, {\ 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 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 }},}

где

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

После анализа с ущербом

B 2 (r, μ) = 1 + 1 2 [2 r μ + 2 - (r + 1) μ + 2 - | г - 1 | μ + 2] 2 (1-2 μ). {\ displaystyle B_ {2} (r, \ mu) = {\ frac {1 + {\ frac {1} {2}} \ left [2r ^ {\ mu +2} - (r + 1) ^ {\ mu +2} - | р-1 | ^ {\ mu +2} \ right]} {2 (1-2 ^ {\ mu})}}.}{\ displaystyle B_ {2} (r, \ mu) = {\ frac {1 + {\ f rac {1} {2}} \ left [2r ^ {\ mu +2} - (r + 1) ^ {\ mu +2} - | р-1 | ^ {\ mu +2} \ right]} { 2 (1-2 ^ {\ mu})}}.}

B3функция с ущербом

B 3 функция с ущерба связывает двухвыборочную дисперсию для времени выборки MT 0 и времени наблюдения Mτ 0 с двухвыборочной дисперсией (дисперсией Аллана) и определяется как

B 3 (N, M, r, μ) = ⟨σ Y 2 (N, M, T, τ)⟩ ⟨σ Y 2 (N, T, τ)⟩, {\ 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 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} },}

где

T = MT 0, {\ displaystyle T = MT_ {0},}{\ displaystyle T = MT_ {0 },}
τ = M τ 0. {\ displaystyle \ tau = M \ tau _ { 0}.}{\displaystyle \tau =M\tau _{0}.}

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

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

B 3 (2, M, r, μ) = 2 M + MF (M r) - ∑ n = 1 M - 1 ( M - n) [2 F (nr) - F ((M + n) r) + F ((M - n) r)] M μ + 2 [F (r) + 2], {\ displaystyle B_ {3 } (2, M, r, \ mu) = {\ frac {2M + MF (Mr) - \ sum _ {n = 1} ^ {M-1} (Mn) \ left [2F (nr) - F { \ big (} (M + n) r {\ big)} + F {\ big (} (Mn) r {\ big)} \ right]} {M ^ {\ mu +2} [F (r) + 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]}},}

где

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

Функция с ущербом τ

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

B τ (τ 1, τ 2, μ) = ⟨σ y 2 ( 2, τ 2, τ 2)⟩ ⟨σ y 2 (2, τ 1, τ 1)⟩. {\ 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)={\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 }}.}

После анализа функция с повреждением становится

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

Преобразование между значениями

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

⟨σ y 2 (N 2, T 2, τ 2)⟩ = (τ 2 τ 1) μ [B 1 (N 2, r 2, μ) B 2 (r 2, μ) В 1 (N 1, р 1, μ) В 2 (r 1, μ)] ⟨σ Y 2 (N 1, T 1, τ 1)⟩, {\ Displaystyle \ влево \ langle \ sigma _ {y} ^ { 2} (N_ {2}, T_ {2}, \ tau _ {2}) \ right \ rangle = \ left ({\ frac {\ tau _ {2}} {\ tau _ {1}}} \ справа) ^ {\ 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 \ 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,}

где

r 1 = T 1 r 1, {\ displaystyle r_ {1} = {\ frac {T_ {1 }} {r_ {1}}},}{\ displaystyle r_ {1} = {\ frac {T_ {1}} {r_ {1}}},}
r 2 = T 2 r 2. {\ displaystyle r_ {2} = {\ frac {T_ {2}} {r_ {2}}}.}{\ displaystyle r_ {2} = {\ frac {T_ {2}} {r_ {2}}}.}

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

⟨σ y 2 (N 2, M 2, T 2, τ 2)⟩ = (τ 2 τ 1) μ [B 3 (N 2, M 2, r 2, μ) B 1 (N 2, r 2, μ) B 2 (r 2, μ) B 3 (N 1, M 1, r 1, μ) B 1 (N 1, r 1, μ) B 2 (r 1, μ)] ⟨σ y 2 (N 1, M 1, T 1, τ 1)⟩. {\ 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.}{\ displaystyle \ left \ langle \ sigma _ {y} ^ {2} (N_ {2}, M_ {2}, T_ {2}, \ tau _ {2}) \ right \ rangle = \ left ({\ frac {\ tau _ {2}} {\ tau _ {1}}} \ справа) ^ {\ 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.}
Проблемы измерения

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

Пределы пропускания измерения

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

τ ≫ 1 2 π f H. {\ displaystyle \ tau \ gg {\ frac {1} {2 \ pi f_ {H}}}.}\tau \gg {\frac {1}{2\pi f_{H}}}.

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

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

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

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

Мертвое время в измерениях

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

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

Измерения, выполняемые с мертвым временем, можно скорректировать с помощью функций с нарушением B 1, B 2 и B 3. Таким образом, мертвое время как таковое не запрещает доступ к дисперсии Аллана, но делает его более проблематичным. Мертвое время должно быть известно, чтобы можно было установить время между выборками T.

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

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

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

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

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

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

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

Тип доминирующего шума

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

Линейный дрейф

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

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

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

Практические измерения

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

Измерение

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

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

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

Постобработка

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

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

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

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

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

История исследований

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

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

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

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

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

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

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

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

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

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

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

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

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

Классическим справочником является монография NBS № 140 от 1974 г., в главе 8 есть «Статистика анализа временных и частотных данных». Это расширенный вариант NBS Technical Note 394, который используется в методах измерения и практическую обработку значений.

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

Стандартные определения физических величин для основного метрологии частоты и времени в стандарте IEEE 1139 выходят за рамки стандартного справочного и образовательного ресурса.

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

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

Использует

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

50-летие

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

См. Также
Ссылки
Внешние ссылки
Последняя правка сделана 2021-06-11 00:38:31
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте