análisis de Fourier -Fourier analysis

Señal de tiempo de guitarra baja de la nota A de cuerda abierta (55 Hz).
Transformada de Fourier de la señal de tiempo del bajo de la nota A de cuerda abierta (55 Hz). El análisis de Fourier revela los componentes oscilatorios de señales y funciones .

En matemáticas , el análisis de Fourier ( / ˈ f ʊr i , - i ər / ) es el estudio de la forma en que las funciones generales pueden representarse o aproximarse mediante sumas de funciones trigonométricas más simples . El análisis de Fourier surgió del estudio de las series de Fourier y lleva el nombre de Joseph Fourier , quien demostró que representar una función como una suma de funciones trigonométricas simplifica enormemente el estudio de la transferencia de calor .

El tema del análisis de Fourier abarca un amplio espectro de las matemáticas. En las ciencias y la ingeniería, el proceso de descomposición de una función en componentes oscilatorios a menudo se denomina análisis de Fourier, mientras que la operación de reconstrucción de la función a partir de estas piezas se conoce como síntesis de Fourier . Por ejemplo, determinar qué frecuencias componentes están presentes en una nota musical implicaría calcular la transformada de Fourier de una nota musical muestreada. Luego, se podría volver a sintetizar el mismo sonido al incluir los componentes de frecuencia como se revela en el análisis de Fourier. En matemáticas, el término análisis de Fourier a menudo se refiere al estudio de ambas operaciones.

El proceso de descomposición en sí mismo se llama transformación de Fourier . Su salida, la transformada de Fourier , a menudo recibe un nombre más específico, que depende del dominio y otras propiedades de la función que se transforma. Además, el concepto original del análisis de Fourier se ha ampliado con el tiempo para aplicarse a situaciones cada vez más abstractas y generales, y el campo general se conoce a menudo como análisis armónico . Cada transformada utilizada para el análisis (consulte la lista de transformadas relacionadas con Fourier ) tiene una transformada inversa correspondiente que se puede utilizar para la síntesis.

Para utilizar el análisis de Fourier, los datos deben estar igualmente espaciados. Se han desarrollado diferentes enfoques para analizar datos espaciados de manera desigual, en particular los métodos de análisis espectral de mínimos cuadrados (LSSA) que utilizan un ajuste de mínimos cuadrados de sinusoides a muestras de datos, similar al análisis de Fourier. El análisis de Fourier, el método espectral más utilizado en la ciencia, generalmente aumenta el ruido de período largo en registros con intervalos prolongados; LSSA mitiga tales problemas.

Aplicaciones

El análisis de Fourier tiene muchas aplicaciones científicas: en física , ecuaciones diferenciales parciales , teoría de números , combinatoria , procesamiento de señales , procesamiento de imágenes digitales , teoría de probabilidad , estadística , ciencia forense , fijación de precios de opciones , criptografía , análisis numérico , acústica , oceanografía , sonar , óptica , difracción . , geometría , análisis de estructura de proteínas y otras áreas.

Esta amplia aplicabilidad se deriva de muchas propiedades útiles de las transformadas:

En medicina forense, los espectrofotómetros infrarrojos de laboratorio utilizan el análisis de transformada de Fourier para medir las longitudes de onda de la luz a las que un material absorberá en el espectro infrarrojo. El método FT se utiliza para decodificar las señales medidas y registrar los datos de longitud de onda. Y mediante el uso de una computadora, estos cálculos de Fourier se llevan a cabo rápidamente, de modo que en cuestión de segundos, un instrumento FT-IR operado por computadora puede producir un patrón de absorción de infrarrojos comparable al de un instrumento de prisma.

La transformación de Fourier también es útil como representación compacta de una señal. Por ejemplo, la compresión JPEG utiliza una variante de la transformación de Fourier ( transformada de coseno discreta ) de pequeñas piezas cuadradas de una imagen digital. Los componentes de Fourier de cada cuadrado se redondean a una precisión aritmética más baja y los componentes débiles se eliminan por completo, de modo que los componentes restantes se pueden almacenar de forma muy compacta. En la reconstrucción de imágenes, cada cuadrado de la imagen se vuelve a ensamblar a partir de los componentes transformados de Fourier aproximados conservados, que luego se transforman inversamente para producir una aproximación de la imagen original.

En el procesamiento de señales , la transformada de Fourier a menudo toma una serie de tiempo o una función de tiempo continuo y la mapea en un espectro de frecuencia . Es decir, lleva una función del dominio del tiempo al dominio de la frecuencia ; es una descomposición de una función en sinusoides de diferentes frecuencias; en el caso de una serie de Fourier o una transformada discreta de Fourier , las sinusoides son armónicos de la frecuencia fundamental de la función que se analiza.

Cuando una función es una función del tiempo y representa una señal física , la transformada tiene una interpretación estándar como el espectro de frecuencia de la señal. La magnitud de la función de valor complejo resultante en la frecuencia representa la amplitud de un componente de frecuencia cuya fase inicial viene dada por el ángulo de (coordenadas polares).

Las transformadas de Fourier no se limitan a funciones de tiempo y frecuencias temporales. Se pueden aplicar igualmente para analizar frecuencias espaciales y, de hecho, para casi cualquier dominio de función. Esto justifica su uso en ramas tan diversas como el procesamiento de imágenes , la conducción de calor y el control automático .

Al procesar señales, como audio , ondas de radio , ondas de luz, ondas sísmicas e incluso imágenes, el análisis de Fourier puede aislar los componentes de banda estrecha de una forma de onda compuesta, concentrándolos para detectarlos o eliminarlos más fácilmente. Una gran familia de técnicas de procesamiento de señales consiste en la transformación de Fourier de una señal, la manipulación de los datos transformados por Fourier de manera sencilla y la inversión de la transformación.

Algunos ejemplos incluyen:

Variantes del análisis de Fourier

Una transformada de Fourier y 3 variaciones causadas por muestreo periódico (en el intervalo T) y/o suma periódica (en el intervalo P) de la función subyacente en el dominio del tiempo. La relativa facilidad de cálculo de la secuencia DFT y la comprensión que brinda de S ( f ) la convierten en una herramienta de análisis popular.

(Continua) Transformada de Fourier

La mayoría de las veces, el término sin calificar transformada de Fourier se refiere a la transformada de funciones de un argumento real continuo , y produce una función continua de frecuencia, conocida como distribución de frecuencia . Una función se transforma en otra y la operación es reversible. Cuando el dominio de la función de entrada (inicial) es el tiempo ( t ), y el dominio de la función de salida (final) es la frecuencia ordinaria , la transformada de la función s ( t ) en la frecuencia f viene dada por el número complejo:

La evaluación de esta cantidad para todos los valores de f produce la función de dominio de frecuencia . Entonces s ( t ) se puede representar como una recombinación de exponenciales complejas de todas las frecuencias posibles:

que es la fórmula de la transformada inversa. El número complejo, S ( f ) , transmite tanto la amplitud como la fase de la frecuencia f .

Consulte la transformada de Fourier para obtener mucha más información, que incluye:

  • Convenciones para la normalización de amplitud y escalado/unidades de frecuencia
  • transformar propiedades
  • transformadas tabuladas de funciones específicas
  • una extensión/generalización para funciones de múltiples dimensiones, como imágenes.

series de Fourier

La transformada de Fourier de una función periódica, s P ( t ) , con período P , se convierte en una función peine de Dirac , modulada por una secuencia de coeficientes complejos :

    (donde P es la integral sobre cualquier intervalo de longitud P ).

La transformada inversa, conocida como serie de Fourier , es una representación de s P ( t ) en términos de una suma de un número potencialmente infinito de sinusoides relacionados armónicamente o funciones exponenciales complejas , cada una con una amplitud y una fase especificadas por uno de los coeficientes:

Cualquier s P ( t ) se puede expresar como una suma periódica de otra función, s ( t ) :

y los coeficientes son proporcionales a muestras de S ( f ) a intervalos discretos de 1/PAG:

Tenga en cuenta que cualquier s ( t ) cuya transformada tenga los mismos valores de muestra discretos se puede usar en la suma periódica. Una condición suficiente para recuperar s ( t ) (y por lo tanto S ( f ) ) solo de estas muestras (es decir, de la serie de Fourier) es que la parte distinta de cero de s ( t ) se limite a un intervalo conocido de duración P , que es el dominio de frecuencia dual del teorema de muestreo de Nyquist-Shannon .

Consulte la serie de Fourier para obtener más información, incluido el desarrollo histórico.

Transformada de Fourier en tiempo discreto (DTFT)

La DTFT es el dual matemático de la serie de Fourier en el dominio del tiempo. Por lo tanto, una suma periódica convergente en el dominio de la frecuencia se puede representar mediante una serie de Fourier, cuyos coeficientes son muestras de una función de tiempo continua relacionada:

que se conoce como DTFT. Por lo tanto, la DTFT de la secuencia s [ n ] es también la transformada de Fourier de la función de peine de Dirac modulada .

Los coeficientes de la serie de Fourier (y la transformada inversa), están definidos por:

El parámetro T corresponde al intervalo de muestreo, y esta serie de Fourier ahora puede reconocerse como una forma de la fórmula de suma de Poisson . Así tenemos el importante resultado de que cuando una secuencia de datos discretos, s [ n ] , es proporcional a las muestras de una función continua subyacente, s ( t ) , se puede observar una suma periódica de la transformada continua de Fourier, S ( f ) . Tenga en cuenta que cualquier s ( t ) con los mismos valores de muestra discretos produce la misma DTFT Pero bajo ciertas condiciones idealizadas, teóricamente se puede recuperar S ( f ) y s ( t ) exactamente. Una condición suficiente para la recuperación perfecta es que la parte distinta de cero de S ( f ) se limite a un intervalo de frecuencia conocido de ancho1/T. Cuando ese intervalo es [−1/2 toneladas,1/2 toneladas] , la fórmula de reconstrucción aplicable es la fórmula de interpolación de Whittaker-Shannon . Esta es una piedra angular en la base del procesamiento de señales digitales .

Otra razón para estar interesado en S 1/ T ( f ) es que a menudo proporciona información sobre la cantidad de alias causado por el proceso de muestreo.

Las aplicaciones de la DTFT no se limitan a funciones muestreadas. Consulte Transformada de Fourier en tiempo discreto para obtener más información sobre este y otros temas, incluidos:

  • unidades de frecuencia normalizadas
  • ventanas (secuencias de longitud finita)
  • transformar propiedades
  • transformadas tabuladas de funciones específicas

Transformada discreta de Fourier (DFT)

Similar a una serie de Fourier, la DTFT de una secuencia periódica, con período , se convierte en una función peine de Dirac, modulada por una secuencia de coeficientes complejos (ver DTFT § Datos periódicos ):

    (donde Σ n es la suma de cualquier secuencia de longitud N ).

La secuencia S [ k ] es lo que habitualmente se conoce como la DFT de un ciclo de sN . También es N -periódico, por lo que nunca es necesario calcular más de N coeficientes. La transformada inversa, también conocida como serie discreta de Fourier , viene dada por:

  donde Σ k es la suma de cualquier secuencia de longitud N .

Cuando s N [ n ] se expresa como una suma periódica de otra función:

  y  

los coeficientes son proporcionales a muestras de S 1/ T ( f ) a intervalos discretos de1/PAG=1/Nuevo Testamento:

Por el contrario, cuando se quiere calcular un número arbitrario ( N ) de muestras discretas de un ciclo de una DTFT continua, S 1/ T ( f ) , se puede hacer calculando la DFT relativamente simple de s N [ n ] , como definido anteriormente. En la mayoría de los casos, N se elige igual a la longitud de la parte distinta de cero de s [ n ] . El aumento de N , conocido como relleno de ceros o interpolación , da como resultado muestras más próximas entre sí de un ciclo de S 1/ T ( f ) . La disminución de N provoca superposición (adición) en el dominio del tiempo (análogo al aliasing ), que corresponde a la aniquilación en el dominio de la frecuencia. (ver Transformada de Fourier de tiempo discreto § L=N×I ) En la mayoría de los casos de interés práctico, la secuencia s [ n ] representa una secuencia más larga que se truncó mediante la aplicación de una función de ventana de longitud finita o una matriz de filtro FIR .

La DFT se puede calcular utilizando un algoritmo de transformada rápida de Fourier (FFT), lo que la convierte en una transformación práctica e importante en las computadoras.

Consulte Transformada discreta de Fourier para obtener mucha más información, que incluye:

  • transformar propiedades
  • aplicaciones
  • transformadas tabuladas de funciones específicas

Resumen

Para funciones periódicas, tanto la transformada de Fourier como la DTFT comprenden solo un conjunto discreto de componentes de frecuencia (serie de Fourier), y las transformadas divergen en esas frecuencias. Una práctica común (no discutida anteriormente) es manejar esa divergencia a través de las funciones delta de Dirac y peine de Dirac . Pero la misma información espectral se puede discernir de un solo ciclo de la función periódica, ya que todos los demás ciclos son idénticos. De manera similar, las funciones de duración finita se pueden representar como una serie de Fourier, sin pérdida real de información, excepto que la periodicidad de la transformada inversa es un mero artefacto.

Es común en la práctica que la duración de s (•) se limite al período, P o N . Pero estas fórmulas no exigen esa condición.

s ( t ) se transforma (tiempo continuo)
Frecuencia continua frecuencias discretas
Transformar
Inverso
s ( nT ) transforma (tiempo discreto)
Frecuencia continua frecuencias discretas
Transformar

Inverso

Propiedades de simetría

Cuando las partes real e imaginaria de una función compleja se descomponen en sus partes pares e impares , hay cuatro componentes, indicados a continuación por los subíndices RE, RO, IE e IO. Y hay un mapeo uno a uno entre los cuatro componentes de una función de tiempo compleja y los cuatro componentes de su transformada de frecuencia compleja:

A partir de esto, varias relaciones son evidentes, por ejemplo:

  • La transformada de una función de valor real ( s RE + s RO ) es la función simétrica par S RE + i S IO . Por el contrario, una transformación incluso simétrica implica un dominio de tiempo de valor real.
  • La transformada de una función con valores imaginarios ( i s IE + i s IO ) es la función simétrica impar S RO + i S IE , y lo contrario es cierto.
  • La transformada de una función incluso simétrica ( s RE + i s IO ) es la función de valor real S RE + S RO , y lo contrario es cierto.
  • La transformada de una función simétrica impar ( s RO + i s IE ) es la función de valor imaginario i S IE + i S IO , y lo contrario es cierto.

Historia

Una forma temprana de series armónicas se remonta a las antiguas matemáticas babilónicas , donde se usaban para calcular efemérides (tablas de posiciones astronómicas).

Los conceptos griegos clásicos de deferente y epiciclo en el sistema ptolemaico de astronomía estaban relacionados con las series de Fourier (ver Deferente y epiciclo § Formalismo matemático ).

En tiempos modernos, Alexis Clairaut utilizó variantes de la transformada discreta de Fourier en 1754 para calcular una órbita, que se ha descrito como la primera fórmula para la DFT, y en 1759 Joseph Louis Lagrange , para calcular los coeficientes de una serie trigonométrica. para una cuerda vibrante. Técnicamente, el trabajo de Clairaut era una serie de solo coseno (una forma de transformada de coseno discreta ), mientras que el trabajo de Lagrange era una serie de solo seno (una forma de transformada de seno discreta ); Gauss utilizó en 1805 una verdadera DFT coseno+seno para la interpolación trigonométrica de las órbitas de los asteroides . Tanto Euler como Lagrange discretizaron el problema de la cuerda vibrante, usando lo que hoy se llamaría muestras.

Un desarrollo moderno temprano hacia el análisis de Fourier fue el artículo de 1770 Réflexions sur la résolution algébrique des équations de Lagrange, que en el método de los resolventes de Lagrange utilizó una descomposición compleja de Fourier para estudiar la solución de un cúbico: Lagrange transformó las raíces x 1 , x 2 , x 3 en los resolventes:

donde ζ es una raíz cúbica de la unidad , que es la DFT de orden 3.

Varios autores, en particular Jean le Rond d'Alembert y Carl Friedrich Gauss , utilizaron series trigonométricas para estudiar la ecuación del calor , pero el avance más importante fue el artículo de 1807 Mémoire sur la propagation de la chaleur dans les corps solides de Joseph Fourier , cuyo La idea crucial fue modelar todas las funciones mediante series trigonométricas, introduciendo la serie de Fourier.

Los historiadores están divididos en cuanto a cuánto darle crédito a Lagrange y a otros por el desarrollo de la teoría de Fourier: Daniel Bernoulli y Leonhard Euler habían introducido representaciones trigonométricas de funciones, y Lagrange había dado la solución de la serie de Fourier a la ecuación de onda, por lo que la contribución de Fourier fue principalmente la audaz afirmación de que una función arbitraria podría ser representada por una serie de Fourier.

El desarrollo subsiguiente del campo se conoce como análisis armónico , y es también un ejemplo temprano de la teoría de la representación .

El primer algoritmo de transformada rápida de Fourier (FFT) para la DFT fue descubierto alrededor de 1805 por Carl Friedrich Gauss al interpolar medidas de la órbita de los asteroides Juno y Pallas , aunque ese algoritmo de FFT en particular se atribuye más a menudo a sus redescubridores modernos Cooley y Tukey .

Transformadas de tiempo-frecuencia

En términos de procesamiento de señales , una función (de tiempo) es una representación de una señal con resolución de tiempo perfecta , pero sin información de frecuencia, mientras que la transformada de Fourier tiene resolución de frecuencia perfecta , pero sin información de tiempo.

Como alternativas a la transformada de Fourier, en el análisis de tiempo-frecuencia , se utilizan transformadas de tiempo-frecuencia para representar señales en una forma que tiene alguna información de tiempo y alguna información de frecuencia; por el principio de incertidumbre , existe un compromiso entre estos. Estas pueden ser generalizaciones de la transformada de Fourier, como la transformada de Fourier de tiempo corto , la transformada de Gabor o la transformada fraccionada de Fourier (FRFT), o pueden usar diferentes funciones para representar señales, como en transformadas wavelet y transformadas chirplet , con la wavelet analógica de la transformada de Fourier (continua) siendo la transformada wavelet continua .

Transformadas de Fourier en grupos topológicos abelianos localmente compactos arbitrarios

Las variantes de Fourier también se pueden generalizar a transformadas de Fourier en grupos topológicos abelianos localmente compactos arbitrarios , que se estudian en el análisis armónico ; allí, la transformada de Fourier lleva funciones en un grupo a funciones en el grupo dual. Este tratamiento también permite una formulación general del teorema de convolución , que relaciona las transformadas de Fourier y las convoluciones . Consulte también la dualidad de Pontryagin para conocer los fundamentos generalizados de la transformada de Fourier.

Más específico, el análisis de Fourier se puede realizar en clases laterales, incluso en clases laterales discretas.

Ver también

notas

Referencias

Otras lecturas

enlaces externos