Docsity
Docsity

Prepara tus exámenes
Prepara tus exámenes

Prepara tus exámenes y mejora tus resultados gracias a la gran cantidad de recursos disponibles en Docsity


Consigue puntos base para descargar
Consigue puntos base para descargar

Gana puntos ayudando a otros estudiantes o consíguelos activando un Plan Premium


Orientación Universidad
Orientación Universidad


Transformadas integrales: Tranformada de fourier, Apuntes de Análisis Matemático

Descripción y características de las transformadas de fourier, comparación con sus series.

Tipo: Apuntes

2014/2015

Subido el 16/01/2015

AntonioRuiz
AntonioRuiz 🇲🇽

5

(1)

1 documento

1 / 16

Toggle sidebar

Esta página no es visible en la vista previa

¡No te pierdas las partes importantes!

bg1
Capítulo 9
La Transformada de Fourier
En este capítulo estudiaremos otra de las transformadas integrales que es ampliamente utilizada
tanto en matemáticas puras y aplicadas como en algunos campos de la ingeniería. Nos referimos a
la transformada de Fourier la cual es, además, una potente herramienta para resolver ecuaciones
en derivadas parciales. La idea básica de la transformada de Fourier en lo referente a las EDPs
es la misma que en el caso de la de Laplace, esto es, transformar un problema complicado en
otro más fácil de resolver y luego obtener la solución del problema original como la transformada
de Fourier inversa de la solución del problema transformado.
Veremos también que la transformada de Fourier es una herramienta básica en el análisis de
señales aperiódicas que tienen energía finita. En este sentido, la transformada de Fourier juega
el mismo papel que las series de Fourier para señales periódicas.
9.1 Definición y Propiedades Básicas
Antes de dar la definición rigurosa de transformada de Fourier, vamos a motivar dicha definición
a partir de las series de Fourier. Sea f:RC,una función dada. Para cada l>0pode-
mos calcular el desarrollo en serie de Fourier de la función fel cual, en su forma compleja y
suponiendo suficiente regularidad en la función f, viene dado por
f(x)= 1
2l
+
X
n=−∞
cnleinπx/l,con cnl =Zl
l
f(y)einπy/l dy.
Con el cambio ξn=nπ
llas expresiones anteriores se reescriben como
f(x)= 1
2π
+
X
n=−∞
cnleiξnxπ
l,c
nl =Zl
l
f(y)eiξnydy. (9.1)
Si suponemos que fdecae muy rápidamente cuando |x|±,entonces
cnl Z+
−∞
f(y)eiξnydy
y a medida que l→∞la primera expresión en (9.1) tiene el aspecto de una suma de Riemann.
De esta forma, haciendo tender l→∞se tiene que
f(x)= 1
2πZ
−∞ b
f(ξ)eiξxdx, con b
f(ξ)=Z+
−∞
f(y)eiξydy.
115
pf3
pf4
pf5
pf8
pf9
pfa
pfd
pfe
pff

Vista previa parcial del texto

¡Descarga Transformadas integrales: Tranformada de fourier y más Apuntes en PDF de Análisis Matemático solo en Docsity!

Capítulo 9

La Transformada de Fourier

En este capítulo estudiaremos otra de las transformadas integrales que es ampliamente utilizada tanto en matemáticas puras y aplicadas como en algunos campos de la ingeniería. Nos referimos a la transformada de Fourier la cual es, además, una potente herramienta para resolver ecuaciones en derivadas parciales. La idea básica de la transformada de Fourier en lo referente a las EDPs es la misma que en el caso de la de Laplace, esto es, transformar un problema complicado en otro más fácil de resolver y luego obtener la solución del problema original como la transformada de Fourier inversa de la solución del problema transformado. Veremos también que la transformada de Fourier es una herramienta básica en el análisis de señales aperiódicas que tienen energía finita. En este sentido, la transformada de Fourier juega el mismo papel que las series de Fourier para señales periódicas.

9.1 Definición y Propiedades Básicas

Antes de dar la definición rigurosa de transformada de Fourier, vamos a motivar dicha definición a partir de las series de Fourier. Sea f : R → C, una función dada. Para cada l > 0 pode- mos calcular el desarrollo en serie de Fourier de la función f el cual, en su forma compleja y suponiendo suficiente regularidad en la función f, viene dado por

f (x) =

2 l

+X∞

n=−∞

cnleinπx/l, con cnl =

Z (^) l

−l

f (y) e−inπy/l^ dy.

Con el cambio ξn = n lπ las expresiones anteriores se reescriben como

f (x) =

2 π

+X∞

n=−∞

cnleiξnx^

π l , cnl =

Z (^) l

−l

f (y) e−iξny^ dy. (9.1)

Si suponemos que f decae muy rápidamente cuando |x| → ±∞, entonces

cnl ≈

Z +∞

−∞

f (y) e−iξny^ dy

y a medida que l → ∞ la primera expresión en (9.1) tiene el aspecto de una suma de Riemann. De esta forma, haciendo tender l → ∞ se tiene que

f (x) =

2 π

Z ∞

−∞

f^ b (ξ) eiξx^ dx, con fb (ξ) =

Z +∞

−∞

f (y) e−iξy^ dy.

116 Capítulo 9. La Transformada de Fourier

Los coeficientes de Fourier cnl pasan a tener la forma fb (ξ) que es lo que enseguida llamaremos transformada de Fourier de la función f. Tenemos pues que la transformada de Fourier es, en el sentido anterior, lo que equivale a los coeficientes de Fourier de una función periódica. Dicho de otra forma, la transformada de Fourier juega el mismo papel para funciones definidas en todo R (o de periodo infinito si se quiere) que las series de Fourier para funciones periódicas o definidas en un intervalo acotado. Todo lo anterior motiva la siguiente definición.

Definición 9.1.1 Dada una función f : R → C, se define formalmente su transformada de Fourier como la función de variable real fb : R → C, definida como

f^ b (ξ) = 1 2 π

Z +∞

−∞

e−iξxf (x) dx, ξ ∈ R. (9.2)

Nota 9.1.1 La definición anterior es formal en el sentido de que la integral que aparece en dicha definición no tiene porque existir para una f cualquiera. La convergencia de dicha integral está garantizada en caso de ser f absolutamente integrable, es decir, si consideramos el espacio

L^1 (R) =

f : R → C :

Z +∞

−∞

|f (x)| dx < ∞

entonces la transformada de Fourier de f está bien definida siempre que f ∈ L^1 (R). Recordemos que L^1 (R) , equipado de la norma

kfk 1 =

Z +∞

−∞

|f (x)| dx

es un espacio de Banach. En realidad, la integral que aparece en (9.2) no se entiende en sentido de Riemann impropio sino en el sentido de Lebesgue. Lo mismo sucede en la definición del espacio L^1 (R). Sin embargo, si f es continua a trozos y si existen los límites y son finitos

lim a→−∞

Z 0

a

|f (x)| dx y lim b→+∞

Z (^) b

0

|f (x)| dx,

es decir, si |f | es integrable en sentido de Riemann impropio, entonces f ∈ L^1 (R) y por tanto existe la transformada de Fourier de f la cual viene dada por (9.2). A efectos de cálculo, en este capítulo entenderemos que las funciones que consideramos son absolutamente integrables en sentido de Riemann impropio.

Nota 9.1.2 Señalemos también que hay una gran diversidad en la bibliografía en lo que hace referencia al término (^21) π que aparece en la definición de la transformada de Fourier. Hay autores que incluyen esta constante y autores que prescinden de ella. Es una cuestión de gustos. Por supuesto en nada afecta ésto a la teoría que desarrollaremos a continuación.

Calcularemos ahora la transformada de Fourier de algunas funciones concretas.

Ejemplo 9.1.1 Consideremos la función característica del intervalo [a, b] , con −∞ < a < b < +∞, esto es, la función

X[a,b] (x) =

1 si a ≤ x ≤ b 0 en caso contrario

118 Capítulo 9. La Transformada de Fourier

Demostraremos ahora que (^) Z +∞

−∞

e−(x+^

iξ 2 )

2 dx =

π.

Para ello partimos de la identidad Z (^) +∞

−∞

e−x

2 dx =

π

la cual es una consecuencia inmediata del hecho de ser Z (^) +∞

−∞

Z +∞

−∞

e−(x (^2) +y (^2) ) dxdy = π.

Esta última identidad se demuestra fácilmente haciendo un cambio a coordenadas polares. Vol-

vamos al cálculo de

R +∞

−∞ e

−(x+ i 2 ξ )^2 dx. Supongamos que ξ > 0. Si ξ < 0 , el razonamiento que

sigue es completamente análogo. Para cada R > 0 consideremos el camino ΓR formado por la unión de los segmentos

Γ^1 R = {z ∈ C : z = x , − R ≤ x ≤ R} , Γ^2 R = {z ∈ C : z = R + iy , 0 ≤ y ≤ ξ/ 2 } , Γ^3 R =

n z ∈ C : z = x + i ξ 2 , − R ≤ x ≤ R

o , Γ^4 R = {z ∈ C : z = −R + iy , 0 ≤ y ≤ ξ/ 2 }

y orientado positivamente, esto es, en sentido contrario al de las agujas del reloj.

  • R (^) x

y

R

C ξ/ 2

ΓR

Dado que la función e−z 2 es holomorfa en todo el plano complejo, por el Teorema de Cauchy se tiene que

0 =

Z

ΓR

e−z 2 dz

Z

Γ^1 R

e−z 2 dz +

Z

Γ^2 R

e−z 2 dz −

Z

Γ^3 R

e−z 2 dz −

Z

Γ^4 R

e−z 2 dz.

Dado que sobre Γ^2 R y sobre Γ^4 R se tiene que ¯ ¯¯e−z^2

¯¯ = e−R^2 +y^2 , 0 ≤ y ≤ ξ/ 2 ,

entonces fácilmente se deduce que

lim R→∞

Z

Γ^2 R

e−z 2 dz = lim R→∞

Z

Γ^4 R

e−z 2 dz = 0.

9.1. Definición y Propiedades Básicas 119

Por otra parte,

lim R→∞

Z

Γ^1 R

e−z 2 dz = lim R→∞

Z R

−R

e−x 2 dx =

π

y

lim R→∞

Z

Γ^3 R

e−z 2 dz = lim R→∞

Z R

−R

e−(x+^

iξ 2 )

2 dx =

Z ∞

−∞

e−(x+^

iξ 2 )

2 dx.

Con todo ello se tiene que (^) Z ∞ −∞

e−(x+^

iξ 2 )

2 dx =

π

y por tanto, (^) h

e−x 2 i∧ (ξ) =

π

e−ξ (^2) / 4 .

Enunciamos a continuación algunas de las propiedades básicas de la transformada de Fourier. Su demostración es un sencillo ejercicio de cálculo.

Proposición 9.1.

(a) Linealidad. Sean f, g ∈ L^1 (R) y α, β ∈ C. Entonces

[αf + βg]∧^ = α fb + βbg.

(b) Translación y cambio de escala. Sean f ∈ L^1 (R) y a, b ∈ R, con a 6 = 0. Entonces la función g (x) = f (ax + b) pertenece a L^1 (R) y además

bg (ξ) =

|a|

e

iξb a (^) fb (ξ/a).

(c) Fórmula de Modulación. Sean f ∈ L^1 (R) y a ∈ R. Entonces la función eiaxf ∈ L^1 (R) y (^) h eiaxf (x)

i∧ (ξ) = fb (ξ − a).

Como consecuencia,

[f (x) cos ax]∧^ (ξ) = fb (ξ − a) + fb (ξ + a) 2 y

[f (x) sin ax]∧^ (ξ) = fb (ξ − a) − fb (ξ + a) 2 i

9.1.1 La Transformada de Fourier y la Derivación

Sea f : R → C una función de clase C^1 a trozos y supongamos que tanto f como f 0 están en L^1 (R). Del teorema fundamental del cálculo y del hecho de estar f 0 en L^1 (R) se deduce que

lim x→+∞ f (x) = f (0) + lim x→+∞

Z (^) x

0

f 0 (t) dt

9.2. El Teorema de Plancherel y la Transformada Inversa de Fourier 121

Z ∞

−∞

μ 1 2 π

Z ∞

−∞

e−iξ(x−y)f (x − y) dx

e−iξyg (y) dy

Z ∞

−∞

μ 1 2 π

Z ∞

−∞

e−iξz^ f (z) dz

e−iξyg (y) dy

= 2 π fb (ξ) bg (ξ).

En resumen: si f, g ∈ L^1 (R) , entonces

(f ∗ g)∧^ (ξ) = 2π fb (ξ) bg (ξ).

En la práctica usaremos este resultado en sentido inverso, es decir, habitualmente nos en- contraremos con fb (ξ) bg (ξ) , donde f y g son conocidas, y estaremos interesados en calcular una función h tal que bh (ξ) = fb (ξ) bg (ξ). El resultado anterior nos dice que esa tal función h no es sino (^21) π (f ∗ g).

9.2 El Teorema de Plancherel y la Transformada Inversa de Fou-

rier

En la sección anterior hemos visto que el dominio natural para definir la transformada de Fourier es el espacio de funciones L^1 (R). Pero, ¿cuál es su imagen o rango? ¿es inyectiva la transformada de Fourier? ¿y sobreyectiva?. En este caso, ¿cómo se calcula la transformada inversa?. De responder a todas estas cuestiones nos ocuparemos en esta sección. Consideremos en primer lugar el espacio

C∞ (R) =

f : R → C , f continua y lim |x|→∞

f (x) = 0

el cual equipado con la norma

kfk∞ = sup {|f (x)| : x ∈ R}

es un espacio de Banach. El resultado que sigue establece que la transformada de Fourier del espacio L^1 (R) está contenida en C∞ (R). Para una demostración véase [8, p. 217].

Lema 9.2.1 (Riemann-Lebesgue) Si f ∈ L^1 (R) , entonces fb ∈ C∞ (R). Además, °° ° fb

2 π

kfk 1.

Nos ocuparemos ahora del problema de recuperar la función f a partir de su transformada de Fourier f.b Para ello, dada f ∈ L^1 (R) se define su transformada inversa de Fourier como

∨ f (x) =

Z ∞

−∞

eiξxf (ξ) dξ

Un primer resultado de inversión es el siguiente Teorema de Inversión de Fourier. Una demos- tración puede encontrarse en [8, p.p. 218-219].

Teorema 9.2.2 Sea f ∈ L^1 (R) y supongamos que fb ∈ L^1 (R). Entonces f =

fb

, es decir,

f (x) =

Z ∞

−∞

eiξx^ fb (ξ) dξ para todo x ∈ R. (9.3)

122 Capítulo 9. La Transformada de Fourier

Como consecuencia de este teorema se tiene que la transformada de Fourier es inyectiva. De manera precisa:

Corolario 9.2.3 Sean f, g ∈ L^1 (R) tales que fb = bg. Entonces f = g.

La fórmula (9.3) expresa una función f como una superposición continua de funciones expo- nenciales. En este sentido, la transformada de Fourier proporciona un análogo para funciones no periódicas de las series de Fourier para funciones periódicas. Son muchas las funciones de L^1 cuya transformada de Fourier está también en L^1. Por ejem- plo, si f ∈ C^2 (R) y si f, f 0 , f 00 ∈ L^1 (R) , entonces (f 00 )∧^ (ξ) = −ξ^2 fb (ξ) es acotada. De esta

forma,

¯ fb^ (ξ)

¯ ≤^ |Cξ| 2 ,^ lo que garantiza que^ fb^ ∈^ L^1 (R)^.^ Sin embargo, es importante señalar tam- bién que la transformada de Fourier de algunas funciones importantes en las aplicaciones, como la considerada en el Ejemplo 9.1.1, no están en L^1 (ver Ejercicio 3). Resulta pues interesante obtener alguna otra versión del teorema de inversión anterior. Presentamos a continuación una de estas versiones que hace uso del concepto de valor principal (de Cauchy) de una integral. Dada una función f : R → C, se llama valor principal de f al número

lim M→∞

Z M

−M

f (x) dx

que denotamos por P.V.

R ∞

−∞ f^ (x)^ dx.^ Es importante señalar que puede existir y ser^ finito P.V.

R ∞

−∞ f^ (x)^ dx^ y sin embargo no existir^

R ∞

−∞ f^ (x)^ dx.^ Sin embargo, si^ f^ ∈^ L

(^1) (R) , enton-

ces existe P.V.

R ∞

−∞ f^ (x)^ dx^ y además Z (^) ∞

−∞

f (x) dx = P.V.

Z ∞

−∞

f (x) dx.

Ya estamos en condiciones de enunciar una segunda versión del Teorema de Inversión de Fourier. Para la demostración véase [8, p.p. 220-221].

Teorema 9.2.4 Sea f ∈ L^1 (R) y supongamos que f es diferenciable a trozos. Entonces

P.V.

Z ∞

−∞

eiξx^ fb (ξ) dξ = f (x+) + f (x−) 2

donde f (x+) = lim h h>→ 00 f (x + h) y f (x−) = lim h h>→ 00 f (x − h).

Nota 9.2.1 Nótese que dentro del teorema anterior si que tienen cabida funciones como la considerada en el Ejemplo 9.1.1.

Hasta ahora hemos desarrollado toda la teoría de la transformada de Fourier en el marco del espacio L^1 (R). Sin embargo, nuestra experiencia con las series de Fourier sugiere que el espacio L^2 (R) también ha de jugar algún papel importante en toda esta teoría de la transformada de Fourier. Esto es efectivamente cierto. Como veremos enseguida, es L^2 (R) y no L^1 (R) el espacio en el que la transformada de Fourier tiene un mejor comportamiento. Veamos, en primer lugar, que la identidad de Parseval de las series de Fourier tiene también su análogo para la transformada de Fourier. Supongamos que f, g, f ,b gb ∈ L^1 (R). Entonces no

124 Capítulo 9. La Transformada de Fourier

9.3.1 Teoría de la Señal

Supongamos que f (t) representa la amplitud de una señal (una onda de sonido o una onda electromagnética, por ejemplo) en el instante de tiempo t. La fórmula

f (t) =

Z

fb (ω) eiωt^ dt , fb (ω) = 1 2 π

Z

f (t) e−iωt^ dt

describe a f como una superposición de las ondas periódicas simples eiωt, donde ω cubre el rango de todas las posibles frecuencias. Esta representación es básica en el análisis de señales en ingeniería eléctrica. En esta sección nos limitaremos a mostrar algunas ideas subyacentes a la Teoría de la Señal y, en particular, veremos como las series de Fourier y la transformada de Fourier son una herramienta básica en este campo de la ingeniería. En primer lugar, la potencia de la señal f (t) es proporcional al cuadrado de su amplitud |f (t)|^2 , y por tanto, la energía total de una señal es proporcional a

R

|f (t)|^2 dt. De esta forma, la condición de ser finita la energía total de la señal f significa simplemente que f ∈ L^2. Vemos, por tanto, que L^2 es un espacio que no sólo tiene muy buenas propiedades desde un punto de vista matemático (es un Hilbert, y eso es mucho) sino que también tiene una realidad física clara: es, en particular, el espacio de las señales que tienen energía finita. En segundo lugar, los sistemas eléctricos pueden ser modelizados matemáticamente como operadores A que transforman una señal de entrada f en una señal de salida Af. Muchos de estos sistemas son (al menos en una primera aproximación) lineales, lo cual significa que el operador A es lineal. Además, la acción de estos sistemas no se ve afectada por el paso del tiempo (para una misma entrada f, la señal de salida es la misma hoy que mañana). Esto significa que el operador A conmuta con las translaciones, es decir, si (Af ) (x) = g (x) , entonces (Af ) (x + s) = g (x + s) para cualquier s ∈ R. Se puede demostrar (ver [8, p.p. 225-226]) que entonces el operador A es del tipo

Af =

h fb

, o lo que es equivalente, Af = H ∗ f (9.5)

donde h : R → C y H es su transformada inversa de Fourier. La función h se denomina función del sistema y H se suele llamar la respuesta impulso. Si escribimos h (ω) en forma polar como h (ω) = A (ω) eiθ(ω), entonces A (ω) y θ (ω) representan, respectivamente, la amplitud y la fase debidos al operador A en la frecuencia ω. Los operadores del tipo (9.5) han sido y son amplia- mente estudiados en matemáticas puras, especialmente a partir de los trabajos de A. Zygmund y A. Calderón quienes crearon a principios del siglo XX en la universidad de Chicago una muy fructífera escuela de matemáticos que se ocupó, entre otras cosas, del estudio de este tipo de operadores. Pero volvamos a las señales. Supongamos que estamos interesados en conocer una señal f a partir del conocimiento de los valores de la amplitud de dicha señal en una sucesión de tiempos t 1 < t 2 < · · ·. ¿Es ésto posible? Por supuesto, para una función arbitraria f (t) , el conocer f (t 1 ) , f (t 2 ) , · · · nos dice absolutamente nada sobre f (t) para un t 6 = t 1 , t 2 , · · ·. Sin embargo, si la señal f involucra frecuencias ω limitadas es un cierto rango, es decir, si ω ∈ [−L, L], para un cierto L > 0 , entonces sí que vamos a poder describir apropiadamente la señal f. Para ello necesitaremos hacer uso de técnicas que involucran a las series de Fourier y a la transformada de Fourier. Antes de seguir adelante, diremos que una señal f es limitada en banda si fb (ω) = 0 para ω ∈/ [−L, L]. Si f ∈ L^2 y es limitada en banda, entonces fb ∈ L^2 (teorema de Plancharel) y además fb (ω) = 0 para ω ∈/ [−L, L]. Por tanto, fb ∈ L^1 y vale el Teorema de Inversión de Fourier.

9.3. Algunas Aplicaciones de la Transformada de Fourier 125

Recordemos también que dada una función 2 L−periódica f, su serie de Fourier asociada se escribe en la forma A 0 2

X^ ∞

n=

An cos nπx L

  • Bn sin nπx L

Usando las identidades trigonométricas

cos

nπx L

ei^ nLπx

  • e−i^ nLπx

2 , sin

nπx L

ei^ nπLx − e−i^ nπLx

2 i

la expresión (9.6) se transforma en X∞

n=−∞

cneinπx/L, (9.7)

donde

c 0 =

A 0

, cn = An − iBn 2

y c−n = An + iBn 2

Con todo esto en mente ya estamos en condiciones de enunciar y probar uno de los resultados más importantes en Teoría de la Señal, el Teorema del Muestreo de Shannon.

Teorema 9.3.1 (del muestreo de Shannon) Supongamos que f es continua a trozos, f ∈ L^2 (R) ∩ L^1 (R) y que fb (ω) = 0 para ω ∈/ [−L, L]. Entonces

f (t) =

X^ ∞

n=−∞

f

³ (^) nπ L

´ (^) sin (Lt − nπ) Lt − nπ

es decir, f está completamente determinada por los valores de f en los puntos n Lπ.

La demostración es muy formativa pues hace intervenir algunos de los resultados más im- portantes del Análisis de Fourier. La esbozamos a continuación. Como fb ∈ L^2 , podemos calcular la serie de Fourier asociada a f ,b la cual en su forma compleja viene dada por

f^ b (ω) =

X^ ∞

n=−∞

c−ne−inπω/L^ (|ω| ≤ L) (9.8)

donde los coeficientes de Fourier c−n están dados por

c−n =

2 L

Z L

−L

f^ b (ω) einπω/L^ dω = 1 2 L

f

³ (^) nπ L

Nótese que la identidad (9.8) es consecuencia del teorema de convergencia puntual para series de Fourier (ya que fb es diferenciable a trozos al ser f continua a trozos, y fb es continua por el Lema de Riemann-Lebesgue) y que (9.9) es consecuencia del Teorema de Inversión de Fourier (ya que f y fb están en L^1 (R)). Usando de nuevo este mismo argumento se tiene que

f (t) =

Z L

−L

f^ b (ω) eiωt^ dω

2 L

Z L

−L

X^ ∞

n=−∞

f

³ (^) nπ L

e−inπω/Leiωt^ dω

2 L

X^ ∞

n=−∞

f

³ (^) nπ L

´ (^) ei(Lt−nπ)ω/L i (Lt − nπ) /L

L

−L

=

X^ ∞

n=−∞

f

³ (^) nπ L

´ (^) sin (Lt − nπ) Lt − nπ

9.4. Ejercicios 127

Usando el teorema de derivación para integrales paramétricas y suponiendo, por ejemplo, que f ∈ L^1 (R) , se puede probar que la función dada en (9.10) es solución clásica de la ecuación del calor. Esto es debido al rápido decaimiento del núcleo de Dirichlet cuando |x| → ±∞. De nuevo aquí nos encontramos con el efecto regularizante de la ecuación del calor. Para obtener unicidad y estabilidad de soluciones es preciso suponer, además, que f está acotada.

El problema de Dirichlet para la ecuación de Laplace en un semiplano

Consideremos ahora el problema de Dirichlet en un semiplano, es decir, se trata de encontrar una función u = u (x, y) , con −∞ < x < ∞ e y > 0 , de clase C^2 y que satisfaga la EDP y la condición de frontera ½ uxx (x, y) + uyy (x, y) = 0, −∞ < x < ∞, y > 0 u (x, 0) = f (x) , −∞ < x < ∞

Para obtener unicidad de soluciones es preciso suponer alguna condición de acotación sobre u. El motivo es sencillo: si u (x, y) es solución del problema anterior, entonces la función u (x, y) + y también es solución de dicho problema. Supondremos, por el momento, que f es acotada e integrable y buscaremos una solución acotada del problema anterior. Al igual que en el caso de la ecuación del calor, si aplicamos la transformada de Fourier a la ecuación de Laplace y a la condición de contorno obtenemos el problema ( −ξ^2 ub (ξ, y) + (ub (ξ, y))yy = 0, y > 0 bu (ξ, 0) = fb (ξ)

La solución de este problema es

bu (ξ, y) = c 1 (ξ) e|ξ|y^ + c 2 (ξ) e−|ξ|y

donde además se satisface que c 1 (ξ) + c 2 (ξ) = fb (ξ)

Con el fin de obtener una solución acotada ponemos c 1 = 0. Como la transformada inversa de e−|ξ|y^ es la función

Py (x) = 2 y x^2 + y^2

que se denomina núcleo de Poisson, la solución de nuestro problema original es

u (x, y) =

2 π

(f ∗ Py) (x) =

Z ∞

−∞

y π (t^2 + y^2 )

f (x − t) dt.

Se puede demostrar que si f es continua y acotada, entonces esta función es la única solución clásica del problema (9.11).

9.4 Ejercicios

  1. Demostrar el Corolario 9.2.3.
  2. Comprobar que la función f (x) = x, x ∈ R, no está en L^1 (R) , y que sin embargo P.V.

R ∞

−∞ x dx^ = 0.

128 Capítulo 9. La Transformada de Fourier

  1. Demuestra que para cada n ∈ N se tiene que Z (^) (n+1)π

¯¯^ sin^ x x

¯¯ dx ≥ 2 (n + 1) π

Deducir de ello que la función f (x) = sinx^ xno está en L^1 (R).

  1. Probar que si f, fb ∈ L^1 (R) , entonces f ∈ L^2 (R).
  2. Demostrar que la transformada de Fourier de la función f (x) = e−|x|^ es fb (ξ) = (^) π^1 (1+ξ^2 )
  1. Sea f ∈ L^2 (R) y supongamos que f (x) = 0 para todo |x| > L > 0. Probar que entonces f ∈ L^1 (R). Indicación: usar la desigualdad de Cauchy-Schwartz para espacios de Hilbert.
  2. Consideremos las funciones

f (x) =

x−^2 /^3 si 0 < x < 1 0 si |x| ≥ 1 y g (x) =

x−^2 /^3 si x > 1 0 si x ≤ 1

Comprueba que f ∈ L^1 (R) , f /∈ L^2 (R) y que g /∈ L^1 (R) , g ∈ L^2 (R).

  1. Usar la transformada de Fourier para resolver el siguiente problema de EDPs para la ecuación del calor   

ut (t, x) = 4uxx (t, x) , t > 0 , − ∞ < x < ∞ u (0, x) =

1 si |x| ≤ 2 0 si |x| > 2

  1. Para cada ε > 0 consideremos la función

fε (x) =

½ (^) π ε si^ −ε^ ≤^ x^ ≤^ ε 0 si |x| > ε

Demostrar que para cada ξ ∈ R, fbε (ξ) converge a 1 cuando ε → 0. Por ello se suele decir que la transformada de Fourier de la delta de Dirac es 1. Sin embargo, esto contradice la conclusión del Lema de Riemann-Lebesgue que afirma que la transformada de Fourier de cualquier función de L^1 se anula en infinito. El motivo de ésto es que la delta de Dirac no es una función de L^1.

  1. Consideremos el siguiente problema de valores iniciales para la ecuación de ondas ½ utt (t, x) = c^2 uxx (t, x) , t > 0 , − ∞ < x < ∞ u (0, x) = f (x) ; ut (0, x) = g (x) , −∞ < x < ∞

Suponiendo que todas las transformadas de Fourier existen, demuestra que

u b (t, ξ) = fb (ξ) cos (ctξ) + gb (ξ) (cξ)−^1 sin (ctξ).

Usando las propiedades de la transformada de Fourier, deduce la fórmula de d’Alembert para la ecuación de ondas, esto es,

u (t, x) =

f (x + ct) + f (x − ct) +

c

Z (^) x+ct

x−ct

g (y) dy

130 Capítulo 9. La Transformada de Fourier

Probar que {em}N m−=0^1 es una base ortogonal de CN^ , y que kemk^2 = N para todo m. Una vez probado este resultado se tiene que para todo a ∈ CN^ ,

a =

N

NX− 1

m=

< a, em > em.

Probar que entonces las componentes del vector a, que denotamos por an, están dadas por

an =

N

NX− 1

m=

e^2 πimn/N^ bam,

donde bam =< a, em >. La fórmula anterior es la fórmula de inversión para la transformada discreta de Fourier.

9.5 Objetivos

  • Estudiar las propiedades básicas de la transformada de Fourier.
  • Mostrar la aplicación de la transformada de Fourier en la resolución de EDPs y en proble- mas relacionados con la Teoría de la Señal.

9.6 Comentarios sobre la Bibliografía

En la exposición de este tema hemos seguido básicamente el libro de Folland [8, Ch. 7]. Una presentación también muy directa y asequible al nivel del alumno medio es la dada en [19, Ch. 3] donde se trabaja en todo momento con funciones continuas a trozos y con la integral de Riemann. Es, por tanto, una excelente referencia. Como se ha mencionado a lo largo del capítulo, la transformada de Fourier hace un uso continuo de las propiedades de la integral de Lebesgue. Por ello, en la exposición del tema no hemos justificado adecuadamente algunos cálculos. El alumno que desee conocer todos estos matices matemáticos ha de acudir necesariamente a libros más especializados como por ejemplo [1, 10, 17]. El libro de Gasquet-Witomski está especialmente enfocado a la aplicación de la transformada de Fourier a la Teoría de la Señal y, aunque el nivel matemático es bastante elevado, puede ser interesante para conectar los aspectos teóricos de la transformada de Fourier con algunas de sus aplicaciones.