Ir al contenido

Diferencia entre revisiones de «Varianza de Allan»

De Wikipedia, la enciclopedia libre
Contenido eliminado Contenido añadido
Sin resumen de edición
Eloygd (discusión · contribs.)
Traducción aún sin concluir
Línea 172: Línea 172:


=== Estimadores <math>\tau</math> de variables superpuestas ===
=== Estimadores <math>\tau</math> de variables superpuestas ===
Una técnica presentada por J. J. Snyder proporcionaba una herramienta mejorada, ya que las mediciones se superponían en n series superpuestas de la serie original. Howe, Allan y Barnes introdujeron el estimador de varianza de Allan superpuesto, que puede demostrarse que es equivalente a promediar las muestras de tiempo o de frecuencia normalizada en bloques de n muestras antes del procesamiento. El predictor resultante es
Una técnica presentada por J. J. Snyder<ref>{{Cita publicación|url=http://dx.doi.org/10.1109/freq.1981.200513|título=An Ultra-High Resolution Frequency Meter|apellidos=Snyder|nombre=J.J.|fecha=1981|publicación=Thirty Fifth Annual Frequency Control Symposium|editorial=IEEE|fechaacceso=2022-02-25|doi=10.1109/freq.1981.200513}}</ref> proporcionaba una herramienta mejorada, ya que las mediciones se superponían en n series superpuestas de la serie original. Howe, Allan y Barnes<ref name=":0">{{Cita publicación|url=http://dx.doi.org/10.1109/freq.1981.200541|título=Properties of Signal Sources and Measurement Methods|apellidos=Howe|nombre=D.A.|apellidos2=Allan|nombre2=D.U.|fecha=1981|publicación=Thirty Fifth Annual Frequency Control Symposium|editorial=IEEE|fechaacceso=2022-02-25|doi=10.1109/freq.1981.200541|apellidos3=Barnes|nombre3=J.A.}}</ref> introdujeron el estimador de varianza de Allan superpuesto, que puede demostrarse que es equivalente a promediar las muestras de tiempo o de frecuencia normalizada en bloques de n muestras antes del procesamiento. El predictor resultante es
: <math>\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, </math>
: <math>\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, </math>
o, para la serie temporal:
o, para la serie temporal:
: <math>\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.</math>
: <math>\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.</math>
Los estimadores solapados tienen un rendimiento muy superior al de los estimadores no solapados a medida que n aumenta y la serie temporal es de longitud moderada. Los estimadores solapados han sido aceptados como los estimadores de varianza de Allan preferidos en las normas IEEE, UIT-T y ETSI para mediciones comparables, como las necesarias para la calificación de las telecomunicaciones.
Los estimadores solapados tienen un rendimiento muy superior al de los estimadores no solapados a medida que n aumenta y la serie temporal es de longitud moderada. Los estimadores solapados han sido aceptados como los estimadores de varianza de Allan preferidos en las normas IEEE<ref name="IEEE11392" />, UIT-T<ref>{{Obra citada|título=DIN EN 61883-7:2003-08, Audio/Video-Geräte der Unterhaltungselektronik_- Digitale Schnittstelle_- Teil_7: Übertragung nach Rec. ITU-R_BO.1294 System_B (IEC_61883-7:2003); Deutsche Fassung EN_61883-7:2003|url=http://dx.doi.org/10.31030/9499058|fechaacceso=2022-02-25|editorial=Beuth Verlag GmbH}}</ref> y ETSI<ref>{{Cita libro|título=ETSI ITS-G5 Standard|url=http://dx.doi.org/10.1007/978-3-319-78262-1_300222|editorial=Springer International Publishing|fecha=2020|fechaacceso=2022-02-25|páginas=442–442}}</ref> para mediciones comparables, como las necesarias para la calificación de las telecomunicaciones.


=== Varianza de Allan modificada ===
=== Varianza de Allan modificada ===
Para solucionar la incapacidad de separar la modulación de fase blanca de la modulación de fase de parpadeo utilizando los estimadores de varianza de Allan tradicionales, un filtrado algorítmico reduce el ancho de banda en n. Este filtrado proporciona una modificación a la definición y a los estimadores y ahora se identifica como una clase de varianza separada llamada varianza Allan modificada. La medida de la varianza de Allan modificada es una medida de estabilidad de la frecuencia, al igual que la varianza de Allan.


=== Estimadores de estabilidad temporal ===
=== Estimadores de estabilidad temporal ===
A partir de la desviación Allan modificada (MDEV) puede calcularse una medida estadística de estabilidad temporal (<math>\sigma_x</math>), que suele denominarse desviación temporal (TDEV). La TDEV se basa en la MDEV en lugar de la desviación de Allan original, porque la MDEV puede discriminar entre la modulación de fase (PM) blanca y la de parpadeo. A continuación se presenta la estimación de la varianza temporal basada en la varianza de Allan modificada:
: <math>\sigma_x^2(\tau) = \frac{\tau^2}{3}\bmod\sigma_y^2(\tau),</math>
: <math>\sigma_x^2(\tau) = \frac{\tau^2}{3}\bmod\sigma_y^2(\tau),</math>
y de forma similar para la desviación modificada de Allan a la desviación temporal:

: <math>\sigma_x(\tau) = \frac{\tau}{\sqrt{3}}\bmod\sigma_y(\tau).</math>
: <math>\sigma_x(\tau) = \frac{\tau}{\sqrt{3}}\bmod\sigma_y(\tau).</math>
La TDEV se normaliza de forma que sea igual a la desviación clásica para el PM blanco para la constante de tiempo <math>\tau=\tau_0</math>. Para entender el factor de escala de normalización entre las medidas estadísticas, la siguiente es la regla estadística pertinente: Para variables aleatorias independientes X e Y, la varianza (<math>\sigma_z^2</math>) de una suma o diferencia (z = x - y) es la suma cuadrada de sus varianzas (<math>\sigma_z^2=\sigma_x^2+\sigma_y^2</math>). La varianza de la suma o diferencia (y = x<sub>2τ</sub> - x<sub>τ</sub>) de dos muestras independientes de una variable aleatoria es el doble de la varianza de la variable aleatoria (<math>\sigma_y^2=2\sigma_x^2</math>). El MDEV es la segunda diferencia de las medidas de fase independientes (x) que tienen una varianza (<math>\sigma_x^2</math>). Como el cálculo es la doble diferencia, que requiere tres mediciones de fase independientes (x<sub>2τ</sub> - 2x<sub>τ</sub> + x), la varianza de Allan modificada (MVAR) es tres veces las varianzas de las mediciones de fase.


=== Otros estimadores ===
=== Otros estimadores ===
Otros desarrollos han producido métodos de estimación mejorados para la misma medida de estabilidad, la varianza/desviación de la frecuencia, pero se conocen con nombres distintos como la [[varianza de Hadamard]], la [[varianza de Hadamard modificada]], la [[varianza total]], la [[varianza total modificada]] y la [[varianza de Theo]]. Estas se distinguen por un mejor uso de la estadística para mejorar los límites de confianza o la capacidad de manejar la deriva lineal de la frecuencia.


== Intervalos de confianza y grados de libertad equivalentes ==
== Intervalos de confianza y grados de libertad equivalentes ==
Los estimadores estadísticos calcularán un valor estimado sobre la serie de muestras utilizada. Las estimaciones pueden desviarse del valor real y el rango de valores que, con cierta probabilidad, contendrá el valor real se denomina [[intervalo de confianza]]. El intervalo de confianza depende del número de observaciones de la serie muestral, del tipo de ruido dominante y del estimador utilizado. La anchura también depende de la certeza estadística para la que los valores del intervalo de confianza forman un rango acotado, por lo que la certeza estadística de que el valor verdadero está dentro de ese rango de valores. En el caso de los estimadores variables-τ, el múltiplo n de τ<sub>0</sub> también es una variable.


=== Intervalo de confianza ===
=== Intervalo de confianza ===
El intervalo de confianza puede establecerse mediante la [[Distribución χ²|distribución chi-cuadrado]] utilizando la distribución de la varianza de la muestra:<ref name="IEEE11392" /><ref name=":0" />
: <math>\chi^2 = \frac{\text{df}\,s^2}{\sigma^2},</math>
: <math>\chi^2 = \frac{\text{df}\,s^2}{\sigma^2},</math>
donde s<sup>2</sup> es la varianza muestral de nuestra estimación, σ<sup>2</sup> es el valor de la varianza verdadera, df son los grados de libertad para el estimador, y χ<sup>2</sup> son los grados de libertad para una determinada probabilidad. Para una probabilidad del 90 %, que cubre el rango del 5 % al 95 % en la curva de probabilidad, los límites superior e inferior se pueden encontrar utilizando la desigualdad

: <math>\chi^2(0.05) \le \frac{\text{df}\,s^2}{\sigma^2} \le \chi^2(0.95),</math>
: <math>\chi^2(0.05) \le \frac{\text{df}\,s^2}{\sigma^2} \le \chi^2(0.95),</math>
que después de reordenar para la varianza verdadera se convierte en

: <math>\frac{\text{df}\,s^2}{\chi^2(0.95)} \le \sigma^2 \le \frac{\text{df}\,s^2}{\chi^2(0.05)}.</math>
: <math>\frac{\text{df}\,s^2}{\chi^2(0.95)} \le \sigma^2 \le \frac{\text{df}\,s^2}{\chi^2(0.05)}.</math>


=== Grados de libertad efectivos ===
=== Grados de libertad efectivos ===
Los grados de libertad representan el número de variables libres capaces de contribuir a la estimación. Dependiendo del estimador y del tipo de ruido, los grados de libertad efectivos varían. Se han encontrado empíricamente fórmulas de estimadores que dependen de N y n:<ref name=":0" />

{| class="wikitable"
{| class="wikitable"
|+ Grados de libertad de la varianza de Allan
|-
!Tipo de ruido
!Tipo de ruido
!Grados de libertad
!Grados de libertad
|-
|-
|modulación en fase blanco (WPM)
|modulación de fase blanca (WPM)
|<math>\text{df} \cong \frac{(N + 1)(N - 2n)}{2(N - n)}</math>
|<math>\text{df} \cong \frac{(N + 1)(N - 2n)}{2(N - n)}</math>
|-
|-
|modulación en fase flicker (FPM)
|modulación de fase de parpadeo (FPM)
|<math>\text{df} \cong \exp\left[\left(\ln \frac{N - 1}{2n} \ln \frac{(2n + 1)(N - 1)}{4}\right)^{-1/2}\right]</math>
|<math>\text{df} \cong \exp\left[\left(\ln \frac{N - 1}{2n} \ln \frac{(2n + 1)(N - 1)}{4}\right)^{-1/2}\right]</math>
|-
|-
|modulación en freciencia blanco (WFM)
|modulación de frecuencia blanca (WFM)
|<math>\text{df} \cong \left[ \frac{3(N - 1)}{2n} - \frac{2(N - 2)}{N}\right] \frac{4n^2}{4n^2 + 5}</math>
|<math>\text{df} \cong \left[ \frac{3(N - 1)}{2n} - \frac{2(N - 2)}{N}\right] \frac{4n^2}{4n^2 + 5}</math>
|-
|-
|modulación en frecuencia flicker (FFM)
|modulación de frecuencia de parpadeo (FFM)
|<math>\text{df} \cong \begin{cases}\frac{2(N - 2)}{2.3N - 4.9} & n = 1 \\ \frac{5N^2}{4n(N + 3n)} & n \ge 2\end{cases}</math>
|<math>\text{df} \cong \begin{cases}\frac{2(N - 2)}{2.3N - 4.9} & n = 1 \\ \frac{5N^2}{4n(N + 3n)} & n \ge 2\end{cases}</math>
|-
|-
|modulación en freciencia aleatoria (RWFM)
|modulación de frecuencia aleatoria (RWFM)
|<math>\text{df} \cong \frac{N - 2}{n}\frac{(N - 1)^2 - 3n(N - 1) + 4n^2}{(N - 3)^2}</math>
|<math>\text{df} \cong \frac{N - 2}{n}\frac{(N - 1)^2 - 3n(N - 1) + 4n^2}{(N - 3)^2}</math>
|}
|}
Línea 219: Línea 229:
== Ruido de potencia ==
== Ruido de potencia ==
{| class="wikitable"
{| class="wikitable"
|+Respuesta de ley de potencia de la varianza de Allan
|+Allan variance power-law response
!Tipo de ruido de potencia
!Power-law noise type
!Pendiente de ruido de fase
!Phase noise slope
!Pendiente de ruido en frecuencia
!Frequency noise slope
!Coeficiente de potencia
!Power coefficient
!Phase noise<br /><br /> <math>S_x(f)</math>
!Ruido de fase<br /><br /> <math>S_x(f)</math>
!Allan variance<br /><br /> <math>\sigma_y^2(\tau)</math>
!Varianza de Allan<br /><br /> <math>\sigma_y^2(\tau)</math>
!Allan deviation<br /><br /> <math>\sigma_y(\tau)</math>
!Desviación de Allan<br /><br /> <math>\sigma_y(\tau)</math>
|-
|-
|white phase modulation (WPM)
|white phase modulation (WPM)

Revisión del 12:01 25 feb 2022

La forma más fácil de examinar un reloj es compararlo con un reloj de referencia mucho más preciso. Durante un intervalo de tiempo τ, medido por el reloj de referencia, el reloj bajo prueba avanza en τy, donde y es la frecuencia media (relativa) del reloj durante ese intervalo. Si medimos dos intervalos consecutivos como se muestra, podemos obtener un valor de (y - y′)2 - un valor menor indica un reloj más estable y preciso. Si repetimos este procedimiento muchas veces, el valor medio de (y - y′)2 es igual al doble de la varianza de Allan (o desviación de Allan al cuadrado) para el tiempo de observación τ.

La varianza de Allan (abreviadamente en inglés, AVAR,), también conocida como varianza de dos muestras, es una medida de estabilidad de la frecuencia en relojes, osciladores y amplificadores. Lleva el nombre de David W. Allan y se expresa matemáticamente como .

La desviación de Allan (abreviado en inglés, ADEV), también conocida como sigma-tau, es la raíz cuadrada de la varianza de Allan, .

La varianza de M muestras es una medida de la estabilidad de la frecuencia utilizando M muestras, el tiempo T entre las mediciones y el tiempo de observación . La varianza de M muestras se expresa como

La varianza de Allan tiene por objeto estimar la estabilidad debida a los procesos de ruido y no la de los errores sistemáticos o las imperfecciones, como la deriva de la frecuencia o los efectos de la temperatura. La varianza Allan y la desviación Allan describen la estabilidad de la frecuencia. Véase también el apartado Interpretación del valor más abajo.

También existen diferentes adaptaciones o alteraciones de la varianza de Allan, especialmente la varianza de Allan modificada MAVAR o MVAR, la varianza total y la varianza de Hadamard. También existen variantes de la estabilidad temporal, como la desviación temporal (TDEV) o la varianza temporal (TVAR). La varianza de Allan y sus variantes han demostrado su utilidad fuera del ámbito de la medición del tiempo y son un conjunto de herramientas estadísticas mejoradas que se utilizan siempre que los procesos de ruido no son incondicionalmente estables, por lo que existe una derivada.

La varianza general de M muestras sigue siendo importante, ya que permite determinar el tiempo muerto (el tiempo después de cada evento durante el cual el sistema no puede registrar otro evento en las mediciones), y las funciones de sesgo permiten la conversión en valores de varianza de Allan. Sin embargo, para la mayoría de las aplicaciones, el caso especial de 2 muestras, o "varianza de Allan" con es el de mayor interés.

Gráfica ejemplo de la desviación de Allan de un reloj. Con un tiempo de observación τ muy corto, la desviación de Allan es alta debido al ruido. Con un tiempo τ más largo, disminuye porque el ruido se promedia. A un tiempo τ aún mayor, la desviación de Allan comienza a aumentar de nuevo, lo que sugiere que la frecuencia del reloj se desvía gradualmente debido a los cambios de temperatura, el envejecimiento de los componentes u otros factores similares. Las barras de error aumentan con τ simplemente porque se tarda mucho en obtener muchos puntos de datos para τ grandes.

Antecedentes

Al investigar la estabilidad de los osciladores de cristal y los relojes atómicos, se descubrió que no tenían un ruido de fase compuesto únicamente por ruido blanco, sino también por ruido de frecuencia de parpadeo. Estas formas de ruido se convierten en un reto para las herramientas estadísticas tradicionales, como la desviación estándar, ya que el estimador no converge. Por tanto, se dice que el ruido es divergente. Los primeros esfuerzos por analizar la estabilidad incluyeron tanto análisis teóricos como mediciones prácticas.[1][2]

Una importante consecuencia colateral de contar con estos tipos de ruido era que, al no coincidir los distintos métodos de medición entre sí, no se podía conseguir el aspecto clave de la repetibilidad de una medición. Esto limita la posibilidad de comparar fuentes y hacer especificaciones significativas para exigir a los proveedores. Esencialmente, todas las formas de uso científico y comercial se limitaron entonces a mediciones dedicadas, que se espera que capten la necesidad de esa aplicación.

Para abordar estos problemas David Allan introdujo la varianza de M muestras e (indirectamente) la varianza de dos muestras.[3]​ Aunque la varianza de dos muestras no permitía distinguir completamente todos los tipos de ruido, proporcionaba un medio para separar de forma significativa muchas formas de ruido para las series temporales de mediciones de fase o frecuencia entre dos o más osciladores. Allan proporcionó un método para convertir cualquier varianza de M muestras en cualquier varianza de N muestras a través de la varianza común de dos muestras, haciendo así comparables todas las varianzas de M muestras. El mecanismo de conversión también demostró que la varianza de M muestras no converge para valores grandes de M, lo que las hace menos útiles. El IEEE identificó posteriormente la varianza de 2 muestras como la medida preferida.[4]

Una de las primeras preocupaciones estaba relacionada con los instrumentos de medición de tiempo y frecuencia que tenían un tiempo muerto entre las mediciones. Estas series de mediciones no formaban una observación continua de la señal y, por tanto, introducían un sesgo sistemático en la medición. Se puso mucho cuidado en la estimación de estos sesgos. La introducción de los contadores de tiempo muerto cero eliminó esta necesidad, pero las herramientas de análisis de los sesgos han resultado útiles.

Otro aspecto que preocupaba al principio estaba relacionado con la forma en que el ancho de banda del instrumento de medición influiría en la medición, de manera que era necesario tomar nota. Más tarde se comprobó que cambiando algorítmicamente la observación sólo los valores bajos de se verían afectados, mientras que los valores más altos no lo harían. El cambio de se hace dejando que sea un múltiplo entero de la base de tiempo de medición :

La física de los osciladores de cristal fue analizada por D. B. Leeson,[2]​ y el resultado se conoce ahora como ecuación de Leeson. La retroalimentación en el oscilador hará que el ruido blanco y el ruido de parpadeo del amplificador de retroalimentación y del cristal se conviertan en los ruidos de color de frecuencia de ruido blanco y de frecuencia de ruido de parpadeo respectivamente. Estas formas de ruido tienen el efecto de que el estimador de varianza estándar no converge al procesar las muestras de error temporal. Esta mecánica de los osciladores de retroalimentación era desconocida cuando se iniciaron los trabajos sobre la estabilidad de los osciladores, pero fue presentada por Leeson al mismo tiempo que el conjunto de herramientas estadísticas fue puesto a disposición por David W. Allan. Para una presentación más completa del efecto Leeson, véase la literatura moderna sobre el ruido de fase.[5]

Interpretación del valor

La varianza de Allan se define como la mitad de la media temporal de los cuadrados de las diferencias entre las lecturas sucesivas de la desviación de la frecuencia muestreada a lo largo del periodo de muestreo. La varianza de Allan depende del período de tiempo utilizado entre las muestras, por lo tanto, es una función del período de muestreo, comúnmente denotado como , al igual que la distribución que se está midiendo, y se muestra como un gráfico en lugar de un solo número. Una varianza Allan baja es una característica de un reloj con buena estabilidad a lo largo del periodo medido.

La desviación de Allan se utiliza ampliamente para las gráficas (convencionalmente en formato log-log) y la presentación de los números. Se prefiere, ya que da la estabilidad de la amplitud relativa, lo que permite una fácil comparación con otras fuentes de error.

Una desviación de Allan de 1,3×10-9 en un tiempo de observación de 1 s (es decir, = 1 s) debe interpretarse como la existencia de una inestabilidad en la frecuencia entre dos observaciones separadas por 1 segundo con un valor medio cuadrático relativo (RMS) de 1,3×10-9. Para un reloj de 10 MHz, esto equivaldría a un movimiento RMS de 13 mHz. Si se necesita la estabilidad de fase de un oscilador, hay que consultar y utilizar las variantes de desviación temporal.

Se puede convertir la varianza de Allan y otras varianzas del dominio del tiempo en medidas del dominio de la frecuencia de la estabilidad del tiempo (fase) y de la frecuencia.[6]

Definiciones

Varianza de M muestras

La varianza de M muestras se define[3]​ (aquí en una forma de notación modernizada) como

donde es la lectura del reloj (en segundos) medida en el momento , o con la serie temporal de frecuencia media fraccionaria

donde es el número de muestras de frecuencia utilizadas en la varianza, es el tiempo entre cada muestra de frecuencia, y es la duración de cada estimación de frecuencia.

Un aspecto importante es que la varianza de muestras puede incluir el tiempo muerto dejando que el tiempo sea diferente al de .

Varianza de Allan

La varianza de Allan se define como

donde denota valor esperado. Esto se puede expresar convenientemente como

donde es el período de observación, e es la enésima frecuencia fraccionaria media a lo largo del tiempo de observación .

Las muestras se toman sin tiempo muerto entre ellas, lo que se consigue dejando

Desviación de Allan

Al igual que la desviación estándar y la varianza, la desviación de Allan se define como la raíz cuadrada de la varianza de Allan:

Definiciones de apoyo

Modelo de oscilador

Se supone que el oscilador analizado sigue el modelo básico de

Se supone que el oscilador tiene una frecuencia nominal de dada en ciclos por segundo (unidad del SI: hertz). La frecuencia angular nominal (en radianes por segundo) viene dada por

La fase total puede separarse en un componente perfectamente cíclico junto con un componente fluctuante :

Error temporal

La función de error temporal x(t) es la diferencia entre el tiempo nominal esperado y el tiempo normal real:

Para los valores medidos se define una serie de errores temporales TE(t) a partir de la función temporal de referencia Tref(t) como

Función de frecuencia

La función de frecuencia es la frecuencia a lo largo del tiempo, definida como

Frecuencia fraccionaria

La frecuencia fraccionaria y(t) es la diferencia normalizada entre la frecuencia y la frecuencia nominal :

Frecuencia fraccionaria media

La frecuencia fraccionaria media se define como

donde la media se toma sobre el tiempo de observación , y(t) es el error de frecuencia fraccional en el tiempo t, y es el tiempo de observación.

Como y(t) es la derivada de x(t), podemos reescribirla sin pérdida de generalidad como

Estimadores

Esta definición se basa en el valor estadístico esperado, que se integra en un tiempo infinito. La situación en el mundo real no permite este tipo de series temporales, en cuyo caso es necesario utilizar un estimador estadístico en su lugar. Se presentarán y discutirán varios estimadores diferentes.

Convenciones

  • El número de muestras de frecuencia en una serie de frecuencia fraccionaria se denota con M.
  • El número de muestras de error temporal en una serie de error temporal se denota con N. La relación entre el número de muestras de frecuencia fraccionaria y la serie de errores temporales se fija en la relación
  • Para las series de muestras de errores temporales, xi se refiere a la muestra i de la función de tiempo continuo x(t) dada por
    donde T es el tiempo entre mediciones. Para la varianza de Allan, el tiempo que se utiliza tiene T ajustado al tiempo de observación τ. La serie de muestras de errores temporales permite que N denote el número de muestras (x0 ... xN-1) en la serie. La convención tradicional lleva el índice desde 1 hasta N en lugar de ir desde 0 hasta "N-1"
  • Para las series muestrales de Frecuencia fraccionaria media, se refiere a la muestra i de la función continua de frecuencia fraccionaria media y(t) dada por
    lo que da
    Para el supuesto de que T sea τ se convierte en
    La serie de muestras de Frecuencias fraccionarias medias deja que M denote el número de muestras () en la serie. La convención tradicional lleva el índice desde 1 hasta M. Como abreviatura, la Frecuencia fraccionaria media suele escribirse sin la barra de media encima. Sin embargo, esto es formalmente incorrecto, ya que la frecuencia fraccionaria y la frecuencia fraccionaria media son dos funciones diferentes. Un instrumento de medición capaz de producir estimaciones de frecuencia sin tiempo muerto entregará en realidad una serie temporal de frecuencia media, que sólo necesita ser convertida en frecuencia fraccionaria media y puede entonces ser utilizada directamente.
  • El tiempo entre las mediciones se denota con T, que es la suma del tiempo de observación τ y el tiempo muerto.

Estimadores de fijo

Un primer estimador sencillo sería traducir directamente la definición en

o para las series temporales:

Sin embargo, estas fórmulas sólo proporcionan el cálculo para el caso . Para calcular para un valor diferente de , es necesario proporcionar una nueva serie temporal.

Estimadores de variables no solapadas

Si se toma la serie temporal y se salta más allá de n - 1 muestras, se obtendría una nueva serie temporal (más corta) con como tiempo entre las muestras adyacentes, para la que se podría calcular la varianza de Allan con los estimadores simples. Éstos podrían modificarse para introducir la nueva variable n, de modo que no habría que generar una nueva serie temporal, sino que podría reutilizarse la serie temporal original para diversos valores de n. Los estimadores se convierten en

con

y para la serie temporal:

con

Estos estimadores tienen un inconveniente importante, ya que dejan caer una cantidad significativa de datos de la muestra, ya que sólo se utiliza 1/n de las muestras disponibles.

Estimadores de variables superpuestas

Una técnica presentada por J. J. Snyder[7]​ proporcionaba una herramienta mejorada, ya que las mediciones se superponían en n series superpuestas de la serie original. Howe, Allan y Barnes[8]​ introdujeron el estimador de varianza de Allan superpuesto, que puede demostrarse que es equivalente a promediar las muestras de tiempo o de frecuencia normalizada en bloques de n muestras antes del procesamiento. El predictor resultante es

o, para la serie temporal:

Los estimadores solapados tienen un rendimiento muy superior al de los estimadores no solapados a medida que n aumenta y la serie temporal es de longitud moderada. Los estimadores solapados han sido aceptados como los estimadores de varianza de Allan preferidos en las normas IEEE[4]​, UIT-T[9]​ y ETSI[10]​ para mediciones comparables, como las necesarias para la calificación de las telecomunicaciones.

Varianza de Allan modificada

Para solucionar la incapacidad de separar la modulación de fase blanca de la modulación de fase de parpadeo utilizando los estimadores de varianza de Allan tradicionales, un filtrado algorítmico reduce el ancho de banda en n. Este filtrado proporciona una modificación a la definición y a los estimadores y ahora se identifica como una clase de varianza separada llamada varianza Allan modificada. La medida de la varianza de Allan modificada es una medida de estabilidad de la frecuencia, al igual que la varianza de Allan.

Estimadores de estabilidad temporal

A partir de la desviación Allan modificada (MDEV) puede calcularse una medida estadística de estabilidad temporal (), que suele denominarse desviación temporal (TDEV). La TDEV se basa en la MDEV en lugar de la desviación de Allan original, porque la MDEV puede discriminar entre la modulación de fase (PM) blanca y la de parpadeo. A continuación se presenta la estimación de la varianza temporal basada en la varianza de Allan modificada:

y de forma similar para la desviación modificada de Allan a la desviación temporal:

La TDEV se normaliza de forma que sea igual a la desviación clásica para el PM blanco para la constante de tiempo . Para entender el factor de escala de normalización entre las medidas estadísticas, la siguiente es la regla estadística pertinente: Para variables aleatorias independientes X e Y, la varianza () de una suma o diferencia (z = x - y) es la suma cuadrada de sus varianzas (). La varianza de la suma o diferencia (y = x - xτ) de dos muestras independientes de una variable aleatoria es el doble de la varianza de la variable aleatoria (). El MDEV es la segunda diferencia de las medidas de fase independientes (x) que tienen una varianza (). Como el cálculo es la doble diferencia, que requiere tres mediciones de fase independientes (x - 2xτ + x), la varianza de Allan modificada (MVAR) es tres veces las varianzas de las mediciones de fase.

Otros estimadores

Otros desarrollos han producido métodos de estimación mejorados para la misma medida de estabilidad, la varianza/desviación de la frecuencia, pero se conocen con nombres distintos como la varianza de Hadamard, la varianza de Hadamard modificada, la varianza total, la varianza total modificada y la varianza de Theo. Estas se distinguen por un mejor uso de la estadística para mejorar los límites de confianza o la capacidad de manejar la deriva lineal de la frecuencia.

Intervalos de confianza y grados de libertad equivalentes

Los estimadores estadísticos calcularán un valor estimado sobre la serie de muestras utilizada. Las estimaciones pueden desviarse del valor real y el rango de valores que, con cierta probabilidad, contendrá el valor real se denomina intervalo de confianza. El intervalo de confianza depende del número de observaciones de la serie muestral, del tipo de ruido dominante y del estimador utilizado. La anchura también depende de la certeza estadística para la que los valores del intervalo de confianza forman un rango acotado, por lo que la certeza estadística de que el valor verdadero está dentro de ese rango de valores. En el caso de los estimadores variables-τ, el múltiplo n de τ0 también es una variable.

Intervalo de confianza

El intervalo de confianza puede establecerse mediante la distribución chi-cuadrado utilizando la distribución de la varianza de la muestra:[4][8]

donde s2 es la varianza muestral de nuestra estimación, σ2 es el valor de la varianza verdadera, df son los grados de libertad para el estimador, y χ2 son los grados de libertad para una determinada probabilidad. Para una probabilidad del 90 %, que cubre el rango del 5 % al 95 % en la curva de probabilidad, los límites superior e inferior se pueden encontrar utilizando la desigualdad

que después de reordenar para la varianza verdadera se convierte en

Grados de libertad efectivos

Los grados de libertad representan el número de variables libres capaces de contribuir a la estimación. Dependiendo del estimador y del tipo de ruido, los grados de libertad efectivos varían. Se han encontrado empíricamente fórmulas de estimadores que dependen de N y n:[8]

Grados de libertad de la varianza de Allan
Tipo de ruido Grados de libertad
modulación de fase blanca (WPM)
modulación de fase de parpadeo (FPM)
modulación de frecuencia blanca (WFM)
modulación de frecuencia de parpadeo (FFM)
modulación de frecuencia aleatoria (RWFM)

Ruido de potencia

Respuesta de ley de potencia de la varianza de Allan
Tipo de ruido de potencia Pendiente de ruido de fase Pendiente de ruido en frecuencia Coeficiente de potencia Ruido de fase

Varianza de Allan

Desviación de Allan

white phase modulation (WPM)
flicker phase modulation (FPM)
white frequency modulation (WFM)
flicker frequency modulation (FFM)
random walk frequency modulation (RWFM)

Mapeo α-μ

Conversión general del ruido de fase

Respuesta lineal

Propiedades de los filtros de tiempo y frecuencia

Funciones de sesgo

Referencias

  1. Cutler, L. S.; Searle, C. L. (febrero de 1966), «Some Aspects of the Theory and Measurements of Frequency Fluctuations in Frequency Standards», Proceedings of the IEEE (en inglés) 54 (2): 136-154, doi:10.1109/proc.1966.4627 .
  2. a b Leeson, D. B (febrero de 1966), «A simple Model of Feedback Oscillator Noise Spectrum», Proceedings of the IEEE (en inglés) 54 (2): 329-330, doi:10.1109/proc.1966.4682, archivado desde el original el 1 de febrero de 2014, consultado el 20 de septiembre de 2012 .
  3. a b Allan, D. Statistics of Atomic Frequency Standards, pages 221–230. Proceedings of the IEEE, Vol. 54, No 2, febrero de 1966.
  4. a b c «Definitions of physical quantities for fundamental frequency and time metrology – Random Instabilities». IEEE STD 1139-1999 (en inglés). 1999. ISBN 978-0-7381-1753-9. doi:10.1109/IEEESTD.1999.90575. 
  5. Rubiola, Enrico (2008), Phase Noise and Frequency Stability in Oscillators, Cambridge university press, ISBN 978-0-521-88677-2 .
  6. http://www.allanstime.com/Publications/DWA/Conversion_from_Allan_variance_to_Spectral_Densities.pdf. Uso incorrecto de la plantilla enlace roto (enlace roto disponible en Internet Archive; véase el historial, la primera versión y la última).
  7. Snyder, J.J. (1981). «An Ultra-High Resolution Frequency Meter». Thirty Fifth Annual Frequency Control Symposium (IEEE). doi:10.1109/freq.1981.200513. Consultado el 25 de febrero de 2022. 
  8. a b c Howe, D.A.; Allan, D.U.; Barnes, J.A. (1981). «Properties of Signal Sources and Measurement Methods». Thirty Fifth Annual Frequency Control Symposium (IEEE). doi:10.1109/freq.1981.200541. Consultado el 25 de febrero de 2022. 
  9. DIN EN 61883-7:2003-08, Audio/Video-Geräte der Unterhaltungselektronik_- Digitale Schnittstelle_- Teil_7: Übertragung nach Rec. ITU-R_BO.1294 System_B (IEC_61883-7:2003); Deutsche Fassung EN_61883-7:2003, Beuth Verlag GmbH, consultado el 25 de febrero de 2022 .
  10. ETSI ITS-G5 Standard. Springer International Publishing. 2020. pp. 442-442. Consultado el 25 de febrero de 2022.