T11 Series Fourier

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 16

Tema 11: SERIES DE FOURIER

11.1 Introducción.
11.2 Nociones previas.
11.3 Serie trigonométrica de Fourier.
11.4 Forma Amplitud-Fase de la serie de Fourier.
11.5 Convergencia de la serie de Fourier.
11.6 Serie de Fourier en senos y cosenos.
11.7 Serie exponencial de Fourier.
11.7.1 Fórmulas de paso.
11.7.2 Espectro de frecuencias.
11.8 Identidad de Parseval.
11.9 Ejercicios propuestos.

Fundamentos de matemáticas II Escuela Politécnica de Cuenca


3

11.1 Introducción

En este tema vamos a introducir las series de Fourier las cuales constituyen la
herramienta matemática básica del análisis de Fourier empleado para analizar fun-
ciones periódicas a través de la descomposición de dicha función en una suma innita
de senos y cosenos.

Fue el matemático francés J.B.J. Fourier (1758- 1830) quien, en diciembre de 1807
escribió Théorie Analytique de la Chaleur (Teoría Analítica del Calor), donde reali-
zó un estudio detallado de las series trigonométricas las cuales utilizó para resolver
varios problemas relacionados con la conducción de calor. Demostró que cualquier
función periódica se puede representar como una serie innita de funciones senos y
cosenos, cada una con frecuencias angulares múltiplos enteros de la frecuencia angu-
lar fundamental. Los argumentos establecidos por Fourier eran imprecisos y en 1829
P.G.L. Dirichlet proporcionó las condiciones precisas para que una señal periódica
pueda ser representada por una serie de Fourier.

Hoy en día las series de Fourier son importantes, tanto práctica como teórica-
mente y siguen siendo objeto de investigación. Es una herramienta usada en muchas
ramas de la ingeniería: análisis vibratorio, acústica, óptica, procesamiento de imá-
genes y señales, y compresión de datos.

11.2 Nociones previas

Denición.- Se dice que una función f : R −→ R (ó C) es periódica de periodo


T si ∀t ∈ R f (t+T ) = f (t) siendo T > 0. El mínimo valor de T es el periodo de f (t).

Denición.- Una función f se dice que es continua a trozos en un intervalo [a, b] si


existe una partición a = x0 < x1 < x2 < · · · < xn−1 < xn = b del intervalo [a, b] de
forma que:
(1) f es continua en cada intervalo (xi , xi+1 ) ∀i = 1, 2, 3, · · · , n.
(2) f tiene límites laterales en los puntos xi ∀i = 1, 2, 3, · · · , n.
Denición.- Diremos que una función f es derivable a trozos en un intervalo [a, b]
si existe una partición a = x0 < x1 < x2 < · · · < xn−1 < xn = b del intervalo [a, b]
de forma que:
(1) f es derivable en cada intervalo (xi , xi+1 ) ∀i = 1, 2, 3, · · · , n.
(2) La función f 0 tiene límites laterales en los puntos xi ∀i = 1, 2, 3, · · · , n.
Proposición. Toda función derivable a trozos en un intervalo también es continua a
trozos en dicho en dicho intervalo. Las funciones continuas a trozos en un intervalo
son integrables en dicho intervalo. Además la integral la podemos calcular como suma
de integrales en cada uno de los intervalos donde la función es continua.

Fundamentos de matemáticas II Escuela Politécnica de Cuenca


4

Ejemplo. Funciones periódicas

(a) f (t) = sen(t) (b) f (t) = cos(t)

(c) Onda triangular (d) Diente de sierra

Figura 1: Funciones periódicas

Denición.-Sea f : R −→ R:
Decimos que la función f es función par si f (−t) = f (t).
Decimos que la función f es función impar si f (−t) = −f (t).
Proposición.
El producto de dos funciones pares es par.
El producto de dos funciones impares es par.
El producto de una función par por otra impar es impar.
Z a Z a
Si f es par ⇒ f (t)dt = 2 f (t)dt.
−a 0
Z a
Si f es impar ⇒ f (t)dt = 0
−a

Denición.- Se dice que las funciones {fk (t)} son ortogonales en el intervalo (a, b)
si dos funciones cualesquiera de dicho conjunto cumple:
Z b
fm (t)fn (t) dt = 0 m 6= n
a

Tema 11: Series de Fourier


5

Ejemplo.
La funciones f1 (t) = t y f2 (t) = t2 son ortogonales en (−1, 1).
La funciones f1 (t) = sen t y f2 (t) = cos t son ortogonales en (−π/2, π/2).

En la teoría de funciones hay muchas familias de funciones ortogonales. Para la


teoría de Fourier hay una familia particularmente importante, que es la familia de
funciones sinusoidales, es decir:
{1, sen ω0 t, sen 2ω0 t, sen 3ω0 t, · · · , cos ω0 t, cos 2ω0 t, cos 3ω0 t, · · · }

Por tanto:
Z T /2 Z T /2 Z T /2
cos(nω0 t) dt = sen(nω0 t) dt = cos(nω0 t) sen(nω0 t) dt = 0 ∀n ∈ N
−T /2 −T /2 −T /2
(1)
y ∀n, m ∈ N, m 6= n
Z T /2 Z T /2 Z T /2
cos(nω0 t) cos(mω0 t) dt = sen(nω0 t) sen(mω0 t) dt = cos(nω0 t) sen(mω0 t) dt = 0
−T /2 −T /2 −T /2
(2)

11.3 Serie trigonométrica de Fourier

Denición.- Sea f función 2π−periódica, integrable en el intervalo [−π, π]. Se dene


desarrollo en serie trigonométrico de Fourier a la serie:

a0 X
f (t) = + (an cos(nt) + bn sen(nt)) (3)
2 n=1

a0 , an y bn reciben el nombre de coecientes de Fourier y quedan completamente


determinados por la función f .

Para determinar la fórmula del coeciente a0 integramos (3) y aplicando (1)


obtenemos:
Z π Z π ∞ Z π ∞ Z π
a0 X X
f (t) dt = dt + an cos(nt) dt + bn sen(nt) dt = πa0
−π −π 2 n=1 −π n=1 −π

Z π
1
⇒ a0 = f (t) dt.
π −π

Determinemos ahora an para n ≥ 1. Para ello multiplicamos ambos miembros


de la expresión (3) por cos(nt), integramos y por (2) obtenemos:
Z π
1
an = f (t) cos(nt) dt n ≥ 1
π −π

Fundamentos de matemáticas II Escuela Politécnica de Cuenca


6

De manera análoga calculamos bn para n ≥ 1. Multiplicando ambos miembros de la


expresión (3) por sen(nt), integramos y por (2) obtenemos:
Z π
1
bn = f (t) sen(nt) dt n ≥ 1
π −π

El intervalo de integración no necesita ser simétrico respecto al origen. Como la


ortogonalidad de las funciones seno y coseno no solo se da en el intervalo [−T /2, T /2],
sino en cualquier intervalo que cubra un periodo completo, las fórmulas anteriores
son válidas si en la integral reemplazamos el intervalo [−π, π] por otro intervalo
cualquiera de longitud T : [t0 , t0 + T ], con t0 arbitrario.

Denición.- Sea f función periódica de periodo T , integrable en el intervalo [t0 , t0 + T ].


Se dene desarrollo en serie trigonométrico de Fourier a la serie:

a0 X
f (t) v + (an cos(nω0 t) + bn sen(nω0 t)) (4)
2 n=1


siendo ω0 = la frecuencia angular.
T
a0 , an y bn reciben el nombre de coecientes de Fourier y se determinan:

Z t0 +T
2
a0 = f (t) dt.
T t0

Z t0 +T
2
an = f (t) cos(nω0 t) dt.
T t0

Z t0 +T
2
bn = f (t) sen(nω0 t) dt.
T t0

11.4 Forma Amplitud-Fase de la serie de Fourier

Podemos escribir el desarrollo en serie (4) con una expresión más simplicada de
la siguiente manera: multiplicamos y dividimos el término an cos(nω0 t)+bn sen(nω0 t)
por an + b2n :
p
2

!
an bn
(5)
p
a2n + b2n p cos(nω0 t) + p sen(nω0 t)
a2n + b2n a2n + b2n

y consideramos el triángulo rectángulo:

Tema 11: Series de Fourier


7

donde:
bn a
sen θn = p
2 2
; cos θn = p 2 n 2 ;
an +  bn an + b n
θn = arctan abnn

Entonces, sustituyendo en (5), obtenemos:

An [cos θn cos(nω0 t) + sen θn sen(nω0 t)]


a
Teniendo en cuenta que cos(α − β) = cos α cos β + sen α sen β y deniendo A0 = 0
2
el desarrollo en serie trigonométrico (4) queda

X
f (t) v A0 + An cos(nω0 t − θn )
n=1

La serie de Fourier así escrita se conoce Forma amplitud-fase de la serie de


Fourier donde

An es la amplitud.

θn es la fase.

A la componente sinusoidal An cos(nω0 t − θn ) se le llama n-ésimo armónico


de f (t).
Al primer armónico (n = 0) se le llama armónico fundamental.

Nota.- La serie de Fourier también puede escribirse de la forma



X
f (t) v A0 + An sen(nω0 t + θn )
n=1

Queda como ejercicio para el estudiante denir correctamente A0 , An y θn para este


caso.

Ejemplo. Si calculamos los coecientes de Fourier de f (t) = t denida sobre


[−π, π], obtenemos

2(−1)n+1 π
a0 = 0; an = 0; bn = ; θn =
n 2
Por lo que
∞ ∞
X 2(−1)n+1 π X 2(−1)n+1
f (t) v cos(nt − ) = sen(nt)
n=1
n 2 n=1
n

Fundamentos de matemáticas II Escuela Politécnica de Cuenca


8

Representamos f y algunas de sus sumas parciales (ver Figura 2). Podemos ver
como según aumentamos los términos en las sumas parciales, la aproximación es
cada vez mejor a la función original. Sin embargo, en los extremos de denición,
que se corresponden con los puntos en que es discontinua la extensión periódica de
la función, parece que la serie de Fourier no aproxima bien a la función original.
Tal fenómeno se conoce como el fenómeno de Gibbs.

Figura 2: Representación de f (t) y f1 , f2 , f5 y f10

Proposición. Si la función f denida en el intervalo [t0 , t0 + T ] es periódica, la


representación de la función en serie de Fourier es válida para toda la recta real

puesto que las funciones sen(nω0 t) y cos(nω0 t) son periódicas de periodo T = .
ω0

11.5 Convergencia de la serie de Fourier

Veamos ahora en que condiciones razonablemente generales la serie de Fourier


de una función converge puntualmente a dicha función.

Teorema. (Condiciones de Dirichlet). Sea f: [t0 , t0 + T ] → R una función


periódica derivable a trozos en su intervalo de denición. Entonces ∀t ∈ [t0 , t0 + T ]
se cumple que:

(1) Si f es continua en t ∈ [t0 , t0 + T ], la serie de Fourier converge a f .


(2) Si t es punto de discontinuidad de f , entonces la serie de Fourier converge a
f (t+ ) + f (t− )
, siendo f (t+ ) y f (t− ) los límites laterales por la derecha y por
2
la izquierda de f en t respectivamente.

Tema 11: Series de Fourier


9

Nota.- Las condiciones anteriores son sucientes, por lo que puede haber funcio-
nes que aun no cumpliendolas sean convergentes. En la práctica, por lo general, se
cumplen.

11.6 Serie de Fourier en senos y cosenos

Veamos cómo se nos simplica el trabajo si la función presenta simetría.

 
−T T
Proposición. Sea f : , −→ R:
2 2

(1) Si f es una función impar. Entonces:

ˆ f (t) cos(nω0 t) es impar ⇒ a0 = an = 0


4 T /2
Z
ˆ f (t) sen(nω0 t) es par ⇒ bn = f (t) sen(nω0 t)dt
T 0
 
−T T
Es decir, la serie de Fourier de una función impar en el intervalo ,
2 2
es una serie en la que solo aparecen senos:


X
f (t) v bn sen(nω0 t)
n=1

(2) Si f es una función par. Entonces:


Z T /2 Z T /2
4 4
ˆ f (t) cos(nω0 t) es par ⇒ a0 = f (t)dt y an = f (t) cos(nω0 t)dt
T 0 T 0

ˆ f (t) sen(nω0 t) es impar ⇒ bn = 0


 
−T T
Es decir, la serie de Fourier de una función par en el intervalo , es
2 2
una serie en la que solo aparecen cosenos:


a0 X
f (t) v + an cos(nω0 t)
2 n=1

 
T
¾Qué hacer si f es una función denida en el intervalo 0, , y se busca obtener
2
un desarrollo en serie de Fourier de solo cosenos (o de solo senos) que la represente?
−T
Para ello extenderemos la función al intervalo , 0 . Veamos cómo.
2

Fundamentos de matemáticas II Escuela Politécnica de Cuenca


10

 
T
Denición.- Sea f : 0, → R, se llama extensión periódica impar de f a la
  2
−T T
función f1 : , −→ R denida como:
2 2
(
−f (−t) − T /2 ≤ t < 0
f1 (t)
f (t) 0 ≤ t ≤ T /2

El desarrollo en serie de Fourier de f1 , dado que es una función impar, verica


que a0 = an = 0, por lo que se le denomina serie de Fourier de senos de f y viene
dada por:

X
f (t) v bn sen(nω0 t)
n=1

 
T
Denición.- Sea f : 0, → R, se llama extensión periódica par de f a la función
  2
−T T
f2 : , −→ R denida como:
2 2
(
f (−t) − T /2 ≤ t < 0
f2 (t)
f (t) 0 ≤ t ≤ T /2

El desarrollo en serie de Fourier de f2 , dado que es una función par, verica que
bn = 0, por lo que se le denomina serie de Fourier de cosenos de f y viene dada por:


a0 X
f (t) v + an cos(nω0 t)
2 n=1

11.7 Serie exponencial de Fourier

En muchas aplicaciones de las series de Fourier es conveniente expresar dicha


serie en términos de la exponencial compleja.

Denición.- Sea f : [t0 , t0 + T ] → R cuya serie trigonométrica de Fourier es:



a0 X
f (t) v + an cos(nω0 t) + bn sen(nω0 t) (6)
2 n=1

Sustituimos en (6) el seno y el coseno por su expresión exponencial:

eiϕ − e−iϕ eiϕ + e−iϕ


sen ϕ = cos ϕ =
2i 2
Tema 11: Series de Fourier
11


a0 X einω0 t + e−inω0 t einω0 t − e−inω0 t
f (t) = + an + bn =
2 n=1
2 2i

Operamos y agrupamos términos:



a0 X 1 1 1 1
= + an einω0 t + an e−inω0 t + bn einω0 t − bn e−inω0 t =
2 n=1
2 2 2i 2i
1
Teniendo en cuenta que = −i:
i

a0 X 1 1
= + [an − ibn ] einω0 t + [an + ibn ] e−inω0 t
2 n=1
2 2
a0
LLamamos C0 = , Cn = 12 [an − ibn ] y entonces 1
2
[an + ibn ] = Cn :
2


X
Cn enω0 ti + Cn e−nω0 ti
 
f (t) = C0 +
n=1

Z t0 +T
1
C0 = f (t) dt.
T t0

1 2 t0 +T 2 t0 +T
 Z Z 
1
Cn = [an − ibn ] = f (t) cos(nω0 t)dt − i f (t) sen(nω0 t)dt =
2 2 T t0 T t0
1 t0 +T 1 t0 +T
Z Z
= f (t) [cos(nω0 t) − i sen(nω0 t)] dt = f (t)e−nω0 ti dt
T t0 T t0

Z t0 +T
1
=⇒ Cn = f (t)e−inω0 t dt
T t0

Z t0 +T
1 1
Cn = [an + ibn ] = f (t)enω0 ti dt = C−n .
2 T t0

Así:
X∞ X∞ X∞ −1
X ∞
X
inω0 t −inω0 t inω0 t inω0 t
f (t) = C0 + Cn e + C−n e = C0 + Cn e + Cn e = Cn einω0 t
n=1 n=1 n=1 n=−∞ n=−∞

Por tanto, la serie exponencial (o compleja) de Fourier será de la forma:


X
f (t) = Cn einω0 t
n=−∞

Fundamentos de matemáticas II Escuela Politécnica de Cuenca


12

Nota.- Puesto que la serie exponencial de Fourier no es más que otra manera de
escribir la serie trigonométrica de Fourier, los teoremas de convergencia se aplican
sin necesidad de ninguna adaptación.

Nota.- Si la señal es periódica su representación en serie exponencial de Fourier es


válida para toda la recta real ya que einω0 t = einω0 (t+T ) para T = 2π/ω0

11.7.1 Fórmulas de paso

Veamos cómo pasar de la serie trigonométrica a la exponencial y viceversa.

Trigonométrica a exponencial
a0
ˆ C0 = .
2
1
ˆ Cn = [an − ibn ]
2
Exponencial a trigonométrica
ˆ a0 = 2C0
ˆ an = Cn + C−n = 2Re(Cn )
ˆ bn = i (Cn − C−n ) = −2Im(Cn )

11.7.2 Espectro de frecuencia compleja

El espectro es la representación de una señal en dominio frecuencia. Hemos visto


que dada una función periódica f (t), le corresponde una y solo una serie de Fourier,
es decir le corresponde un conjunto único de coecientes Cn . Estos coecientes son
números complejos, por lo que su expresión en forma exponencial sería:
Cn = |Cn |eφn i
1
donde, teniendo en cuenta que Cn = [an − bn i],
2
 
1p 2 bn
|Cn | = an + b2n and φn = arctan − ∀n 6= 0
2 an
a0
Para n = 0, C0 es un número real, con lo que |C0 | = | | y φ0 = 0.
2
Por ello, la representación de los coecientes complejos Cn frente a la variable
nω0 especica la función periódica f en el dominio frecuencia al igual que la repre-
sentación (t, f (t)) especica la función en el dominio tiempo.

Denición.- El espectro de amplitud de la función periódica f es la gráca de los


puntos (nω0 , |Cn |), siendo |Cn | la magnitud del coeciente complejo Cn .

Tema 11: Series de Fourier


13

Denición.- Se le llama espectro de fases de f (t) a la representación del ángulo


fase φn de Cn frente a nω0 ((nω0 , φn )) .

Nota.- Recordar que θn = arctan


 
bn
an
, y dado que la función arcotangente es im-
par, φn = −θn .

Dado que n solo toma valores enteros, nω0 es una variable discreta y por tanto
el espectro es una gráca discreta formada por líneas verticales (componentes espec-
trales) regularmente espaciadas.

Proposición. El desplazamiento en el tiempo de una función periódica no tiene


efecto sobre el espectro de amplitud pero modica el espectro de fase en una cantidad
de −nω0 K radianes para la componente de frecuencia nω0 si el desplazamiento en
el tiempo es K .

Nota.- Si la señal es real |Cn | = |C−n | y θn = −θ−n

Nota.- Si para obtener los espectros de frecuencia utilizamos la forma Amplitud-


Fase, es decir, An y θn (en lugar de Cn y φn ), se obtienen los espectros llamados de
una sola cara (solo para n positivos).

11.8 Identidad de Parseval

Denición.- El contenido de potencia de una señal periódica f (t) en un periodo T


se dene como el error cuadrático medio:
Z T /2
1
f 2 (t)dt
T −T /2

y representa la energía disipada por una resistencia de un ohm cuando se aplica una
onda f (t) de voltaje o corriente.

Si f (t) es periódica también lo será f 2 (t) y el promedio en un periodo será el


promedio en cualquier periodo.

El siguiente teorema nos permite calcular la integral de f 2 (t) mediante los coe-
ciente de Fourier o mediante los coecientes complejos.
Teorema. (Identidad de Parseval)
En términos de los coecientes de Fourier:
Z T /2 ∞
1 1 1X 2
f (t)dt = a20 +
2
(an + b2n )
T −T /2 4 2 n=1

Fundamentos de matemáticas II Escuela Politécnica de Cuenca


14

O en términos de los coecientes complejos:


Z T /2 ∞
1 X
2
f (t)dt = |Cn |2
T −T /2 n=−∞

11.9 Ejercicios propuestos

1. Representa y desarrolla en serie trigonométrica de Fourier la señal periódica


denida por:
(
0 si −1<t<0
f (t) =
1 si 0<t<1

2. (a) Hallar el desarrollo en serie trigonométrica de Fourier la señal periódica


denida por:
(
−π/2 si −π <t<0
f (t) =
π/2 si 0<t<π

(b) Dene la serie en los puntos de discontinuidad de la función.

3. Representar y clasicar, por pariedad, las funciones periódicas siguientes:


(
−π/2 si −π <t<0
a) f (t + 2π) = f (t) =
π/2 si 0<t<π
(
cos t si 0 < t < π
b) f (t + 2π) = f (t) =
0 si π < t < 2π
(
sen t si 0 < t < π
c) f (t + 2π) = f (t) =
0 si π < t < 2π
d) f (t + 4) = f (t) = 4t − t2 0<t<4

4. Calcular el desarrollo en serie trigonométrico de Fourier de la señal:


t
f (t) = −π ≤t≤π
2
Dene la serie en los puntos de dicontinuidad de la función.
5. Con los resultados obtenidos de los ejercicios 2 y 4 desarrolla en serie:
 −π t
 − si −π <t<0
f (t) = π2 t 2
 − si 0<t<π
2 2
Tema 11: Series de Fourier
15

6. Dada la señal :

f (t) = f (t + 2) = |t| −1≤t≤1

Se pide:

a) Escribir la serie trigonométrica de Fourier.


b) Escribir la serie exponencial de Fourier.
c) Utilizando la Identidad de Parseval, calcular las sumas:
∞ ∞
X 1 X 1
a) b)
n=1
(2n − 1)4 n=1
n4

d) Dibujar el espectro de amplitud.


7. Calcular el desarrollo en serie trigonométrico de Fourier de la señal:

f (t) = f (t + 2π) = |t| −π <t<π

Del resultado obtenido calcular:


∞ ∞
X 1 X 1
a) b)
n=1
(2n − 1)2 n=1
n2

Pista: Tener en cuenta el valor de la función en t = 0


8. Hallar la serie exponencial de Fourier de la función diente de sierra denida
por:
A
f (t) = f (t + T ) = t 0≤t≤T (T = 1)
T
Con lo obtenido en el apartado anterior calcula los coecientes de Fou-rier.
9. Desarrollar en serie exponencial de Fourier la señal denida por:

f (t) = f (t + 2π) = et 0 ≤ t ≤ 2π

10. Desarrollar en serie exponencial de Fourier la señal denida por:

f (t) = f (t + T ) = sen2 t 0 ≤ t ≤ 2π

Ayuda : En este caso es conveniente sustituir la función seno por su expresión


exponencial.
11. Representa el espectro de amplitud de la señal:
(
A sen t si 0 < t < T /2
f (t) = f (t + 2π) =
0 si T /2 < t < T

Fundamentos de matemáticas II Escuela Politécnica de Cuenca


16

12. Calcula la serie de Fourier seno de f (t) = 1 para t ∈ [0, π]


13. Calcula la serie de Fourier coseno de f (t) = t para t ∈ [0, π]
14. Sea f (t) = t(1 − t) 0 ≤ t ≤ 1. Escribe la extensión impar de f de periodo 2.
a) Calcula la serie de Fourier seno de f .

1
b) Escribe la identidad de Parseval y utilizala para calcular
X

n=1
(2n − 1)6

15. Sea la señal f (t) = f (t + 10) = 2t 0 ≤ t ≤ 10. Se pide:

a) Calcular la serie exponencial de Fourier.


b) Calcula los coecientes de Fourier.
c) Espectro de amplitud.
16. Desarrolla la serie trigonométrica de Fourier de
(
0 si 0<t<π
f (t) = f (t + 2π) =
π si π < t < 2π

Dene sus puntos de discontinuidad.


17. a) Calcula la serie trigonométrica de Fourier para la función:
(
0 si −5<t<0
f (t) = f (t + 10) =
3 si 0<t<5

b) ¾Qué valor tiene la señal en los puntos de discontinuidad?.

c) Utilizando las fórmulas de paso, expresa la serie en forma exponencial.


18. a) Calcula la serie trigonométrica de Fourier para la función:
(
1 si −5<t≤0
f (t) = f (t + 10) =
1+t si 0<t<5

b) Calcular el valor de la serie en los puntos de discontinuidad.


19. Sea la señal periódica f (t + 2π) = f (t) = t2 − π ≤ t ≤ π . Se pide:

a) El desarrollo en serie de Fourier.



(−1)n−1
b) Calcular la suma .
X

n=1
n2

Tema 11: Series de Fourier

También podría gustarte