Transformada de Fourier

23 de mayo de 2019
La Teoría
∞
2 2𝑗𝑘𝜋
𝑥 𝑡 = 0.8 + ෍ sin 0.4𝑘𝜋 𝑒 𝑇 𝑡
𝑘𝜋
• Con la serie de Fourier, obtenemos coeficientes 𝑘=1
espectrales relacionados a exponenciales
complejas de frecuencias (fundamental y
armónicas) 𝑘𝜔0
• Podríamos representar los coeficientes con
respecto a una recta 𝜔 tal que en el punto en la
recta equivalente a 𝜔0 se grafique un punto
equivalente en magnitud a 𝑎1 , en el punto 2𝜔0
se grafique 𝑎2 , y así sucesivamente.
• La curva resultante sería equivalente a una
suma infinita de impulsos ponderados
respectivamente por cada coeficiente espectral.
 • Esta sumatoria de impulsos va a ser un ejemplo de una función definida en el
dominio de la frecuencia.
• (Notación: el nombre de la función normalmente irá en mayúsculas, y su dominio va
a ser j𝜔, no solo 𝜔. Ejemplos: 𝑋 𝑗𝜔 , 𝑌 𝑗𝜔 , 𝐻(𝑗𝜔))
• Continuando la analogía del tiempo de las diapositivas anteriores, estas son
funciones definidas no con respecto al tiempo, sino que con respecto a ciclos o
armonías. Su posición abscisa no va a ser el instante relativo en que ocurre, sino la
frecuencia asociada a ese armónico específico.
Dominio de • Dicho esto, el análisis espectral de una señal normalmente irá enfrentado al tiempo,
por medio de una cascada de magnitudes. En un espectrograma, por ejemplo, el eje
frecuencia de abscisas representa el tiempo t, el eje de ordenadas las frecuencias armónicas 𝜔, y
por código de colores se representa la intensidad.

Espectrograma de la parte opereta de Bohemian Rhapsody, de Queen (1975), capturado con Audacity desde .mp3
(I see a little silhouette-o of a man…)
 Espectrograma de la parte opereta de Bohemian Rhapsody, de Queen (1975), capturado con Audacity desde .mp3
(I see a little silhouette-o of a man…)

• Considerando por ejemplo este espectrograma, donde la intensidad se representa en una escala
de azul a blanco pasando por el rojo.
Dominio de • Empieza en la parte opereta (minuto 3:03 aprox.) de Bohemian Rhapsody de Queen (1975), en el

frecuencia piano staccato: el acorde de piano constante como metrónomo. En el espectrograma, esto se ve
como pilares de rojo, un sonido de media intensidad, contrastando con pilares azules:
esencialmente, silencios.
• Cuando Freddy empieza a cantar, el timbre brillante de su voz y las armonías en falsetto colorean
las frecuencias más altas de rojo por la duración de cada frase.
• Asimismo, cuando el coro se apodera del verso “thunderbolts and lightning/ very very
frightening” con su vigor característico, se ve un aumento a blanco en la intensidad del sonido
en frecuencias que, se puede asumir, son vecindades de las fundamentales y primeras
armónicas de las notas que componen los acordes cantados.
• Y, así, se puede continuar analizando este espectrograma comparándolo con la canción original.
La transformada de Fourier

• A medida que aumentamos el periodo fundamental de una función periódica, su

frecuencia fundamental disminuye. En dominio de frecuencia, los impulsos
asociados a su serie de Fourier se van acercando cada vez más.
• Eventualmente, para un periodo infinitesimal, la distancia entre cada impulso es
infinitesimalmente pequeña y despreciable.
• Pero, ¡una función periódica que se repite luego de una infinidad es una función
aperiódica! Y, una suma de elementos infinitesimalmente cercanos es una integral.
• Entonces, la ecuación de análisis de la serie, ahora
∞
transformada, de Fourier, será:
𝑋 𝑗𝜔 = lim 𝑇 ⋅ 𝑎𝑘 = න 𝑥 𝑡 𝑒 −𝑗𝜔𝑡 𝑑𝑡
𝑇→∞
Δ𝑘→0 −∞
• Esta es una transformación que permite representar una función definida en el
tiempo como una función definida en frecuencia.
La transformada inversa de fourier
• Reemplazando en la ecuación de síntesis esta expresión para la
transformada de Fourier en lugar de los coeficientes espectrales, el periodo
fundamental por una frecuencia fundamental según 𝑇 = 2𝜋/𝜔0 , y
sabiendo que esta frecuencia fundamental se volverá diferencialmente
pequeña en la integral, obtenemos una transformación que permite
representar una función definida en frecuencia como una función definida
en el tiempo.
• En otras palabras, la ecuación de síntesis se convierte en la transformada
inversa de Fourier. Con estos reemplazos algebraicos, esta transformada se
ve así:
1 ∞
𝑥 𝑡 = න 𝑋 𝑗𝜔 𝑒 𝑗𝜔𝑡 𝑑𝜔
2𝜋 −∞
Resumiendo las definiciones…
Serie de Fourier Transformada de Fourier

Ecuación de análisis: Transformada de Fourier:

∞
Para expresar una función 1
𝑎𝑘 = න 𝑥 𝑡 𝑒 −𝑗𝑘𝜔0𝑡 𝑑𝑡 𝑋 𝑗𝜔 = න 𝑥 𝑡 𝑒 −𝑗𝜔𝑡 𝑑𝑡
en el tiempo según 𝑇 𝑇 −∞
frecuencia (Otorga los coeficientes de Fourier, o (Otorga la transformada de
coeficientes espectrales, de x) Fourier, o el espectro, de x)

Para expresar una función Ecuación de síntesis: Transformada inversa:

∞
en frecuencia según el 1 ∞
𝑥 𝑡 = ෍ 𝑎𝑘 𝑒 𝑗𝑘𝜔0 𝑡 𝑥 𝑡 = න 𝑋 𝑗𝜔 𝑒 𝑗𝜔𝑡 𝑑𝜔
tiempo 2𝜋 −∞
𝑘=−∞

Para un 𝑋(𝑗𝜔) equivalente a la Para coeficientes espectrales {𝑎𝑘 }

Para expresar una serie de ∞
transformada de un periodo de x(t),
Fourier como función en 1 𝑋 𝑗𝜔 = ෍ 2𝜋𝑎𝑘 𝛿(𝜔 − 𝑘𝜔0 )
frecuencia y viceversa. 𝑎𝑘 = 𝑋 𝑗𝜔 ቚ
𝑇 𝜔=𝑘𝜔0 𝑘=−∞
Problemas
• Encontrar la transformada de Fourier de:
3𝜋 3𝜋
𝑥 𝑡 = cos 𝑡 − sen 𝑡
8 4

𝑥 𝑡 = 𝑒 𝑎𝑡 𝑢(𝑡)
 3𝜋 3𝜋
𝑥 𝑡 = cos 𝑡 − sen 𝑡
8 4
Se ve que esta es una suma trigonométrica periódica, por lo que, de
resolver algebraicamente, resulta conveniente llegar a una expresión en
serie de Fourier antes de pasar a la transformada.
1 −3𝜋𝑗𝑡 3𝜋
𝑗𝑡
3𝜋 3𝜋
𝑥 𝑡 = (𝑒 8 + 𝑒 8 − 𝑗𝑒 4 + 𝑗𝑒 4 𝑗𝑡 )
− 𝑗𝑡
2
Convirtiendo esta serie por tren de impulsos, donde todo lo que viene
antes del t en la exponencial es su desfase, se obtiene la siguiente
función en frecuencia:
𝑋 𝑗𝜔
3𝜋 3𝜋 3𝜋 3𝜋
= 𝜋𝛿 𝜔 − + 𝜋𝑗𝛿 𝜔 − − 𝜋𝑗𝛿 𝜔 + + 𝜋𝑗𝛿 𝜔 +
8 4 4 8
 𝑥 𝑡 = 𝑒 −𝑎𝑡 𝑢 𝑡 , 𝑎 > 0, 𝑎∈ℝ
Esta señal es aperiódica, por lo que su espectro se obtendría con la integral de
Fourier. Se desarrolla ésta:
∞
𝑋 𝑗𝜔 = න 𝑥 𝑡 𝑒 −𝑗𝜔t 𝑑𝑡
∞ −∞ ∞
𝑋 𝑗𝜔 = න 𝑒 −𝑎𝑡 𝑒 −𝑗𝜔𝑡 𝑑𝑡 = න 𝑒 −(𝑎+𝑗𝜔)𝑡 𝑑𝑡
0 0
1 ∞ 1
𝑋 𝑗𝜔 = − 𝑒 − 𝑎+𝑗𝜔 𝑡
ቚ = −0 − −
𝑎 + 𝑗𝜔 0 𝑎 + 𝑗𝜔
1
𝑋 𝑗𝜔 =
𝑎 + 𝑗𝜔
¡Somos ingenieros en informática y
telecomunicaciones, no matemáticas!
• ¿Podemos calcular las integrales de Fourier para cada caso? Sí.
¿Deberíamos? No. Otros ya lo han hecho por nosotros.
Señal: 𝒙 𝒕 = ⋯ Transformada de Fourier Coeficientes de Fourier
∞ ∞ 𝑎𝑘
෍ 𝑎𝑘 𝑒 𝑗𝑘𝜔0 𝑡 2𝜋 ෍ 𝑎𝑘 𝛿(𝜔 − 𝑘𝜔0 )
𝑘=−∞ 𝑥=−∞
∞
𝑇 2 sin 𝑘𝜔0 𝑇1 𝜔0 𝑇1 𝑘𝜔0 𝑇1 𝑘𝜔0 𝑇1
1, 𝑡 < 𝑇1 < ෍ 𝛿(𝜔 − 𝑘𝜔0 ) senc = sen
2 𝑘 𝜋 𝜋 𝑘
𝑘=−∞
𝑇
0, 𝑇1 < 𝑡 <
2
𝑥 𝑡 + 𝑇 = 𝑥(𝑡)
∞ ∞ 1
2𝜋 2𝑘𝜋 𝑎𝑘 =
෍ 𝛿(𝑡 − 𝑛𝑇) ෍ 𝛿 𝜔− 𝑇
𝑇 𝑇
𝑘=−∞ 𝑘=−∞
𝛿(𝑡 − 𝑡0 ) 𝑒 −𝑗𝜔𝑡0
𝑢(𝑡) 1
+ 𝜋𝛿 𝜔
𝑗𝜔
𝑒 −𝑎𝑡 𝑢 𝑡 1
ℛℯ 𝑎 > 0 𝑎 + 𝑗𝜔
𝑡 𝑛−1 1
𝑒 −𝑎𝑡 𝑢 𝑡 𝑛
𝑛−1 ! 𝑎 + 𝑗𝜔
ℛℯ 𝑎 > 0, 𝑛∈ℕ
Problemas
• Determinar transformada de Fourier de:
2, 2≤𝑡<3
𝑥 𝑡 = ቐ3, 3≤𝑡<5
0, 𝑒. 𝑜. 𝑐
𝑥 𝑡 = 𝑡 2 𝑒 −4𝑡 𝑢(𝑡)
• Determinar transformada inversa de Fourier de:
21 + 16𝑗𝜔 − 4𝜔2
𝑋 𝑗𝜔 =
27 + 45𝑗𝜔 − 24𝜔 2 − 4𝑗𝜔 3
 2, 2≤𝑡<3
𝑥 𝑡 = ቐ3, 3≤𝑡<5
0, 𝑒. 𝑜. 𝑐
Esta expresión es equivalente a
𝑥 𝑡 = 2𝑢 𝑡 − 2 + 𝑢 𝑡 − 3 − 3𝑢(𝑡 − 5)
Calculemos esta transformada con la integral correspondiente
∞ 3 5
𝑋 𝑗𝜔 = න 𝑥 𝑡 𝑒 −𝑗𝜔t 𝑑𝑡 = 2 න 𝑒 −𝑗𝜔t 𝑑𝑡 + 3 න 𝑒 −𝑗𝜔t 𝑑𝑡
−∞ 2 3
1 3 5
𝑋 𝑗𝜔 = − 2𝑒 −𝑗𝜔t ቚ
+ 3𝑒 −𝑗𝜔t ቚ
𝑗𝜔 2 3
1
𝑋 𝑗𝜔 = − 2𝑒 −3𝑗𝜔 − 2𝑒 −2𝑗𝜔 + 3𝑒 −5𝑗𝜔 − 3𝑒 −3𝑗𝜔
𝑗𝜔
1
𝑋 𝑗𝜔 = 2𝑒 −2𝑗𝜔 + 𝑒 −3𝑗𝜔 − 3𝑒 −5𝑗𝜔
𝑗𝜔
 𝑥 𝑡 = 2𝑢 𝑡 − 2 + 𝑢 𝑡 − 3 − 3𝑢(𝑡 − 5)
1
𝑋 𝑗𝜔 = 2𝑒 −2𝑗𝜔 + 𝑒 −3𝑗𝜔 − 3𝑒 −5𝑗𝜔
𝑗𝜔
Con este resultado se pueden comprobar dos cosas:
1. La transformada de Fourier de la señal se parece en linealidad a la señal expresada en escalones,
donde
2. Cada desfase en tiempo pareciera traducirse en un factor exponencial complejo.
Asumiendo que esto es cierto, se puede recalcular esta transformada con datos de tabla:
𝑋 𝑗𝜔 = 2ℱ 𝑢 𝑡 − 2 + ℱ 𝑢 𝑡 − 3 − 3ℱ 𝑢 𝑡 − 5
1 1 1
𝑋 𝑗𝜔 = 2 + 𝜋𝛿 𝜔 𝑒 −2𝑗𝜔 + + 𝜋𝛿 𝜔 𝑒 −3𝑗𝜔 − 3 + 𝜋𝛿 𝜔 𝑒 −5𝑗𝜔
𝑗𝜔 𝑗𝜔 𝑗𝜔
1
𝑋 𝑗𝜔 = 2𝑒 −2𝑗𝜔 + 𝑒 −3𝑗𝜔 − 3𝑒 −5𝑗𝜔 + 𝜋𝛿 𝜔 2𝑒 −2𝑗𝜔 + 𝑒 −3𝑗𝜔 − 3𝑒 −5𝑗𝜔
𝑗𝜔
Se nota que como el impulso va a ser distinto a 0 solo en w=0, se pueden evaluar las exponenciales
complejas:
1
𝑋 𝑗𝜔 = 2𝑒 −2𝑗𝜔 + 𝑒 −3𝑗𝜔 − 3𝑒 −5𝑗𝜔 + 𝜋𝛿 𝜔 2𝑒 −2𝑗⋅0 + 𝑒 −3𝑗⋅0 − 3𝑒 −5𝑗⋅0
𝑗𝜔
1
𝑋 𝑗𝜔 = 2𝑒 −2𝑗𝜔 + 𝑒 −3𝑗𝜔 − 3𝑒 −5𝑗𝜔 + 𝜋𝛿 𝜔 2 + 1 − 3
𝑗𝜔
1
𝑋 𝑗𝜔 = 2𝑒 −2𝑗𝜔 + 𝑒 −3𝑗𝜔 − 3𝑒 −5𝑗𝜔
𝑗𝜔
Finalmente, esta expresión es idéntica a la obtenida por la integral, comprobándose estas propiedades
de linealidad y desfase. (Éstas y otras propiedades se van a poder ver mejor en una tabla en la siguiente
presentación.)
 𝑥 𝑡 = 𝑡 2 𝑒 −4𝑡 𝑢(𝑡)
Viendo en la tabla de pares básicos, se comprueba que esta señal
corresponde a:
𝑡 (𝑛−1) −𝑎𝑡
𝑥 𝑡 =2⋅ 𝑒 𝑢 𝑡 , 𝑛 − 1 = 2, 𝑎 = 4
(𝑛 − 1)!
Entonces, la transformada de Fourier para esta señal será:
2
𝑋 𝑗𝜔 =
4 + 𝑗𝜔 3
(Acá, gracias a la linealidad comprobada en el ejercicio anterior,
podemos tener la certeza que el factor de 2 de la señal original se
traspasa a su transformada también como factor escalar.)
 21 + 16𝑗𝜔 − 4𝜔2
𝑋 𝑗𝜔 =
27 + 45𝑗𝜔 − 24𝜔 2 − 4𝑗𝜔 3
Al igual que uno sabe calcular transformadas, uno debe saber calcular transformadas inversas.
Como tenemos una transformada que es una fracción de polinomios, nos puede convenir
factorizar el denominador y separar esta expresión en sumas parciales, y ver si así se logra
parecer a alguna expresión de la segunda columna de la tabla de pares básicos.
Se prueba a dividir el denominador por un binomio de raíz divisora de 27. Divisores
potenciales serían 1, 3, 9, o 27. Se prueba con el 3:
−4𝑗𝜔3 − 24𝜔2 + 45𝑗𝜔 + 27: 𝑗𝜔 + 3 = −4𝜔2 + 12𝑗𝜔 + 9
−4𝑗𝜔3 − 12𝜔2
−12𝜔2
−12𝜔2 + 36𝑗𝜔
9𝑗𝜔
9𝑗𝜔 + 27
0
Se ve que el cociente polinómico obtenido es un cuadrado perfecto de 2𝑗𝜔 + 3 2 = −4𝜔2 +
12𝑗𝜔 + 9
 21 + 16𝑗𝜔 − 4𝜔2 𝐴 𝐵𝑗𝜔 + 𝐶
𝑋 𝑗𝜔 = 2 3
= +
27 + 45𝑗𝜔 − 24𝜔 − 4𝑗𝜔 𝑗𝜔 + 3 9 + 12𝑗𝜔 − 4𝜔 2
De la descomposición en fracciones parciales se formula el siguiente sistema de
ecuaciones:
21 = 9𝐴 + 3𝐶 ⇔ 7 = 3𝐴 + 𝐶
ቐ 16 = 12𝐴 + 3𝐵 + 𝐶
−4 = −4𝐴 − 𝐵 ⇔ 4 = 4𝐴 + 𝐵
Sumando el inverso de la primera ecuación simplificada con la segunda y el triple de
la tercera nos permite encontrar A
16 − 7 − 12 = 12𝐴 + 3𝐵 + 𝐶 − 3𝐴 − 𝐶 − 12𝐴 − 3𝐵
−3 = −3𝐴 ⇔ 𝐴 = 1
Insertando este valor en la primera y tercera ecuación nos permite determinar B y C
21 − 9
21 = 9 + 3𝐶 ⇔ 𝐶 = =4
3
4=4+𝐵 ⇔ 𝐵 =0
 21 + 16𝑗𝜔 − 4𝜔2 1 4
𝑋 𝑗𝜔 = 2 3
= +
27 + 45𝑗𝜔 − 24𝜔 − 4𝑗𝜔 𝑗𝜔 + 3 9 + 12𝑗𝜔 − 4𝜔 2
Con esta descomposición precisa en fracciones parciales, se trata de
ajustar los parámetros para que encajen con algún par básico conocido.
Para ello, se multiplica un 1 al denominador de la segunda fracción:
1 4
𝑋 𝑗𝜔 = +
3 + 𝑗𝜔 4 ⋅ 1 ⋅ 3 + 2𝑗𝜔 2
4
1 1
𝑋 𝑗𝜔 = +
3 + 𝑗𝜔 1.5 + 𝑗𝜔 2
Esta expresión corresponde exactamente a una exponencial con escalón,
y una exponencial con rampa. Se reemplazan estas expresiones para
explicitar esta transformada inversa y terminar el ejercicio.
𝑥 𝑡 = 𝑒 −3𝑡 𝑢 𝑡 + 𝑡𝑒 −1.5𝑡 𝑢(𝑡)