Apunte Fourier
Apunte Fourier
Departamento de Matemática
Notas mat – 023
Series de Fourier
Segundo Semestre 2012
1. Funciones Ortogonales
Definición 1.1 El conjunto de funciones seccionalmente continuas {φn (x) : n = 0, 1, . . . } en un intervalo finito
[α, β] o infinito, se dice ortogonal en [α, β] con respecto a la función no negativa r(x) si
Z β
r(x)φm (x)φ̄n (x) dx = 0 ∀ m 6= n
α
(1.1) hφm , φn i =
β
Z
r(x)|φn (x)|2 dx 6= 0
∀n
α
Z β
r(x)|φn (x)|2 dx = 1 ∀ n.
α
Ası́, las funciones ortonormales poseen las mismas propiedades que las funciones ortogonales, pero además, han
sido normalizadas, i.e., cada función φn (x) del conjunto ortogonal ha sido dividida por su propia norma, que se
define como
!1/2
Z β
2
(1.2) kφn k = r(x)|φn (x)| dx .
α
JCM/jcm. [email protected].
1
es ortogonal sobre 0 < x < π con r(x) = 1. En efecto, para m 6= n, tenemos
Z π
1 π
Z
cos mx cos nx dx = [cos(m − n)x + cos(m + n)x] dx
0 2 0
π
1 sen(m − n)x sen(m + n)x
= +
2 m−n m+n
0
1 sen(m − n)π sen(m + n)π
= −0+ − 0 = 0.
2 m−n m+n
Ahora, puesto que
π 1/2
√
Z
kφ0 k = dx = π
0
Z π 1/2 Z π 1/2 r
2 1 + cos 2nx π
kφn k = cos nx dx dx = dx = ,
0 0 2 2
sigue que el conjunto ( )
r
1 2
√ , cos nx, n = 1, 2, . . . ,
π π
es ortonormal en 0 < x < π con r(x) = 1.
Ejemplo 1.2 El conjunto
1 1 1
(1.3) √ , √ cos nx, √ sen nx, n = 1, 2, . . .
2π π π
es ortonormal sobre −π < x < π con r(x) = 1. Para esto, primero necesitamos chequear que
Z π 2 Z π 2 Z π 2
1 1 1
√ dx = 1, √ cos nx dx = 1, √ sen nx dx = 1.
−π 2π −π π −π π
La primera integral es evidente. En tanto, las siguientes se calculan usando el ejemplo anterior más la propiedad
Rp Rp
de integración de funciones pares en intervalos simétricos: −p h(z)dz = 2 0 h(z)dz. En segundo lugar, debemos
chequear que
π π
1 1 1 1
Z Z
√ × √ cos nx dx = 0, √ × √ sen nx dx = 0,
Z π−π 2π π −π 2π π
π
1 1 1 1
Z
√ cos nx √ cos mx dx = 0, √ sen nx √ sen mx dx = 0, m 6= n
−π π π −π π π
y
π
1 1
Z
√ cos nx √ sen mx dx = 0 m 6= n.
−π π π
Es fácil chequear esto recordando que todas las funciones involucradas son trigonométricas integradas en el
perı́odo (−π, π).
Ejemplo 1.3 El conjunto de funciones a valores complejos
1
(1.4) √ einx , n ∈ Z ,
2π
es ortonormal en el intervalo −π < x < π con la función peso r(x) = 1. En efecto
Z π Z π
1 1
kφn k2 = einx e−inx dx = dx = 1.
2π −π 2π −π
2
Además, usando la fórmula de Euler si m 6= n:
Z π Z π
1 1
hφm , φn i = eimx e−inx dx = ei(m−n)x dx
2π −π 2π −π
π
1 1 n o
ei(m−n)x = ei(m−n)π − e−i(m−n)π
=
2π(m − n)i 2π(m − n)i
−π
1
= {cos(m − n)π + i sen(m − n)π − cos(m − n)π + i sen(m − n)π} = 0.
2π(m − n)i
Ejemplo 1.4 Los polinomios de Chebyshev de primer orden
son ortogonales en el intervalo −1 < x < 1 con la función peso r(x) = (1 − x2 )−1/2 . En efecto
Z 1
1
hTm , Tn i = √ cos(m arc cos x) cos(n arc cos x) dx
−1 1 − x2
Z π
= cos mθ cos nθ dθ (x = cos θ)
0
0, m 6= n
= π/2, m = n 6= 0
π,
m = n = 0.
2. Series de Fourier
A menudo necesitaremos expandir una función dada en términos de otras funciones con una exactitud es-
pecı́fica de tal forma de realizar operaciones en la práctica. En esta sección mostraremos que los conjuntos de
funciones que hemos estudiado en la sección anterior pueden ser utilizados efectivamente como la base en la
expansión de funciones más generales.
La base {e1 , e2 , . . . , en } (ek es un vector unitario) de Rn posee una importante caracterı́stica: para cualquier
Pn
u ∈ Rn existe una elección única de constantes α1 , . . . , αn para las cuales u = i
i=1 αi e . Más aún, por la
ortonormalidad de los vectores ei , 1 < i < n podemos determinar estas constantes αi , 1 < i < n como sigue
n n
!
X X
j
u·e = αi e · ej =
i
αi (ei · ej ) = αj , 1 ≤ j ≤ n.
i=1 i=1
Ası́, bajo condiciones adecuadas de convergencia los coeficientes constantes cn , n = 0, 1, . . . están dados por la
fórmula
,
Z β
(2.2) cn = r(x)φ̄n (x)f (x) dx kφn k2 .
α
Sin embargo, si el conjunto {φn (x)} es ortonormal, es decir kφn k = 1, entonces tenemos
Z β
(2.3) cn = r(x)φ̄n (x)f (x) dx.
α
P∞
Si la serie n=0 cn φn (x) converge uniformemente a f (x) en [α, β], entonces el procedimiento previo está justi-
ficado, y entonces los coeficientes cn están dados por 2.2.
Los coeficientes cn obtenidos en 2.2 son llamados coeficientes de Fourier de la función f (x) con respecto al
P∞
conjunto ortogonal {φn (x)} y la serie n=0 cn φn (x) con coeficientes 2.2 es llamada la Serie de Fourier de f (x).
Vamos a escribir
∞
X
f (x) ∼ cn φn (x)
n=0
P∞
lo cual, en general, es sólo una correspondencia, i.e., a menudo f (x) 6= n=0 cn φn (x), a menos que se pruebe.
A continuación, listamos algunos de los casos con más frecuencia utilizados.
2.1. Serie de Fourier de Cosenos. En la sección pasada hemos probado que el conjunto de funciones
( r )
1 2
φ0 (x) = √ , φn (x) = cos nx, n = 1, 2, . . . ,
π π
es ortonormal en 0 < x < π. Ası́, para cualquier función seccionalmente continua f (x) definida sobre 0 < x < π,
∞ ∞ r
X c0 X 2
f (x) ∼ c0 φ0 (x) + cn φn (x) = √ + cn cos nx,
n=1
π n=1 π
donde
π π
r
1 2
Z Z
c0 = f (t) √ dt, cn = f (t) cos nt dt, n ≥ 1.
0 π 0 π
Luego,
π ∞ Z π
1 2 2
Z X
f (x) ∼ × f (t) dt + f (t) cos nt dt cos nx,
2 π 0 n=1
π 0
4
lo que puede ser escrito como
∞
a0 X
(2.4) f (x) ∼ + an cos nx,
2 n=1
donde
π
2
Z
(2.5) an = f (t) cos nt dt, n ≥ 0.
π 0
Esta serie es conocida como serie de Fourier de cosenos.
Ejemplo 2.1 Vamos a encontrar la serie de Fourier de cosenos de la función f (x) = x, 0 < x < π. Claramente
2 π 2 π
Z Z
a0 = t cos 0t dt = t dt = π,
π 0 π 0
π Z π
2 π
2 sen nt sen nt
Z
an = t cos nt dt = t − dt
π 0 π n 0 0 n
2 (−1)n − 1
2 0 − 0 cos nπ
= + = , n ≥ 1.
π n n2 π n2
Ası́, de 2.4, tenemos
∞
π X 2 (−1)n − 1
x∼ + cos nx, 0 < x < π.
2 n=1 π n2
En la figura se muestra la gráfica de la función y la aproximación por la serie truncada para diversos valores.
3.0
2.5
2.0
1.5
1.0
0.5
0.0
0.0 0.5 1.0 1.5 2.0 2.5 3.0
Figura 2.1 Gráfico de f (x) = x (púrpura) y de las aproximaciones con 1 coeficiente (rojo), 4 coeficientes (naranja) y 10 coeficientes
(verde). A medida que se toman más coeficientes de la expansión, la gráfica de la serie se ajusta a la función, como se esperaba.
donde
π
2
Z
(2.7) bn = f (t) sen nt dt, n ≥ 1.
π 0
2.3. Serie de Fourier Trigonométrica. En la sección pasada, hemos verificado que el conjunto
1 1 1
φ0 (x) = √ , φ2n−1 (x) = √ cos nx, φ2n (x) = √ sen nx, n = 1, 2, . . .
2π π π
es ortonormal sobre −π < x < π. Ası́, para cualquier función seccionalmente continua f (x) definida sobre
−π < x < π,
∞
c0 X cos nx sen nx
f (x) ∼ √ + c2n−1 √ + c2n √ ,
2π n=1 π π
donde
π
1
Z
c0 = f (t) √ dt
−π 2π
Z π
1
c2n−1 = f (t) √ cos nt dt, n≥1
−π π
Z π
1
c2n = f (t) √ sen nt dt, n ≥ 1.
−π π
Al reescribir las relaciones de arriba, tenemos la serie de Fourier trigonométrica
∞
a0 X
(2.8) f (x) ∼ + (an cos nx + bn sen nx) ,
2 n=1
donde
π
1
Z
(2.9a) an = f (t) cos nt dt, n ≥ 0,
π −π
π
1
Z
(2.9b) bn = f (t) sen nt dt, n ≥ 1.
π −π
6
Para nuestra discusión posterior notamos las siguientes relaciones
√ √
r
π
(2.10) c0 = a0 , c2n−1 = πan , c2n = πbn .
2
Ejemplo 2.3 Vamos a encontrar la serie de Fourier trigonométrica de la función
(
1, si −π < x < 0
f (x) =
x, si 0 < x < π.
De 2.9a-2.9b, tenemos
π 0 π
1 1 1
Z Z Z
an = f (t) cos nt dt = 1 · cos nt dt + t cos nt dt
π −π π −π π 0
π
1 t sen nt cos nt
= 0+ +
π n n2 0
1 (−1)n − 1
= , n ≥ 1.
π n2
0 π
1 1 π
Z Z
a0 = 1 · 1 dt + t · 1 dt = 1+ ,
π −π π 0 2
0 π
1 1 −1 + (1 − π)(−1)n
Z Z
bn = 1 · sen nt dt + t sen nt dt = , n ≥ 1.
π −π π 0 nπ
3.0
2.5
2.0
1.5
1.0
0.5
0.0
-3 -2 -1 0 1 2 3
Figura 2.2 Función f (x) y su serie de Fourier trigonométrica para 2 coeficientes (rojo), 6 (naranja) y 100 coeficientes (verde).
Note la oscilación de la serie en los extremos y en los puntos de discontinuidad.
Es claro que las constantes an , bn para la serie de Fourier trigonométrica son diferentes de aquellos obtenidos
para la serie de Fourier de cosenos y senos. Sin embargo, si la función f (x) es impar, entonces
1 π
Z
an = f (t) cos nt dt = 0, n≥0
π −π
1 π 2 π
Z Z
bn = f (t) sen nt dt = f (t) sen nt dt, n≥1
π −π π 0
7
y la serie de Fourier trigonométrica se reduce a la serie de Fourier de senos. Ası́, concluimos que si f (x) es impar,
o definida sólo sobre (0, π) y hacemos su extensión impar entonces la serie de Fourier de senos es válida sobre
(−π, π). Exactamente en el mismo sentido, si f (x) es par, o definida sólo sobre (0, π) y hacemos su extensión
par entonces
π π
1 2
Z Z
an = f (t) cos nt dt = f (t) cos nt dt, n≥0
π −π π 0
π
1
Z
bn = f (t) sen nt dt = 0, n≥1
π −π
π
−(x + π) si −π ≤ x ≤ − 2
f (x) = x si − π2 ≤ x ≤ π
2
π
π−x si < x ≤ π,
2
Π
Puesto que la función es impar, sigue que a0 = an = 0, 2
+ π sen nx dx . Figura 2.3 La función f (x) del ejemplo. Note que es impar.
π/2
R sen nz z cos nz
Las dos primeras integrales se calculan por partes y se obtiene z sen nz dz = − . Luego
n2 n
2 sen nx π/2 x cos nx π/2 x cos nx π sen nx π π cos nx π
bn = − + − −
π n2 0 n n n2 π/2 n
0 π/2 π/2
2 2 sen nπ π cos nπ
π cos nπ
π cos nπ
2 2 π cos nπ 2 π cos nπ 2
= − + − − +
π n2 2n n 2n n n
4 nπ
= 2
sen .
πn 2
8
Estos coeficientes son distintos de cero si n es impar, o sea, n = 2k − 1, k ∈ N. Utilizando la identidad
trigonométrica del seno de la diferencia de dos ángulos, sigue que sen(2k − 1) π2 = − cos kπ = (−1)k+1 . Entonces
4 (−1)k+1
bk = ; y
π (2k − 1)2
∞
4 X (−1)k+1
f (x) ∼ sen(2k − 1)x .
π (2k − 1)2
k=1
Finalmente, hacemos notar que si la función seccionalmente continua f (x) está definida sobre −L < x < L,
entonces su serie de Fourier trigonométrica 2.8 toma la forma
∞
a0 X nπx nπx
(2.11) f (x) ∼ + an cos + bn sen ,
2 n=1
L L
donde
L
1 nπt
Z
(2.12a) an = f (t) cos dt, n≥0
L −L L
L
1 nπt
Z
(2.12b) bn = f (t) sen dt, n ≥ 1.
L −L L
2.3.1. Expansión de Medio Rango. Suponga que una función f está definida sobre 0 < x < π y que deseamos
expandirla en una serie de Fourier trigonométrica. Dado que los coeficientes an y bn involucran integrales desde
−π a π, debemos de alguna manera extender la definición de f al intervalo (−π, π). Podemos hacer esto de varias
maneras a nuestro beneficio. Una forma, es extender f tal que sea una función par sobre el intervalo (−π, π)
(figura 2.4). Dado que es una función par, se tiene que bn = 0, n = 1, 2, . . . y la serie de Fourier trigonométrica
sólo posee términos de cosenos, es decir, coincide con la serie de Fourier de cosenos antes estudiada. Esta
expansión es denominada expansión de medio rango.
También, una función definida en el intervalo (0, π) puede también ser extendida al intervalo (−π, π) como una
función impar. La figura 2.5 muestra tal extensión, y notamos que se introducen discontinuidades en x = −π, 0, π,
a menos que f (−π) = f (0) = f (π) = 0. esto no nos preocupa, puesto que como se verá en la siguiente sección,
los teoremas de convergencia para series de Fourier son válidos para funciones seccionalmente continuas. Dado
que esta expansión es impar, se tiene que an = 0, n = 0, 1, 2, . . . , la serie resultante es una serie de Fourier de
senos.
y y
x
-Π Π
x
-Π Π
Figura 2.4 Extensión par de f (x) al intervalo (−π, π). 2.5 Extensión impar de f (x) al intervalo (−π, π).
9 Figura
Ejemplo 2.5 Vamos a expandir de forma par la función
(
0, si 0 < x < π/2
f (x) =
1, si π/2 < x < π.
1
El gráfico de la extensión par se muestra en la fi-
gura 2.6. Puesto que es una expansión par, sigue que
bn = 0, n = 1, 2, . . . Además
2 π 2 π
Z Z
a0 = f (t) dt = dt = 1,
π 0 π π/2
y
2 π 2 π
Z Z
ak = f (t) cos kt dt = cos kt dt
π 0 π π/2
π
2 2 kπ
= =− sen
kπ π/2 kπ 2 -Π -
Π Π Π
2 2
Los coeficientes son cero, salvo los impares, i.e., si Figura 2.6 Extensión par de f (x).
k = 2n − 1, n = 1, 2, . . . Usando la identidad
π
sen nπ − = − cos nπ = (−1)n+1 , obtenemos
2
0, si k = 2n, n = 1, 2, . . .
ak = 2 (−1)n
, si k = 2n − 1, n = 1, 2, . . .
π 2n − 1
Por lo tanto,
∞
1 2 X (−1)n
f (x) ∼ + cos nx, −π < x < π.
2 π n=1 2n − 1
Se deja como ejercicio realizar la expansión impar.
Si deseamos expandir una función seccionalmente continua definida sobre el intervalo [c − π, c + π], formamos
su extensión periódica f0 como sigue
(
f (x), si c − π < x < c + π,
f0 (x) =
f0 (x − 2π), para cualquier otro x.
No importa cuan suave sea la función f , estaremos introduciendo discontinuidades en x = c − π y x = c + π
al menos que f (c − π) = f (c + π). Dado que tanto f0 (x), cos nx y sen nx son todas funciones periódicas con
periodo 2π, podemos reemplazar el intervalo [−π, π] por [c − π, c + π], sobre el cual f0 coincide con f . Entonces,
los coeficientes de Fourier están dados por
c+π
1
Z
(2.13a) an = f (t) cos nt dt,
π c−π
c+π
1
Z
(2.13b) bn = f (t) sen nt dt.
π c−π
La serie de Fourier formada con los coeficientes en 2.13a-2.13b entregan una expansión que es válida sobre el
intervalo c − π < x < c + π.
10
Ejemplo 2.6 Vamos a expandir la función f (x) = x sobre 0 < x < 2π. El gráfico de la extensión periódica se
muestra en la figura 2.7. Calculamos
2 Π
2π
1
Z
a0 = t dt = 2π,
π 0
2π 2π Π
1 1 t sen nt cos nt
Z
an = t cos nt dt = + = 0,
π 0 π n n2 0
2π
1 2π
1 t cos nt sen nt
Z
bn = t sen nt dt = − +
π 0 π n n2 0
2 -2 Π -Π Π 2 Π 4 Π
= − .
n Figura 2.7 Extensión periódica de f (x).
Por lo tanto,
∞
X sen nx
f (x) ∼ π − 2 , 0 < x < 2π.
n=1
n
Observación 2.1 Por supuesto, podrı́amos haber hecho los cálculos sobre el intervalo [−π, π] usando la fun-
ción extendida f0 . Sin embargo, en tal caso, las integrales que determinan los coeficientes de Fourier deberı́an
separarse en dos, lo que hace más engorrosa su determinación, aun cuando el resultado serı́a el mismo.
2.4. Serie de Fourier Compleja. En la sección anterior hemos visto que el conjunto de funciones
1 inx
φn (x) = √ e , n ∈ Z ,
2π
es ortonormal en el intervalo −π < x < π con la función peso r(x) = 1. Ası́, para cualquier función seccional-
mente continua f (x) definida sobre −π < x < π, (hemos abusado de la notación para estar de acuerdo con la
literatura existente)
∞
X
f (x) ∼ c̃n φn (x),
n=−∞
donde
π
1
Z
c̃n = f (t) √ e−int dt.
−π 2π
Ası́
∞ Z π
X 1 1
f (x) ∼ √ f (t)e−int dt √ einx
n=−∞ 2π −π 2π
∞ Z π
X 1
= f (t)e−int dt einx .
n=−∞
2π −π
Ejemplo 2.7 Vamos a encontrar la serie de Fourier exponencial de la función f (x) = ex definida sobre
−π < x < π. Tenemos que
π Z π
1 1
Z
t −int
cn = ee dt = e(1−in)t dt
2π −π 2π −π
1 1 h (1−in)π i
= e − e−(1−in)π
2π 1 − in
1 1 π −inπ
− e−π einπ .
= e e
2π 1 − in
Usando la fórmula de Euler y la definición de seno hiperbólico
1 1 π
cn = e cos nπ − e−π cos nπ
2π 1 − in
senh π (−1)n
= .
π 1 − in
Luego
∞
senh π X (−1)n inx
f (x) ∼ e , −π < x < π.
π n=−∞ 1 − in
En esta sección vamos a examinar la convergencia de la serie de Fourier de la función f (x). Por esto, para
hacer el análisis ampliamente aplicable, asumiremos que las funciones φn (x), n = 0, 1, . . . y f (x) son sólo
PN
seccionalemente continuas sobre [α, β]. Sea la suma de los primeros N + 1 términos de n=0 cn φn (x) que
denotaremos por SN (x). Consideraremos la diferencia |SN (x) − f (x)| para varios valores de N y x. Si para una
constante arbitraria ǫ > 0 existe un entero N (ǫ) > 0 tal que |SN (x) − f (x)| < ǫ, entonces la serie de Fourier
converge uniformemente a f (x) para todo x ∈ [α, β]. Por otro lado, si N depende tanto de x como de ǫ, entonces
la serie de Fourier converge puntualmente a f (x). Sin embargo, por el momento ambos tipos de convergencia
son demasiado exigentes, y nos conformaremos con algo menos. Para este fin, estudiaremos lo siguiente.
N
X N
X
= kφn k2 (d2n − 2dn cn + c2n ) − kφn k2 c2n + kf k2
n=0 n=0
N
X N
X
(3.3) = kφn k2 (dn − cn )2 − kφn k2 c2n + kf k2.
n=0 n=0
Teorema 3.1. Para cualquier entero no negativo N , la mejor aproximación en media a la función f (x) por
PN
una expresión de la forma n=0 dn φn (x) es obtenida cuando los coeficientes dn son los coeficientes de Fourier
de f (x).
0 ≤ kSN +1 − f k ≤ kSN − f k.
Ası́, la serie {kSN − f k, N = 0, 1, . . . } es no creciente y acotada inferiormente por cero, y por lo tanto debe
converger. Si converge a cero, entonces la serie de Fourier de f (x) converge en media a f (x). Además, de 3.4,
tenemos la desigualdad
N
X
kφn k2 c2n ≤ kf k2 .
n=0
PN
Dado que la serie {CN , N = 0, 1, . . . }, donde CN = n=0 kφn k2 c2n es no decreciente y acotada superiormente
por kf k2 , debe converger. Por lo tanto, tenemos
∞
X
(3.7) kφn k2 c2n ≤ kf k2 .
n=0
Luego, de 3.4 vemos que la serie de Fourier de f (x) converge en media a f (x) si y sólo si
∞
X
(3.8) kf k2 = kφn k2 c2n .
n=0
Para el caso en que φn (x), n = 0, 1, . . . sean ortonormales, 3.7 se reduce a la desigualdad de Bessel
∞
X
(3.9) c2n ≤ kf k2
n=0
Teorema 3.3. Sean {φn (x), n = 0, 1, . . . un conjunto ortonormal de funciones y sean cn los coeficientes de
Fourier de f (x) dados por 2.3. Entonces son válidas las siguientes
2
i La serie ∞
P
n=0 cn converge, y por lo tanto
Z β
lı́m cn = lı́m r(x)φn (x)f (x) dx = 0,
n→∞ n→∞ α
ii La desigualdad de Bessel 3.9 es verdadera,
iii La serie de Fourier de f (x) converge en media a f (x) si y sólo si la identidad de Parseval 3.10 se
cumple.
Recordando las relaciones 2.10, para la serie de Fourier trigonométrica la identidad de Parseval puede ser
escrita como 2
r ∞
√ √ π
π X Z
( πan )2 + ( πbn )2 = f 2 (x) dx,
a0 +
2 n=1 −π
o bien
∞ π
a20 X 2 1
Z
(3.11) + (an + b2n ) = f 2 (x) dx.
2 n=1
π −π
14
Ejemplo 3.1 Encontraremos la serie de Fourier trigonométrica de la función |x| sobre el intervalo −π < x < π.
Puesto que es una función par, sigue que bk = 0, ∀k ≥ 1. Además
π
2
Z
a0 = t dt = π
π 0
π π
1 π
2 2 t sen kt
Z Z
ak = t cos kt dt = − sen kt dt
π 0 π k 0 k 0
π
2 cos kt 2
= 0+ = 2 [(−1)k − 1] k≥1
kπ k 0 k π
0 si k = 2n, n = 1, 2, . . .
= 4
− si k = 2n − 1, n = 1, 2, . . .
π(2n − 1)2
Luego
∞
π X −4
|x| ∼ + cos(2n − 1)x, −π < x < π.
2 n=1 π(2n − 1)2
Evidentemente, |x| es una función seccionalmente continua sobre [−π, π]. Además, comparando su serie de
Fourier trigonométrica tenemos que
a0 π 4
= , a2n = 0, a2n−1 = − , bn = 0.
2 2 π(2n − 1)2
Y luego
∞
X 1 π4
4
= .
n=1
(2n − 1) 96
Teorema 3.4. Sean f (x) y f ′ (x) seccionalmente continuas en el intervalo [−π, π]. Entonces, la serie de Fourier
trigonométrica de f (x) converge a [f (x+ ) + f (x− )]/2 en cada punto del intervalo abierto (−π, π) y en x = ±π
la serie converge a [f (−π + ) + f (π − )]/2.
15
(
0, si x ∈ [−π, 0)
Ejemplo 3.2 Vamos a encontrar la serie de Fourier trigonométrica de la función f (x) =
1 si x ∈ [0, π].
Tenemos que
π
1
Z
a0 = 1 dt = 1.
π 0
π π
1 1
Z
ak = 1 cos kt = sen kt = 0, k ≥ 1,
π 0 kπ 0
π
1 π −1
Z
bk = 1 sen kt = cos kt = 0
π 0 kπ 0
−1 0, si k = 2n, n = 1, 2, . . .
= ((−1)k − 1) = 2
kπ , si k = 2n − 1, n = 1, 2, . . .
(2n − 1)π
Luego
∞
1 2X 1
f (x) ∼ + sen(2n − 1)x = S(x), digamos, −π < x < π.
2 π n=1 (2n − 1)
1.0
Del teorema anterior, la igualdad S(x) = f (x) es
válida en cada punto de los intervalos abiertos (−π, 0) 0.8
0.0
-3 -2 -1 0 1 2 3
Ejemplo 3.3 Vamos a obtener la serie de Fourier de f (x) = x + x2 /4, −π < x < π para mostrar que
P∞ 2 2
n=1 1/n = π /6. En primer lugar tenemos que
π
t2 1 π t2 π2
1
Z Z
a0 = t+ dt = dt = .
π −π 4 π −π 4 6
16
En segundo lugar, notando que ∀n ≥ 1, t cos nt y t2 sen nt son funciones impares sigue que
1 π
Z π
t2
Z
1
an = t+ cos nt dt = t2 cos nt dt
π −π 4 4π −π
" π #
2 π 2 π
Z π
1 t2 sen nt
1 1
Z Z
= −n t sen nt dt = 0− t sen nt dt = − t sen nt dt
4π n −π −π 4π n −π 2nπ −π
" π #
1 π
1 t cos nt 1 π cos nπ + π cos(−nπ)
Z
= − − − cos nt dt = + 0
2nπ n −π n −π 2nπ n
(−1)n
= , n ≥ 1,
n2
Z π
t2 1 π
1
Z
bn = t+ sen nt dt = t sen nt dt
π −π 4 π −π
" π #
1 π
1 −t cos nt 1 π cos nπ + π cos(−nπ)
Z
= +n cos nt dt = − −0
π n −π −π π n
2
= − (−1)n , n ≥ 1.
n
Ası́
∞
π 2 X (−1)n
f (x) ∼ + (cos nx − 2n sen nx) , −π < x < π.
12 n=1 n2
Del teorema 3.4, al evaluar en x = π, encontramos que
∞
π 2 X (−1)n f (π − ) + f (−π + )
+ 2
(cos nπ − 2n sen nπ) =
12 n=1 n 2
y luego
∞
π 2 X (−1)n n (π + π 2 /4) + (−π + π 2 /4) π2
+ ((−1) − 0) = = ,
12 n=1 n2 2 4
o ∞
X 1 π2 π2 π2
= − = .
n=1
n2 4 12 6
Para las series de Fourier de cosenos y senos también tenemos resultados similares, esta vez sobre el intervalo
[0, π].
Teorema 3.5. Sean f (x) y f ′ (x) seccionalmente continuas en el intervalo [0, π]. Entonces, la serie de Fourier
de cosenos de f (x) converge a [f (x+ ) + f (x− )]/2 en cada punto del intervalo abierto (0, π) y en x = 0 y x = π
la serie converge a f (0+ ) y f (π − ), respectivamente.
Teorema 3.6. Sean f (x) y f ′ (x) seccionalmente continuas en el intervalo [0, π]. Entonces, la serie de Fourier
de senos de f (x) converge a [f (x+ ) + f (x− )]/2 en cada punto del intervalo abierto (0, π). En x = 0 y x = π la
serie converge a 0.
Para derivación e integración de series de Fourier trigonométricas tenemos los siguientes resultados
Teorema 3.7. Suponga que f (x) es continua en (−π, π), f (−π) = f (π) y que f ′ (x) es seccionalmente continua.
Entonces,
∞
X
f ′ (x) = −nan sen nx + nbn cos nx,
n=1
en cada punto x ∈ (−π, π) donde f ′′ (x) exista.
Ejemplo 3.5 En el ejemplo 3.1 mostramos que la serie de Fourier trigonométrica de la función f (x) = |x|
está dada por
∞
π 4 X cos(2n − 1)x
f (x) ∼ − , −π < x < π.
2 π n=1 (2n − 1)2
Claramente, f (x) es continua sobre (−π, π), f (−π) = f (π) = π. Además,
−1, si −π < x < 0
′
f (x) = 1, si 0 < x < π
∄ si x = 0.
Luego, f ′ (x) es seccionalmente continua en (−π, π). Por último, f ′′ (x) = 0 si x ∈ (−π, 0) ∪ (0, π) y f ′′ (x) no
existe si x = 0. Ası́, el teorema 3.7 es aplicable y tenemos
−1, si −π < x < 0 ∞ ∞
4 X −(2n − 1) sen(2n − 1)x 4 X sen(2n − 1)x
1, si 0 < x < π = f ′ (x) ∼ − = .
π n=1 (2n − 1)2 π n=1 2n − 1
∄ si x = 0.
Note que si bien la igualdad es válida para todo x ∈ (−π, 0) ∪ (0, π), en x = 0 la derivada no existe, pero la
serie de Fourier vale 0.
Claramente, la función f (x) = x2 cumple todas las condiciones del teorema 3.8. Ası́, integrando desde −π a x,
obtenemos ∞
x
x3 π3 π2 1 (−1)n
Z X
t2 dt = + = (x + π) + 4 sen nx ,
−π 3 3 3 n=1
n n2
o
∞
x3 π2 x X (−1)n
= +4 sen nx
3 3 n=1
n3
y luego
∞
x(x2 − π 2 ) X (−1)n
p(x) = = sen nx, −π ≤ x ≤ π.
12 n=1
n3
Referencias
[1] Ordinary and Partial Differential Equations. R. P. Agarwal & D. O’Regan. Universitext, Springer Verlag, 2009.
[2] Modern Mathematical Analysis. Murray H. Protter & Charles B. Morrey. Addison–Wesley Publishing Company, 1964.
19
Ejercicios
Apéndice A. Ejercicios Resueltos
√
1. a) Pruebe que las funciones φ1 (x) = 1, φ2 (x) = 3(2x − 1) son ortonormales en el intervalo 0 < x < 1
con la función peso r(x) = 1.
√ 2
b) Encontrar la serie de Fourier de f (x) = 3x en términos de φ1 , φ2 en el intervalo 0 < x < 1.
c) Sea
Z 1 √ 2
F (a, b) = 3x2 − aφ1 (x) − bφ2 (x) dx.
0
Encontrar a, b tal que F (a, b) alcanza su mı́nimo. Además, encuentre el mı́nimo valor de F (a, b).
a) Solución. Tenemos que
√ Z 1 √ 1
hφ1 , φ2 i = 3 (2x − 1) dx = 3 [x2 − x]0 = 0,
0
Z 1
kφ1 k2 = hφ1 , φ1 i = 12 dx = 1,
0
1
(2x − 1)3
Z
kφ2 k2 = hφ2 , φ2 i = 3 (2x − 1)2 dx = = 1.
0 2
Luego, las funciones φ1 (x), φ2 (x) son ortonormales en el intervalo 0 < x < 1.
b) Los coeficientes de Fourier están dados por 2.3, luego
√
√ Z 1 2 3
c1 = hf, φ1 i = 3 x dx = ,
0 3
Z 1
1
c2 = hf, φ2 i = 3 x2 (2x − 1) dx = .
0 2
Ası́, la serie de Fourier de f (x) está dada por
√ √
3 3
f (x) ∼ + (2x − 1), 0 < x < 1.
3 2
c) Primero, notamos que F (a, b) ≥ 0. Ahora, si desarrollamos
Z 1h √ i
F (a, b) = 3x2 − 2 3x2 (aφ1 (x) + bφ2 (x)) + (aφ1 (x) + bφ2 (x))2 dx
0
Z 1 √ Z 1 √ 2 Z 1 √ Z 1
= ( 3x2 )2 dx − 2a 3x φ1 (x) dx − 2b 3x2 φ2 (x) dx + (aφ1 (x) + bφ2 (x))2 dx
0 0 0 0
= kf k2 − 2ac1 − 2bc2 + a2 + b2 ,
donde los últimos dos términos se han obtenido usando la ortonormalidad de φ1 , φ2 . Ahora derivamos
a F con respecto a a y b para obtener
∂F
= 2a − 2c1
∂a
∂F
= 2a − 2c2 .
∂b
20
Al igualar a cero ambas ecuaciones obtenemos los valores a = c1 , b = c2 . Puesto que F (a, b) ≥ 0
sigue que estos valores minimizan F . Además, este valor mı́nimo es
3 1 1 1
F (c1 , c2 ) = kf k2 − c21 − c22 =
− − = .
5 3 4 60
Note que este valor no es cero. Esto es natural, es imposible obtener una aproximación exacta al
aproximar una función cuadrática usando polinomios(de grado cero y uno.
0, si −π < x < 0,
2. Encontrar la serie de Fourier trigonométrica de f (x) = Utilice este resultado
1, si 0 < x < π.
P∞ 1
para calcular la suma n=1 (2n−1) 2.
n2 −a2 .
Solución. Dejamos propuesto como ejercicio el caso a ∈ Z. Para los otros casos, dado que cos ax es una
función par, sigue que bn = 0, n ≥ 1 y
2 π 1
Z
an = cos ax cos nx dx = [cos(n − a)x + cos(n + a)x] dx
π 0 π
π
1 sen(n − a)x sen(n + a)x
= +
π n−a n+a
0
1 sen aπ sen aπ 1 2a
= (−1)n+1 − = (−1)n+1 sen aπ 2 , n ≥ 0.
π n−a n+a π n − a2
Luego, sigue que
" ∞
#
2a 1 X (−1)n+1
cos ax ∼ sen aπ + cos nx , −π < x < π.
π 2a2 n=1 n2 − a2
21
Al evaluar en el punto de continuidad x = 0 y usando el teorema 3.4 tenemos
" ∞
#
2a 1 X (−1)n+1
sen aπ + = cos 0 = 1
π 2a2 n=1 n2 − a2
y luego
∞
X (−1)n+1 1 h aπ i
2 2
= 2 −1 .
n=1
n −a 2a sen aπ
(
0 si −π < x < 0
4. Considere a la función de periodo 2π definida por f (x) =
sen x si 0<x<π
a) Pruebe que la serie de Fourier S(x) de f (x) es
∞
1 sen x 2 X cos 2nx
S(x) = + − ,
π 2 π n=1 4n2 − 1
1. Encuentre la serie de Fourier de cosenos sobre el intervalo 0 < x < π para cada una de las siguientes
funciones:
(
0, si 0 < x < π/2
(i) f (x) = x2 (ii) f (x) = (iii) f (x) = cos x.
1, si π/2 < x < π
23
2. Encuentre la serie de Fourier de senos sobre el intervalo 0 < x < π para cada una de las siguientes
funciones:
Además, para cada serie grafique la suma mostrando algunos términos gráficamente, comparando con la
función expandida.
4. Encuentre la serie de Fourier trigonométrica de las siguientes funciones:
0, si −2 < x < −1
(i) f (x) = x2 , 0 < x < 2π (ii) f (x) = 1, si −1 < x < 1
0, si
1 < x < 2.
5. Sea f (x) una función periódica de perı́odo 2π tal que sus coeficientes de Fourier existan. Suponga que
f (π − x) = f (x), ∀x. Muestre que an = 0 si n es impar y que bn = 0 si n es par. Encuentre explicitamente
c0 , c2 y c4 .
6. Encontrar, si existe, una serie de Fourier de perı́odo fundamental 2π para la función f (x) = sen3 x, ∀x ∈ R.
Justifique brevemente su respuesta.
7. Sea f una función arbitraria, seccionalmente continua en [0, π], simétrica con respecto a la recta x = π/2.
Demostrar que los únicos términos no nulos en el desarrollo de Fourier senoidal de f , con perı́do 2π, son
los términos b2k−1 , k = 1, 2, . . .
8. Desarrolle la función:
f (x) = x − 1, 1≤x≤2
en serie de Fourier, de tal modo que esa expansión tenga sólo términos pares.
9. Demostrar o refutar la identidad
∞
π 2 X cos 2nx
x(π − x) = − , 0 ≤ x ≤ π.
6 n=1
n2
(
1, si 0 < x < π/2 2
q
10. Considere la función f (x) = Es sabido que las funciones φn (x) =
π sen nx, n =
0, si π/2 < x < π.
1, 2, 3 forman un conjunto ortonormal en el intervalo 0 < x < π. Encuentre los valores de a, b para los
cuales la integral
Z π
F (a, b) = (f (x) − aφ1 (x) − b(φ2 (x) + φ3 (x)))2 dx
0
alcanza su mı́nimo.
24
(
0, si x ∈ [−1, 0)
11. Sea f (x) = Muestre que
1, si x ∈ [0, 1].
Z 1 2 Z 1
1 3 2
f (x) − − x dx ≤ f (x) − d0 − d1 x − d2 x2 dx
−1 2 4 −1
25