Explicación de PSD (densidad espectral de potencia)


14

Estoy tratando de entender cómo se calcula el PSD. He buscado en algunos de mis libros de texto de Ingeniería de la Comunicación, pero fue en vano. También he buscado en línea. Wikipedia parece tener la mejor explicación; sin embargo, me pierdo en la parte donde deciden hacer la CDF (función de distribución acumulativa) y luego, por alguna razón, decido relacionar eso con la función de autocorrelación.

Supongo que lo que no entiendo es, ¿cómo la autocorrelación tiene algo que ver con el cálculo de la PSD? Pensé que la PSD simple sería la Transformada de Fourier de (donde es la potencia de la señal con respecto al tiempo).P ( t )PAG(t)PAG(t)


¿Cómo se define ? P(t)
Phonon

Realmente no lo defino como algo. Es solo una señal de energía. Supongo que si tuviera que definirlo, sería ... Supongo que el punto es que el PSD no es F { P ( t ) } y tiene algo que ver con la autocorrelación y no conseguir lo ...P(t)=v(t)i(t)F{P(t)}
user968243

Realmente no se puede definir un poder así para señales arbitrarias. No hay conceptos de voltaje y corriente. La potencia en este caso se define como la potencia de una onda (electromagnética si lo desea). Entonces es , y es un número único, no una cantidad variable en el tiempo. 1T0 0TX2(t)ret
Phonon

1
Lea sobre el teorema de Wiener-Khinchin . Se niega a comprender lo que Phonon le señala que el límite que está calculando es una constante y, por lo tanto, su transformación de Fourier es solo un impulso en en el dominio de la frecuencia. Si eso flota su bote, hágalo, pero no es la densidad espectral de potencia como todos los demás lo entienden. f=0
Dilip Sarwate

1
He leído sobre ese teorema ... Y entiendo cómo se relaciona la transformada de Fourier con la autocorrelación. Y no me niego a entender lo que dijo Phonon ... Entiendo exactamente lo que dijo @Phonon. Lo que no entiendo es por qué se usa la fórmula de autocorrelación y tampoco entiendo por qué se usa la forma de transformación de Fourier (para obtener el PSD, puede tomar la transformación de Fourier, tomar su magnitud, cuadrarla, etc.) ... No tengo idea de por qué hacerlo daría un PSD y no he podido encontrar una derivación decente.
user968243

Respuestas:


17

Tienes razón, PSD tiene que ver con el cálculo de la Transformada de Fourier de la potencia de la señal y adivina qué ..... hace. Pero primero veamos la relación matemática entre el PSD y la función de autocorrelación.

  1. Anotaciones:

    • Transformada de Fourier:
      F[x(t)]=X(ω)=x(t)ejωtdt
    • (Tiempo) Función de autocorrelación:
      R(τ)=x(τ)x(τ)=x(t)x(t+τ)dt
  2. Probemos que la Transformada de Fourier de la función de Auto-Correlación es, de hecho, igual a la densidad de Potencia Espectral de nuestra señal de señal estocástica .X(t)

= - - x ( t ) x ( t + τ ) e - j ω τ d t d τ = - x ( t )

F[R(τ)]=-R(τ)mi-jωτreτ
=--X(t)X(t+τ)mi-jωτretreτ
=-X(t)-X(t+τ)mi-jωτreτF[X(t+τ)]=X(ω)mijωtret
=X(ω)-X(t)mijωtret

=X(ω)X(ω)=El |X(ω)El |2

Que significa todo esto? Nota: esta explicación es un poco "hacky". Pero aquí va

F[X(t)]

¿Qué sucede si tomas el valor esperado de la transformada de Fourier entonces? Esto no funcionaria. Tomemos una señal media cero, por ejemplo.

mi{F[X(t)]}=F[mi{X(t)}]=0 0

En cambio, ¿qué pasa si tomas la transformada de Fourier del cuadrado de la señal?

mi{F[X2(t)]}=F[mi{X2(t)}AV. Poder de la señal]

La función de autocorrelación es esencialmente la PAG(t) a lo que te referías

Referencias

[1] Comunicaciones 1, PL. Dragotti, Imperial College London

[2] Ruido blanco y estimación, F. Tobar [Informe no publicado]


¡Fantástica explicación! Una pequeña pregunta de cálculo: ¿puede intercambiar elret y el reτdentro de las integrales dobles, solo porque sus límites son ambos de: a +?
Spacey

si eso es correcto
ssk08

De acuerdo, creo que lo entiendo. Puedo ver cómo la transformación de Fourier está relacionada con la autocorrelación. Sin embargo, no entiendo realmente cuál es el problema con la transformación de Fourier deX(t) o X2(t). Realmente no veo por qué es necesario tomar el valor esperado (sé que lo promedia, pero no sé por qué es necesario) y realmente no entiendo lo que quieres decir con 'para cada realización de lo aleatorio proceso, tendrás diferentes expresiones para '. Si pudieras elaborar un poco, ¡sería genial! ¡Gracias por tu tiempo!
user968243

1
@ user968243 En cuanto a la parte "para cada realización", piense de esta manera: su señal original, digamos longitud norte, para el que desea encontrar el PSD, es un vector aleatorio. Entonces es un vector connortecomponentes. Ahora, dado que este es un vector aleatorio, cada vez que 'tira el dado', obtiene diferentes valores para sus componentes. Una posibilidad podría ser [3 4 1 9 ...]. Otra posibilidad podría ser [2.9 4.2 1.1 9.02 ...]. Esto es lo que quiere decir cuando dice: "Por cada realización de un proceso aleatorio, (su vector), obtiene diferentes expresiones para" (la transformación de Fourier. ¿Tiene sentido?)
Spacey

@Mohammad lo resumió perfectamente.
ssk08

6

Buena derivación, pero creo que puedes hacerlo aún más fácilmente.

Correlación automática r(t)=X(t)X(-t), es la convolución de la señal con su tiempo invertido.

La convolución en el dominio del tiempo es la multiplicación en el dominio de la frecuencia.

El cambio de tiempo en el dominio del tiempo es "conjugado complejo" en el dominio de la frecuencia.

De ahí obtenemos

R(ω)=F{r(t)}=F{X(t)}F{X(-t)}=X(ω)X(ω)=El |X(ω)El |2=PAGSre

¿No es la autocorrelación la convolución de la señal con su complejo conjugado, auto-invertido?
Jim Clay

Creo que está asumiendo que la señal es real.
ssk08

1
@ Jim y ssk08: ambos tienen razón, por supuesto. Gracias por limpiar las ecuaciones.
Hilmar
Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.