BMAI U2 Contenido
BMAI U2 Contenido
BMAI U2 Contenido
Series de Fourier
Programa de la asignatura:
U2 Series de Fourier
Índice
Presentación de la Unidad
En esta Unidad se presenta una introducción a las series de Fourier, la primera sección
comienza con el concepto de producto interno de funciones, a partir de éste se define la
ortogonalidad de funciones e inmediatamente se presenta la definición de función
periódica.
La segunda parte presenta cómo se obtiene la serie de Fourier de una función periódica
junto con las propiedades que posee la misma. Finalmente, se presentan los conceptos
básicos de la aproximación para las series de Fourier.
Propósitos
Competencia específica
2. Series de Fourier
Solución: Para resolver este ejercicio solo basta aplicar la definición de producto
interno del siguiente modo:
b 2 2 2
1
∫ ∫ ( x )( x= ) dx ∫−=
3 2
f,g
= f ( x) g =
( x)dx x 5 dx x 6
a −1 1 6 −1
1 1 32 1 21
= ( 2 ) − ( −1) =
6 6
− =
6 6 3 6 2
21
Por lo tanto f,g = .
2
En el siguiente paso se aplica la técnica de integración por partes, lo que permite obtener:
1 1
∫ x sen(2 x)dx =
− x cos(2 x) + sen(2 x)
2 4
2 2 2
x3 x 4
f , g = ∫ ( x ) ( 2 x − 1) dx = ∫ ( 2 x − x ) dx = − +
2 3 2
−1 −1 3 2 −1
( 2 )3 ( 2 )4 ( −1)3 ( −1)4
= − + − − +
3 2 3 2
9
=
2
9
Por lo tanto, se tiene que f,g =
2
( f + g ) ( x) =f ( x) + g ( x) .
b b
=f,g ∫=
a
f ( x) g ( x)dx ∫=
a
g ( x) f ( x)dx g, f
= c f , h + g, h
A partir de la relación anterior se puede definir una norma sobre C[ a, b ] de forma similar
al caso en 3 :
∫ [ f ( x)]
2
=f f, f
= dx .
a
b 2
x) ] dx
∫ [ f (= ∫ [3sen(π x)]
2 2
=f f, f
= dx
a o
2 2
1 1 x 1
= 9 ∫ − cos(2π x)=
dx 3 − sen(2π x)
o
2 2 2 4 0
(2) 1 (0) 1
= 3 − sen ( 2π (2) ) − − sen ( 2π (0) )
2 4 2 4
= 3 1− 0 =3
1 1
Por lo tanto f = 3 y su forma normalizada es f ( x) = sen(π x) .
f 3
Utilizando la identidad:
1
sen A cos
= B sen ( A − B ) + sen ( A + B )
2
En consecuencia:
1 1
1
∫ cos(
−1
= π x) sen(2π x)dx ∫ 2 ( sen(π x + 2π x) + sen(2π x − π x) )dx
−1
1 1
1 1 1
∫−1 2 ( sen(3π x) + sen(π x) )dx =
= − 6π cos(3π x) − 2π cos(π x)
−1
1 1 1 1
= − cos(3π ) − cos(π ) − − cos(−3π ) − cos(−π )
6π 2π 6π 2π
=0
2 2 2
1
∫ ( x ) ( x= ) dx ∫−=
2
f,g
= x 3 dx x 4
−2 2 4 −2
1 1
( 2 ) − ( −2 ) = 4 − 4 = 0
4 4
=
4 4
Definición: Sea {φn }n∈I un subconjunto de C[ a, b ] , se dice que {φn }n∈I es un conjunto
Solución: Por comodidad, sea φn ( x) = cos(nπ x) , con n ∈ , hay que mostrar que
φn es ortogonal a φm para m ≠ n , esto se tiene aplicando la definición de producto interno
de funciones del siguiente modo:
1 1
=φn , φm ∫=
φ ( x)φ ( x)dx ∫ cos(nπ x) cos(mπ x)dx
0
n m
0
Se tiene lo siguiente:
1 1
1
∫ cos(nπ x) cos(mπ=
0
x)dx ∫ 2 cos ( nπ x − mπ x ) + cos ( nπ x − mπ x ) dx
0
1
1
= ∫ 2 cos π ( n − m ) x + cos π ( n + m ) x dx
0
1
1 1 1
= sen π ( n − m ) x + sen π ( n + m ) x
2π n − m n+m 0
1 1 1
= sen π ( n − m ) + sen π ( n + m ) − 0
2π n − m n+m
=0
Cuando m = n ≠ 0 , es decir:
1 1
∫ [φ ( x=
) ] dx ∫ cos
2 2
φn
= φ=
n , φn n
2
(nπ x)dx
0 0
Utilizando la relación:
1 1
cos 2 A= + cos ( 2 A )
2 2
Se tiene lo siguiente:
1 1 1
1 1 x 1
∫0 2 + 2 cos ( 2nπ x ) =
∫0 cos (nπ x)dx = 2 + 4π n sen ( 2nπ x )
2
1 1 1
= + sen ( 2nπ ) − 0 =
2 4π n 2
1
Por consiguiente φn = . Finalmente, normalizando se tiene lo siguiente:
2
1 1
= φn ( x) = cos(nπ x) 2 cos(nπ x)
φn 1
2
{
Por lo tanto, el conjunto 1, 2 cos(nπ x)} n∈ \{0}
es ortonormal sobre [ 0,1] .
Considera que tienes una hoja de papel, en ella traza un sistema de ejes cartesianos, uno
horizontal y otro vertical.
Enrolla dicha hoja de tal forma en que el eje horizontal coincida consigo mismo, formando
un cilindro circular.
Luego, sobre la superficie del cilindro, observando la posición del eje vertical marca el
punto que está del otro lado del círculo que forma el eje horizontal con respecto al origen
de coordenadas.
Dibuja una curva que se marque sobre todas las capas del cilindro:
Lo que se obtiene en una figura que se repite a partir de un intervalo igual a la longitud del
círculo base del cilindro que formaste, una función periódica tiene este comportamiento, la
diferencia radica en el hecho de que en vez de considerar un hoja, hay que considerar un
plano infinito. Ahora se presenta la definición formal del proceso anterior.
Esta definición afirma que basta conocer cómo se comporta la función en un intervalo de
longitud T para tener el comportamiento de toda la función, por simplicidad se toma dicho
T T
intervalo alrededor del 0 , es decir, el intervalo base o fundamental es −, , aquí se
2 2
considerará el intervalo base indistintamente abierto o cerrado. De forma gráfica se tiene
lo siguiente:
Por lo tanto las funciones sen( x) y cos x son periódicas y de periodo 2π . Gráficamente
se tiene lo siguiente:
=
función h( x) sen ( 3x ) + sen ( 4 x ) es periódica y encontrar dicho periodo.
=
Implica que h( x) sen ( 3x ) + sen ( 4 x ) es periódica. Para encontrar el periodo primero hay
que observar que el denominador común 6 por consiguiente:
2 4 1 3
=T1 = π = π y T2 π= π
3 6 2 6
Ahora, hay que calcular el mínimo común múltiplo entre los numeradores, en este caso es
12 . Por consiguiente:
4 3
= π y 4T2 4 π
3T1 3=
6 6
12
=
Por lo tanto, h( x) sen ( 3x ) + sen ( 4 x ) es periódica y de periodo T = π = 2π .
6
x x
Ejercicio: Dadas las función f ( x) = sen y g ( x) = sen , muestre que la
3 2
x x
=
función h( x) sen + sen es periódica y encontrar dicho periodo.
3 2
x
respectivamente. Gráficamente, para la función f ( x) = sen :
3
x
Para la función g ( x) = sen :
2
=
Esto implica que h( x) sen ( 5 x ) + sen ( 2 x ) es periódica. Para encontrar el periodo, hay
que observar que los periodos son múltiplos enteros de π , en consecuencia, solo hay
que calcular el mínimo común múltiplo entre los coeficientes, en este caso es 12 . Por lo
=
tanto, h( x) sen ( 3x ) + sen ( 4 x ) es periódica y de periodo T = 12π . Gráficamente, se
tiene lo siguiente:
A continuación se presentan algunas relaciones que hay entre las funciones periódicas y
sus integrales: Sean α , β ∈ , con α < β , entonces hay que analizar el comportamiento
de la integral.
β +T
∫
α +T
f ( x)dx
f ( y ) = f ( y − T + T ) = f ( x + T ) = f ( x) :
β β +T β +T
∫
α
f ( x)dx= ∫
α +T
f ( y − T )dy=
α +T
∫ f ( y )dy
Lema: Sea f una función periódica de periodo T , dados α , β ∈ , con α < β , entonces
β +T β
∫
α +T
f ( x)dx = ∫ f ( x)dx
α
∫
T
f ( x)dx = ∫ f ( x)dx
0
T
∫ f ( x)dx
= ∫ T
f ( x)dx + ∫T
f ( x)dx
a− a− −
2 2 2
∫ f ( x)dx
=
T
∫ f ( x)dx
=
T T
∫ f ( x)dx
a− a + −T a+
2 2 2
∫ T
f ( x)dx = ∫ T
f ( x)dx + ∫ T
f ( x)dx = ∫ T
f ( x)dx + ∫ T
f ( x)dx
a− a− − a+ −
2 2 2 2 2
T T T
a+
2 2 2
= ∫T
f ( x)dx + ∫ T
f ( x)dx = ∫ T
f ( x)dx
− a+ −
2 2 2
Corolario: Sea f una función periódica de periodo T entonces para todo a∈ se tiene
que:
T T
a+
2 2
∫ T
f ( x)dx = ∫ T
f ( x)dx
a− −
2 2
∫
0
f ( x)dx = ∫
T
f ( x)dx .
T
Aplicando el corolario anterior al caso particular cuando a = se tiene que:
2
T
T 2
∫ f ( x)dx = ∫
0 T
f ( x)dx .
−
2
Como es sabido, las funciones senos y cosenos son funciones continuas y periódicas, la
esencia detrás de la serie de Fourier es utilizar estas funciones para poder aproximarse a
cualquier función periódicas.
Para poder expresar una función periódica en términos de seno y cosenos, primero hay
que considerar el siguiente resultado, que muestra que este conjunto es ortogonal el
intervalo de longitud 2π centrado en el origen:
Proposición: El conjunto
{1,sen( x),cos( x),sen(2 x),cos(2 x),sen(3x),cos(3x),…}
Es ortogonal en ( −π , π ) .
Es decir, φn es ortogonal a 1 .
π
1
ϕn ,1
= ∫π cos(
−
= nx)dx
n
=sen(nx) 0
Es decir, φn es ortogonal a 1 .
En resumen se tiene:
0 si m ≠ 0
φn ,φm = .
π si m = 0
1
2 cos ( m − n ) x + cos ( m + n ) x si m ≠ n
cos ( mx ) cos ( nx ) =
1
2
(1 + cos ( 2nx ) ) si m = n
Es ortogonal en ( −π , π ) .
Ahora se tienen todos los ingredientes para poder presentar la serie de Fourier de una
función periódica, lo que se realiza a continuación: Sea f una función periódica de
periodo 2π y supóngase que tal función puede escribirse de la siguiente manera:
a0 ∞
f ( x) = + ∑ an cos ( nx ) + bn sen ( nx )
2 n =1
Partiendo de una compatibilidad entre la integral y la serie, se puede integrar tal función
de −π a π para obtener lo siguiente:
π π
a0 ∞
∫
−π
∫−π 2 ∑
f ( x)dx =+
n =1
an cos ( nx ) + bn sen ( nx ) dx
a0 π ∞ π
2 −∫π
=+ d x ∑ ∫ ( an cos ( nx ) + bn sen ( nx ) ) dx
n =1 −π
= a0π
π
1
Es decir, se tiene que a0 =
π ∫
−π
f ( x)dx . Por otra parte para k ∈ se multiplica la función
π π
a0 ∞
∫
−π
∫−π 2 ∑
f ( x) cos(kx)dx =+
n =1
an cos ( nx ) + bn sen ( nx ) cos(kx)dx
π
a0 ∞
∫−π 2 sen(kx) + ∑
=
n =1
an cos ( nx ) cos(kx) + bn sen ( nx ) cos(kx) dx
a0 π ∞ π
2 −∫π
= cos( kx ) dx + ∑ ∫ ( an cos ( nx ) cos(kx) + bn sen ( nx ) cos(kx) ) dx
n =1 −π
π
=
−
∫π a k cos 2 (kx)dx = ak π
π
1
Por consiguiente ak =
π ∫
−π
f ( x)cos(kx)dx . Si el proceso anterior se realiza con sen(kx) se
tiene lo siguiente:
π π
a0 ∞
∫
−π
∫−π 2 ∑
f ( x)sen(kx)dx =+
n =1
an cos ( nx ) + bn sen ( nx ) sen(kx)dx
π
a0 ∞
∫−π 2
= sen( kx ) + ∑
n =1
an cos ( nx ) sen(kx) + bn sen ( nx ) sen(kx) dx
a0 π ∞ π
sen(kx)dx + ∑ ∫ ( an cos ( nx ) sen(kx) + bn sen ( nx ) sen(kx) ) dx
2 −∫π
=
n =1 −π
π
=
−
∫π b k sen 2 (kx)dx = bk π
π
1
Por consiguiente bk =
π ∫
−π
f ( x)sen(kx)dx .
Cabe mencionar que cuando no hay peligro de confusión usualmente se dice serie de
Fourier en vez de serie trigonométrica de Fourier.
Solución: Para tener una idea visual observa que f en el intervalo fundamental
tiene la siguiente gráfica:
π π
1
0
1
an = ∫ f ( x ) cos( nx ) dx = ∫ ( − 2 ) cos( nx ) dx + ∫ x cos(nx)dx
π −π π −π 0
π
1 2 cos(nx) x sen(nx)
0
= − sen(nx) + +
π n −π n
2
n 0
1 cos(nπ ) 1 1 ( −1) 1
n
= 0 + − 2
= − 2
π n 2
n π n 2
n
2
Finalmente se tiene que a2 k −1 = − .
π ( 2k − 1)
2
1 2 2 n
π ( −1)n
= − ( −1) + −
π n n n
2 − 2 ( −1) − π ( −1)
n n
=
nπ
2 − 2 ( −1) − π ( −1)
n n
Solución: Para tener una idea visual observa que f en el intervalo base tiene la
siguiente gráfica:
π π
1
0
1
an = ∫ f ( x ) cos( nx ) dx = ∫ ( − x ) cos( nx ) dx + ∫ x cos(nx)dx
π −π π −π 0
π
1 cos(nx) x sen(nx) cos(nx) x sen(nx)
0
= − − + +
π n 2
n −π n
2
n 0
cos(nπ ) 1 2 ( −1) 1
n
1 1 cos(nπ )
0
= 2− + + − 2
= − 2
π n n −π n 2
n π n 2
n
4
Finalmente se tiene que a2 k −1 = − .
π ( 2k − 1)
2
= − + + − +
π n n 2
−π n n 2 0
1 π ( −1) π ( −1)
n n
= = +− 0
π n n
π ∞ 2
= −∑ cos ( (2n − 1) x )
2 n =1 π ( 2n − 1)2
Donde
π π π
1 1 1
=a0 = ∫ g ( x)dx an
π −π
π ∫ g ( x)=
−π
cos(nx)dx b n
π ∫ g ( x)sen(nx)dx
−π
2π
Sustituyendo x x se tiene que:
T
2π a0 ∞ 2π n 2π n
f ( x) =
g x = + ∑ an cos x + bn sen x
T 2 n =1 T T
Además
T T
π π
1 1 T 1 2 2π 2 2
=a0 = ∫ g ( x)dx
π −π
π ∫
−π
f=
x
2π
dx
π ∫T
f
= (u )
T
du
T ∫T
f (u )du
− −
2 2
Luego
π π
1 1 T
=an =
π ∫ g ( x) cos(nx)dx
−π
π ∫π f 2π x cos ( nx ) dx
−
T
1 2π 2π
2
=
π ∫ f ( u ) cos T nu T du
T
−
2
T
2 2π n
2
=
T ∫ f ( u ) cos T u du
T
−
2
2 2π n
2
bn =
T ∫ f ( u ) sen T u du
T
−
2
Esto permite obtener la serie de Fourier de una función periódica de periodo arbitrario
como se presenta a continuación:
2 2 2 2π n 2 2 2π n
2
=a0 = ∫ f ( x)dx an ∫ f ( x ) cos = x dx bn ∫ f ( x ) sen x dx .
T T T T T T T T
− − −
2 2 2
1 4cos ( nπ ) 4cos ( − nπ )
= = − 0
2 n 2π 2 n 2π 2
1 2 ( 2 ) cos ( nπ ) 2 ( −2 ) x cos ( nπ )
= − +
2 nπ nπ
4 ( −1)
n +1
=
nπ
4 ( −1)
n +1
3 ( −1) − 1
n
∫ x sen (π nx ) dx + ∫ x sen (π nx ) dx
2
=
−| 0
=
2
−
( 2 − n 2π 2 ) cos(− nπ ) cos(nπ )
+ −
− 0
n π nπ nπ
3 3 3 3
=
2
−
( 2 − n π ) ( −1) − ( −1)
2 2 n
2 − ( 2 − n 2π 2 ) ( −1) − n 2π 2 ( −1)
=
n n n
=
n3π 3
4 n es impar
En consecuencia
0 n es par
bn = 4
π 3 n3 n es impar
4
Finalmente se tiene que b2 k −1 = .
π ( 2k − 1)
3 3
El proceso anterior, se puede generalizar del siguiente modo: Sea f una función
continua a pedazos definida sobre el intervalo ( a, b ) . Dado {φn }n∈ un conjunto ortonormal
Dado k ∈ con k fijo, suponiendo que la integral y la serie pueden conmutar, entonces
se tiene las siguientes relaciones:
b
∞ ∞
=f ,φk ∫ ∑ anφn ( x),φk ( x)
∑ anφn ( x) φk ( x)dx
=
=n 1=a n 1
b ∞ ∞ b
= n n k ∑ a φ ( x)φ ( x)dx ∑ ∫ a φ ( x)φ ( x)dx
∫= n n k
=a n 1= n 1 a
∞ b ∞
=
=n 1= n 1
∑=
an ∫ φn ( x)φk ( x)dx ∑a n φn ,φk
a
= ak
Por lo tanto se tiene que:
∞ b
f ( x) = ∑ anφn ( x) donde
= an f ,φn
= ∫ f ( x)φ ( x)dx .
n
n =1 a
Esta última relación se conoce como la serie de Fourier generalizada con respecto al
conjunto ortonormal {φn }n∈ de la función f en el intervalo ( a, b ) , finalmente los
=
n n
n 0=
0 f ( x=
n 1
) ∑ a φ ( x=) a + ∑ an 2 cos ( nπ x )
1 1 1
Donde a0 = ∫ =
f ( x)dx y an ∫=
f ( x)φn ( x)dx 2 ∫ f ( x) cos(nπ x)dx .
0 0 0
Como observarás a medida que se k crece la discrepancia que hay entre las gráficas de
las funciones f y f k es menor, este proceso toma el nombre de convergencia de
funciones.
De forma gráfica, una función a suave a pedazos debe de tener un número finito de
discontinuidades y un número finito de picos.
Cabe mencionar que estas condiciones son de suficiencia, esto quiere decir, que si se
cumplen se garantiza la convergencia de la serie de Fourier, sin embargo cuando una
función no se cumpla con dichas condiciones su serie de Fourier puede o no converger.
Aplicando este teorema al ejemplo presentando anteriormente, debes observa que hay
una discontinuidad en x0 = 0 . Luego f ( x) → −2 cuando x → 0− y además f ( x) → 0
cuando x → 0+ , entonces la serie de Fourier converge en x0 = 0 a
1 1
2
( )
f ( x0+ ) + f ( x0− ) = ( 0 + ( −2 ) ) =−1
2
Teorema: Sea f una función suave a pedazos y periódica de periodo T tal que f ' es
suave a pedazos, entonces
∞
2π n 2π 2π
∑
f '( x) =
n =1 T
−an sen
T
nx + bn cos
T
nx .
Donde
T T T
2 2 2π n 2 2π n
2 2 2
=α0
T ∫T
= f '( x)dx α n
T ∫ f ' ( x ) cos=
T
T
x dx β n
T ∫ f ' ( x ) sen T x dx
T
− − −
2 2 2
2 2
2 T T T T
α0 =
T ∫ T
f '( x)dx = f −
T 2
f − = 0 ya que f = f − .
2 2 2
−
2
Por otra parte, por medio de la integración por partes se tienen las siguientes relaciones:
T
2 2π n
2
αn =
T ∫ f ' ( x ) cos T x dx
T
−
2
T T
2 2 π n 2 2π n 2
2π n
f ( x ) sen
T ∫T
= f ( x) cos x + x dx
T T − T T
2 −
2
T2
2π n 2 2π n 2π n
f ( x ) sen
T T ∫T
= = x dx bn
T T
− 2
2 2π n
2
βn =
T ∫ f ' ( x ) sen T x dx
T
−
2
T T
2 2π n 2 2π n 2 2π n
f ( x ) cos
T ∫T
= f ( x)sen x − x dx
T T − T T
2 −
2
T
2π n 2 2 2π n 2π n
f ( x ) cos
T T ∫T
= − x dx =− an
T T
−
2
Solución: Antes que nada hay que observar que la función f tiene periodo T = 6 ,
Los coeficientes de la serie de Fourier, se calculan de forma similar a los ejemplos de la
sección anterior:
=
(
3 ( −1) e3 − ( −1) e −3
n n
) = 3( −1) ( e
n 3
− e −3 )
9 + n 2π 2 9 + n 2π 2
3 ( −1) ( e3 − e −3 )
n
=
( −nπ ( −1) e n 3
− nπ ( −1) e −3
n
) = − nπ ( −1) ( en 3
+ e −3 )
9 + n 2π 2 9 + n 2π 2
nπ ( −1) ( e3 + e −3 )
n
+
3 ( −1)
( 2)
(e 3
− e −3 ) π ( 2 ) π ( −1) ( e + e )
cos ( 2 ) x −
3 −3( 2)
π
sen ( 2 ) x
9 + ( 2) π 2 9 + ( 2) π
2 2
3 2
3
Utilice los primeros 7 términos diferentes de cero de la serie de Fourier para aproximar a
la función f .
= ( 0 + ( −1) ) + ( (1) − 0 ) =
[ − x ]−1 + [ x ]0 =
0 1
0
=
πn
4 n es impar
En consecuencia:
0 n es par
bn = 4
π n n es impar
4
Finalmente se tiene que b2 k −1 = .
π ( 2k − 1)
Como habrás visto, hay funciones periódicas que tienen serie de Fourier con
exclusivamente cosenos o exclusivamente senos, en esta sección se presenta las
condiciones para que esto resulte. Para esto se comienza con la siguiente definición:
Definición: Sea f una función definida en algún intervalo I centrado en cero, entonces
se dice que f es par si y solos si f (− x) =f ( x) para todo x ∈ I . Por otro lado se dice que
f es impar si y solo si f (− x) =− f ( x) para todo x ∈ I .
Gráficamente, una función par es simétrica con respecto al eje vertical, como lo muestra
la siguiente figura:
De forma similar, una función impar es simétrica con respecto al origen de coordenadas,
como lo muestra la siguiente figura:
× P I
P P I
I P I
Ejemplo: La tabla anterior a firma que la función cos(mx) cos(nx) es par, que la función
cos(mx)sen(nx) es impar.
Ahora se presenta las propiedades que tienen las funciones pares e impares con respecto
a la integral sobre un intervalo centrado en el origen.
∫
−a
f ( x)dx = 2 ∫ f ( x)dx
0
∫=
−a
f ( x)dx ∫ f ( x)dx + ∫ f ( x)dx
−a 0
∫
−a
∫ f (−u ) ( −du ) =
f ( x)dx =
a
− ∫ f (u )du =
a
∫ f (u )du
0
Por consiguiente:
0 a a a a
Por lo tanto:
a a
∫
−a
f ( x)dx = 2 ∫ f ( x)dx .
0
Este resultado se puede obtener de forma gráfica, observando la simetría con respecto al
eje vertical, como lo muestra la siguiente figura:
−a
∫ f ( x)dx = 0
∫=
−a
f ( x)dx ∫
−a
f ( x)dx + ∫ f ( x)dx
0
∫
−a
f ( x)dx =∫ f (−u ) ( −du ) =−∫ ( − f (u ) ) du =−∫ f (u )du
a a 0
Por consiguiente:
0 a a a
∫
−a
f ( x)dx + ∫ f ( x)dx =
0
− ∫ f ( x)dx + ∫ f ( x)dx =
0 0
0
Por lo tanto
a
∫ f ( x)dx = 0 .
−a
Ahora se aplican las dos proposiciones anteriores para el cálculo de los coeficientes de
funciones pares e impares.
2π
Demostración: supóngase que f es par, dado que la función f ( x)sen nx es impar
T
se tiene que:
T
2 2 2π
=bn = ∫
T T
f ( x)sen
T
nx dx 0
−
2
2 2
T ∫T
=a0 = f ( x)dx 0
−
2
2π
Además la función f ( x)sen nx es impar, lo que implica que:
T
T
2 2π
2
=an =∫
T T
f ( x) cos
T
nx dx 0 .
−
2
1 π
2 2
2 2 2π 2 2π
=bn =∫
T T
f ( x )sen
T
nx
dx = ∫
4 −2
f ( x )sen
4
nx
dx 2
2 0 ∫ f ( x)sen nx dx
2
−
2
2 2
π 4 2 x nπ x 4 nπ x
= ∫ ( − x + 2 ) sen nx dx = − + cos − 2 2 sen
0 2 nπ nπ 2 nπ 2 0
4 2 ( 2 ) 4
= − + cos ( nπ ) − −
nπ nπ nπ
4
=
nπ
4
Es decir bn = . Por consiguiente la serie de Fourier de f es:
nπ
a0 ∞ π π
f ( x) = + ∑ an cos nx + bn sen nx
2 n =1 2 2
∞
4 π
=∑ sen nx
n =1 nπ 2
Además, f es par ya que es simétrica con respecto al eje vertical, por consiguiente
bn = 0 para todo n ≥ 1 . Así, solo resta calcular los valores de an para n ≥ 0 , tomando
1
2 2 2
2 2
2
a0 =
T ∫ T
f ( x)dx = ∫
4 −2
f ( x ) dx = 2
2 0 ∫ f ( x)dx =
∫ ( − x + 2 ) dx
0
−
2
2
1 2 1 2
− 2 x + 2 x =
= − ( 2) + 2 ( 2) =
2
0 2
2π
Luego, la función f ( x)cos nx es par, lo cual implica lo siguiente:
T
T
1 π
2 2
2 2 2π 2 2π
=an =∫
T T
f ( x )sen
T
nx
dx = ∫
4 −2
f ( x ) cos
4
nx
dx ∫ f ( x) cos nx dx
2
2 0 2
−
2
2 2
π 4 nπ x 4 2 x nπ x
∫0 ( − x + 2 ) cos 2 nx dx = − n2π 2 cos 2 + nπ − nπ sen 2 0
=
4
= − 2 2 cos ( nπ ) − − 2 =
nπ
4
nπ nπ
2
4
2 2
− ( −1) + 1
n
( )
La expresión:
0 n es par
1 − ( −1) =
n
2 n es impar
En consecuencia
0 n es par
an = 8
π 2 n 2 n es impar
8
Finalmente se tiene que a2 k −1 = . Por consiguiente la serie de Fourier de f es:
( 2k − 1)
2
π 2
a0 ∞ π π
f ( x) = + ∑ an cos nx + bn sen nx
2 n =1 2 2
∞
8 π
1+ ∑ 2
= cos ( 2n − 1) x
n =1 π ( 2n − 1)
2
2
Para comenzar esta sección, se presenta la definición de distancia entre dos funciones, la
cual es análoga a la distancia entre dos vectores tridimensionales. Para vectores en 3 la
distancia entre punto α y el punto β se define por:
d (α , β ) = α − β
∫ ( f ( x) − g ( x) )
2
d ( f , g) = f − g = dx .
a
Como se vio en la sección anterior, una función periódica pude ser aproximada por un
número finito de elementos de su serie de Fourier, en esta sección se presenta una forma
de calcular dicha aproximación.
Como se vio en la sección anterior, para una función continua a pedazos f de periodo T
su serie de Fourier es:
a0 ∞ 2π 2π
f ( x) = + ∑ an cos nx + bn sen nx
2 n =1 T T
Donde
T T
π
2 2 2π 2 2π
2 2
=a0 = ∫ f ( x)dx an ∫ f ( x ) cos = nx dx bn ∫ f ( x ) sen T nx dx
T −π T T T T T
− −
2 2
a0 N 2π 2π
Φ N ( x) = + ∑ an cos nx + bn sen nx
2 n =1 T T
2 1 N 2π 2π 2π 2π
2
N ( x)
Φ= ∫ f (t ) + ∑ cos nt cos nx + sen nt sen nx dt
T
−
T 2 n =1 T T T T
2
T
2 2 1 N 2π
= ∫ f (t ) + ∑ cos n ( t − x ) dt
T T
−
2 n =1 T
2
Considerando la relación:
A− B A+ B
sen ( A ) − sen ( B ) =
2sen cos
2 2
Es decir
2π 1 2π 1
sen n + ( t − x ) − sen n − ( t − x )
2π T 2 T 2
cos n ( t − x ) =
T t − x
2sen
2
Esta es una manera de calcular la suma parcial Φ n de una serie de Fourier a través de
una integral, esta integra toma el nombre de Integral de Dirichlet. En la práctica no es
fácil calcular dicha integral, pero existen métodos como la regla de los trapecios o la regla
de Simpson que permiten aproximar el valor de dicha integral adecuadamente.
Ahora hay que calcular la serie de Fourier de f , para esto observa que f es impar lo
que implica que an = 0 para toda n ∈ . Para los coeficientes bn se tiene lo siguiente:
1 π
2 2
2 2π
=bn ∫ f ( x )sen nx
= dx ∫ ( −2 ) sen nx dx
2
4 −2 4 2 0 2
2
4 π 4 4
= cos = nx cos (π n ) −
π n 2 0 π n πn
=
4
πn
( −1) − 1
n
( )
Luego se tiene que:
a0 8 2π 2π
Φ 8 ( x) = + ∑ an cos nx + bn sen nx
2 n =1 T T
4
( ) π
8
= ∑ ( −1) − 1 sen nx
n
n =1 π n 2
8 πx 8 3π x 8 5π x 8 7π x
= − sen − sen − sen − sen
π 2 3π 2 5π 2 7π 2
Tienes que observar que también se toman en cuenta los coeficientes iguales a cero.
Finalmente, se tiene que
1
2 sen 4π + ( t − x )
1 2 dx = − 8 sen π x − 8 sen 3π x
∫
4 −2
f (t )
t−x π
2 3π
2
sen
2
8 5π x 8 7π x
− sen − sen
5π 2 7π 2
Por consiguiente:
1
2 sen 4π + ( t − x )
2 dx = − 32 sen π x + 1 sen 3π x + 1 sen 5π x + 1 sen 7π x
∫ f (t )
t−x π 2 3 2 5 2 7
2
−2 sen
2
1
1 sen 3π + ( t − x )
2 dx .
∫−1 f (t ) t−x
sen
2
Es similar a la integral:
2π 1
T
sen N + ( t − x )
2
T 2 dt
∫ f (t )
t−x
−
T
sen
2
2
Toca el turno de calcular la serie de Fourier de f , para esto observa que bn = 0 para toda
n ∈ ya que f es par. Para el coeficiente a0 se tiene lo siguiente:
1 1 1
2 1 2 1
a0
= ∫
2 −1
x)dx 2 ∫=
f (=
0
xdx 2 =x 2=
2 0
2 1
Al igual que en el ejemplo anterior, también se toman en cuenta los coeficientes iguales a
cero. Finalmente, se tiene que:
1
1 sen 3π + ( t − x )
2 dx = 1 − 8 cos π x − 8 cos 3π v .
∫−1 f (t ) ( ) ( )
t−x π2 9π 2
sen
2
Dada una función f continua a pedazos definida sobre el intervalo base ( a, b ) , dado un
conjunto ortogonal de funciones continuas a pedazos {φn }n∈ definidas sobre ( a, b ) , con
Donde a0 ,…, an ∈ .
El siguiente número:
b
∫ [ f ( x) − Φ ( x) ] dx
2 2
f − ΦN = n
a
Observa que:
[ f ( x) − Φ N ( x)] =[ f ( x)] − 2 f ( x )Φ N ( x ) + [ Φ N ( x ) ] .
2 2 2
∫ [ f ( x) − Φ N ( x)] dx
2 2
EN =f − Φ N =
a
b
∫ ([ f ( x)] − 2 f ( x)Φ N ( x) + [ Φ N ( x) ] dx )
2 2
=
a
b b b
[Φ N ( x)] dx ∫ ( γ φ ( x) + γ φ ( x) + + γ φN ( x) ) dx
2 2
ΦN ,ΦN
= ∫=
a a
0 0 1 1 N
b
∫a (γ 0φ0 ( x) ) + + (γ N φN ( x) ) + 2m∑≠ n γ mγ nφm ( x)φn ( x) dx
2 2
=
= γ 02 + + γ N 2
Además:
b b
=f ,ΦN ∫=
a
f ( x)Φ N ( x)dx ∫ f ( x ) [γ φ ( x ) + γ φ ( x ) + + γ
a
0 0 1 1 N φN ( x)] dx
b b
= γ 0 ∫ f ( x)φ0 ( x)dx + + γ N ∫ f ( x)φN ( x)dx
a a
= γ 0 a0 + + γ N an
+ ( a0 2 + + aN 2 ) − 2 ( γ 0 a0 + + γ N an ) + ( γ 0 2 + + γ N 2 ) − ( a0 2 + + aN 2 )
2
= f
+ ( a0 2 − 2γ 0 a0 + γ 0 2 ) + + ( an 2 − 2γ n an + γ n 2 ) − ( a0 2 + + aN 2 )
2
= f
+ ( a0 − γ 0 ) + + ( an − γ n ) − ( a0 2 + + aN 2 )
2 2 2
= f
2
Finalmente, como EN = f − Φ N ≥ 0 trae como consecuencia que:
− ( a0 2 + + aN 2 ) ≥ 0
2
f
2
f ≥ a0 2 + + aN 2
Como la relación anterior es independiente del valor N , N puede ser tan grande como se
desee, lo que implica que:
∑ an 2 ≤ f
2
a0 2 + + aN 2 + =
n =0
∑a
2
n
2
≤ f .
n=0
∞
Como f es finito, la serie ∑a
n=0
n
2
es convergente, así an → 0 cuando n → ∞ , es decir:
b
lim f ,φn lim
= =
n →∞ ∫ f ( x)φn ( x)dx 0
n →∞
a
Como se vio en la sección anterior, para una función continua a pedazos f de periodo T
su serie de Fourier es:
a0 ∞ 2π 2π
f ( x) = + ∑ an cos nx + bn sen nx
2 n =1 T T
Donde:
T T
π
2 2 2π 2 2π
2 2
=a0 =∫
T −π
f ( x)dx an
T ∫ f ( x ) cos =
T T
nx dx b
n
T ∫ f ( x ) sen T nx dx
T
− −
2 2
a0 N 2π 2π
Φ N ( x) = + ∑ an cos nx + bn sen nx
2 n =1 T T
T T
Donde x ∈ − , . Luego:
2 2
T
2
a0 N 2π 2π
2
ΦN =∫T 2 + ∑
n =1
an cos
T
nx + bn sen
T
nx dx
−
2
T
2
a0 2 N 2 2π 2π 2π 2π
=∫T 4 + ∑n =1
an cos 2
T
nx + bn 2 sen 2
T
nx + 2 ∑ cos
n≠m T
nx sen
T
mx dx
−
2
T
N T T
a0 2 2
2
2 2π
2
2 2π
dx + ∑ an ∫ cos
4 ∫T
= 2
nx dx + bn ∫ sen
2
nx dx
n =1 T T
−
2
−
T
2
−
T
2
2
T T a
2
a N
T N
= 0 (T ) + ∑ an 2 + bn 2 = 0 + ∑ ( an 2 + bn 2 )
= 4 n 1= 2 2 2 2 n 1
T 2π 2π
El elemento 1,cos
se obtiene porque el conjunto nx ,sen nx no es
2 T T
normalizado, por consiguiente el error medio cuadrático:
2 a 2 N
− 0 + ∑ ( an 2 + bn 2 )
2
EN = f
T 2 n =1
a02 N
2 2 2
+ ∑ ( an 2 + bn 2 ) ≤ [ f ( x)] dx
2 2
2 n =1 T
f = ∫
T T
−
2
Y la función Φ que está formada por los primeros seis términos no nulos de la serie de
Fourier.
−2 si − π < x < 0
f ( x) = x si 0 ≤ x ≤ π
f ( x + 2π ) en otro caso
Donde:
2 − 2 ( −1) − π ( −1)
n n
π 2
a=
0 −2 a2 n −1 = − bn =
( 2n − 1) π
2
2 nπ
En consecuencia:
1 4 4
(π − 4 )
2
a02
= a12 = a32 =
4 π 2
81π 2
(4 + π ) (4 + π )
2 2
2 2 1 2
b
1 = b =
2 b 3 =
π2 4 9π 2
π −∫π
= f
T π −π 0 π 3 0
1 π3 π2
= = 4π + 4 +
π 3
3
2 a2 N
− 0 + ∑ ( an2 + bn2 )
2
E= f
T 2 n =1
π 1 1 328 10 ( 4 + π )
2
2
( )
2
= 4 + − + π − 4 + +
3 4 8 81π 2 9π 2
≈ 0.7957
Y la función Φ que está formada por los primeros siete términos no nulos de la serie de
Fourier.
La figura muestra que f es impar, lo que implica que an = 0 para toda n ∈ . Luego:
1 1 1
bn= ∫ f ( x)sen (π nx ) dx= 2 ∫ f ( x)sen (π nx ) dx= 2 ∫ ( − x + 1) sen (π nx ) dx
−1 0 0
1
1 x 1
2 −
= cos ( nπ x ) + cos ( nπ x ) − 2 2 sen ( nπ x )
nπ nπ nπ 0
1 1 1 2
− nπ cos ( nπ ) + nπ cos ( nπ ) − − nπ nπ
2=
En consecuencia:
4 1 4 1 4 1 4
= b12 = b22 = b32 = b42 = b52 = b62 = b72
2
π π 2
9π 2
4π 2
25π 2
9π 2
49π 2
2 a2 N
− 0 + ∑ ( an2 + bn2 )
2
E= f
T 2 n =1
2 4 1 4 1 4 1 4
= − 2 + 2 + 2 + 2 + + 2+
3 π π 9π 4π 25π 2
9π 49π 2
266681
= ≈ 0.05395
44100π 2
∑ an 2
2
= f
= si y solo si lim EN 0 .
N →∞
n=0
Corolario: Sea f una función continua a pedazos y periódica de periodo T con serie de
Fourier
a0 ∞ 2π 2π
+ ∑ an cos nx + bn sen nx
2 n =1 T T
a0 2 ∞ 2
+ ∑ ( an 2 + bn 2 ) =f
2
Entonces si y solo si lim EN = 0 .
2 n =1 T N →∞
Corolario: Supóngase que f que satisface las hipótesis del corolario anterior. Si f es
par, entonces:
T
2
a0 ∞
2 4 2
+ ∑ an 2 = [ f ( x)] dx .
2 2
2 n =1 T
f = ∫
T 0
T ∫0
f =
2 n =1 T
π π
1 1 1
=a0 ∫π nx ) dx
f ( x) cos ( = ∫π x + 4 x
2
cos ( nx ) dx
π −
π −
1 2n ( 2 + x ) cos ( nx ) + ( −2 + n x ( 4 + x ) ) sen ( nx )
π
2
=
π 4n 3
−π
1 2n ( 2 + π ) cos ( nπ ) 2n ( 2 − π ) cos ( nπ )
= −
π 4n 3 4n 3
1 2n ( 2 + π )( −1) 2n ( 2 − π )( −1)n ( −1)
n n
= − =
π 4n 3
4n3
n2
π π
1 1 1
=b0 ∫π f ( x)sen ( =
nx ) dx ∫π x + 4 x
2
sen ( nx ) dx
π −
π −
1 ( 2 − n x ( 4 + x ) ) cos ( nx ) + 2n ( 2 + x ) sen ( nx )
π
2
=
π 4n 3
−π
1 ( 2 − n π ( 4 + π ) ) cos ( nπ ) ( 2 + n π ( 4 − π ) ) cos ( − nπ )
2 2
= −
π 4n 3 4n 3
= − =−
π 4n 3 4n 3 n
π 2 ∞ ( −1) 2 ( −1)
n n
f (=
x) + ∑ 2 cos ( nx ) − sen ( nx )
12 n =1 n n
2 n =1 T
∞ −
2 2
( 1)
1π2
2 n
2 ( −1)n 1
π
1 2
2
=
∑ + +−
π −∫π
x + x dx
2 6 n =1 n 2
n
4
π
π4 1 4 1 x3 x 4 x5
∞
+ ∑ 4 + 2 = + +
72 n =1 n n π3 8 80 −π
∞
1 + 4n 2 2π 2 π 4 π 4
∑ n4
= + −
n =1 3 40 72
∞
1 + 4n 2 π 2
Por lo tanto ∑ 4
= ( 60 + π 2 ) .
n =1 n 90
−1 0 3 0 3
Además:
1 1
=an ∫ f= ∫ x cos ( nπ x ) dx
( x) cos ( nπ x ) dx 2= 2
−1 0
= 2
n3π 3
0
= 2=
n3π 3 nπ
2 2
a0 2 ∞ 2 4 2
+ ∑ an = [ f ( x)] dx
2
2 n =1 T ∫0
2
∞
4 ( −1)
2 n 1
12 4 2
∑ + = ∫ x 2 dx
2 3 n =1 n 2π 2
20
1
2 ∞ 16 1
+∑ 4 4 = 2 x5
9 n =1 n π 5 0
2 16 ∞
1 2
+ 4∑ 4 =
9 π n =1 n 5
∞
1 π4 2 2
∑
=
n =1 n
4 −
16 5 9
Por lo tanto:
∞
1 π4
∑
n =1 n
4
=
90
.
2
1 nπ x 2
2 2
2 nπ x nπ x
=bn = ∫ f ( x )sen dx =
2 ∫ (1) sen dx − cos
4 −2 2 2 0 2 nπ 2 0
= =
2
−
nπ
cos (π n ) +
n
2 2
π
nπ
1 − ( −1)
n
( )
Por consiguiente:
4
b2 n −1 =
( 2n − 1) π
n =1 T0
2
∞ 4 42 2
∑ = ∫ [1] dx
n =1 ( 2n − 1) π 40
∞
16
∑ = [ x ]0
2
( 2n − 1) π
2 2
n =1
16 ∞
1
∑ =2
( 2n − 1)
2
π 2
n =1
∞
1 π2
∑ =
( 2n − 1)
2
n =1 18
Actividades
Autorreflexiones
Cierre de la Unidad
En esta Unidad aprendiste la definición y las propiedades del producto interno y como
este caracteriza la ortogonalidad de funciones, luego estudiaste funciones periódicas y
sus propiedades fundamentales. Luego estudiaste el concepto la serie de Fourier
trigonométrica y su generalización de una función periódica. Finalmente, aprendiste a
calcular el error medio cuadrático de una aproximación junto con la desigualdad de Bessel
y la identidad de Parseval-Liapunov.
Para más tener ejemplos de series de Fourier de funciones periódicas puedes consultar
las siguientes páginas:
• http://www.fourier-series.com/
• https://math24.net/fourier-series-definition-typical-examples.html
Fuentes de consulta
• Churchill, R., Brown, J. (2011). Fourier series and boundary valued problems. 8
edición, USA: Mc Graw Hill.