Análisis de Fourier
En matemáticas, el análisis de Fourier es el estudio de la forma general en que las funciones pueden ser representados o aproximadas por sumas de funciones trigonométricas 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.[1]
Hoy, 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 reveló 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í se llama transformación de Fourier. Su producto resultado, 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 está transformando. Además, el concepto original del análisis de Fourier se ha extendido a lo largo del tiempo para aplicarse a situaciones cada vez más abstractas y generales, y el campo general a menudo se conoce 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.
Usos
El análisis de Fourier tiene muchos usos científicos - en la física, ecuaciones diferenciales parciales, teoría de números, combinatoria, procesamiento de señales, procesamiento digital de imágenes, teoría de la probabilidad, estadística, análisis forense, valoración de opciones, la criptografía, análisis numérico, acústica, oceanografía, el sonar, óptica, la difracción, geometría, análisis de estructuras de proteínas y otras áreas.
Esta amplia aplicabilidad se debe a muchas propiedades útiles de las transformadas:
- Las transformadas son aplicaciones lineales y, con la normalización adecuada, también son unitarias (una propiedad conocida como teorema de Parseval o, más generalmente, como el teorema de Plancherel, y más generalmente a través de la dualidad de Pontryagin).[2]
- Las transformadas suelen ser invertibles.
- Las funciones exponenciales son funciones propias de diferenciación, lo que significa que esta representación transforma ecuaciones diferenciales lineales con coeficientes constantes en algebraicas ordinarias.[3] Por lo tanto, el comportamiento de un sistema lineal invariante en el tiempo se puede analizar en cada frecuencia de forma independiente.
- Según el teorema de convolución, las transformadas de Fourier convierten la complicada operación de convolución en una multiplicación simple, lo que significa que proporcionan una forma eficiente de calcular operaciones basadas en convolución, como la multiplicación de polinomios y la multiplicación de números grandes.[4]
- La versión discreta de la transformada de Fourier (ver más abajo) se puede evaluar rápidamente en computadoras usando algoritmos de transformada rápida de Fourier (FFT).[5]
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 al usar 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.[6]
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 para reducir la precisión aritmética 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 conservados aproximadamente transformados de Fourier, que luego se transforman a la inversa para producir una aproximación de la imagen original.
Uso en procesamiento de señales
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 componentes de banda estrecha de una forma de onda compuesta, concentrándolos para una detección o eliminación más fácil. Una gran familia de técnicas de procesamiento de señales consiste en transformar una señal de Fourier, manipular los datos transformados de Fourier de una manera simple e invertir la transformación.[7]
Algunos ejemplos incluyen:
- Ecualización de grabaciones de audio con una serie de filtros de paso de banda ;
- Recepción de radio digital sin un circuito superheterodino, como en un teléfono celular moderno o un escáner de radio;
- Procesamiento de imágenes para eliminar artefactos periódicos o anisotrópicos tales como irregularidades de video entrelazado, artefactos de bandas de fotografías aéreas de bandas o patrones de ondas de interferencias de radiofrecuencia en una cámara digital;
- Correlación cruzada de imágenes similares para la co-alineación;
- Cristalografía de rayos X para reconstruir una estructura cristalina a partir de su patrón de difracción;
- Espectrometría de resonancia ciclotrónica por transformada de Fourier para determinar la masa de iones a partir de la frecuencia del movimiento del ciclotrón en un campo magnético;
- Muchas otras formas de espectroscopia, incluyendo infrarrojos y de resonancia magnética nuclear espectroscopias;
- Generación de espectrogramas de sonido utilizados para analizar sonidos;
- Sonar pasivo utilizado para clasificar objetivos según el ruido de la maquinaria.
Variantes del análisis de Fourier
Transformada de Fourier (continua)
En la mayoría de los casos, el término sin calificar transformación de Fourier se refiere a la transformación de funciones de un argumento continuo real, y produce una función continua de frecuencia, conocida como distribución de frecuencias. 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) a la frecuencia f viene dada por el número complejo:
Evaluando esta cantidad para todos los valores de f se obtiene la función dominio de la frecuencia. Entonces s(t) puede representarse como una recombinación de exponenciales complejas de todas las frecuencias posibles:
que es la fórmula de la transformación inversa. El número complejo, S(f), transmite tanto la amplitud como la fase de la frecuencia f.
Ver Transformación de Fourier para mucha más información, incluyendo:
- convenciones para la normalización de la amplitud y el escalado de la frecuencia/unidades.
- propiedades de la transformada
- transformaciones tabuladas de funciones específicas
- una extensión/generalización para funciones de múltiples dimensiones, como las imágenes.
Series de Fourier
La transformada de Fourier de una función periódica, sP(t), con periodo 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 sP(t) en términos de una suma de un número potencialmente infinito de sinusoides o funciones exponenciales complejas relacionadas armónicamente, cada una con una amplitud y una fase especificadas por uno de los coeficientes:
Cualquier sP(t) puede expresarse como una suma periódica de otra función, s(t):
y los coeficientes son proporcionales a las muestras de S(f) en intervalos discretos de 1P:
- .A
donde A =
Obsérvese que cualquier s(t) cuya transformada tenga los mismos valores muestrales discretos puede utilizarse en el sumatorio periódico. Una condición suficiente para recuperar s(t) (y por lo tanto S(f)) a partir de sólo estas muestras (es decir, de la serie de Fourier) es que la porción no nula de s(t) esté confinada a un intervalo conocido de duración P, que es el dual del dominio de la frecuencia del teorema de muestreo de Nyquist-Shannon.
Transformada de Fourier en tiempo discreto
La DTFT es el dual matemático de la serie de Fourier en el dominio del tiempo. Así, una suma periódica convergente en el dominio de la frecuencia puede representarse mediante una serie de Fourier, cuyos coeficientes son muestras de una función temporal continua relacionada:
que se conoce como la DTFT. Así, la DTFT de la secuencia s[n] es también la transformada de Fourier de la función peine de Dirac modulada.
También se puede señalar que:
En consecuencia, una práctica común es modelar el "muestreo" como una multiplicación por la función peine de Dirac, que por supuesto sólo es "posible" en un sentido puramente matemático.
Los coeficientes de la serie de Fourier (y la transformada inversa), se definen por:
El parámetro T corresponde al intervalo de muestreo, y esta serie de Fourier puede reconocerse ahora 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 un sumatorio periódico de la transformada continua de Fourier, S(f). Nótese que cualquier s(t) con los mismos valores discretos de la muestra produce la misma DTFT Pero bajo ciertas condiciones idealizadas uno puede recuperar teóricamente S(f) y s(t) exactamente. Una condición suficiente para la recuperación perfecta es que la porción no nula de S(f) esté confinada a un intervalo de frecuencia conocido de ancho {sfrac}}. Cuando ese intervalo es [−12T, 12T], la fórmula de reconstrucción aplicable es la Fórmula de Interpolación de Whittaker-Shannon. Esta es una piedra angular en los fundamentos del procesamiento digital de señales.
Otra razón para estar interesado en SPlantilla:Sub(f) es que a menudo proporciona una visión de la cantidad de aliasing causado por el proceso de muestreo.
Las aplicaciones de la DTFT no se limitan a las funciones muestreadas. Ver Transformada de Fourier en tiempo discreto para más información sobre este y otros temas, incluyendo:
- unidades de frecuencia normalizadas
- Ventana (secuencias de longitud finita)
- propiedades de la transformada
- transformaciones tabuladas de funciones específicas
Transformada discreta de Fourier (DFT)
De forma similar a una serie de Fourier, la DTFT de una secuencia periódica, sN[n], con periodo N, se convierte en una función de peine de Dirac, modulada por una secuencia de coeficientes complejos:
- (donde Σn es la suma sobre cualquier secuencia de longitud N).
La secuencia S[k] es lo que se conoce habitualmente 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 sobre cualquier secuencia de longitud N.
Cuando sN[n] se expresa como una suma periódica de otra función:
- y Nótese que esta definición difiere intencionalmente de la sección DTFT por un factor de T. Esto facilita la tabla de transformaciones "".
Alternativamente, puede definirse como en cuyo caso
los coeficientes son proporcionales a las muestras de S1/T(f) a intervalos disretos de 1P = 1NT:
Por el contrario, cuando se quiere calcular un número arbitrario (N) de muestras discretas de un ciclo de una DTFT continua, S1/T(f) , se puede hacer calculando la DFT relativamente simple de S1/N(f), como se ha 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 con ceros o interpolación, da como resultado muestras más próximas entre sí de un ciclo de S1/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 puede calcularse mediante un algoritmo de transformada rápida de Fourier (FFT), lo que la convierte en una transformación práctica e importante en los ordenadores.
Ver Transformada discreta de Fourier para mucha más información, incluyendo:
- propiedades de la transformada
- aplicaciones
- transformaciones tabuladas de funciones específicas
Resumen
Para funciones periódicas, tanto la transformada de Fourier como la DTFT comprenden sólo un conjunto discreto de componentes de frecuencia (series 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 puede discernirse a partir de un solo ciclo de la función periódica, ya que todos los demás ciclos son idénticos. Del mismo modo, las funciones de duración finita pueden representarse como una serie de Fourier, sin pérdida real de información, salvo que la periodicidad de la transformada inversa es un mero artefacto.
Es común en la práctica que la duración de s(*) esté limitada al período, P o N. Pero estas fórmulas no requieren esa condición.
Frecuencia continua | Frecuencia discreta | |
---|---|---|
Transformada | ||
Inversa |
Frecuencia continua | Frecuencia discreta | |
---|---|---|
Transformada |
| |
Inversa |
|
|
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, denotadas a continuación por los subíndices RE, RO, IE e IO. Y existe un mapeo uno a uno entre los cuatro componentes de una función temporal compleja y los cuatro componentes de su transformada de frecuencia compleja:[8]
De ello se desprenden varias relaciones, por ejemplo:
- La transformada de una función de valor real sRE + sRO es la función SRE + i S IO. A la inversa, una transformación par-simétrica implica un dominio de tiempo de valor real.
- La transformada de una función de valor imaginario i sIE + i sIO es la función SRO + i S IE, y lo contrario es cierto.
- La transformada de una función par-simétrica sRE + i sIO es la función de valor real SRE + SRO, y la inversa es cierta.
- La transformada de una función impar-simétrica SRO + i sIE es la función de valor imaginario i sIE + i S IO, y la inversa es cierta.
Historia
Una forma temprana de series armónicas se remonta a las antiguas matemáticas babilónicas, donde se utilizaban para calcular las efemérides (tablas de posiciones astronómicas).[9][10][11][12] Los conceptos griegos clásicos de deferente y epiciclo en el sistema ptolemaico de astronomía estaban relacionados con las series de Fourier.
En tiempos modernos, variantes de la transformada discreta de Fourier fueron utilizadas por Alexis Clairaut en 1754 para calcular una órbita,[13] que se ha descrito como la primera fórmula de la DFT,[14] y en 1759 por Joseph Louis Lagrange, al calcular los coeficientes de una serie trigonométrica para una cuerda vibrante.[15] Técnicamente, el trabajo de Clairaut era una serie sólo de coseno (una forma de Transformada de coseno discreta), mientras que el trabajo de Lagrange era una serie sólo de seno (una forma de transformada discreta del seno); una verdadera DFT de coseno+seno fue utilizada por Gauss en 1805 para la interpolación trigonométrica de las órbitas de asteroides.[15] Tanto Euler como Lagrange discretizaron el problema de la cuerda vibrante, utilizando lo que hoy se llamaría muestras.[14]
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 resolventes de Lagrange utilizó una descomposición compleja de Fourier para estudiar la solución de una cúbica:[16] Lagrange transformó las raíces x1, x2, x3 en los resolventes:
donde ζ es una raíz de la unidad cúbica, que es la DFT de orden 3.
Varios autores, especialmente Jean le Rond d'Alembert, y Carl Friedrich Gauss utilizaron series trigonométricas para estudiar la ecuación del calor,[17] pero el avance decisivo fue el artículo de 1807 Memoria sobre la propagación del calor en los cuerpos sólidos 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 al crédito que hay que dar 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 las funciones, y Lagrange había dado la solución en 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 podía representarse mediante una serie de Fourier.[14]
El desarrollo posterior del campo se conoce como análisis armónico, y es también una instancia temprana de la teoría de la representación.
El primer algoritmo de la transformada rápida de Fourier (FFT) para la DFT fue descubierto alrededor de 1805 por Carl Friedrich Gauss al interpolar las mediciones 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.[15][13]
Transformaciones 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 perfecta resolución de tiempo, pero sin información de frecuencia, mientras que la transformada de Fourier tiene perfecta resolución de frecuencia, pero sin información de tiempo.
Como alternativas a la transformada de Fourier, en el análisis tiempo-frecuencia, se utilizan las transformadas de tiempo-frecuencia para representar las señales en una forma que tiene algo de información de tiempo y algo de información de frecuencia - por el principio de incertidumbre, hay un compromiso entre estos. Pueden ser generalizaciones de la transformada de Fourier, como la Transformada de Fourier de Tiempo Reducido, la transformada de Gabor o la transformada de Fourier fraccional (FRFT), o pueden utilizar diferentes funciones para representar las señales, como en las transformada ondícula y las transformadas chirplet, siendo el análogo wavelet de la transformada de Fourier (continua) la ondícula ontínua.
Véase también
Referencias
- Elias M. Stein, Rami Shakarchi. Fourier Analysis: An Introduction (Princeton Lectures in Analysis). (2003) 328 pág. ISBN 069111384X ISBN 978-0691113845
- Rudin, Walter (1990). Fourier Analysis on Groups. Wiley-Interscience. ISBN 978-0-471-52364-2.
- Evans, L. (1998). Partial Differential Equations. American Mathematical Society. ISBN 978-3-540-76124-2.
- Knuth, Donald E. (1997). The Art of Computer Programming Volume 2: Seminumerical Algorithms (3rd edición). Addison-Wesley Professional. Section 4.3.3.C: Discrete Fourier transforms, pg.305. ISBN 978-0-201-89684-8.
- Conte, S. D.; de Boor, Carl (1980). Elementary Numerical Analysis (Third edición). New York: McGraw Hill, Inc. ISBN 978-0-07-066228-5.
- Saferstein, Richard (2013). Criminalistics: An Introduction to Forensic Science.
- Rabiner, Lawrence R.; Gold, Bernard (1975). Theory and Application of Digital Signal Processing. Englewood Cliffs, NJ. (requiere registro).
- Proakis, John G.; Manolakis, Dimitri G. (1996). Digital Signal Processing: Principles, Algorithms and Applications (en inglés) (3 edición). New Jersey: Prentice-Hall International. p. 291. ISBN 9780133942897.
- Prestini, Elena (2004). The Evolution of Applied Harmonic Analysis: Models of the Real World. Birkhäuser. p. 62. ISBN 978-0-8176-4125-2.
- Rota, Gian-Carlo; Palombi, Fabrizio (1997). Indiscrete Thoughts. Birkhäuser. p. 11. ISBN 978-0-8176-3866-5.
- Neugebauer, Otto (1969). «The Exact Sciences in Antiquity». Acta Historica Scientiarum Naturalium et Medicinalium 9 (2nd edición) (Dover Publications). pp. 1-191. ISBN 978-0-486-22332-2. PMID 14884919.
- Brack-Bernsen, Lis; Brack, Matthias (2004). «Analyzing shell structure from Babylonian and modern times». International Journal of Modern Physics E 13 (1): 247. Bibcode:2004IJMPE..13..247B. arXiv:physics/0310126.
- Terras, Audrey (1999). Fourier Analysis on Finite Groups and Applications. Cambridge University Press. pp. 30-32. ISBN 978-0-521-45718-7.
- Briggs, William L.; Henson, Van Emden (1995). The DFT: An Owner's Manual for the Discrete Fourier Transform. SIAM. pp. 2-4. ISBN 978-0-89871-342-8.
- Heideman, M.T.; Johnson, D. H.; Burrus, C. S. (1984). «Gauss and the history of the fast Fourier transform». IEEE ASSP Magazine 1 (4): 14-21. S2CID 10032502. doi:10.1109/MASSP.1984.1162257.
- Knapp, Anthony W. (2006). Basic Algebra. Springer. p. 501. ISBN 978-0-8176-3248-9.
- Narasimhan, T.N. (February 1999). «Fourier's heat conduction equation: History, influence, and connections». Reviews of Geophysics 37 (1): 151-172. Bibcode:1999RvGeo..37..151N. ISSN 1944-9208. OCLC 5156426043.
Bibliografía
- Howell, Kenneth B. (2001). Principles of Fourier Analysis. CRC Press. ISBN 978-0-8493-8275-8.
- Kamen, E.W.; Heck, B.S. (2 de marzo de 2000). Fundamentals of Signals and Systems Using the Web and Matlab (2 edición). Prentiss-Hall. ISBN 978-0-13-017293-8. (requiere registro).
- Müller, Meinard (2015). The Fourier Transform in a Nutshell. Springer. In Fundamentals of Music Processing, Section 2.1, pp. 40–56. ISBN 978-3-319-21944-8. S2CID 8691186. doi:10.1007/978-3-319-21945-5.
- Polyanin, A. D.; Manzhirov, A. V. (1998). Handbook of Integral Equations. Boca Raton: CRC Press. ISBN 978-0-8493-2876-3.
- Smith, Steven W. (1999). The Scientist and Engineer's Guide to Digital Signal Processing (Second edición). San Diego: California Technical Publishing. ISBN 978-0-9660176-3-2.
- Stein, E. M.; Weiss, G. (1971). Introduction to Fourier Analysis on Euclidean Spaces. Princeton University Press. ISBN 978-0-691-08078-9. (requiere registro).
Enlaces externos
- Tables of Integral Transforms at EqWorld: The World of Mathematical Equations.
- An Intuitive Explanation of Fourier Theory by Steven Lehar.
- Lectures on Image Processing: A collection of 18 lectures in pdf format from Vanderbilt University. Lecture 6 is on the 1- and 2-D Fourier Transform. Lectures 7–15 make use of it., by Alan Peters
- Moriarty, Philip; Bowley, Roger (2009). «Σ Summation (and Fourier Analysis)». Sixty Symbols. Brady Haran for the University of Nottingham.