En matemáticas , el análisis de Fourier ( / ˈfʊrieɪ , -iər / ) [ 1] 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 recibe su 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 se suele denominar análisis de Fourier, mientras que la operación de reconstruir la función a partir de estos fragmentos 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. A continuación, se podría volver a sintetizar el mismo sonido incluyendo los componentes de frecuencia tal como se revelan en el análisis de Fourier. En matemáticas, el término análisis de Fourier suele referirse al estudio de ambas operaciones.
El proceso de descomposición en sí se denomina transformación de Fourier . Su resultado, la transformada de Fourier , suele recibir 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 suele conocerse 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 espaciados de manera uniforme. 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 las sinusoides a las muestras de datos, similar al análisis de Fourier. [2] [3] 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 espacios largos; LSSA mitiga estos problemas. [4]
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 transformaciones :
En el campo 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 infrarroja comparable al de un instrumento de prisma. [9]
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ños fragmentos cuadrados de una imagen digital. Los componentes de Fourier de cada cuadrado se redondean a una precisión aritmética menor y se eliminan los componentes débiles, 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 aproximados conservados de la transformada de Fourier, que luego se transforman de forma inversa para producir una aproximación de la imagen original.
En el procesamiento de señales , la transformada de Fourier a menudo toma una serie temporal o una función de tiempo continuo y la convierte en un espectro de frecuencias . Es decir, toma 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 de Fourier discreta , 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 transformación tiene una interpretación estándar como el espectro de frecuencia de la señal. La magnitud de la función compleja resultante en la frecuencia representa la amplitud de un componente de frecuencia cuya fase inicial está dada por el ángulo de (coordenadas polares).
Las transformadas de Fourier no se limitan a funciones de tiempo y frecuencias temporales. También se pueden aplicar para analizar frecuencias espaciales y, de hecho, para casi cualquier dominio de funciones. 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 una detección o eliminación más sencilla. Una gran familia de técnicas de procesamiento de señales consiste en transformar una señal mediante Fourier, manipular los datos transformados mediante Fourier de una manera sencilla e invertir la transformación. [10]
Algunos ejemplos incluyen :
En la mayoría de los casos, el término no calificado 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 ( ), y el dominio de la función de salida (final) es la frecuencia ordinaria , la transformada de la función en la frecuencia viene dada por el número complejo :
Al evaluar esta cantidad para todos los valores de se obtiene la función en el dominio de la frecuencia . Luego se puede representar como una recombinación de exponenciales complejos de todas las frecuencias posibles :
que es la fórmula de transformación inversa. El número complejo transmite tanto la amplitud como la fase de la frecuencia.
Consulte la transformada de Fourier para obtener mucha más información, incluyendo :
La transformada de Fourier de una función periódica, con período, se convierte en una función peine de Dirac , modulada por una secuencia de coeficientes complejos :
La transformada inversa, conocida como serie de Fourier , es una representación en términos de una suma de un número potencialmente infinito de senos armónicamente relacionados o funciones exponenciales complejas , cada una con una amplitud y fase especificadas por uno de los coeficientes :
Cualquiera puede expresarse como una suma periódica de otra función :
y los coeficientes son proporcionales a muestras de en intervalos discretos de :
Nótese que cualquier cuya transformada tenga los mismos valores de muestra discretos puede usarse en la suma periódica. Una condición suficiente para recuperar (y por lo tanto ) solo a partir de estas muestras (es decir, de la serie de Fourier) es que la porción distinta de cero de esté confinada a un intervalo conocido de duración que es el dual del dominio de frecuencia del teorema de muestreo de Nyquist-Shannon .
Consulte la serie de Fourier para obtener más información, incluido el desarrollo histórico.
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 es también la transformada de Fourier de la función peine de Dirac modulada . [B]
Los coeficientes de la serie de Fourier (y la transformada inversa), se definen por :
El parámetro corresponde al intervalo de muestreo, y esta serie de Fourier ahora puede reconocerse como una forma de la fórmula de suma de Poisson . Por lo tanto, tenemos el resultado importante de que cuando una secuencia de datos discretos es proporcional a las muestras de una función continua subyacente, se puede observar una suma periódica de la transformada de Fourier continua. Nótese que cualquier con los mismos valores de muestra discretos produce la misma DTFT. Pero bajo ciertas condiciones idealizadas, uno puede recuperar teóricamente y exactamente. Una condición suficiente para una recuperación perfecta es que la porción distinta de cero de esté confinada a un intervalo de frecuencia conocido de ancho. Cuando ese intervalo es 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 .
Otro motivo de interés 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 las funciones muestreadas. Consulte Transformada de Fourier de tiempo discreto para obtener más información sobre este y otros temas, incluidos :
De manera 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 ) :
La secuencia se conoce habitualmente como la DFT de un ciclo de También es periódica, por lo que nunca es necesario calcular más de coeficientes. La transformada inversa, también conocida como serie de Fourier discreta , viene dada por :
Cuando se expresa como suma periódica de otra función :
Los coeficientes son muestras de intervalos discretos de :
Por el contrario, cuando se desea calcular un número arbitrario de muestras discretas de un ciclo de una DTFT continua, se puede hacer calculando la DFT relativamente simple de como se definió anteriormente. En la mayoría de los casos, se elige igual a la longitud de la porción distinta de cero de Aumentar, conocido como relleno de ceros o interpolación , da como resultado muestras más espaciadas de un ciclo de Decreciente provoca superposición (adición) en el dominio del tiempo (análogo al aliasing ), que corresponde a la decimació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 representa una secuencia más larga que fue truncada por la aplicación de una función de ventana de longitud finita o una matriz de filtros 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 de Fourier discreta para obtener mucha más información, incluyendo :
En el caso de las funciones periódicas, tanto la transformada de Fourier como la DTFT comprenden únicamente un conjunto discreto de componentes de frecuencia (series de Fourier), y las transformadas divergen en esas frecuencias. Una práctica común (no analizada anteriormente) es manejar esa divergencia mediante las funciones delta de Dirac y peine de Dirac . Pero la misma información espectral se puede discernir a partir 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.
En la práctica, es habitual que la duración de s (•) se limite al período P o N. Pero estas fórmulas no requieren esa condición.
Cuando las partes reales e imaginarias de una función compleja se descomponen en sus partes pares e impares , hay cuatro componentes, denotados a continuación por los subíndices RE, RO, IE e IO. Y hay una correspondencia biunívoca entre los cuatro componentes de una función de tiempo compleja y los cuatro componentes de su transformada de frecuencia compleja : [11]
De esto se desprenden diversas relaciones, por ejemplo :
Una forma temprana de series armónicas se remonta a las matemáticas de la antigua Babilonia , donde se utilizaban para calcular efemérides (tablas de posiciones astronómicas). [12] [13] [14] [15]
Los conceptos griegos clásicos de deferente y epiciclo en el sistema astronómico ptolemaico estaban relacionados con las series de Fourier (ver Deferente y epiciclo § Formalismo matemático ).
En tiempos modernos, variantes de la transformada de Fourier discreta fueron utilizadas por Alexis Clairaut en 1754 para calcular una órbita, [16] que ha sido descrita como la primera fórmula para la DFT, [17] y en 1759 por Joseph Louis Lagrange , para calcular los coeficientes de una serie trigonométrica para una cuerda vibrante. [17] Técnicamente, el trabajo de Clairaut fue una serie de solo coseno (una forma de transformada de coseno discreta ), mientras que el trabajo de Lagrange fue una serie de solo seno (una forma de transformada de seno discreta ); Gauss utilizó una verdadera DFT de coseno + seno en 1805 para la interpolación trigonométrica de órbitas de asteroides . [18] Euler y Lagrange discretizaron el problema de la cuerda vibrante, utilizando lo que hoy se llamarían muestras. [17]
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 de Fourier compleja para estudiar la solución de una ecuación cúbica : [19] Lagrange transformó las raíces 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 , [20] pero el avance decisivo fue el artículo de 1807 Mémoire sur la propagation de la chaleur dans les corps solides de Joseph Fourier , cuya 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 crédito se le debe a Lagrange y 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 afirmación audaz de que una función arbitraria podía representarse mediante una serie de Fourier. [17]
El desarrollo posterior del campo se conoce como análisis armónico y también es 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 mediciones de la órbita de los asteroides Juno y Pallas , aunque ese algoritmo FFT en particular se atribuye más a menudo a sus redescubridores modernos Cooley y Tukey . [18] [16]
En términos de procesamiento de señales , una función (de tiempo) es una representación de una señal con una resolución temporal perfecta , pero sin información de frecuencia, mientras que la transformada de Fourier tiene una 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 equilibrio 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 de Fourier fraccionaria (FRFT), o pueden usar diferentes funciones para representar señales, como en las transformadas wavelet y las transformadas chirplet , siendo el análogo wavelet de la transformada de Fourier (continua) la transformada wavelet continua .
Las variantes de Fourier también pueden generalizarse 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 convierte las funciones de un grupo en funciones del 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 . Véase también la dualidad de Pontryagin para los fundamentos generalizados de la transformada de Fourier.
Más específicamente, el análisis de Fourier se puede realizar en clases laterales, [21] incluso clases laterales discretas.