Densidad espectral
El espectro de potencia Sxx()f){displaystyle S_{xx}(f)} de una serie de tiempo x()t){displaystyle x(t)} describe la distribución de energía en componentes de frecuencia que componen esa señal. Según el análisis de Fourier, cualquier señal física se puede descomponer en una serie de frecuencias discretas, o un espectro de frecuencias sobre un rango continuo. El promedio estadístico de una determinada señal o tipo de señal (incluyendo el ruido) analizado en términos de su contenido de frecuencia, se llama su espectro.
Cuando la energía de la señal se concentra alrededor de un intervalo de tiempo finito, especialmente si su energía total es finita, se puede calcular energía densidad espectral. Más comúnmente utilizado es el potencia densidad espectral (o simplemente espectro de potencia), que se aplica a las señales existentes sobre Todos tiempo, o durante un período de tiempo lo suficientemente grande (especialmente en relación con la duración de una medición) que también podría haber sido durante un intervalo de tiempo infinito. La densidad espectral de potencia (PSD) se refiere a la distribución de energía espectral que se encontraría por unidad de tiempo, ya que la energía total de tal señal a lo largo de todo el tiempo sería generalmente infinita. Summation or integration of the spectral components yields the total power (for a physical process) or variation (in a statistical process), identical to what would be obtained by integrating x2()t){displaystyle x^{2}(t)} sobre el dominio del tiempo, como lo dicta el teorema de Parseval.
El espectro de un proceso físico x()t){displaystyle x(t)} a menudo contiene información esencial sobre la naturaleza x{displaystyle x}. Por ejemplo, el tono y el timbre de un instrumento musical se determinan inmediatamente a partir de un análisis espectral. El color de una fuente de luz está determinado por el espectro del campo eléctrico de la onda electromagnética E()t){displaystyle E(t)} fluctúa a una frecuencia extremadamente alta. Obtener un espectro de series temporales como estas implica la transformación Fourier, y generalizaciones basadas en el análisis Fourier. En muchos casos el dominio del tiempo no se emplea específicamente en la práctica, como cuando se utiliza un prisma dispersivo para obtener un espectro de luz en un espectrógrafo, o cuando se percibe un sonido a través de su efecto en los receptores auditivos del oído interno, cada uno de los cuales es sensible a una frecuencia particular.
Sin embargo, este artículo se concentra en situaciones en las que la serie temporal se conoce (al menos en un sentido estadístico) o se mide directamente (como por ejemplo mediante un micrófono muestreado por una computadora). El espectro de potencia es importante en el procesamiento estadístico de señales y en el estudio estadístico de procesos estocásticos, así como en muchas otras ramas de la física y la ingeniería. Por lo general, el proceso es una función del tiempo, pero de manera similar se pueden analizar los datos en el dominio espacial que se descomponen en términos de frecuencia espacial.
Unidades
En física, la señal puede ser una onda, como una onda electromagnética, una onda acústica o la vibración de un mecanismo. La densidad espectral de potencia (PSD) de la señal describe la potencia presente en la señal en función de la frecuencia, por unidad de frecuencia. La densidad espectral de potencia se expresa comúnmente en vatios por hercio (W/Hz).
Cuando una señal se define únicamente en términos de voltaje, por ejemplo, no hay una potencia única asociada con la amplitud establecida. En este caso, "poder" se calcula simplemente en términos del cuadrado de la señal, ya que esto siempre sería proporcional a la potencia real entregada por esa señal en una impedancia dada. Entonces, uno podría usar unidades de V2 Hz−1 para la PSD. La densidad espectral de energía (ESD) tendría unidades que serían V2 s Hz−1, ya que la energía tiene unidades de potencia multiplicadas por el tiempo (por ejemplo,, vatio-hora).
En el caso general, las unidades de PSD serán la relación de unidades de varianza por unidad de frecuencia; así, por ejemplo, una serie de valores de desplazamiento (en metros) a lo largo del tiempo (en segundos) tendrá PSD en unidades de metros cuadrados por hercio, m2/Hz. En el análisis de vibraciones aleatorias, las unidades de g2 Hz−1 se utilizan con frecuencia para la PSD de aceleración, donde g denota la fuerza g.
Matemáticamente, no es necesario asignar dimensiones físicas a la señal ni a la variable independiente. En la siguiente discusión, el significado de x(t) permanecerá sin especificar, pero se supondrá que la variable independiente es la del tiempo.
Definición
Densidad espectral de energía
La densidad espectral de energía describe cómo la energía de una señal o una serie de tiempo se distribuye con frecuencia. Aquí, el término energía se utiliza en el sentido generalizado del procesamiento de señales; es decir, la energía E{displaystyle E} de una señal x()t){displaystyle x(t)} es:
- E≜ ≜ ∫ ∫ − − JUEGO JUEGO JUEGO JUEGO Silenciox()t)Silencio2dt.{displaystyle Etriangleq int _{-infty }left quietx(t)right WordPress^{2} dt.}
La densidad espectral de energía es más adecuada para transitorios, es decir, señales similares a pulsos, que tienen una energía total finita. Finito o no, el teorema de Parseval (o el teorema de Plancherel) nos da una expresión alternativa para la energía de la señal:
- ∫ ∫ − − JUEGO JUEGO JUEGO JUEGO Silenciox()t)Silencio2dt=∫ ∫ − − JUEGO JUEGO JUEGO JUEGO Silenciox^ ^ ()f)Silencio2df,{displaystyle int _{-infty}{infty } ¿Qué?
donde:
- x^ ^ ()f)≜ ≜ ∫ ∫ − − JUEGO JUEGO JUEGO JUEGO e− − i2π π ftx()t)dt{displaystyle {hat {x}(f)triangleq int _{-infty }{infty }e^{-i2pi ft}x(t) dt}
es el valor de la transformación de Fourier x()t){displaystyle x(t)} a frecuencia f{displaystyle f} (en Hz). El teorema también es cierto en los casos discretos. Puesto que el integral del lado izquierdo es la energía de la señal, el valor deSilenciox^ ^ ()f)Silencio2df{fnMicrosoft Sans Serif} puede ser interpretado como una función de densidad multiplicada por un intervalo de frecuencia infinitamente pequeño, describiendo la energía contenida en la señal a frecuencia f{displaystyle f} en el intervalo de frecuencia f+df{displaystyle f+df}.
Por tanto, el energía densidad espectral de x()t){displaystyle x(t)} se define como:
S̄ ̄ xx()f)≜ ≜ Silenciox^ ^ ()f)Silencio2{displaystyle {bar {S}_{xx}(f)triangleq left sometida{hat {x}(f)right sobre la vida {2}} | ()Eq.1) |
La función S̄ ̄ xx()f){displaystyle {bar {S}_{xx}(f)} y la autocorrelación de x()t){displaystyle x(t)} formar un par de transformación Fourier, un resultado también conocido como el teorema Wiener-Khinchin (ver también Periodograma).
Como ejemplo físico de cómo se puede medir la densidad espectral de energía de una señal, supongamos V()t){displaystyle V(t)} representa el potencial (en voltios) de un pulso eléctrico propagando a lo largo de una línea de transmisión de impedancia Z{displaystyle Z}, y suponer que la línea se termina con un resistor emparejado (para que toda la energía del pulso se entrega al resistor y ninguno se refleja de nuevo). Por la ley de Ohm, el poder entregado al resistor a la vez t{displaystyle t} es igual a V()t)2/Z{displaystyle V(t)^{2}/Z}, por lo que la energía total se encuentra mediante la integración V()t)2/Z{displaystyle V(t)^{2}/Z} con respecto al tiempo durante la duración del pulso. Para encontrar el valor de la densidad espectral de energía S̄ ̄ xx()f){displaystyle {bar {S}_{xx}(f)} a frecuencia f{displaystyle f}, se puede insertar entre la línea de transmisión y el resistor un filtro de bandpass que pasa sólo una gama estrecha de frecuencias (Δ Δ f{displaystyle Delta f}, digamos) cerca de la frecuencia de interés y luego mide la energía total E()f){displaystyle E(f)} disipado a través de la resistencia. El valor de la densidad espectral de energía f{displaystyle f} se calcula que E()f)/Δ Δ f{displaystyle E(f)/Delta f}. En este ejemplo, desde el poder V()t)2/Z{displaystyle V(t)^{2}/Z} tiene unidades de V2 Ω−1, la energía E()f){displaystyle E(f)} tiene unidades de V2SΩ−1= J, y por consiguiente la estimación E()f)/Δ Δ f{displaystyle E(f)/Delta f} de la densidad espectral de energía tiene unidades de J Hz−1, según sea necesario. En muchas situaciones, es común olvidar el paso de la división por Z{displaystyle Z} para que la densidad espectral de energía tenga unidades de V2Hz−2.
Esta definición generaliza de manera directa a una señal discreta con un número contablemente infinito de valores xn{displaystyle x_{n} como una señal muestrada en tiempos discretos tn=t0+()nΔ Δ t){displaystyle t_{n}=t_{0}+(n,Delta t)}:
- S̄ ̄ xx()f)=limN→ → JUEGO JUEGO ()Δ Δ t)2Silencio.. n=− − NNxne− − i2π π fnΔ Δ tSilencio2⏟ ⏟ Silenciox^ ^ d()f)Silencio2,{displaystyle {bar {S}_{xx}(f)=lim _{Nto infty }(Delta t)^{2}underbrace {left durablesum ¿Por qué? - ¿Por qué?
Donde x^ ^ d()f){fnMicrosoft Sans Serif} es la transformación discreta de Fourier xn.{displaystyle x_{n}El intervalo de muestreo Δ Δ t{displaystyle Delta t} es necesario para mantener las unidades físicas correctas y para asegurar que recuperamos el caso continuo en el límite Δ Δ t→ → 0.{displaystyle Delta tto 0.}Pero en las ciencias matemáticas el intervalo se establece a menudo a 1, que simplifica los resultados a expensas de la generalidad. (también ver frecuencia normalizada)
Densidad espectral de potencia
La definición anterior de densidad espectral de energía es adecuada para los transitorios (señales similares al púlsculo) cuya energía se concentra alrededor de una ventana de tiempo; entonces los Fourier transforma de las señales generalmente existen. Para señales continuas a lo largo de todo el tiempo, uno debe definir el potencia densidad espectral (PSD) que existe para procesos estacionarios; esto describe cómo el poder de una serie de señal o tiempo se distribuye a través de la frecuencia, como en el simple ejemplo dado anteriormente. Aquí, el poder puede ser el poder físico real, o más a menudo, para comodidad con señales abstractas, se identifica simplemente con el valor cuadrado de la señal. Por ejemplo, los estadísticos estudian la diferencia de una función con el tiempo x()t){displaystyle x(t)} (o sobre otra variable independiente), y utilizando una analogía con señales eléctricas (entre otros procesos físicos), es habitual referirse a ella como la espectro de potencia incluso cuando no hay poder físico involucrado. Si se creara una fuente de tensión física que siguió x()t){displaystyle x(t)} y lo aplicó a los terminales de un oohm resistor, entonces de hecho el poder instantáneo disipado en ese resistor sería dado por x()t)2{displaystyle x(t)^{2} Watts.
El poder promedio P{displaystyle P} de una señal x()t){displaystyle x(t)} en todo el tiempo se da por el promedio del tiempo siguiente, donde el período T{displaystyle T} está centrado en algún tiempo arbitrario t=t0{displaystyle t=t_{0}:
- P=limT→ → JUEGO JUEGO 1T∫ ∫ t0− − T/2t0+T/2Silenciox()t)Silencio2dt{displaystyle P=lim _{Tto infty}{frac {1}{T}int _{t_{0}-T/2} {t_{0}
Sin embargo, en aras de tratar con las matemáticas que siguen, es más conveniente tratar con los límites de tiempo en la propia señal en lugar de los límites de tiempo en los límites de la integral. Como tal, tenemos una representación alternativa del poder promedio, donde xT()t)=x()t)wT()t){displaystyle x_{T}(t)=x(t)w_{T}(t)} y wT()t){displaystyle w_{T}(t)} es unidad dentro del período arbitrario y cero en otros lugares.
- P=limT→ → JUEGO JUEGO 1T∫ ∫ − − JUEGO JUEGO JUEGO JUEGO SilencioxT()t)Silencio2dt{displaystyle P=lim _{Tto infty}{frac {1}{T}int - No. - ¿Qué?
Claramente, en los casos en que la expresión anterior para P no es cero (incluso cuando T crece sin límite), la integral misma también debe crecer sin límite. Esa es la razón por la que no podemos usar la densidad espectral de energía en sí misma, que es esa integral divergente, en tales casos.
Al analizar el contenido de frecuencia de la señal x()t){displaystyle x(t)}, a uno le gustaría computar la transformación ordinaria Fourier x^ ^ ()f){displaystyle {hat {x}(f)}; sin embargo, para muchas señales de interés la transformación Fourier no existe formalmente. Independientemente, El teorema de Parseval nos dice que podemos reescribir el poder promedio como sigue.
- P=limT→ → JUEGO JUEGO 1T∫ ∫ − − JUEGO JUEGO JUEGO JUEGO Silenciox^ ^ T()f)Silencio2df{displaystyle P=lim _{Tto infty}{frac {1}{T}int ¿Por qué? T}(f) forever^{2},df}
Entonces, la densidad espectral de potencia se define simplemente como el integrando anterior.
Sxx()f)=limT→ → JUEGO JUEGO 1TSilenciox^ ^ T()f)Silencio2{displaystyle S_{xx}(f)=lim _{Tto infty }{frac {1}{fn} {x}_{T}(f) habit^{2},} | ()Eq.2) |
Desde aquí también podemos ver Silenciox^ ^ T()f)Silencio2{fnMicrosoft Sans Serif} como la transformación Fourier del tiempo convolución xTAlternativa Alternativa ()− − t){displaystyle ¿Qué? y xT()t){displaystyle x_{T}(t)}
- Silenciox^ ^ T()f)Silencio2=F{}xTAlternativa Alternativa ()− − t)Alternativa Alternativa xT()t)}=∫ ∫ − − JUEGO JUEGO JUEGO JUEGO [∫ ∫ − − JUEGO JUEGO JUEGO JUEGO xTAlternativa Alternativa ()t− − τ τ )xT()t)dt]e− − i2π π fτ τ dτ τ {fnMicrosoft Sans Serif}= {fnMitcal} [F]left{x_}{*}(-t)mathbin {mathbf {*} x_{T}(t)right=int - No.. - No. }x_{T}{*}(t-tau)x_{T}(t)dtright]e^{-i2pi ftau } dtau }
Ahora, si dividimos la convolución del tiempo arriba por el período T{displaystyle T} y tomar el límite como T→ → JUEGO JUEGO {displaystyle Trightarrow infty }, se convierte en la función de autocorrelación de la señal no ventana x()t){displaystyle x(t)}, que se denota como Rxx()τ τ ){displaystyle R_{xx}(tau)}, siempre que x()t){displaystyle x(t)} es ergodic, que es cierto en la mayoría, pero no todos, casos prácticos..
- limT→ → JUEGO JUEGO 1TSilenciox^ ^ T()f)Silencio2=∫ ∫ − − JUEGO JUEGO JUEGO JUEGO [limT→ → JUEGO JUEGO 1T∫ ∫ − − JUEGO JUEGO JUEGO JUEGO xTAlternativa Alternativa ()t− − τ τ )xT()t)dt]e− − i2π π fτ τ dτ τ =∫ ∫ − − JUEGO JUEGO JUEGO JUEGO Rxx()τ τ )e− − i2π π fτ τ dτ τ {displaystyle lim _{Tto infty}{frac {1} {T}fn} {x}_{T}(f)right perpetua{2}=int ¿Por qué? {1}{T}int - No. }x_{T}{*}(t-tau)x_{T}(t)dtright]e^{-i2pi ftau } ♪ ♪ ♪ ♪ ♪ ¿Qué?
Desde aquí vemos, otra vez asumiendo la ergodicidad de x()t){displaystyle x(t)}, que la densidad espectral de potencia se puede encontrar como la transformación Fourier de la función de autocorrelación (Teorema de Windows–Khinchin).
Sxx()f)=∫ ∫ − − JUEGO JUEGO JUEGO JUEGO Rxx()τ τ )e− − i2π π fτ τ dτ τ =R^ ^ xx()f){displaystyle S_{xx}(f)=int _{-infty }R_{xx}(tau)e^{-i2pi ftau },dtau ={hat {R}_{x}(f)}} | ()Eq.3) |
Muchos autores utilizan esta igualdad para definir la densidad espectral de potencia.
El poder de la señal en una banda de frecuencia dada [f1,f2]{displaystyle [f_{1},f_{2}}, donde <math alttext="{displaystyle 0<f_{1}0.f1.f2{displaystyle 0 0 0 0 0<img alt="{displaystyle 0<f_{1}, se puede calcular mediante la integración sobre frecuencia. Desde Sxx()− − f)=Sxx()f){displaystyle S_{xx}(-f)=S_{xx}(f)}, una cantidad igual de poder se puede atribuir a bandas de frecuencia positiva y negativa, lo que representa el factor de 2 en la siguiente forma (estos factores triviales dependen de las convenciones utilizadas):
- Pbandlimited=2∫ ∫ f1f2Sxx()f)df{displaystyle P_{textsf {bandlimited}=2int ¿Qué?
En términos más generales, se pueden utilizar técnicas similares para estimar una densidad espectral que varía el tiempo. En este caso el intervalo de tiempo T{displaystyle T} es finito en lugar de acercarse al infinito. Esto da lugar a una disminución de la cobertura espectral y la resolución desde las frecuencias inferiores a las 1/T{displaystyle 1/T} no se muestra, y resultados en frecuencias que no son un número entero de 1/T{displaystyle 1/T} no son independientes. Sólo utilizando una serie de tiempo, el espectro de potencia estimado será muy "noisy"; sin embargo esto se puede aliviar si es posible evaluar el valor esperado (en la ecuación anterior) utilizando un gran (o infinito) número de espectros de corto plazo correspondientes a conjuntos estadísticos de realizaciones de realizaciones de x()t){displaystyle x(t)} evaluado en la ventana de tiempo especificada.
Al igual que con la densidad espectral de energía, la definición de la densidad espectral de potencia se puede generalizar para variables de tiempo discretas xn{displaystyle x_{n}. Como antes, podemos considerar una ventana − − N≤ ≤ n≤ ≤ N{displaystyle -Nleq nleq N} con la señal mostrada en tiempos discretos xn=x0+()nΔ Δ t){displaystyle x_{n}=x_{0}+(n,Delta t)} para un período de medición total T=()2N+1)Δ Δ t{displaystyle T=(2N+1),Delta t}.
- Sxx()f)=limN→ → JUEGO JUEGO ()Δ Δ t)2TSilencio.. n=− − NNxne− − i2π π fnΔ Δ tSilencio2{displaystyle S_{xx}(f)=lim _{Nto infty }{frac {Delta t)^{2}}{T}}left permanentlysum ¿Por qué?
Tenga en cuenta que una sola estimación de la PSD se puede obtener a través de un número finito de muestreos. Como antes, la PSD real se logra cuando N{displaystyle N} (y así T{displaystyle T}) aborda el infinito y el valor esperado se aplica formalmente. En una aplicación del mundo real, uno normalmente promediaría una PSD de medición finita sobre muchos ensayos para obtener una estimación más exacta de la PSD teórica del proceso físico subyacente a las mediciones individuales. Esta PSD computada a veces se llama periodograma. Este periodograma converge al verdadero PSD como el número de estimaciones, así como el intervalo de tiempo promedio T{displaystyle T} infinidad de enfoque (Brown & Hwang).
Si dos señales poseen densidades espectrales de potencia, entonces la densidad espectral cruzada se puede calcular de manera similar; como la PSD está relacionada con la autocorrelación, también lo está la densidad espectral cruzada con la correlación cruzada.
Propiedades de la densidad espectral de potencia
Algunas propiedades del PSD incluyen:
- El espectro de potencia es siempre real y no negativo, y el espectro de un proceso real valorado es también una función uniforme de la frecuencia: Sxx()− − f)=Sxx()f){displaystyle S_{xx}(-f)=S_{xx}(f)}.
- Para un proceso estocástico continuo x(t), la función de autocorrelación Rxx()t) puede ser reconstruido de su espectro de potencia Sxxf) mediante la transformación inversa de Fourier
- Utilizando el teorema de Parseval, se puede calcular la varianza (poder de promedio) de un proceso integrando el espectro de potencia sobre toda frecuencia:
- P=Var ()x)=∫ ∫ − − JUEGO JUEGO JUEGO JUEGO Sxx()f)df{displaystyle P=operatorname {Var} (x)=int _{-infty }infty }!S_{x}(f),df}
- Para un proceso real x()t) con densidad espectral de potencia Sxx()f){displaystyle S_{xx}(f)}, uno puede calcular espectro integrado o distribución espectral F()f){displaystyle F(f)}, que especifica el promedio bandlimited potencia contenida en frecuencias de DC a f utilizando:
- F()f)=2∫ ∫ 0fSxx()f.)df..{displaystyle F(f)=2int _{0}S_{xx}(f'),df'}
Tenga en cuenta que la expresión anterior para el poder total (variasignal) es un caso especial en el que ô → ∞.
Densidad espectral de potencia cruzada
Dadas dos señales x()t){displaystyle x(t)} y Sí.()t){displaystyle y(t)}, cada uno de los cuales poseen densidades espectral Sxx()f){displaystyle S_{xx}(f)} y SSí.Sí.()f){displaystyle S_{yy}(f)}, es posible definir un fuerza transversal densidad espectral ()CPSD) o transversal densidad espectral ()CSD). Para comenzar, consideremos el poder promedio de tal señal combinada.
- P=limT→ → JUEGO JUEGO 1T∫ ∫ − − JUEGO JUEGO JUEGO JUEGO [xT()t)+Sí.T()t)]Alternativa Alternativa [xT()t)+Sí.T()t)]dt=limT→ → JUEGO JUEGO 1T∫ ∫ − − JUEGO JUEGO JUEGO JUEGO SilencioxT()t)Silencio2+xTAlternativa Alternativa ()t)Sí.T()t)+Sí.TAlternativa Alternativa ()t)xT()t)+SilencioSí.T()t)Silencio2dt{displaystyle {begin{aligned}P pl=lim ¿Qué? {1}{T}int ¿Por qué? ¿Qué? {1}{T}int - No. }Princex_{T}(t) Anterior{2}+x_{T}{*}(t)y_{T}(t)+y_{T} {*}(t)x_{T}(t)+ Solicioso_{T}(t) Sober_{2}dt\\\end{aligned}}}}}}}}
Usando la misma notación y los mismos métodos que se usaron para la derivación de la densidad espectral de potencia, explotamos el teorema de Parseval y obtenemos
- SxSí.()f)=limT→ → JUEGO JUEGO 1T[x^ ^ TAlternativa Alternativa ()f)Sí.^ ^ T()f)]SSí.x()f)=limT→ → JUEGO JUEGO 1T[Sí.^ ^ TAlternativa Alternativa ()f)x^ ^ T()f)]{displaystyle {begin{aligned}S_{xy}(f) limit=lim _{Tto infty }{frac {1} {fn}left[{hat] {x}_{*}(f){hat {y}_{T}(f)right] {1} {fn}left[{hat] {fnMicrosoft Sans Serif}
donde, de nuevo, las contribuciones Sxx()f){displaystyle S_{xx}(f)} y SSí.Sí.()f){displaystyle S_{yy}(f)} ya se entienden. Note que SxSí.Alternativa Alternativa ()f)=SSí.x()f){displaystyle S_{xy} {}(f)=S_{yx}(f)}, por lo que la contribución completa al poder de la cruz es, en general, del doble de la parte real de cada individuo CPSD. Como antes, desde aquí retransmitimos estos productos como la transformación Fourier de una convolución temporal, que cuando se divide por el período y se lleva al límite T→ → JUEGO JUEGO {displaystyle Tto infty } se convierte en la transformación de Fourier de una función de corelación cruzada.
- SxSí.()f)=∫ ∫ − − JUEGO JUEGO JUEGO JUEGO [limT→ → JUEGO JUEGO 1T∫ ∫ − − JUEGO JUEGO JUEGO JUEGO xTAlternativa Alternativa ()t− − τ τ )Sí.T()t)dt]e− − i2π π fτ τ dτ τ =∫ ∫ − − JUEGO JUEGO JUEGO JUEGO RxSí.()τ τ )e− − i2π π fτ τ dτ τ SSí.x()f)=∫ ∫ − − JUEGO JUEGO JUEGO JUEGO [limT→ → JUEGO JUEGO 1T∫ ∫ − − JUEGO JUEGO JUEGO JUEGO Sí.TAlternativa Alternativa ()t− − τ τ )xT()t)dt]e− − i2π π fτ τ dτ τ =∫ ∫ − − JUEGO JUEGO JUEGO JUEGO RSí.x()τ τ )e− − i2π π fτ τ dτ τ {displaystyle {begin{aligned}S_{xy}(f) ¿Por qué? {1}{T}int - No. }x_{T} {*}(t-tau)y_{T}(t)dtright)e^{-i2pi ftau }dtau =int _{-infty } {infty }R_{xy}(tau)e^{-i2pi ftau }dtautau ¿Por qué? {1}{T}int - No. }y_{T} {*}(t-tau)x_{T}(t)dtright)e^{-i2pi ftau }dtau =int _{-infty }R_{yx}(tau)e^{-i2pi ftau)d {tau}d}dtauend}
Donde RxSí.()τ τ ){displaystyle R_{xy}(tau)} es la cruz-correlación de x()t){displaystyle x(t)} con Sí.()t){displaystyle y(t)} y RSí.x()τ τ ){displaystyle R_{yx}(tau)} es la cruz-correlación de Sí.()t){displaystyle y(t)} con x()t){displaystyle x(t)}. A la luz de esto, se considera que el PSD es un caso especial de la CDS para x()t)=Sí.()t){displaystyle x(t)=y(t)}. Si x()t){displaystyle x(t)} y Sí.()t){displaystyle y(t)} son señales reales (por ejemplo, tensión o corriente), sus transformaciones Fourier x^ ^ ()f){displaystyle {hat {x}(f)} y Sí.^ ^ ()f){displaystyle {hat {y}(f)} suelen limitarse a frecuencias positivas por convención. Por lo tanto, en el procesamiento de señal típica, el completo CPSD es sólo uno de los CPSDs escalada por un factor de dos.
- CPSDTotal=2SxSí.()f)=2SSí.x()f){displaystyle operatorname {CPSD} _{text{Full}=2S_{xy}(f)=2S_{yx}(f)}
Para señales discretas xn e yn, la relación entre la densidad espectral cruzada y la la covarianza cruzada es
- SxSí.()f)=.. n=− − JUEGO JUEGO JUEGO JUEGO RxSí.()τ τ n)e− − i2π π fτ τ nΔ Δ τ τ {displaystyle S_{xy}(f)=sum _{n=-infty }{infty }R_{xy}(tau _{n})e^{-i2pi ftau _{n},Delta tau }
Estimación
El objetivo de la estimación de densidad espectral es estimar la densidad espectral de una señal aleatoria a partir de una secuencia de muestras de tiempo. Dependiendo de lo que se sepa sobre la señal, las técnicas de estimación pueden implicar enfoques paramétricos o no paramétricos, y pueden basarse en análisis en el dominio del tiempo o en el dominio de la frecuencia. Por ejemplo, una técnica paramétrica común consiste en ajustar las observaciones a un modelo autorregresivo. Una técnica no paramétrica común es el periodograma.
La densidad espectral generalmente se estima mediante métodos de transformada de Fourier (como el método de Welch), pero también se pueden utilizar otras técnicas, como el método de máxima entropía.
Conceptos relacionados
- El centroide espectralide de una señal es el punto medio de su función de densidad espectral, es decir, la frecuencia que divide la distribución en dos partes iguales.
- El frecuencia espectral ()SEF), generalmente expresado como "SEF x", representa la frecuencia debajo de la cual x por ciento del poder total de una señal dada se encuentran; típicamente, x está en el rango 75 a 95. Es más particularmente una medida popular utilizada en el monitoreo de la EEG, en cuyo caso SEF se ha utilizado de manera diferente para estimar la profundidad de la anestesia y las etapas del sueño.
- A sobre espectral es la curva de sobre de la densidad del espectro. Describe un punto en el tiempo (una ventana, para ser precisa). Por ejemplo, en la teleobservación utilizando un espectrómetro, el sobre espectral de una característica es el límite de sus propiedades espectrales, tal como se define por la gama de niveles de brillo en cada una de las bandas espectrales de interés.
- La densidad espectral es una función de frecuencia, no una función del tiempo. Sin embargo, la densidad espectral de una pequeña ventana de una señal más larga puede ser calculada, y trazada contra el tiempo asociado con la ventana. Tal gráfico se llama espectrograma. Esta es la base de una serie de técnicas de análisis espectral tales como la transformación y ondas Fourier de corto plazo.
- Un "espectro" generalmente significa la densidad espectral de potencia, como se discutió anteriormente, que representa la distribución del contenido de señal sobre la frecuencia. Para las funciones de transferencia (por ejemplo, Bode plot, chirp) la respuesta de frecuencia completa puede ser graficada en dos partes: potencia versus frecuencia y fase versus frecuencia; la fase de densidad espectral, espectro de fase, o fase espectral. Menos comúnmente, las dos partes pueden ser las partes reales e imaginarias de la función de transferencia. Esto no debe confundirse con el frecuencia de respuesta de una función de transferencia, que también incluye una fase (o equivalente, una parte real e imaginaria) como función de frecuencia. La respuesta del impulso de dominio del tiempo h()t){displaystyle h(t)} generalmente no se puede recuperar únicamente de la densidad espectral de potencia sola sin la parte de fase. Aunque también son pares de transformación Fourier, no hay simetría (como hay para la autocorrelación) que obligue a la transformación Fourier a ser real-valorada. Ver Pulso Ultrashort#Fase espacial, ruido de fase, retraso de grupo.
- A veces uno encuentra un densidad espectral ()ASD), que es la raíz cuadrada de la PSD; la ASD de una señal de tensión tiene unidades de V Hz1/2−. Esto es útil cuando el forma del espectro es bastante constante, ya que las variaciones en la ASD serán proporcionales a las variaciones en el nivel de tensión de la señal en sí. Pero es matemáticamente preferido utilizar el PSD, ya que sólo en ese caso es el área bajo la curva significativa en términos de potencia real sobre toda frecuencia o sobre un ancho de banda especificado.
Aplicaciones
Cualquier señal que pueda representarse como una variable que varía en el tiempo tiene un espectro de frecuencia correspondiente. Esto incluye entidades familiares como la luz visible (percibida como color), las notas musicales (percibidas como tono), la radio/TV (especificadas por su frecuencia o, a veces, por su longitud de onda) e incluso la rotación regular de la tierra. Cuando estas señales se ven en forma de espectro de frecuencia, se revelan ciertos aspectos de las señales recibidas o los procesos subyacentes que las producen. En algunos casos, el espectro de frecuencias puede incluir un pico distinto correspondiente a un componente de onda sinusoidal. Y adicionalmente puede haber picos correspondientes a armónicos de un pico fundamental, indicando una señal periódica que no es simplemente sinusoidal. O un espectro continuo puede mostrar intervalos de frecuencia estrechos que están fuertemente realzados correspondientes a resonancias, o intervalos de frecuencia que contienen una potencia casi nula, como los produciría un filtro de muesca.
Ingeniería eléctrica
El concepto y el uso del espectro de potencia de una señal es fundamental en la ingeniería eléctrica, especialmente en los sistemas de comunicación electrónica, incluidas las comunicaciones por radio, los radares y los sistemas relacionados, además de la tecnología de detección remota pasiva. Los instrumentos electrónicos llamados analizadores de espectro se utilizan para observar y medir los espectros de potencia de las señales.
El analizador de espectro mide la magnitud de la transformada de Fourier de tiempo corto (STFT) de una señal de entrada. Si la señal que se analiza puede considerarse un proceso estacionario, la STFT es una buena estimación suavizada de su densidad espectral de potencia.
Cosmología
Las fluctuaciones primordiales, las variaciones de densidad en el universo primitivo, se cuantifican mediante un espectro de potencia que da la potencia de las variaciones en función de la escala espacial.
Ciencias del clima
El análisis espectral de potencia se ha utilizado para examinar las estructuras espaciales para la investigación climática. Estos resultados sugieren que la turbulencia atmosférica vincula el cambio climático con una mayor volatilidad regional local en las condiciones climáticas.
Contenido relacionado
Conducto atmosférico
Telururo
Instituto de prospectiva