Guía 2 - Ejercicios 15-19: Transformada de Fourier
Práctica 2: Análisis de Fourier - Resoluciones (Ej. 15 a 19)
La Transformada de Fourier
Marco Teórico Previo
Pasamos del dominio de series (espectro discreto para señales periódicas) al dominio de transformadas (espectro continuo para señales aperiódicas).
Transformada de Fourier (Convención Estándar en Matemática/Modelado):
F{f(x)}(ξ)=f^(ξ)=∫−∞∞f(x)e−2πixξdx
Fórmula de Inversión:
f(x)=∫−∞∞f^(ξ)e2πixξdξ
Propiedades Esenciales de Derivación:
Transformada de la derivada: F{f′(x)}=(2πiξ)f^(ξ)
Derivada de la transformada: dξdf^(ξ)=F{−2πixf(x)}⟹F{xf(x)}=2πidξdf^(ξ)
Ejercicio 15: Integral Gaussiana de e−πx2 en (0,∞)
Intuición Didáctica
La campana de Gauss e−πx2 es la función reina del análisis armónico. Su integral no tiene primitiva expresable mediante funciones elementales. El "truco" de elevarla al cuadrado para pasar al plano R2 y usar Fubini para convertirla en coordenadas polares es uno de los recursos analíticos más elegantes de la matemática. Vamos a demostrarlo por dos vías distintas.
Resolución Rigurosa
Llamemos I=∫0∞e−πx2dx. Como el integrando es siempre positivo, sabemos que I>0.
Método 1: Coordenadas Polares (El Clásico)
Multiplicamos la integral por sí misma, cambiando la variable muda (variable de integración) para no mezclar:
I2=(∫0∞e−πx2dx)(∫0∞e−πy2dy)
Por el Teorema de Fubini (podemos aplicarlo ya que el integrando e−π(x2+y2) es medible y absolutamente integrable en el primer cuadrante R+×R+), reescribimos el producto como una integral doble:
I2=∫0∞∫0∞e−π(x2+y2)dxdy
Realizamos un cambio a coordenadas polares. El dominio es el primer cuadrante, por lo que:
x=rcosθ,y=rsinθ
x2+y2=r2
El diferencial de área (Jacobiano) es dxdy=rdrdθ
Los límites: r∈[0,∞) y θ∈[0,π/2].
Sustituyendo:
I2=∫0π/2∫0∞e−πr2rdrdθ
Como los límites son constantes y el integrando es separable, resolvemos primero la integral en θ:
I2=(∫0π/2dθ)(∫0∞re−πr2dr)=2π∫0∞re−πr2dr
Ahora aplicamos sustitución simple en r: sea u=πr2⟹du=2πrdr⟹rdr=2πdu.
Límites: si r=0→u=0; si r→∞→u→∞.
I2=2π∫0∞e−u2πdu=41[−e−u]0∞
Evaluando los límites, en el infinito e−u→0 y en cero vale 1:
I2=41(0−(−1))=41
Como I>0, tomamos la raíz positiva: I=1/2.
(Nota: La integral en (−∞,∞) vale el doble, o sea 1, demostrando que e−πx2 tiene área unitaria).
Método 2: Teorema de Fubini sobre una función auxiliar f(x,y)=ye−π(x2+1)y2
Consideremos la integral doble en el cuadrante positivo:
J=∫0∞∫0∞ye−π(x2+1)y2dxdy
Integraremos de dos formas distintas intercambiando el orden (válido por Fubini/Tonelli).
Forma A (Integramos primero respecto a y):
J=∫0∞(∫0∞ye−π(x2+1)y2dy)dx
Sea la sustitución u=π(x2+1)y2⟹du=2π(x2+1)ydy⟹ydy=2π(x2+1)du.
Hacemos la sustitución t=yx⟹dt=ydx⟹dx=dt/y. (Notar que para x=0,t=0; para x→∞,t→∞).
∫0∞e−π(yx)2dx=∫0∞e−πt2ydt=y1∫0∞e−πt2dt=y1I
Sustituyendo esto en la integral exterior:
J=∫0∞ye−πy2(y1I)dy=I∫0∞e−πy2dy=I⋅I=I2
Igualando los resultados de ambas formas:
I2=1/4⟹I=1/2.
[!infobox] Ejercicio 15: Integral Gaussiana de 0 a ∞
Contexto: Cálculo de área bajo la campana de Gauss mediante Teorema de Fubini.
Demostración Rigurosa:
Por la simetría par, la integral en todo R es el doble que en [0, ∞).
Método 1 (Polares): I² = ∫∫ e^(-π(x²+y²)) dx dy. Se pasa a polares: Jacobiano es r dr dθ. El dominio en el primer cuadrante da θ ∈ [0, π/2]. La integral r*e^(-πr²) sale por sustitución u=πr². I² = 1/4 => I = 1/2.
Método 2 (Función Auxiliar): Usando J = ∫∫ y e^(-π(x²+1)y²) dx dy.
Si integramos dy primero, sale por sustitución u=π(x²+1)y², dejando un arctan(x) que evalúa a 1/4.
Si integramos dx primero, por sustitución t=yx aparece un 1/y que cancela la y externa, dejando el producto de dos integrales gaussianas separadas => I².
Igualando: I² = 1/4 => I=1/2.
Advertencia/Clave: Es vital no olvidar el Jacobiano 'r' al pasar a polares, es lo que "salva" a la integral permitiendo la sustitución que baja el grado del exponente.
Ejercicio 16: Transformada de Fourier de f(x)=e−πx2
Intuición Didáctica
La Gaussiana estándar e−πx2 es su propia Transformada de Fourier. Es el punto fijo (autofunción) del operador de Fourier. Demostrar esto es un hito monumental en la carrera. Haremos la demostración primero con expansiones formales y luego resolviendo un Problema de Valor Inicial (PVI) de variable compleja.
Resolución Rigurosa
Queremos calcular f^(ξ)=∫−∞∞e−πx2e−2πixξdx.
Método 1: Serie de Taylor e Integración Término a Término
Usamos la expansión en serie de Maclaurin de la exponencial: e−2πixξ=∑n=0∞n!(−2πixξ)n.
Sustituyendo en la integral:
f^(ξ)=∫−∞∞e−πx2n=0∑∞n!(−2πixξ)ndx
Intercambiamos suma con integral (justificado por el Teorema de Convergencia Dominada o convergencia uniforme en compactos y decaimiento gaussiano abrumador):
f^(ξ)=n=0∑∞n!(−2πiξ)n∫−∞∞xne−πx2dx
Analizamos la integral: Mn=∫−∞∞xne−πx2dx.
Si n es impar, el integrando xne−πx2 es una función impar en un dominio simétrico, por lo que Mimpar=0. ¡Desaparecen las potencias impares y la unidad imaginaria 'i'!
Si n es par (n=2k), el integrando es par: M2k=2∫0∞x2ke−πx2dx.
Calculamos M2k por recurrencia mediante integración por partes.
Para M2: ∫x2e−πx2dx. Elegimos u=x⟹du=dx, y dv=xe−πx2dx⟹v=2π−1e−πx2.
El término límite uv evaluado en ±∞ da 0 porque la gaussiana mata cualquier polinomio.
Nos queda ∫vdu=2π1∫e−πx2dx=2π1(1)=2π1.
En general, la recurrencia para la integral de Gauss de grado superior (o usando la Función Gamma) da:
M2k=k!(4π)k(2k)!.
Sustituyendo en la sumatoria (sabiendo que solo sobreviven los n=2k):
f^(ξ)=k=0∑∞(2k)!(−2πiξ)2kk!(4π)k(2k)!
Simplificando el factoriales (2k)!:
f^(ξ)=k=0∑∞k!4kπk(−1)k(2π)2k(i2)kξ2k
Recordando que i2=−1, el producto (−1)k(i2)k=(−1)k(−1)k=1k=1.
Agrupando términos: 4kπk(2π)2k=4kπk4kπ2k=πk.
f^(ξ)=k=0∑∞k!πkξ2k=k=0∑∞k!(πξ2)k
Identificamos esta serie como la expansión de Taylor de la exponencial alternante, pero ojo que nos olvidamos de distribuir el signo original del (−2πiξ)n. Revisemos el álgebra de signos:
El término era (−2πiξ)2k=(−1)2k(2)2k(π)2k(i)2k(ξ)2k=1⋅4k⋅π2k⋅(−1)k⋅ξ2k.
Aha! Nos queda el (−1)k vivo.
Por otro lado, calculemos la derivada del integrando pero respecto a x:
∂x∂e−π(x+iξ)2=e−π(x+iξ)2(−2π(x+iξ))
Notamos que el integrando de h′(ξ) es exactamente la derivada respecto a x multiplicada por i:
h′(ξ)=i∫−∞∞[∂x∂e−π(x+iξ)2]dx
Aplicando el Teorema Fundamental del Cálculo:
h′(ξ)=i[e−π(x+iξ)2]x→−∞x→∞
Como Re(−π(x+iξ)2)=−π(x2−ξ2), cuando x→±∞, la parte real del exponente se va a −∞, haciendo que la exponencial decaiga a 0.
Por ende, h′(ξ)=0⟹h(ξ)=C (una constante!).
Para hallar C, evaluamos en ξ=0:
h(0)=∫−∞∞e−π(x+0)2dx=∫−∞∞e−πx2dx=1 (Demostrado en Ej 15).
Por lo tanto, h(ξ)=1 para todo ξ.
Sustituyendo de vuelta:
f^(ξ)=e−πξ2⋅1=e−πξ2
(De nuevo, demostración completada de forma magistral).
[!infobox] Ejercicio 16: La Gaussiana es Autofunción de Fourier
Contexto: F(e^(-πx²)) = e^(-πξ²).
Demostración Rigurosa:
Método de completar cuadrados: Al agrupar el núcleo complejo con el argumento de la Gaussiana, queda e^(-πξ²) * h(ξ).
Para demostrar que h(ξ)=1, se deriva bajo el signo integral respecto a ξ, y se comprueba que equivale a derivar respecto a x (regla de la cadena).
Al aplicar TFC, se evalúa la Gaussiana en infinito, lo que da cero.
Por ende h'(ξ)=0 => h(ξ) es constante.
Al evaluarla en ξ=0, se recupera la integral estándar gaussiana que vale 1.
Advertencia/Clave: El paso crítico es demostrar analíticamente que integrar en una "línea compleja" paralela al eje real (x + iξ) es equivalente a integrar sobre el eje real. La derivación bajo la integral evita tener que usar Análisis Complejo (Teorema de los Residuos y contornos rectangulares).
Ejercicio 17: Fórmula de Inversión para una función Real
Intuición Didáctica
La fórmula de inversión pura de Fourier se expresa en el plano complejo y expone un número imaginario i. Sin embargo, si sabemos a priori que f(x) es una función estrictamente real en el dominio del tiempo, la matemática dictamina que las componentes imaginarias deben colapsar a cero.
Resolución Rigurosa
Arrancamos con la fórmula de inversión de la Transformada de Fourier:
f(x)=∫−∞∞f^(ξ)e2πixξdξ
Sabemos a su vez, por definición, que f^(ξ)=∫−∞∞f(y)e−2πiyξdy (cambiamos x por y para que sea variable muda interna).
Sustituyendo esto dentro de la inversión:
f(x)=∫−∞∞(∫−∞∞f(y)e−2πiyξdy)e2πixξdξ
Agrupando las exponenciales mediante propiedad de exponentes de misma base:
f(x)=∫−∞∞∫−∞∞f(y)e2πiξ(x−y)dydξ
Aplicamos la Fórmula de Euler para expandir la exponencial compleja: eiθ=cosθ+isinθ.
Dado que el enunciado asegura que f(x) es una función real, el lado izquierdo de esta ecuación es real. Por consiguiente, el lado derecho debe ser real también. Esto obliga axiomáticamente a que la parte imaginaria sea exactamente cero.
∫−∞∞∫−∞∞f(y)sin(2πξ(x−y))dydξ=0
(Nota: Otra forma analítica de probar que esto es 0, sin el axioma, es observar que la función seno es IMPAR respecto a la variable de integración exterior ξ. Al integrarla en un dominio simétrico (−∞,∞), toda la integral doble colapsa a cero independientemente de f(y)).
Habiendo descartado la parte imaginaria, nos queda:
f(x)=∫−∞∞∫−∞∞f(y)cos(2πξ(x−y))dydξ
Ahora, observemos el integrando respecto a ξ. El coseno es una función par: cos(θ)=cos(−θ).
Por ende, cos(2πξ(x−y)) es PAR respecto de la variable ξ.
La integral de una función par en (−∞,∞) es equivalente a 2 veces la integral en (0,∞).
Aplicando esta propiedad sobre la integral exterior:
f(x)=2∫0∞∫−∞∞f(y)cos(2π(x−y)ξ)dydξ
(Queda rigurosamente demostrado el postulado).
[!infobox] Ejercicio 17: Inversión en Dominio Real
Contexto: Simplificación de la Antitransformada para funciones reales.
Demostración Rigurosa:
Al agrupar la transformada dentro de la antitransformada, queda la doble integral de f(y) * e^(2πiξ(x-y)).
Al abrir por Euler en coseno + i*seno:
1. La parte imaginaria da 0 porque el seno es impar respecto a ξ y se integra sobre todo R.
2. La parte real sobrevive. El coseno es par respecto a ξ, permitiendo reducir la integral exterior de (-∞, ∞) a 2*∫[0, ∞).
Advertencia/Clave: Es crucial dominar la paridad del seno y coseno respecto a las variables de integración mudas. Lo que es impar en ξ no necesariamente es impar en y, pero la integración exterior mata todo el término.
Ejercicio 18: Condiciones para que f^ sea puramente real o imaginaria
Intuición Didáctica
La Transformada de Fourier mezcla las partes reales e imaginarias de una función mediante una rotación de fase (Euler). Si queremos que la salida (espectro) no tenga parte imaginaria, debe haber un balance o simetría perfecta en la entrada que cancele las contribuciones de los senos cruzados.
Resolución Rigurosa
Por definición: f^(ξ)=∫−∞∞f(x)e−2πixξdx.
Calculamos el complejo conjugado del espectro f^(ξ):
f^(ξ)=∫−∞∞f(x)e−2πixξdx=∫−∞∞f(x)e2πixξdx
Hacemos un cambio de variable u=−x⟹du=−dx. Límites se invierten de ∞ a −∞, y el negativo del diferencial los vuelve a acomodar.
f^(ξ)=∫−∞∞f(−u)e−2πiuξdu=∫−∞∞f(−x)e−2πixξdx
(Renombramos u a x porque es variable muda).
Caso 1: Condición para que f^ sea puramente real.
Un número complejo z es real si y solo si z=z.
Requerimos f^(ξ)=f^(ξ).
∫−∞∞f(x)e−2πixξdx=∫−∞∞f(−x)e−2πixξdx
Como esto debe cumplirse para todo ξ (y el núcleo de Fourier forma una base completa), los integrandos deben ser idénticos en casi todo punto:
f(x)=f(−x)
A esta condición se la llama Simetría Hermítica par.
Corolario: Si la función original f(x) ya sabíamos que era estrictamente real (es decir, f(x)=f(x)), entonces la condición se reduce simplemente a f(x)=f(−x). Es decir, toda función real y par tiene una Transformada puramente real.
Caso 2: Condición para que f^ sea imaginaria pura.
Un número complejo z es puramente imaginario si y solo si z=−z.
Requerimos f^(ξ)=−f^(ξ).
∫−∞∞f(x)e−2πixξdx=−∫−∞∞f(−x)e−2πixξdx
Al igual que en el caso anterior, igualando los núcleos:
f(x)=−f(−x)
A esta condición se la llama Simetría Hermítica impar.
Corolario: Si asumimos que f(x) toma exclusivamente valores reales, entonces la condición se reduce a f(x)=−f(−x). O sea, toda función real e impar tiene una Transformada puramente imaginaria.
[!infobox] Ejercicio 18: Simetrías del Espectro
Contexto: Relación de paridad entre el dominio del tiempo y la fase espectral.
Demostración Rigurosa:
Al operar con el conjugado de la definición de F[f] y aplicando un cambio de variable u = -x, se obtiene que el conjugado de la transformada es la transformada de f*(-x).
Para que F[f] ∈ R => z = z* => f(x) = f*(-x) (Simetría Hermítica Par). Si f es real, requiere que f sea PAR.
Para que F[f] ∈ I => z = -z* => f(x) = -f*(-x) (Simetría Hermítica Impar). Si f es real, requiere que f sea IMPAR.
Advertencia/Clave: No asuman que "f" siempre es real. El análisis de señales de telecomunicaciones a menudo trabaja en banda base compleja (I/Q), donde la simetría conjugada f*(-x) es lo que realmente rige el universo, no la simple paridad.
Ejercicio 19: Cálculo Analítico de Transformadas
Intuición Didáctica
Calcular transformadas por definición (resolviendo integrales de la forma xe−x2e−2πixξ) es engorroso y propenso a errores. Como ingenieros y analistas, usamos el catálogo de propiedades. Particularmente, el teorema multiplicativo de la derivada: multiplicar en el tiempo por "x" equivale (salvo constantes) a derivar el espectro en la frecuencia.
Recordatorio del teorema multiplicativo:
Si F{f(x)}=f^(ξ), entonces F{x⋅f(x)}=2πidξdf^(ξ).
Resolución Rigurosa
(a) Transformada de xe−∣x∣
Sea g(x)=e−∣x∣. Su transformada tabular (que sale fácil integrando por partes o sabiendo que es un pulso laplaciano) es:
g^(ξ)=1+(2πξ)22=1+4π2ξ22.
La función pedida es f(x)=x⋅g(x).
Aplicando la propiedad de derivación frecuencial:
F{xe−∣x∣}=2πidξd(1+4π2ξ22)
Derivamos usando regla de la cadena y del cociente (u=2,v=1+4π2ξ2⟹v′=8π2ξ):
(Nota de control de cordura: la original xe−∣x∣ es impar y real. El resultado nos dio imaginario puro e impar. ¡Cumple perfectamente lo demostrado en el Ejercicio 18!).
(b) Transformada de xe−x2
Sea g(x)=e−x2. Necesitamos su transformada base. Ojo que no es la Gaussiana estándar (e−πx2).
Aplicando propiedad de escalado F{f(ax)}=∣a∣1f^(aξ):
Si h(x)=e−πx2, sabemos que h^(ξ)=e−πξ2.
Nuestra g(x)=e−π(πx)2=h(πx). Entonces a=π1.
g^(ξ)=1/π1h^(πξ)=πe−π(πξ)2=πe−π2ξ2
Ahora, la función pedida es x⋅g(x). Aplicamos nuevamente el teorema multiplicativo:
F{xe−x2}=2πidξd(πe−π2ξ2)
Derivamos la exponencial:
dξd(πe−π2ξ2)=πe−π2ξ2(−2π2ξ)=−2π5/2ξe−π2ξ2
Sustituyendo:
f^(ξ)=2πi(−2π5/2ξe−π2ξ2)=−iπ3/2ξe−π2ξ2
(De nuevo: original impar y real -> transformada imaginaria e impar. Check!).
(c) Transformada de (4x2−2)e−x2
Esta expresión tiene términos cuadráticos. Recordemos que multiplicar por x2 equivale a aplicar la propiedad derivada espectral dos veces:
F{x2g(x)}=2πidξd(2πig^′(ξ))=−4π21dξ2d2g^(ξ).
Aprovechemos el resultado del ítem anterior, donde ya calculamos la primera derivada:
F{xe−x2}=−iπ3/2ξe−π2ξ2.
Calculemos ahora la transformada de x2e−x2. Volvemos a aplicar el operador derivación 2πidξd sobre el espectro de xe−x2:
Ese es el espectro de x2e−x2.
Nuestra función total pedida es f(x)=4x2e−x2−2e−x2. Por linealidad de Fourier, la transformada total es:
f^(ξ)=4⋅F{x2e−x2}−2⋅F{e−x2}
Reemplazamos con lo calculado y con g^(ξ) de base (ítem b):
f^(ξ)=4(2π(1−2π2ξ2)e−π2ξ2)−2(πe−π2ξ2)
Desarrollando:
f^(ξ)=2π(1−2π2ξ2)e−π2ξ2−2πe−π2ξ2
Factor común:
f^(ξ)=2πe−π2ξ2[1−2π2ξ2−1]=2πe−π2ξ2[−2π2ξ2]
Resultado final:
f^(ξ)=−4π5/2ξ2e−π2ξ2
[!infobox] Ejercicio 19: Transformadas por Derivación Espectral
Contexto: Uso del Teorema Multiplicativo: F[x f(x)] = (i / 2π) * d/dξ F[f(x)].
Demostración Rigurosa:
(a) Partiendo de F[e^(-|x|)] = 2 / (1 + 4π²ξ²), se aplica el operador derivada, obteniendo -8πiξ / (1 + 4π²ξ²)².
(b) Partiendo de F[e^(-x²)] = √π e^(-π²ξ²) (consecuencia de escala sobre la estándar), la derivada da -i π^(3/2) ξ e^(-π²ξ²).
(c) Para (4x²-2)e^(-x²), se usa la segunda derivada espectral para hallar la transformada de x²e^(-x²), y luego se compone linealmente. El álgebra produce cancelaciones milagrosas dando -4π^(5/2) ξ² e^(-π²ξ²).
Advertencia/Clave: ¡Las reglas de escala y derivada son tus mejores amigas! Ni se te ocurra meter (4x² - 2) en la integral e intentar resolverla a mano por partes; te vas a llenar de variables y perderás valioso tiempo de examen.