Densidad espectral


Densidad espectral

Densidad espectral

DEP de una señal aleatoria en la banda de 0kHz-20kHz estimada mediante el método de Welch

En matemáticas y en física, la Densidad Espectral (Spectral Density) de una señal es una función matemática que nos informa de cómo está distribuida la potencia o la energía (según el caso) de dicha señal sobre las distintas frecuencias de las que está formada, es decir, su espectro.

La definición matemática de la Densidad Espectral (DE) es diferente dependiendo de si se trata de señales definidas en energía, en cuyo caso hablamos de Densidad Espectral de Energía (DEE), o en potencia, en cuyo caso hablamos de Densidad Espectral de Potencia (DEP).

Aunque la densidad espectral no es exactamente lo mismo que el espectro de una señal, a veces ambos términos se usan indistintamente, lo cual, en rigor, es incorrecto.

Contenido

Definición matemática

Para el caso de señales definidas en energía

Una señal x(t) es definida en energía si su energía media es finita, i.e,  0 < E_x < \infty y por tanto, su potencia media es cero. Otra forma de decir lo mismo es si la integral de su valor absoluto al cuadrado existe y es finita.

Su DEE es

S_{xx}(f) = \left|X(f)\right|^2 expresado en [J/Hz]

donde X(f) es la Transformada de Fourier de x(t), la integral de esta función en todo el eje f es el valor de la energía total de la señal x(t)

E=\int_{-\infty}^{+\infty} S_{xx}(f) df

Para el caso de señales definidas en potencia

Una señal x(t) es definida en potencia si su potencia media es finita, i.e,  0 < P_x < \infty y por tanto, su energía media es infinita,  E_x=\infty.

La DEP se calcula usando el teorema de Wiener-Khinchin el cual relaciona la DEP con la transformada de Fourier de la función de autocorrelación

S_{xx}(f)=\mbox{TF}\left\{R_{xx}(\tau)\right\}=\int_{-\infty}^{\infty}\,R_{xx}(\tau)\,e^{-2\pi\,if\tau}\,d \tau expresado en [W/Hz]

donde TF significa Transformada de Fourier y Rxx es la función de autocorrelación de x(t).

El valor Sxx(0) es la potencia de la componente continua (DC) de la señal. La integral de esta función en todo el eje f es el valor de la potencia total de la señal x(t)

P=\int_{-\infty}^{+\infty} S_{xx}(f) df

Usando el concepto de correlación cruzada es posible definir también la densidad espectral cruzada.

S_{xy}(f)=\mbox{TF}\left\{R_{xy}(\tau)\right\}


Nota: En realidad, la definición de la DEP sirve también para las señales definidas en energía, que serían un caso particular. En este caso la transformada de Fourier de la autocorrelación sería simplemente la transformada de Fourier al cuadrado, es decir, la DEE.

Propiedades

Relativas a los sistemas lineales e invariantes con el tiempo

En un sistema lineal e invariante con el tiempo en el que x(t) es la entrada, h(t) la respuesta al impulso e y(t) la salida del sistema. Tenemos las siguientes propiedades:

m_y = m_x \cdot H(0)
S_{yy}(\omega)=S_{xx}(\omega) \cdot |H(\omega)|^2
S_{xy}(\omega)=S_{xx}(\omega) \cdot H^*(\omega)
S_{yx}(\omega)=S_{xy}^*(\omega) = S_{xx}(\omega) \cdot H(\omega)

donde my es la media de y(t) y Sxy es la densidad espectral cruzada entre x(t) e y(t)

Suma de procesos

En general, la Densidad Espectral de la suma NO es suma de Densidades Espectrales. Esto sólo es cierto si ambos procesos no están correlacionados. En general, si tenemos:

z(t) = x(t) + y(t)

donde x(t) e y(t) son conjuntamente estacionarios, entonces

R_{zz}(\tau) = R_{xx}(\tau) + R_{yy}(\tau) + R_{xy}(\tau) + R_{yx}(\tau) \,
S_{zz}(\omega) = S_{xx}(\omega) + S_{yy}(\omega) + S_{xy}(\omega) + S_{yx}(\omega) \,
S_{zz}(\omega) = S_{xx}(\omega) + S_{yy}(\omega) + S_{xy}(\omega) + S_{xy}^*(\omega) \,
S_{zz}(\omega) = S_{xx}(\omega) + S_{yy}(\omega) + 2\mbox{Re}\left[S_{xy}(\omega)\right]

Estimación de la densidad espectral

Un problema muy común y con grandes aplicaciones prácticas en procesado de señal es el de estimar la densidad espectral de potencia de una señal aleatoria estacionaria. Decimos "estimar" puesto que, como la señal es un proceso estocástico (estacionario) dada la naturaleza estocástica del mismo no es posible determinar con absoluta precisión su DEP a no ser que dispongamos de un registro de señal infinito, lo cual no es posible.

Las técnicas de estimación se dividen en dos grandes grupos:

  • No Paramétricas. Están basadas siempre de una u otra forma en el cálculo del periodograma. Calcular la transformada de fourier (en un ordenador es la DFT) de un registro de señal para estimar su espectro es un ejemplo de técnica no paramétrica.
  • Paramétricas. Consisten en suponer un determinado modelo para el proceso estocástico (modelos AR, MA, ARMA, etc) y en la estimación de los parámetros de estos modelos mediante técnicas de predicción lineal (filtrado lineal óptimo) u otros métodos.

Acerca de los procesos estocásticos no estacionarios

La DE sólo está matemáticamente bien definida en el caso de señales con una función de autocorrelación estacionaria, i.e, que no dependa de la posición de las variables aleatorias que componen el proceso sino sólo de la distancia entre ellas. Es decir, la DE sólo está bien definida para el caso de señales deterministas y señales aleatorias estacionarias.

Un proceso aleatorio no estacionario que es estacionario a trozos se llama cuasi-estacionario y es posible definir la DEP en cada uno de estos trozos. Para estimar la DEP en este tipo de procesos lo normal es usar un método de estimación espectral paramétrico adaptativo (por ejemplo mediante un modelo AR y el algoritmo LMS para identificar el modelo AR).

Aplicaciones

  • Intuitivamente, la Densidad Espectral sirve para identificar periodicidades escondidas en una función de variable continua o de variable discreta (secuencia de números)
  • Estimar la entropía de un proceso aleatorio (las señales deterministas obviamente no tienen entropía). Cuanto más plana es la DEP de una señal aleatoria, más entropía contiene. Una señal aleatoria cuya DEP sea perfectamente plana se llama ruido blanco, no contiene redundancia y por tanto no puede ser comprimida (sin pérdidas).
  • Una vez conocida su entropía, comprimir con o sin pérdidas una señal de audio o vídeo (codecs FLAC/MP3/OGG, DIVX/THEORA/etc) o restaurar sus propiedades
  • Proporciona información muy valiosa sobre de la dinámica interna de muchos sistemas físicos. Sirve para identificar elementos o compuestos químicos (espectroscopía). También sirve para la identificación de modelos matemáticos lineales en teoría de control

Referencias

  • "Tratamiento digital de señales. Principios, algoritmos y aplicaciones". John G. Proakis, Dimitris G. Manolakis. Prentice Hall, ISBN 0-13-373762-4
Obtenido de "Densidad espectral"

Wikimedia foundation. 2010.

Mira otros diccionarios:

  • Manchado espectral — Cuando se calcula la transformada discreta de Fourier (DFT) de una señal (por ejemplo, para obtener la densidad espectral de potencia), el manchado espectral o leakage ocurrirá si alguna de las componentes frecuenciales cae entre frecuencias de… …   Wikipedia Español

  • Línea espectral — Las líneas espectrales se detectan como líneas de absorción (A) o líneas de emisión (B) dependiendo de las posiciones del detector, el gas y la fuente luminosa. Una línea espectral es una línea oscura o brillante en un espectro uniforme y… …   Wikipedia Español

  • Ruido de cuantificación — Saltar a navegación, búsqueda Figura 1: Procesos de la conversión A/D. Se define como error de cuantificación o ruido de cuantificación a la señal en tiempo discreto y amplitud continua introducida por el proceso de cuantificación (uno de los… …   Wikipedia Español

  • Ruido blanco — Saltar a navegación, búsqueda Ejemplo de forma de onda de un ruido blanco …   Wikipedia Español

  • Ruido de Johnson-Nyquist — Saltar a navegación, búsqueda El ruido de Johnson–Nyquist (ruido térmico, ruido de Johnson, o ruido de Nyquist) se genera por la agitación térmica de los portadores de carga (generalmente electrones dentro de un conductor) en equilibrio, lo que… …   Wikipedia Español

  • Ruido de color — Saltar a navegación, búsqueda 10 segundos de ruido blanco …   Wikipedia Español

  • Espectro de frecuencias — Saltar a navegación, búsqueda Espectro de frecuencias de la luz emitida por átomos de hierro libres en la región visible del espectro electromagnético. El espectro de frecuencia de un fenómeno ondulatorio (sonoro, luminoso o electromagnético),… …   Wikipedia Español

  • Cosmología física — Saltar a navegación, búsqueda Para otros usos de este término, véase Cosmología. Imagen de la radiación de fondo de microondas por la sonda WMAP. La Cosmología física, es la rama de la astrofísica, que estudia la est …   Wikipedia Español

  • Modelo autorregresivo de media móvil — En estadística, los modelos autorregresivos de media móvil (en inglés AutoRegressive Moving Average models, abreviados ARMA), también llamados Modelos Box Jenkins, se aplican a series temporales de datos. Dada una serie temporal de datos Xt, el… …   Wikipedia Español

  • Teorema de Shannon-Hartley — En teoría de la información, el teorema de Shannon Hartley es una aplicación del teorema de codificación para canales con ruido. Un caso muy frecuente es el de un canal de comunicación analógico continuo en el tiempo que presenta un ruido… …   Wikipedia Español