Transformada de Fourier

Nivel AvanzadoDificultad ★★★★★Concepto⌖ Ver en el mapa

¿Qué es?

La versión continua para señales no periódicas: f^(ξ)=∫f(t) e−2πiξt dt\hat f(\xi) = \int f(t)\,e^{-2\pi i\xi t}\,\dd t da la cantidad de cada frecuencia ξ\xi. Convierte la convolución en multiplicación y la derivación en multiplicar por 2πiξ2\pi i\xi; por eso filtrar, comprimir y resolver EDP lineales es más fácil en el espacio de frecuencias.

¿Por qué existe?

Los sistemas lineales e invariantes en el tiempo (filtros, lentes, salas, líneas de transmisión, la ecuación del calor) no mezclan frecuencias: entra una sinusoide y sale la misma sinusoide, escalada y desplazada. En el dominio de la frecuencia un sistema así es solo una multiplicación, así que la transformada de Fourier lo diagonaliza.

Intuición

Enrolla la señal alrededor de una circunferencia a frecuencia ξ\xi (multiplica por e−2πiξte^{-2\pi i\xi t}) y calcula el centro de masas. Si la señal tiene una componente en ξ\xi, el enrollado la acumula en un lado y el centro se aleja de cero; si no, todo se compensa. f^(ξ)\hat f(\xi) es ese centro de masas. Una señal estrecha necesita muchas frecuencias (espectro ancho); un tono puro tiene un único pico: el principio de incertidumbre tiempo–frecuencia.

Definición formal

f^(ξ)=∫−∞∞f(t) e−2πiξt dt,f(t)=∫−∞∞f^(ξ) e2πiξt dξ.\hat f(\xi) = \int_{-\infty}^{\infty} f(t)\,e^{-2\pi i\xi t}\,\dd t, \qquad f(t) = \int_{-\infty}^{\infty}\hat f(\xi)\,e^{2\pi i\xi t}\,\dd\xi.

Propiedades clave: f∗g^=f^ g^\widehat{f * g} = \hat f\,\hat g (teorema de convolución), f′^(ξ)=2πiξ f^(ξ)\widehat{f'}(\xi) = 2\pi i\xi\,\hat f(\xi), ∫∣f∣2=∫∣f^∣2\int|f|^2 = \int|\hat f|^2 (Plancherel). La transformada discreta de NN muestras, Xk=∑n=0N−1xne−2πikn/NX_k = \sum_{n=0}^{N-1} x_n e^{-2\pi i kn/N}, se calcula en O(Nlog⁡N)O(N\log N) con la FFT.

Fórmulas

f^(ξ)=∫−∞∞f(t) e−2πiξt dt\hat f(\xi) = \int_{-\infty}^{\infty} f(t)\,e^{-2\pi i\xi t}\,\dd t
f∗g^=f^⋅g^\widehat{f * g} = \hat f\cdot\hat g
teorema de convolución
Xk=∑n=0N−1xn e−2πikn/NX_k = \sum_{n=0}^{N-1} x_n\,e^{-2\pi i k n/N}
transformada discreta de Fourier (DFT)
e−πt2^=e−πξ2\widehat{e^{-\pi t^2}} = e^{-\pi\xi^2}
la gaussiana es su propia transformada

Ejemplo

Una grabación de 3 segundos de un la 440 con algo de ruido: su DFT muestra un pico nítido en 440 Hz (y en 880, 1320… si el instrumento tiene armónicos). Un afinador de guitarra es literalmente este cálculo seguido de un argmax.

¿Por qué importa?

Los códecs de audio, la compresión de imagen, el wifi y el 4G/5G (OFDM), la reconstrucción de resonancias magnéticas, el radar, la espectroscopia, la multiplicación rápida de polinomios y enteros, las convoluciones grandes: todo funciona con transformadas de Fourier calculadas por la FFT, uno de los algoritmos más importantes jamás escritos.

Aplicaciones en informática

  • Procesamiento de señales★★★★★fundamentalSeñales, multimedia y visión

    El análisis espectral, el filtrado y la modulación se definen en el dominio de la frecuencia.

  • Transformada rápida de Fourier (FFT)★★★★★fundamentalSeñales, multimedia y visión

    La FFT calcula la DFT en O(Nlog⁡N)O(N\log N) en vez de O(N2)O(N^2).

  • Teorema de muestreo (Nyquist–Shannon)★★★★★fundamentalSeñales, multimedia y visión

    Muestrear replica el espectro; la condición de Nyquist evita que las copias se solapen.

  • Telecomunicaciones (modulación, OFDM)★★★★★fundamentalSeñales, multimedia y visión

    OFDM (wifi, 4G/5G, ADSL) envía datos en muchas subportadoras ortogonales usando FFT inversa y FFT.

  • Principio de incertidumbre★★★★★fundamentalComputación y física cuántica

    Las funciones de onda de posición y momento son pares de Fourier; el principio de incertidumbre es una desigualdad de Fourier.

  • Procesamiento de imagen y visión artificial★★★★★frecuenteSeñales, multimedia y visión

    Filtrado en el dominio de la frecuencia, eliminación de desenfoque y registro de imágenes (correlación de fase).

Dónde aparece en IA

  • Redes convolucionales (CNN)★★★★★avanzadaIA y machine learning

    Las convoluciones grandes se pueden calcular con FFT; los Fourier neural operators aprenden directamente en el espacio de frecuencias.

¿Dónde se utiliza?

Temas de informática a los que se llega desde aquí, con la cadena de ideas que lleva a ellos:

Qué depende de él

Ejercicios

1Informática

Convolucionar directamente dos señales de longitud N=106N = 10^6 cuesta ≈N2\approx N^2 operaciones. Estima el coste con FFT.

Solución

Tres FFT de tamaño ~2N2N más un producto punto a punto: unas 3⋅2Nlog⁡2(2N)≈1,3⋅1083\cdot 2N\log_2(2N) \approx 1{,}3\cdot10^8 frente a 101210^{12}, unas 8000 veces más rápido.

2Interpretación gráfica

Una señal es un chasquido corto (un pulso estrecho). ¿Cómo es su espectro? ¿Y el de un tono puro largo?

Solución

Un pulso estrecho tiene un espectro muy ancho y plano (todas las frecuencias); un tono puro largo, un pico estrecho. El producto anchura temporal × anchura en frecuencia está acotado inferiormente.

↑ ↓ para navegar · ↵ · Esc