Promedio de dominio de tiempo FFT vs promedio de frecuencia bin

12

Tengo múltiples ensayos de datos fisiológicos. Estoy haciendo un análisis basado en la frecuencia para analizar la potencia (amplitud) en ciertas frecuencias de interés. ¿Promediar múltiples pruebas de igual longitud y luego tomar una sola FFT de la señal promediada versus calcular FFT para cada prueba y luego promediar los intervalos de frecuencia de la misma manera? En la práctica, creo que este no es el caso.

Específicamente, la señal naturalmente tiene un fuerte componente 1 / f y esto se enfatiza si calculo el FFT de cada prueba individual y luego promedio las amplitudes (parte real) de cada bin de frecuencia. ¿Son los dos equivalentes? ¿Hay una manera correcta de hacer las cosas? ¿o bajo qué condiciones de principio debería hacerse la elección entre el promedio del dominio del tiempo frente al promedio del bin de frecuencia?


fuente

Respuestas:

8

Déjame aclarar.

  • La transformada de Fourier no representa el histograma de la señal. La transformación de Fourier es una transformación lineal que toma la señal del dominio del tiempo (función compleja) al dominio de frecuencia (otra función compleja). Lleva una función compleja a otra función compleja.
  • La transformada de Fourier es lineal como se señala en el póster anterior.
  • Fase en sus muestras importa como se señaló anteriormente. Si los datos de prueba por prueba varían en fase, entonces no desea promediar antes de hacer una transformación de Fourier, pero tampoco desea promediar después de la transformación de Fourier. Desea promediar después de la transformación y la norma de Fourier. A continuación detallaré exactamente lo que hay que hacer.

El problema principal aquí es que la pregunta se plantea mal. No es "debería tomar la transformada de Fourier antes de promediar o después de promediar". Porque no hace diferencia debido a la linealidad de la transformada de Fourier.

La pregunta correcta es "debería tomar la amplitud de la transformada de Fourier antes de promediar o después de promediar". Para esta pregunta, la respuesta es antes.

Aquí están los detalles.

Suponga que sus datos muestreados están representados por las secuencias:

d1=d1[n1],d1[n2],...d1[nN]

d2=d2[n1],d2[n2],...d2[nN]

d3=d3[n1],d3[n2],...d3[nN]

...

dM=dM[n1],dM[n2],...dM[nN]

donde son datos de M ensayos y son puntos de tiempo muestreados, luego:d1,...dMn1,...nN

F1=j=1M|F{dj}||F{j=1Mdj}|=F2

Entonces, mientras que la transformación es lineal,no es.F|F|

Además, mientras que es real para todo , no lo es, peroes.dj[ni]i,jF{dj}|F{dj}|

En cuanto a lo que debe hacer, debe tomar la transformada de Fourier de las pruebas individuales (a través de FFT), obtener la amplitud de las pruebas individuales y promediarlas juntas.

Finalmente, qué es . es un término corto para el espectro de frecuencia de las señales "naturales" (generalmente las personas piensan en imágenes).1/f1/f

Cuando la gente dice que hay un componente grande , significa que la amplitud en función de la frecuencia se parece a . Es totalmente ondulado ... probablemente proveniente de un biólogo: p1/f1/f

La transformada inversa de Fourier de es alguna función de signo, pero eso es inútil. ¡Es una función de signo imaginario! Las funciones reales generan transformada simétrica de Fourier.1/f

De hecho, decir que el espectro es , le dice algo sobre la señal, pero no le permite recuperar la señal. Todo lo que sabes es que. Esto no le permite determinar de forma única porque toda la información de fase se ha ido, y sabemos que la estructura de una señal depende en gran medida de su fase .1/f|F{x(t)}|=|1/f|x(t)

¿Qué te dice ? Simplemente que contiene mucha frecuencia baja y poca frecuencia alta.1/f

Una pregunta tan importante, ¿qué te compra el promedio? y lo más importante es cómo interpretar el resultado? Sintonice mañana para una discusión más profunda: p

thang
fuente
1
+1 Gracias por aclarar. Siento que entendí mal el problema subyacente detrás de la pregunta, y creo que esto se vuelve más importante.
jstarr
2

Primero, el FFT es un algoritmo. ¡La transformación se llama Transformada de Fourier! Representa el histograma de las señales. En el caso discreto, una lectura alta en los dominios de frecuencia significa mucha energía a esa frecuencia.

No debe promediar los datos antes de la FFT ya que la información de fase causará cambios significativos en los datos.

Imagine 2 muestras, cada una de las cuales consiste en un coseno puro. En el mundo real, nunca capturará este coseno exactamente en el mismo punto de partida. Un coseno se desplazará reletivo a otro (o ambos tienen diferentes turnos reletivos al comienzo. Matemáticamente esto significa y1 = cos (wt-A) y2 = cos (wt-B) donde A y B son desplazamientos. En su modelo estos dos aparecen mejor como lo mismo. Con un poco de matemática puedo elegir estos valores para que y2-y1 = 0. El promedio de cero es cero y no es lo que quieres. Este es el problema de la fase.

Si su objetivo es encontrar el espectro promedio que debe promediar a través de los espectros, ¡no promedie las señales!

Mikhail
fuente
Gracias por eso. Ciertamente fue una preocupación que si hay ruido con la diferencia de fase correcta de un ensayo a otro, cancelará todas las señales o las frecuencias de interés. Todavía no estoy claro por qué el 1 / f se acentúa por el promedio del bin de frecuencia en lugar del promedio del dominio del tiempo seguido de FFT.
@ user1487551 ¿Qué significa un componente 1 / f fuerte? La transformación inversa de Fourier de 1 / f es una función de signo y puede indicar que sus datos incluyen una gran parte cuando el sistema se está estabilizando. Probablemente debería mostrar una trama o algunos datos.
Mikhail
1
Tenga en cuenta que puede hacer un promedio de tiempo si tiene una referencia de fase con la que puede sincronizar la captura de datos, y esto puede ser muy efectivo para reducir el ruido de fondo, pero de lo contrario es cierto que necesita hacer un promedio de conjunto en el dominio de frecuencia.
Paul R
2

A menos que esté completamente fuera de lugar o malinterprete su pregunta, la respuesta es : por la linealidad de la DFT, promediar las señales a tiempo y luego tomar la DFT del promedio es equivalente a promediar las DFT de las señales.

Para mostrar esto, definamos algunas variables:

  • xn[] : -muestra de dominio de tiempo de prueba en el tiempothn
  • Xk[] : -muestra de dominio de frecuencia auricular a frecuenciathk

La señal "promedio" en el dominio del tiempo viene dada por . Tomando su DFT, tenemos1L=0Lxn[]

n=0N11LLxn[]ei2πkn/N .

Cambiando el orden de las sumas, podemos escribir

1L=0Ln=0N1xn[]ei2πkn/N,

pero esto es lo mismo que

1L=0LXk[l]

que es lo mismo que promediar los DFT de cada trival. Esto es lo que queríamos mostrar.

jstarr
fuente
Para promediar en el dominio del tiempo, necesita algún tipo de referencia de fase a la que pueda sincronizar. Cuando esto está disponible, aunque a veces se pueden obtener mejores resultados que con el promedio de conjunto en el dominio de la frecuencia, debido principalmente a la reducción de los efectos de la cuantificación, etc.
Paul R