Series de Fourier

"Series de Fourier, Transformadas de Fourier y Aplicaciones",

Genaro González
16/11/2022 MSc. RAUL P. CASTRO VIDAL 1
UNMSM-FIEE
La primera serie de Fourier de la historia
Euler 1744 escribe en una carta a un amigo:

 t 
sen(nt )
f (t )   
2 n 1 n
sen(2t ) sen(3t )
sen(t )    ...
2 3
¿Es cierto?

Observemos que en t = 0
hay problemas → π/2 = 0 ¡¡
La clave está en el concepto de función periódica.
16/11/2022 MSc. RAUL P. CASTRO VIDAL 2
UNMSM-FIEE
 Funciones Periódicas
Una función periódica f(t) cumple que para todo
valor de t:
f(t) = f(t + T).
Al valor mínimo, mayor que cero, de la constante
T que cumple lo anterior se le llama el periodo
fundamental (o simplemente periodo) de la
función.
Observa que:
f(t) = f(t + nT), donde n = 0, 1,  2, 3,...

16/11/2022 MSc. RAUL P. CASTRO VIDAL 3

Cuestión: ¿Es f(t) = cte. una función periódica?
UNMSM-FIEE
Ejemplo: ¿Cuál es el periodo de la función

f(t)  cos ( 3t )  cos ( t

4 )?
Si f(t) es periódica se debe cumplir:
f(t  T)  cos ( t T
3 )  cos ( t T
4 )  f(t)  cos ( 3t )  cos ( t
4 )

Como cos(t + 2k) = cos(t) para cualquier entero k,

entonces, para que se cumpla la igualdad, se requiere
que:
T/3 = 2k1 y T/4 = 2k2.
Es decir:
T = 6k1 = 8k2
con k1 y k2 enteros.
El valor mínimo de T se obtiene con k1= 4, k2= 3, es
decir, T = 24. MSc. RAUL P. CASTRO VIDAL 4
16/11/2022 UNMSM-FIEE
 Gráfica de la función f(t)  cos ( 3t )  cos ( t
4 )
3
T f(t)=cos(t/3)+cos(t/4)
2

1
f(t)

-1

-2
24
-3
0 50 100 150 200
t

16/11/2022 MSc. RAUL P. CASTRO VIDAL 5

UNMSM-FIEE
¿Es la suma de dos funciones
periódicas una función periódica?
Depende. Consideremos la función:

f(t) = cos(w1t) + cos(w2t).

Para que sea periódica se requiere encontrar dos

enteros m, n tales que:
w1T = 2 m y w2T = 2 n.
Es decir, que cumplan:
w1 m
T = m/ (2 w1) = n/ (2 w2) 
16/11/2022 MSc. RAUL P. CASTRO VIDAL w2 n 6
UNMSM-FIEE
Ejemplo: para la función cos(3t) + cos((+3)t)
tenemos que w 3 1

w2 3 
¿Es periódica? f(t)=cos(3t)+cos((3+π)t)

1
f(t)

-1

-2
0 5 10 15 20 25 30
16/11/2022 t VIDAL
MSc. RAUL P. CASTRO 7
UNMSM-FIEE
Para que exista periodicidad w1/ w2 debe ser
un número racional (n/m).
Ejercicios: Encontrar el periodo de las
siguientes funciones, si es que son periódicas:

1) f(t) = sen(nt), donde n es un entero.

2) f(t) = sen2(2t)
3) f(t) = sen(t) + sen(t + /2)
4) f(t) = sen(w1t) + cos(w2t)
5) f(t) = sen(2 t)
16/11/2022 MSc. RAUL P. CASTRO VIDAL 8
UNMSM-FIEE
 Si f1(t) tiene periodo T1 y f2(t) tiene periodo T2,
¿es posible que f1(t) + f2(t) tenga periodo
T < min(T1,T2)?

T1 = 5

T2 = 5

T = 2,5

16/11/2022 MSc. RAUL P. CASTRO VIDAL 9

UNMSM-FIEE
Podemos construir incluso un ejemplo de dos funciones de
igual periodo, cuya suma puede tener un periodo tan
pequeño como queramos. Sea N un entero, y definamos:
 1  1
0t 
sen(2 Nt ), 0  t  N 0, N
f1 (t )   f 2 (t )  
1 1
 0,  t 1 sen(2 Nt ),  t 1
 N  N
extendida periódicamente con T = 1: extendida periódicamente con T = 1:
f1 (t )  f1 (t  1),    t   f 2 (t )  f 2 (t  1),    t  
 sen(2 Nt ) , 0  t 1
f1 (t )  f 2 (t )  
 f1 (t  1)  f 2 (t  1),    t  
2 2 1
T  
16/11/2022
w 2 N N 10
MSc. RAUL P. CASTRO VIDAL
UNMSM-FIEE
 ¿Puede una función f(t) cumplir la condición
f(t) = f(t + T) para todo t y no tener un periodo
fundamental?

1 si t es un entero
f1 (t )  
0 si t no es un entero
1 si t y t  T son enteros
f1 (t )  f1 (t  T )  
0 si t y t  T no son enteros
 T 1
16/11/2022 MSc. RAUL P. CASTRO VIDAL 11
UNMSM-FIEE
 1 si t es racional pero no un entero
f 2 (t )  
0 si t es irracional o es un entero
1 si t y t  T son racionales pero no enteros
f 2 (t )  f 2 (t  T )  
0 si t y t  T son irracionales o enteros
 T 1

 1 si t es racional
f1 (t )  f 2 (t )  
 0 si t es irracional
T=?
16/11/2022 MSc. RAUL P. CASTRO VIDAL 12
UNMSM-FIEE
Volvamos al resultado  t sen(2t ) sen(3t )
 sen t    ...
de Euler: 2 2 3
S (t )  eit  ei 2t  ei 3t  ...
¿Cómo lo alcanzó?  it
e S (t )  ei 2t  ei 3t  ...
Utilizando la fórmula de eit 1 1 sen t
Euler para cada término: S (t )    i
1 e it
2 2 1  cos t
S (t )  e  e  e  ... 
it i 2t i 3t

cos t  cos(2t )  cos(3t )  ...  i sen t  sen(2t )  sen(3t )  ...

1
 sen(2t ) sen(3t ) 1
Integrando
término a término:
2
sen t    ...   t  C
2 3 2
Particularizamos t  1 1 1  
para encontrar C: t   1     ...    C ; C 
2 3 5 7 4 2
16/11/2022  13
4
MSc. RAUL P. CASTRO VIDAL UNMSM-
  t sen(2t ) sen(3t )
 sen t    ...
2 2 3
 t sen(2t ) sen(3t )
 sen( t )    ...
2 2 3
t  sen(2t ) sen(3t )
   sen(t )    ...
2 2 2 3
Fourier series java applet
(http://www.falstad.com/fourier/)
16/11/2022 MSc. RAUL P. CASTRO VIDAL UNMSM- 14
FIEE
(1) La función de Euler es periódica de periodo T = 2π.

(2) La serie es una función impar.

No es sorprendente, pues se trata de suma de senos de
periodos enteros.

(3) En el intervalo 0 < t < 2π, la serie aproxima a (π-t)/2.

Pero no fuera del intervalo...

(4) Da saltos bruscos entre valores positivos y negativos.

(5) La aproximación no es buena en "los extremos"...

Ninguna de estas dos últimas cuestiones era conocida o
sospechada ni por Euler, ni por Fourier... 15
16/11/2022 MSc. RAUL P. CASTRO VIDAL UNMSM-
FIEE
 Leonhard Euler
Jean 1707-1783
d'Alembert
1717-1783

Daniel
Lagrange
Bernouilli
1700-1782

16/11/2022 16
MSc. RAUL P. CASTRO VIDAL UNMSM-
FIEE
Se necesita también como condición inicial u(0,x)=f(x) para 0<x<1.
Euler en 1749 demostró la misma solución. Pero difería con D'Alambert en el posible
tipo de f(x) inicial. De hecho, este es el inicio del problema de la "definición" de una
función. Para Euler era posible una función en partes: cualquier gráfica era una 17 función
Para D'Alambert necesariamente: expresión analítica compacta.
Se necesita también como condición inicial
u(0,x)=f(x) para 0<x<1.
Euler en 1749 demostró la misma solución.
Pero difería con D'Alambert en el posible
tipo de f(x) inicial. De hecho, este es el inicio
del problema de la "definición" de una
función. Para Euler era posible una función en
partes: cualquier gráfica era una función.
Para D'Alambert necesariamente: expresión
analítica compacta.
16/11/2022 MSc. RAUL P. CASTRO VIDAL 18
UNMSM-FIEE
16/11/2022 MSc. RAUL P. CASTRO VIDAL 19
UNMSM-FIEE
En realidad la forma de solucionar el problema por parte
de Daniel Bernoulli en 1753 fue completamente distinta.
Se basó en la superposición de ondas y tomó como
solución:

un(x,t) = sin(nx) cos(nt)

donde para cada t fijo cada sin(nx) se anula en n-1 puntos

o nodos.


u( x ,t )   an sen( nx ) cos( nt )
n 1

Pero recordemos que u(x,0) = f(x)...

16/11/2022 MSc. RAUL P. CASTRO VIDAL 20
UNMSM-FIEE
 Resolvamos por variables separadas: u(x,t) = X(x) T(t)

 u( x ,t )  u( x ,t )
2 2
 ; c.i . y c.c.
t 2
x 2

X ' ' ( x )  X ( x )  0 , x  ( 0 ,1 ), X ( 0 )  X ( 1 )  0
T ' ' ( t )  T ( t )  0, t  0.

Por eso Bernouilli optó por tomar f(x) como:


f ( x )  u( x ,0 )   an sen( nx )
n 1

con una adecuada elección de los coeficientes an...

16/11/2022 MSc. RAUL P. CASTRO VIDAL 21
UNMSM-FIEE
 Joseph Fourier
En diciembre de 1807 Joseph
Fourier presentó un sorprendente
artículo a la Academia de Ciencias
en París. En él afirmaba que
cualquier función puede escribirse
en forma de serie trigonométrica
semejante al ejemplo de Euler. Jean Baptiste Joseph Fourier
1768-1830

Polémica: Joseph-Louis Lagrange (1736-1813)

era uno de los muchos que opinaba que algo así
era simplemente imposible...
16/11/2022 MSc. RAUL P. CASTRO VIDAL 22
UNMSM-FIEE
Fourier fue nombrado por Napoleón secretario permanente del Instituto Egipcio.
Contrajo una enfermedad de Tiroides (mixedema).

16/11/2022 MSc. RAUL P. CASTRO VIDAL 23

UNMSM-FIEE
Fourier basó su trabajo en el estudio físico de la
ecuación del calor o de difusión:

 u 1 u
2

x 2
k t

Describe cómo el calor o una gota de tinta se

difunden en un medio.

Lord Kelvin (1736-1813): electricidad por los cables

trasatlánticos, edad de la Tierra,...
16/11/2022 MSc. RAUL P. CASTRO VIDAL 24
UNMSM-FIEE
 2u( x ,t ) 1 u( x ,t ) u( x ,t )  X ( x )T ( t )

x 2
k t X ( x )T ' ( t )  X ' ' ( x )T ( t )
u( 0 ,t )  u( ,t )  0; t  0
con X ( 0 )  X (  )  0
u( x ,0 )  f ( x ); 0  x  

Dividiendo entre X(x)T(t):

T' ( t ) X ' ' ( x )
 A , A  cte.
T( t ) X( x )
T ' ( t )  AT ( t ); T ( t )  C0 e At

X ' ' ( x )  AX ( x ); X ( x )  C1 cos(  A x )  C2 sen(  A x )

C1=0, C0=C2=1, A=-n2 con n = 1, 2, 3, ...

 n 2t
16/11/2022 u ( x ,t )  e
MSc. RAUL P. CASTRO n
VIDAL sen( nx ) 25
UNMSM-FIEE
  n 2t
La combinación lineal de soluciones u n ( x ,t )  e sen( nx )
será también solución:


u( x ,t )   anun ( x ,t )
n 1

Llegando al mismo resultado que Bernoulli, pero pudiendo calcular los

coeficientes an.

16/11/2022 MSc. RAUL P. CASTRO VIDAL 26

UNMSM-FIEE
 Serie trigonométrica de Fourier
Algunas funciones periódicas f(t) de periodo
T pueden expresarse por la siguiente serie,
llamada serie trigonométrica de Fourier
f (t )  12 a0  a1 cos(w0t )  a2 cos(2w0t )  a3 cos(3w0t )  ...
...  b1sen(w0t )  b2 sen(2w0t )  b3 sen(3w0t )  ...

Donde w0 = 2/T se denomina frecuencia

fundamental.

f (t )  12 a0   [an cos(nw0t )  bn sen(nw0t )]
n 1
16/11/2022 MSc. RAUL P. CASTRO VIDAL 27
UNMSM-FIEE
  t sen(2t ) sen(3t )
 sen t    ...
2 2 3

f (t )  12 a0   [an cos(nw0t )  bn sen(nw0t )]
n 1

a0 = 0, a1 = 0, a2 = 0 ...

b1 = 1, b2 = 1/2, b3 = 1/3,...
16/11/2022 MSc. RAUL P. CASTRO VIDAL 28
UNMSM-FIEE
¿Cómo calcular los coeficientes de la serie?
Dada una función periódica f(t), ¿cómo se
obtiene su serie de Fourier?

f(t)  12 a0   [an cos (nω0t)  bn sen(nω0t)]
n 1

Necesitamos calcular los coeficientes

a0,a1,a2,...,b1,b2,...

Lo haremos gracias a la ortogonalidad de

las funciones seno y coseno.
16/11/2022 MSc. RAUL P. CASTRO VIDAL 29
UNMSM-FIEE
 Ortogonalidad

Se dice que las funciones del conjunto {fk(t)}

son ortogonales en el intervalo a < t < b si
dos funciones cualesquiera fm(t), fn(t) de
dicho conjunto cumplen:

b
0 para m  n

a
f m (t)f n(t)dt  
 rn para m  n

16/11/2022 MSc. RAUL P. CASTRO VIDAL 30

UNMSM-FIEE
Ejemplo: las funciones t y t2 son ortogonales en el
intervalo –1 < t < 1, ya que:
1 1 4 1
t
1t t dt  1t dt  4 0
2 3

1

Ejemplo: Las funciones sen t y cos t son

ortogonales en el intervalo – < t <, ya que
π 2 π
sen t
π sent cos tdt  2 π
0

¿Falta algo para demostrar en ambos casos la ortogonalidad?

16/11/2022 MSc. RAUL P. CASTRO VIDAL 31
UNMSM-FIEE
 Ortogonalidad de senos y cosenos
Aunque los ejemplos anteriores se limitaron a un
par de funciones, el siguiente es un conjunto de
una infinidad de funciones ortogonales en el
intervalo -T/2< t < T/2:

{1, cos(w0t), cos(2w0t), cos(3w0t),...,

sen(w0t), sen2w0t, sen3w0t,...}

con w0= 2/T.

16/11/2022 MSc. RAUL P. CASTRO VIDAL 32

UNMSM-FIEE
Vamos a verificarlo probándolo a pares:

1.- f(t) = 1 vs. cos(mw0t): w0= 2/T

T/ 2 T/ 2
sen(mω0t)
 1
T/ 2
cos (mω0 t)dt 
mω0 T/ 2


2sen(mω0T/ 2 ) 2sen(mπ )
  0
mω0 mω0

Ya que m es un entero.

16/11/2022 MSc. RAUL P. CASTRO VIDAL 33

UNMSM-FIEE
2.- f(t) = 1 vs. sen(mw0t): w0= 2/T

 cos (mω0t)
T/ 2 T/ 2

 1 sen(mω0t)dt 
T/ 2
mω0 T/ 2


1
 [ cos (mω0T/ 2 )- cos (mω0T/ 2 )]  0
mω0

3.- cos(mw0t) vs. cos(nw0t): cos A cos B = ½[cos(A+B)+cos(A-B)]

cos2q = ½ (1+cos2q)

T /2
 0 para m  n
 cos(mw0 t)cos(nw0 t)dt  
T / 2 T / 2 para m  n  0
16/11/2022 MSc. RAUL P. CASTRO VIDAL 34
UNMSM-FIEE
 sen A sen B = ½[-cos(A+B)+cos(A-B)]
sen2 A =½ (1-cos2q)
4.- sen(mw0t) vs. sen(nw0t):

T/ 2
 0 para m  n
 sen(mω0t)sen(nω0t)dt  
T/ 2 T/ 2 para m  n  0

5.- sen(mw0t) vs. cos(nw0t):

sen A cos B = ½[sen(A+B)+sen(A-B)]
T/ 2

 sen(mω t) cos (nω t)dt  0

T/ 2
0 0 para cualquier m,n

16/11/2022 MSc. RAUL P. CASTRO VIDAL 35

UNMSM-FIEE
¿Cómo calcular los coeficientes de la serie?

Vamos a aprovechar la ortoganilidad que

acabamos de demostrar del conjunto de
funciones: {1, cos(w0t), cos(2w0t), cos(3w0t),...,
sen(w0t), sen2w0t, sen3w0t,...}
con w0= 2/T, en el intervalo -T/2< t < T/2 ,
para calcular los coeficientes a0,a1,a2,... ,
b1,b2,... de la serie de Fourier:

f(t)  a0   [an cos (nω0t)  bn sen(nω0t)]
1
2
n 1
16/11/2022 MSc. RAUL P. CASTRO VIDAL 36
UNMSM-FIEE
Multiplicando ambos miembros de la igualdad por
cos(mw0t) e integrando de –T/2 a T/2, obtenemos:
T /2 T /2 0, si m ≠ 0

T / 2
f (t ) cos(mw0t )dt  12 a0  cos (mω t)dt 
T / 2
0

 T /2 0, si m ≠ 0
 a  cos (nω t) cos (mω t)dt 
n 1
n 0 0 T/2, si m = n
T / 2
 T /2 0

 b  sen(nω t) cos (mω t)dt

n 1
n 0 0
T / 2
T /2
am  2
T  f (t ) cos( m w t ) dt
T / 2
0 m  1, 2, 3,...
16/11/2022 37
MSc. RAUL P. CASTRO VIDAL
UNMSM-FIEE
Observa que el caso anterior no incluye a a0, m = 0
que debemos tratar a parte:
T /2 T /2 T, si m = 0


T / 2
f (t ) cos(mw0t )dt  12 a0  cos (mω t)dt 
T / 2
0

 T /2 0, si m ≠ 0
 a  cos (nω t) cos (mω t)dt 
n 1
n 0 0
T/2, si m = n
T / 2
0
 T /2

 b  sen(nω t) cos (mω t)dt 

n 1
n 0 0
T / 2

1 T /2
2

a0T
2 a0  f (t ) dt
16/11/2022
T T / 2 38
MSc. RAUL P. CASTRO VIDAL
Similarmente, multiplicando por sen(mw0t) e
integrando de –T/2 a T/2, obtenemos:
T /2 T /2 0


T / 2
f (t ) sen(mω0t) dt  12 a0  sen(mω t)dt 
T / 2
0

0
 T /2

 a  cos (nω t)sen(mω t)dt 

n 1
n 0 0
T / 2
 T /2 0, si m ≠ 0
 b  sen(nω t)sen(mω t)dt
n 1
n 0 0
T/2, si m = n
T / 2
T /2
bm  2
T  f (t ) sen ( mw t ) dt
T / 2
0 m  1, 2, 3,...
16/11/2022 MSc. RAUL P. CASTRO VIDAL 39
UNMSM-FIEE
Un ejemplo históricamente importante:
Encontrar la serie de Fourier para la función
de onda cuadrada de periodo T:
f(t)
1

t
... -T/ 0 T/ T ...
2 2

-1

La expresión para f(t) en –T/2< t < T/2 es:

 1 para  T2  t  0
f (t )   w0= 2/T
 1 para 0  t  T
2
16/11/2022 MSc. RAUL P. CASTRO VIDAL 40
UNMSM-FIEE
Coeficiente a0:
T /2

 1 para   t  0
f (t )  
T
2
a0  T1  f (t )dt
T / 2
 1 para 0  t  T
2

 0 T /2
  0 T /2

a0  T2    dt   dt   T2  t t 0
T / 2 0   T / 2 0 

16/11/2022 41
MSc. RAUL P. CASTRO VIDAL
UNMSM-FIEE
Coeficientes an:

T /2
 1 para   t  0
T
f (t )   2

  T
an  T2  f (t ) cos(nw t )dt
0
 1 para 0 t 2 T / 2

 0 T /2

an  T    1 cos(nw0t )dt   1 cos(nw0t )dt 
2

T / 2 0 
 1 0
1 T /2

 T 
2
sen(nw0t )  sen(nw0t )   0
 nw0 T / 2 nw0 0  
para n  0

16/11/2022 42
MSc. RAUL P. CASTRO VIDAL
UNMSM-FIEE
Coeficientes bn:
T /2
 1 para  T2  t  0
f (t )  
bn  T2  f (t )sen(nw t )dt
0

 1 para 0  t  T
2
T / 2

 0 T /2

bn  T    sen(nw0t )dt   sen(nw0t )dt  
2

T / 2 0 
 1 0
1 T /2

T
2
cos(nw0t )  cos(nw0t ) 
 nw0 T / 2 nw 0 0  


1
(1  cos(n ))  (cos(n )  1)
n

2
n

1  (1) n ) para n  0
43
16/11/2022 MSc. RAUL P. CASTRO VIDAL
UNMSM-FIEE
Finalmente, la serie de Fourier queda como

f (t ) 
4
sen(w0t )  13 sen(3w0t )  15 sen(5w0t )  ...


sen((2n  1)w0t ) )
4 1
f (t )  
 n 1 2n  1
En la siguiente figura se muestran: la
componente fundamental y los armónicos 3,
5 y 7, así como la suma parcial de estos
primeros cuatro términos de la serie para
w0 =  (w0= 2/T), es decir, T = 2:
44
16/11/2022 MSc. RAUL P. CASTRO VIDAL
UNMSM-FIEE
 f (t ) 
4
sen(w0t )  13 sen(3w0t )  15 sen(5w0t )  ...

Componentes de la Serie de Fourier
1.5

1
Componentes

0.5

-0.5
Suma
fundamental
-1 tercer armónico
quinto armónico
séptimo armónico
-1.5
-1 -0.5 0 t 0.5 1
16/11/2022 Fourier series java applet (http://www.falstad.com/fourier/) 45
MSc. RAUL P. CASTRO VIDAL
Nota:
Para expresarse como serie de Fourier
f(t), no necesita estar centrada en el origen.
Simplemente debemos tomar el intervalo,
donde está definida, como el periodo de la
serie.
La ortogonalidad de las funciones seno y
coseno no sólo se da en el intervalo de
–T/2 a T/2, sino en cualquier intervalo que
cubra un periodo completo:
de t0 a t0 + T, con t0 arbitrario,
con el mismo resultado.
16/11/2022 MSc. RAUL P. CASTRO VIDAL 46
UNMSM-FIEE
 Habíamos calculado f(t)
1
los coeficientes para:
t
... -T/ 0 T/ T ...
 1 para  T / 2  t  0 2 2
f (t )  
1 para 0  t  T / 2 -1

Si los calculamos para la misma función desplazada

tienen que ser los mismos: f(t)
1
 1 para 0  t  T / 2
f (t )   t
 1 para T / 2  t  T
... -T/ 0 T/ T ...
2 2

-1
Repite los cálculos y compruébalo.
47
16/11/2022 MSc. RAUL P. CASTRO VIDAL
UNMSM-FIEE
 f(t)
De hecho si repetimos 1
para cualquier intervalo
t
de longitud el periodo
T de la función, será lo
-1
mismo: ... t0 t0 +T ...

T /2 T t 0 T

a0  T1 
T / 2
f (t )dt  T2  f (t )dt  T2
0

t0
f (t )dt  T2  f (t )dt
T

T /2
an  T2 
T / 2
f (t ) cos(nw0t )dt  ...  T2  f (t ) cos(nw0t )dt
T
T /2
bn  T2
16/11/2022

T / 2
f (t ) sen(nw0t )dt  ...  T2  f (t ) sen(nw0t )dt
MSc. RAUL P. CASTRO VIDAL
T 48
UNMSM-FIEE
Ejercicio: encontrar la serie de Fourier para

 t
f (t ) 
2
la función con la que empezamos el tema.
O sea, demostrar que Euler tenía razón.

16/11/2022 MSc. RAUL P. CASTRO VIDAL 49

UNMSM-FIEE
 Calcula la serie de Fourier de la función periódica:
2
f (t )  1  cos(3t ) de periodo T 
3
2
3
2 3
a0   f (t )dt   (1  cos(3t ))dt  2
T T  0
2

2 3 3
1, si n  1
an   f (t ) cos(nw0t )dt  0 (1  cos(3t )) cos(nw0t )dt  0, si n  1
T T 
2
3
2 3
bn   f (t ) sen(nw0t )dt   (1  cos(3t ))sen(nw t )dt  0 para todo n

0
T T 0

en definitiva
 
La serie
f (t )  1   an cos(nw0t )   bn sen(nw0t )  1  cos(3t ) es la
propia
n 1 n 1
función...
16/11/2022 MSc. RAUL P. CASTRO VIDAL 50
UNMSM-FIEE
Nota: a partir de ahora entenderemos que f(t) está definida
sólo en el intervalo que especifiquemos. Y que la serie de
Fourier la extiende periódicamente, con periodo T igual al
intervalo de definición. En muchos libros se habla de
extender de forma par o impar una función. La serie de
Fourier extenderá periódicamente los patrones siguientes:

t
Extensión par

t
16/11/2022 Extensión impar 51
MSc. RAUL P. CASTRO VIDAL
UNMSM-FIEE
 Funciones Pares e Impares

Una función (periódica o no) se dice

función par (o con simetría par) si su
gráfica es simétrica respecto al eje vertical,
es decir, la función f(t) es par si
f(t) = f(-t)
f(t)

t
2   2

16/11/2022 MSc. RAUL P. CASTRO VIDAL 52

UNMSM-FIEE
En forma similar, una función f(t) se dice
función impar (o con simetría impar), si su
gráfica es simétrica respecto al origen, es
decir, si cumple lo siguiente: -f(t) = f(-t)

f(t)

t
2   2

16/11/2022 MSc. RAUL P. CASTRO VIDAL 53

UNMSM-FIEE
Ejemplo: ¿Las siguientes funciones son
pares o impares?
f(t) = t + 1/t ,
g(t) = 1/(t2+1).

Solución:
Como f(-t) = -t - 1/t = - f(t), por lo tanto f(t) es
función impar.
Como g(-t) = 1/((-t)2+1) = 1/(t2+1) = g(t), por
lo tanto g(t) es función par.

16/11/2022 MSc. RAUL P. CASTRO VIDAL 54

UNMSM-FIEE
Ejemplo: ¿La función h(t) = f(1+t2) es par o
impar? (f es una función arbitraria).

Solución:
Sea g(t) = 1 + t2. Entonces h(t) = f(g(t)).
Por lo tanto h(-t) = f(g(-t)).
Pero g(-t) = 1+(-t)2 = 1 + t2 = g(t),
finalmente h(-t) = f(g(t)) = h(t), de modo que
h(t) es función par, sin importar como sea
f(t).

16/11/2022 MSc. RAUL P. CASTRO VIDAL 55

UNMSM-FIEE
Ejemplo: De acuerdo al ejemplo anterior,
todas las funciones siguientes son pares:

h(t) = sen (1+t2)

h(t) = exp(1+t2) + 5/ (1+t2)
h(t) = cos (2+t2) + 1
h(t) = (10+t2) - (1+t2)1/2
etc...
Ya que todas tienen la forma f(1+t2).

16/11/2022 MSc. RAUL P. CASTRO VIDAL 56

UNMSM-FIEE
 • Si f (x) es par:
a a

 f ( x)dx
a
 2 f ( x)dx
0

a a

 f ( x)dx  f ( x)dx
0
a

-a a

16/11/2022 MSc. RAUL P. CASTRO VIDAL 57

UNMSM-FIEE
 • Si f (x) es impar:
a

 f ( x)dx  0
a

 f ( x)dx
a

-a a

16/11/2022 MSc. RAUL P. CASTRO VIDAL 58

UNMSM-FIEE
 Como la función sen(nw0t) es una función
impar para todo n y la función cos(nw0t) es
una función par para todo n, es de esperar
que:

• Si f(t) es par, su serie de Fourier no

contendrá términos seno, por lo tanto
bn= 0 para todo n.

• Si f(t) es impar, su serie de Fourier no

contendrá términos coseno, por lo tanto
an= 0 para todo n.
16/11/2022 MSc. RAUL P. CASTRO VIDAL 59
UNMSM-FIEE
Por ejemplo, la señal cuadrada, que hemos
analizado:
f(t)
1

t
... -T/ 0 T/ T ...
2 2

-1

Es una función impar, por ello su serie de

Fourier no contiene términos coseno:
f (t )  sen(w0t )  13 sen(3w0t )  15 sen(5w0t )  ...
4

16/11/2022 MSc. RAUL P. CASTRO VIDAL 60
UNMSM-FIEE
P2. Septiembre 2005

a) Obtener el desarrollo en serie de Fourier de las funciones

f ( x)  sin x y g ( x)  cos x en    x  

Respuesta.


f ( x)    an cos(nx)  bn sin( nx)
a0
2 n 1

f(x) = |sen(x)|, x є [-π,π], 2π periódica

Función par → desarrollo en cosenos, bn = 0

16/11/2022 MSc. RAUL P. CASTRO VIDAL 61
UNMSM-FIEE
 1  2 
an  
 
f ( x) cos(nx)dx  
  0
sin x cos(nx)dx 

sin(1  n) x  sin(1  n) xdx 

1 

  0


1 2
cos(n  1)  1
 n 1
2

4 4
a0  ; a n  , n par; an  0, n impar
  (n  1)
2

16/11/2022 MSc. RAUL P. CASTRO VIDAL 62

UNMSM-FIEE
 
2 4 cos(2nx)
sin x   
 n 1  4n  1
2

f(x) = |cos(x)|, x є [-π,π], 2π periódica

Función par → desarrollo en cosenos, bn = 0

1  4  /2
an  
 
g ( x) cos( nx )dx  
  0
cos x cos( nx )dx 

cos(n  1) x  cos(n  1) xdx

2  /2

  0

16/11/2022 MSc. RAUL P. CASTRO VIDAL 63

UNMSM-FIEE
 4 4
a0  ; an  , n par; a  0, n impar
  (n  1)
2 n

2 4 (1) cos(2nx)
 n
cos x   
 n 1  4n  1
2

16/11/2022 MSc. RAUL P. CASTRO VIDAL 64

UNMSM-FIEE
 Onda triangular
(Triangle Wave)

 4  cos x cos 3 x cos 5 x 

  2    ⋯
2  1 3 2
5 2


16/11/2022 MSc. RAUL P. CASTRO VIDAL 65

UNMSM-FIEE
 Right Triangular Wave

 sin x sin 2 x sin 3 x 

2    ⋯
 1 2 3 
16/11/2022 MSc. RAUL P. CASTRO VIDAL 66
UNMSM-FIEE
 Saw Tooth Wave

 sin x sin 2 x sin 3 x 

  2    ⋯
 1 2 3 
16/11/2022 MSc. RAUL P. CASTRO VIDAL 67
UNMSM-FIEE
Ejercicio: demostrar que la serie de Fourier para

f (t )  cos(t ),    t  
con periodo T = 2π (frecuencia fundamental
w0 = 1) y  un número real no entero, es:

sen(  )  1 (1)  n

cos( t )    2  2 cos( n t ) 
  n 1   n
2


16/11/2022 MSc. RAUL P. CASTRO VIDAL 68

UNMSM-FIEE
 sen(  )  1 (1) 
 n
cos( t )    2  2 cos(n t ) 
  n 1   n
2

  t  
Observa que si tomamos t = 0 entonces:

 1 (1)  n
  2  2
sen(  )  n 1   n
2

y con  = 1/2.

(1) n
(1) n
  2  2  4
n 1 (1 / 2)  n n 1 1  4n
2 2 2

16/11/2022 MSc. RAUL P. CASTRO VIDAL 69

UNMSM-FIEE
 sen(  )  1 (1) 
 n
cos( t )    2  2 cos(n t ) 
  n 1   n
2

  t  
O que si tomamos t = π entonces: cos( t )  (1) n

sen(  )  1 
1 
cos(  )    2  2 2 
  n 1   n 

 1 
1
  2  2
tan(  )  n 1   n
2

¿Es correcto el resultado?

16/11/2022 MSc. RAUL P. CASTRO VIDAL 70
UNMSM-FIEE
 Convergencia uniforme
Que la integral traspase los sumatorios en la
deducción de las fórmulas para los coeficientes
de la serie de Fourier, equivale a asumir que la
serie converge uniformemente... Recordemos
qué es convergencia uniforme.

Sea la serie infinita: S ( x)   un ( x)
n 1

y definamos sus sumas parciales como:

k
S k ( x)   un ( x)
16/11/2022 MSc. RAUL P. CASTRO VIDAL n 1 71
UNMSM-FIEE
Diremos que S converge a f(x) en un intervalo si
 > 0 existe para todo x del intervalo un N > 0 tq.:

S k ( x)  f ( x)   siempre que k  N
Observemos que en general N dependerá de  y
del punto x (convergencia puntual).
Si N solo depende de , pero no de x, decimos que
la convergencia es uniforme.

Que la serie sea uniformemente convergente es

"bueno" porque:

16/11/2022 MSc. RAUL P. CASTRO VIDAL 72

UNMSM-FIEE
(1) Si cada término un(x) de una serie es
continuo en (a, b) y la serie es uniformemente
convergente a f(x), entonces:

(a) f(x) es también continua en (a, b).

b  

  u ( x)dx    u ( x)dx
b
(b)
n n
a a
n 1 n 1

(2) Si cada término un(x) de una serie posee

derivada en (a, b) y la serie derivada es
uniformemente convergente, entonces:
d  
d
16/11/2022

dx n 1
un ( x)   un ( x)
n 1 dx 73

MSc. RAUL P. CASTRO VIDAL

¿Cómo probar la convergencia uniforme de una serie?

(1) Encontrar una expresión "cerrada" para Sk(x) y

aplicar la definición o

(2) utilizar la prueba M de Weierstrass:

Si existe {Mn}n = 1, 2,... tq. |un(x)|  Mn y además

 

M
n 1
n converge   u ( x) converge
n 1
n uniformemente

16/11/2022 MSc. RAUL P. CASTRO VIDAL 74

UNMSM-FIEE
 Ejemplo:

sen(nx)
S ( x)   2
en ( ,  )
n 1 n
1 sen(nx) 1
Mn  2  2
 2
n n n

1  2


n 1 n
2

6
 S converge uniformemente

16/11/2022 MSc. RAUL P. CASTRO VIDAL 75

UNMSM-FIEE
Condiciones de Dirichlet

Condiciones de convergencia de la serie de Fourier

de f(x), suficientes pero no necesarias.

(1) f(x) tiene un número finito de discontinuidades

en un periodo.

(2) f(x) tiene un número finito de máximos y mínimos

en un periodo.

(3)

T
f ( x) dx  
16/11/2022 MSc. RAUL P. CASTRO VIDAL 76
UNMSM-FIEE
Si se cumplen las condiciones de Dirichlet, entonces
la serie de Fourier converge a f(x) si x es un punto
de continuidad y a:
1
2
(  
f (x )  f (x ) )
si x es un punto de discontinuidad.

16/11/2022 MSc. RAUL P. CASTRO VIDAL 77

UNMSM-FIEE
 Desarrolla 0,   x  0
en serie de Fourier: f ( x)  
  x, 0 x 

T  2
2 
a0  
2 
f ( x) dx

1 0 


 
 
0 dx  0
(  x ) dx


1 x  2
  x   
 2 0 2

16/11/2022 MSc. RAUL P. CASTRO VIDAL 78

UNMSM-FIEE
  1 0 
1 
an 
 

f ( x) cos nx dx 
  
0 dx  0 (  x ) cos nx dx

 
1  sin nx 1   1 cos nx
 (  x)   sin nx dx   
  n 0 n 0
 n n 0

 cos n  1 1  (1) n
 
n
2
n 2
1  1
 0
bn  (  x) sin nxdx 
n

 1  (1) n

1 
f ( x)     2 cos nx  sin nx 
4 n1  n  n 
16/11/2022 MSc. RAUL P. CASTRO VIDAL 79
UNMSM-FIEE
La función f es continua en (−, ) excepto en x = 0. Así su
serie de Fourier converge en x = 0 a:
f ( 0  )  f ( 0 )   0 
 
2 2 2

La serie es una extensión periódica de la

función f. Las discontinuidades en x = 0,
f (0  )  f (0 ) 
2, 4, … convergen a: 
2 2

16/11/2022 MSc. RAUL P. CASTRO VIDAL 80

UNMSM-FIEE
Secuencia de sumas parciales y su representación gráfica

  2  2 1
S1  , S 2   cos x  sin x, S3   cos x  sin x  sin 2 x
4 4  4  2

16/11/2022 81
MSc. RAUL P. CASTRO VIDAL
16/11/2022 MSc. RAUL P. CASTRO VIDAL 82
UNMSM-FIEE
16/11/2022 83
MSc. RAUL P. CASTRO VIDAL
UNMSM-FIEE
16/11/2022 MSc. RAUL P. CASTRO VIDAL 84
UNMSM-FIEE
16/11/2022 MSc. RAUL P. CASTRO VIDAL 85
UNMSM-FIEE
16/11/2022 MSc. RAUL P. CASTRO VIDAL 86
UNMSM-FIEE
16/11/2022 MSc. RAUL P. CASTRO VIDAL 87
UNMSM-FIEE
16/11/2022 MSc. RAUL P. CASTRO VIDAL 88
UNMSM-FIEE
16/11/2022 MSc. RAUL P. CASTRO VIDAL 89
UNMSM-FIEE
16/11/2022 MSc. RAUL P. CASTRO VIDAL 90
UNMSM-FIEE
16/11/2022 MSc. RAUL P. CASTRO VIDAL 91
UNMSM-FIEE
16/11/2022 MSc. RAUL P. CASTRO VIDAL 92
UNMSM-FIEE
16/11/2022 MSc. RAUL P. CASTRO VIDAL 93
UNMSM-FIEE
16/11/2022 MSc. RAUL P. CASTRO VIDAL 94
UNMSM-FIEE
16/11/2022 MSc. RAUL P. CASTRO VIDAL 95
UNMSM-FIEE
16/11/2022 MSc. RAUL P. CASTRO VIDAL 96
UNMSM-FIEE
16/11/2022 MSc. RAUL P. CASTRO VIDAL 97
UNMSM-FIEE
16/11/2022 MSc. RAUL P. CASTRO VIDAL 98
UNMSM-FIEE
16/11/2022 MSc. RAUL P. CASTRO VIDAL 99
UNMSM-FIEE
16/11/2022 MSc. RAUL P. CASTRO VIDAL 100
UNMSM-FIEE
16/11/2022 MSc. RAUL P. CASTRO VIDAL 101
UNMSM-FIEE
Ejercicio de examen: Obtener el desarrollo en serie de Fourier de la
función

f (t )  1  t , t  0,1
2

de modo que converja uniformemente a f(t) en [0,1].

Respuesta.

Para que el desarrollo de Fourier se pueda definir debe ser 2L-

periódica.

Para que converja uniformemente, se debe extender f(t) a

~ de
f (t )
~
modo que: 1. f (t ) sea continua en [-L,L].
~
2. f (t ) sea continua a trozos en [-L,L].

16/11/2022 MSc. RAUL P. CASTRO VIDAL 102

UNMSM-FIEE
 La continuidad se consigue con la extensión par de f (f´ = -2t es
continua en [-L,L] ) con L = 1.

Im (z)

-1 1 Re (z)

~ a0  
 n   n  
f (t )     an cos t   bn sin  t  
2 n 1   L   L 
bn  0 por ser función par
16/11/2022 MSc. RAUL P. CASTRO VIDAL 103
UNMSM-FIEE
 1 1
an   (1  t ) cos(nt )dt  2 (1  t ) cos(nt )dt 
2 2
1 ~ 0
f par

4(1) n

(n ) 2

1 1 2 4
a0   (1  t )dt  2 (1  t )dt  2 
2 2
1 0 3 3

~ (1)
 n
cos(nt )
2 4
f (t )   2
3 

n 1 n
2

~
f (t )  f (t ) 
t0 ,1
16/11/2022 104
MSc. RAUL P. CASTRO VIDAL
UNMSM-FIEE
P2. Septiembre 2018

a) (4 puntos)

1. Obtener el desarrollo en serie de Fourier de la función

f(x) = x2 -π ≤ x ≤ π, con f(x) = f(x + 2π)

2. Estudiar si el desarrollo obtenido converge uniformemente a

f(x) en [-π,π]

3. Basándose en los resultados obtenidos, calcular la suma de la


1
serie numérica 
k 1 k 4

4. A partir del desarrollo de Fourier de la función f(x), obtener el

desarrollo en serie de Fourier de la función

g(x) = x(x2 – π2) -π ≤ x ≤ π, con g(x) = g(x + 2π)

16/11/2022 MSc. RAUL P. CASTRO VIDAL 105

UNMSM-FIEE
Respuesta.

1. f(x) = x2, x є [-π,π], 2π periódica

Función par → desarrollo en cosenos, bn = 0:


a0
f ( x)    an cos(nx)
2 n 1
2  2 2 2
a0   x dx  
 0 3
1  2 2  2
an   x cos(nx)dx   x cos(nx)dx 
   0

16/11/2022 MSc. RAUL P. CASTRO VIDAL 106

UNMSM-FIEE
   

2 1 2 2 2 
  x sin(nx)  2 x cos(nx)  3 sin(nx)  
  n 0 n 0 n 0
2 2 4
 ( 1) n
a  ( 1) n

 n 2 n
n 2

 2
(1)
 n
f ( x)   4 2 cos(nx)
3 n 1 n

16/11/2022 MSc. RAUL P. CASTRO VIDAL 107

UNMSM-FIEE
2.

f continua en -  ,   
hay convergencia uniforme
f  continua en (-  ,  )

3. Por convergencia uniforme, se aplica la identidad de Parseval:

( )

2 
1
  f ( x )  2
dx 
a
  a 20
 b 2


n n
  2 n 1

 1 5 2 5
 ( x ) dx  5 x  
2 2

 5
16/11/2022 MSc. RAUL P. CASTRO VIDAL 108
UNMSM-FIEE
 2 
2 4 2 2 1 1
     16  4
5 3  2 n 1 n


1  2


n 1 n 4

90

4.
( 2 2
)
g ( x)  x x   , x    ,  , 2 periódica
( 1) n 
g ( x)  3x    3 f ( x)    12 2 cos(nx)
 2 2

n1 n

( 1) n 
Por convergenc ia uniforme : g ( x )  12 3 sin( nx )
n 1 n
16/11/2022 109
MSc. RAUL P. CASTRO VIDAL
UNMSM-FIEE
 Fenómeno de Gibbs
Si la serie de Fourier para una función f(t) se
trunca para lograr una aproximación en suma
finita de senos y cosenos, es natural pensar que a
medida que agreguemos más armónicos, el
sumatorio se aproximará más a f(t).
Esto se cumple excepto en las discontinuidades
de f(t), en donde el error de la suma finita no
tiende a cero a medida que agregamos armónicos.
Por ejemplo, consideremos el tren de pulsos u
onda cuadrada:

f (t ) 
4
sen(w0t )  13 sen(3w0t )  15 sen(5w0t )  ...
16/11/2022  MSc. RAUL P. CASTRO VIDAL 110
UNMSM-FIEE
 f (t ) 
4
sen(w0t )

1.5
Serie con 1 arm ónico

0.5

-0.5

-1

-1.5
-1 -0.5 0 0.5 1
16/11/2022 MSc. RAUL P. CASTRO VIDAL 111
UNMSM-FIEE
 f (t ) 
4
sen(w0t )  13 sen(3w0t )  15 sen(5w0t )

1.5
Serie con 3 arm ónicos

0.5

-0.5

-1

-1.5
-1 -0.5 0 0.5 1
16/11/2022 MSc. RAUL P. CASTRO VIDAL 112
UNMSM-FIEE
 1.5
Serie con 5 arm ónicos

0.5

-0.5

-1

-1.5
-1 -0.5 0 0.5 1
16/11/2022 MSc. RAUL P. CASTRO VIDAL 113
UNMSM-FIEE
 1.5
Serie con 7 arm ónicos

0.5

-0.5

-1

-1.5
-1 -0.5 0 0.5 1
16/11/2022 MSc. RAUL P. CASTRO VIDAL 114
UNMSM-FIEE
 1.5
Serie con 13 arm ónicos

0.5

-0.5

-1

-1.5
-1 -0.5 0 0.5 1
16/11/2022 MSc. RAUL P. CASTRO VIDAL 115
UNMSM-FIEE
 Fenómeno de Gibbs

1.5
Serie con 50 arm ónicos

0.5

-0.5

-1

-1.5
-1 -0.5 0 0.5 1
16/11/2022 MSc. RAUL P. CASTRO VIDAL 116
UNMSM-FIEE
 Fenómeno de Gibbs

1.5
Serie con 100 arm ónicos

0.5

-0.5

-1

-1.5
-1 -0.5 0 0.5 1
16/11/2022 MSc. RAUL P. CASTRO VIDAL 117
UNMSM-FIEE
16/11/2022 MSc. RAUL P. CASTRO VIDAL 118
UNMSM-FIEE
Forma compleja de la serie de Fourier

Consideremos la serie de Fourier para una

función periódica f(t), con periodo T = 2/w0.

f (t )  12 a0   [an cos(nw0t )  bn sen(nw0t )]
n 1

Es posible obtener una forma alternativa

usando las fórmulas de Euler:

inw0t  inw0t
cos(nw0t )  (e 1
2 e )
inw0t inw0t
sen(nw0t )  1
2i (e e )
16/11/2022 MSc. RAUL P. CASTRO VIDAL 119
UNMSM-FIEE
Sustituyendo:

f (t )  12 a0   [an 12 (e inw0t  e inw0t )  bn 1
2i (e inw0t  e inw0t )]
n 1

Y usando el hecho de que 1/i = -i:


f (t )  a0   [ (an  ibn )e
1
2
1
2
inw0t
 (an  ibn )e
1
2
inw0t
]
n 1

Y definiendo: c0  12 a0, cn  12 (an ibn), cn  12 (an ibn)


f (t )  c e
n  
n
inw0t
w0 
2
T
16/11/2022 MSc. RAUL P. CASTRO VIDAL 120
UNMSM-FIEE
A la expresión obtenida 
f (t )  c e
n  
n
inw0t

se le llama forma compleja de la serie de

Fourier y sus coeficientes cn pueden obtenerse a
partir de los coeficientes an, bn como ya se dijo, o
bien:
T
 inw0t 


inw0t ¿Forma e
cn  1
T f (t )e dt n  
un conjunto
0
ortogonal?
Para n = 0, 1, 2, 3, ...
Demostrarlo.
16/11/2022 MSc. RAUL P. CASTRO VIDAL 121
UNMSM-FIEE
Ejemplo. Encontrar la forma compleja de la
serie de Fourier para la función ya tratada:
f(t)
1

t
... -T/ 0 T/ T ...
2 2

-1
Solución 1. Como ya se calcularon los
coeficientes de la forma trigonométrica (an y
bn), que eran an= 0 para todo n y
2
bn  [1  (1) ]
n
para todo n
n
16/11/2022 MSc. RAUL P. CASTRO VIDAL 122
UNMSM-FIEE
Podemos calcular los coeficientes cn:

cn  [an  ibn ]  i
1
2
1 2
2 n [1  (1) ] n

cn  i 1
n [1  (1) ] n

Entonces la serie compleja de Fourier

queda:
 i 5w0t  i 3w0 t  iw 0 t
f (t )   i (...  e
2 1
5  e
1
3 e
iw 0 t i 3w0t i 5w0t
e  e 1
3  e1
5  ...)
16/11/2022 123
MSc. RAUL P. CASTRO VIDAL
UNMSM-FIEE
Solución 2. También podemos calcular los
coeficientes cn mediante la integral:
T

 f (t )e
inw0t
cn  1
T
dt
0

1  T /2 T

   e inw0t
dt    e inw0t
dt 
T 0 T /2 
1  1 inw0t 
T /2 T
 e  in1wo e  inw0 t 
T  o
 inw
0

T /2 


1
 inwoT
(e 
inw0T / 2
 1)  (e inw0T
e inw0T / 2
) 
16/11/2022 MSc. RAUL P. CASTRO VIDAL 124
UNMSM-FIEE
Como w0T = 2 y además:
 iq
e  cos q  isenq
cn  1
 inwoT [(1) n  1)  (1  (1) n )]

 i nwoT [1  (1) ]
2 n

 i n1 [1  (1) n ]

que coincide con el resultado ya obtenido.

16/11/2022 125
MSc. RAUL P. CASTRO VIDAL
UNMSM-FIEE
Calcular la serie de Fourier de la función de Heaviside:

0 ,  1  x  0 
H ( x)   H (x )  n
c e inx

 1, 0  x  1 n  

1
1 1 inx 
1 1
1 inx 1 inx
cn   e H ( x)dx   e dx   e 
2 1 20 2   in 0
cn  e  1  cos(n )  isen(n )  1 
1 i in i
2 n 2n


0 ; si n es par
i
cos(n )  1    i ; si n es impar n ≠ 0
2 n  n
16/11/2022 MSc. RAUL P. CASTRO VIDAL 126
UNMSM-FIEE
  
 i inx 1
H (x ) 
1
n
c e inx
  
2 0 n   n
  i inx 
e    2 Re  e 
n   2 n0  n 
n impar n impar

1 -ii0x
1
1

1 0 ; si n es par
1
al0   e c0 H(x)dx
1
 ; cn  dx
lx
i 
2 1 2 2 0 ; si n2es impar
 n
1 2   i (cos(nx)  isen(nx) ) 
H ( x)    Re  
2  n 0  n
n impar
sen(nx)
H (x )   
1 2
2  n 0 n
16/11/2022 MSc. RAUL P. CASTRO VIDAL 127
UNMSM-FIEE n impar
16/11/2022 MSc. RAUL P. CASTRO VIDAL 128
UNMSM-FIEE
16/11/2022 MSc. RAUL P. CASTRO VIDAL 129
UNMSM-FIEE
16/11/2022 MSc. RAUL P. CASTRO VIDAL 130
UNMSM-FIEE
La función impulso o (t)
delta de Dirac
 if t  0
 (t )  
 0 if t  0 t

Se trata de una "función generalizada". Podemos pensar en

la delta de Dirac como el límite de una serie de funciones:
m  (mt ) 2
f m(t)  e
(t) 
f3(t)
f2(t)
f1(t)

16/11/2022 MSc. RAUL P. CASTRO VIDAL t 131

UNMSM-FIEE
 (t)
Propiedades de la función 
t


  (t ) dt  1

 

  (t  a ) f (t ) dt    (t  a ) f ( a ) dt 
 
f (a )

 exp(  iw t ) dt  2  (w )



 exp[  i (w  w ')t ]

16/11/2022
dt  2  (w  w ')
MSc. RAUL P. CASTRO VIDAL 132
UNMSM-FIEE
Calcular la serie de Fourier de (x):
 1
 (x )  c n
inx
e
1 inx
 cn   e  ( x)dx 
1
n   2 1 2

 ( x )   e    (e  e )
1 inx 1 1 inx inx

2 n   2 2 n 0
1
   cos(nx)
2 n 0 1
 (x )    cos(nx)
2 n 0
16/11/2022 MSc. RAUL P. CASTRO VIDAL 133
UNMSM-FIEE
 12

10

-2

-4
-4 -3 -2 -1 0 1 2 3 4
x

16/11/2022 MSc. RAUL P. CASTRO VIDAL 134

UNMSM-FIEE
 12

10

-2

-4
-4 -3 -2 -1 0 1 2 3 4
x

16/11/2022 MSc. RAUL P. CASTRO VIDAL 135

UNMSM-FIEE
 12

10

-2

-4
-4 -3 -2 -1 0 1 2 3 4
x

16/11/2022 MSc. RAUL P. CASTRO VIDAL 136

UNMSM-FIEE
 12

10

-2

-4
-4 -3 -2 -1 0 1 2 3 4
x

16/11/2022 MSc. RAUL P. CASTRO VIDAL 137

UNMSM-FIEE
 12

10

-2

-4
-4 -3 -2 -1 0 1 2 3 4
x

16/11/2022 MSc. RAUL P. CASTRO VIDAL 138

UNMSM-FIEE
 12

10

-2

-4
-4 -3 -2 -1 0 1 2 3 4
x

16/11/2022 MSc. RAUL P. CASTRO VIDAL 139

UNMSM-FIEE
 12

10

-2

-4
-4 -3 -2 -1 0 1 2 3 4
x

16/11/2022 MSc. RAUL P. CASTRO VIDAL 140

UNMSM-FIEE
 12

10

-2

-4
-4 -3 -2 -1 0 1 2 3 4
x

16/11/2022 MSc. RAUL P. CASTRO VIDAL 141

UNMSM-FIEE
 12

10

-2

-4
-4 -3 -2 -1 0 1 2 3 4
x

16/11/2022 MSc. RAUL P. CASTRO VIDAL 142

UNMSM-FIEE
 12

10

-2

-4
-4 -3 -2 -1 0 1 2 3 4
x

16/11/2022 MSc. RAUL P. CASTRO VIDAL 143

UNMSM-FIEE
Los coeficientes cn son números complejos,
y también se pueden escribir en forma polar:
i n
cn  cn e
Observemos que,
 i n
c n  c  cn e
*
n

 bn 
Donde cn  a  b , n  arctan  
1 2 2
2 n n
 an 
para todo n  0.
Y para n = 0, c0 es un número real: c0  a0
1
2
16/11/2022 MSc. RAUL P. CASTRO VIDAL 144
UNMSM-FIEE
 Espectros de frecuencia discreta

Dada una función periódica f(t), le

corresponde una y sólo una serie de
Fourier, es decir, le corresponde un conjunto
único de coeficientes cn.

Por ello, los coeficientes cn especifican a f(t)

en el dominio de la frecuencia de la
misma manera que f(t) especifica la función
en el dominio del tiempo.

16/11/2022 MSc. RAUL P. CASTRO VIDAL 145

UNMSM-FIEE
 Espectros de frecuencia discreta

Ejemplo. Para la función ya analizada:

f(t)
1

t
... -T/ 0 T/ T ...
2 2

-1

Encontramos que: cn  i n1 [1  (1) n ]

1
Por lo tanto: cn  [1  (1) ]
n

n

16/11/2022 MSc. RAUL P. CASTRO VIDAL 146

UNMSM-FIEE
A la gráfica de la magnitud de los
coeficientes cn contra la frecuencia angular
w de la componente correspondiente se le
llama el espectro de amplitud de f(t).

A la gráfica del ángulo de fase n de los

coeficientes cn contra w, se le llama el
espectro de fase de f(t).

Como n sólo toma valores enteros, la

frecuencia angular w = nw0 es una variable
discreta y los espectros mencionados son
gráficas discretas.
16/11/2022 MSc. RAUL P. CASTRO VIDAL 147
UNMSM-FIEE
El espectro de amplitud se muestra a continuación
0.7
Espectro de Amplitud de f(t)
0.6
0.5
Cn 

0.4
0.3
0.2
0.1
0
-30 -20 -10 0 n 10 20 30

Frecuencia negativa (?) Frecuencia

Observación: El eje horizontal es un eje de
frecuencia, (n = número de armónico = múltiplo
de w0).
16/11/2022 MSc. RAUL P. CASTRO VIDAL 148
UNMSM-FIEE
El espectro de magnitud de una f(t) real, es una
función PAR por lo que la gráfica para n  0
contiene toda la información acerca de f(t) y se le
conoce como espectro unilateral de magnitud.

El espectro de fase de una f(t) real, es una

función IMPAR por lo que la gráfica para n  0
contiene toda la información acerca de f(t) y se le
conoce como espectro unilateral de fase.

16/11/2022 MSc. RAUL P. CASTRO VIDAL 149

UNMSM-FIEE
Podemos expresar de una manera ligeramente
diferente la serie de Fourier. Cada par de
términos:
ancos(nw0t) + bnsen(nw0t)

se pueden expresar como:

 
2
sen( nw0t ) 
a b
an  bn
2 n
cos( nw0t )  n
 a2  b2 a 2
 b 2 
 n n n n 

Donde lo único que hemos hecho es multiplicar

y dividir por: a 2
n  bn
2

16/11/2022 150
MSc. RAUL P. CASTRO VIDAL
UNMSM-FIEE
  
2
sen( nw0t ) 
a b
an  bn
2 n
cos( nw0t )  n
 a2  b2 a 2
 b 2 
 n n n n 
 an
 2  cos q n
 an  bn  bn 
2
Cn  an2 bn2
bn  q n  arctan  
qn  bn  senq n  an 
 a2  b2
an  n n

Y la suma puede expresarse, por ejemplo, solo en

función del coseno:
C n cos q n cos( nw0t )  senq n sen (nw0t )
 C n cos( nw0t  q n )
16/11/2022 MSc. RAUL P. CASTRO VIDAL 151
UNMSM-FIEE
 Si además definimos C0 = a0/2, la serie de Fourier
se puede escribir como:

f (t )  C0   Cn cos( nw0t  q n )
n 1

 bn 
Con: Cn  a b 2
n
2
n q n  arctan  
 an 
Ejercicio: Definir adecuadamente los coeficientes
C0, Cn y qn, de manera que la serie de Fourier
pueda escribirse como: 
f (t )  C0   C n sen ( nw0t  q n )
n 1
16/11/2022 152
MSc. RAUL P. CASTRO VIDAL
UNMSM-FIEE
 Componentes y armónicos
Hemos visto que, bajo ciertas condiciones, una
función f(t) puede escribirse como la suma de
componentes sinusoidales de diferentes
frecuencias: wn = nw0.
A la componente sinusoidal de frecuencia nw0:
cn cos(nw0t + qn) se le llama el enésimo armónico
de f(t).
Al primer armónico (n = 1) se le llama la
componente fundamental y su periodo es el
mismo que el de f(t).
A la frecuencia w0= 2 f0 = 2 / T se le llama
16/11/2022 153
frecuencia angular fundamental.
MSc. RAUL P. CASTRO VIDAL
Ejemplo: La función f(t)  cos ( 3t )  cos ( 4t )
Como vimos, tiene un periodo T = 24, por lo tanto su
frecuencia fundamental es w0 = 2/T  1/12 rad/s.
O como w0= 2f0, f0 = 1/T = 1/ 24 Hz.
Su componente fundamental (n = 1) será:
c0 cos(w0t + q0) = 0 cos(t/12).
3
f(t)=cos(t/3)+cos(t/4)
2
Tercer armónico:
1
cos(3t/12) = cos(t/4)
f(t)

0
Cuarto armónico:
-1
cos(4t/12) = cos(t/3)
-2
24
16/11/2022 -3 154
0 50 100 150 200
MSc. RAUL P. CASTRO VIDAL t
 Sea f (t ) una señal periódica con periodo T expresada en términos
de la serie compleja de Fourier siguiente :

f (t )  n
c e inw0t

n  

Derivando f (t ) respecto a t :

d
f ' (t )  f (t )   inw0cn e inw0t

dt n  

f ' (t )   n
d e
n  
inw0t

donde los coeficientes d n vienen dados por

d n  inw0cn
en consecuencia, f ' (t ) también es periódica y está representada
por16/11/2022
una serie de Fourier que es función del desarrollo en serie de 155
f (t ).
MSc. RAUL P. CASTRO VIDAL
UNMSM-FIEE
Ejercicio:
f(t)
20
T0 = 10

-10 -5 5 10 t

f '(t) T0 = 10
4

-10 -5 5 10 t
-4

f ''(t) T0 = 10
8
-10 10
-5 5 t
16/11/2022 -8 156
MSc. RAUL P. CASTRO VIDAL
 Potencia y Teorema de Parseval

El promedio o valor medio de una señal

cualquiera f(t) en un periodo dado T se
puede calcular como la altura de un
rectángulo que tenga la misma área que el
área bajo la curva de f(t)
T
Area   f ( t )dt
f(t) 0
1
h = Altura
Area = T h promedio

t
T
16/11/2022 MSc. RAUL P. CASTRO VIDAL 157
UNMSM-FIEE
De acuerdo a lo anterior, si la función
periódica f(t) representa una señal de voltaje
o corriente, la potencia promedio
entregada a una carga resistiva de 1 ohm
en un periodo está dada por:
T /2


1 2
T [ f (t )] dt
T / 2

Si f(t) es periódica, también lo será [f(t)]2 y el

promedio en un periodo será el promedio en
cualquier otro periodo.
16/11/2022 MSc. RAUL P. CASTRO VIDAL 158
UNMSM-FIEE
El teorema de Parseval nos permite calcular
la integral de [f(t)]2 mediante los coeficientes
complejos cn de Fourier de la función
periódica f(t):
T /2 

 c
2
1
T [ f (t )] 2
dt  n
T / 2 n  

O bien, en términos de los coeficientes an,

b n: T /2 
1
T  [ f (t )] dt 
2
a  1
4
2
0
1
2  (a
n 1
2
n b )
2
n
T / 2
16/11/2022 MSc. RAUL P. CASTRO VIDAL 159
UNMSM-FIEE
 Teorema o identidad de Parseval
T /2 
1 1 2 1
 [
T T / 2
f (t )]2
dt 
4
a 0  
2 n 1
( a 2
n  bn
2
)


f (t )  12 a0   [an cos(nw0t )  bn sen(nw0t )]
n 1

1 

T /2 T /2
1
T  f (t ) f (t )dt  T1  f (t ) 2 a0   [an cos(nw0t )  bn sen(nw0t )]dt 
T / 2 T / 2  n 1 
T /2  T /2  T /2
a0 an bn
 f (t )dt    f (t ) cos(nw0t )dt    f (t )sen(nw0t )dt 
T T / 2 n 1 T T / 2 n 1 T T / 2

a02 1  2
(
  an  bn
4 2 n 1
2
)
16/11/2022 160
MSc. RAUL P. CASTRO VIDAL
UNMSM-FIEE
Ejemplo. Calcular el valor cuadrático medio
de la función f(t): f(t)
1
t
... -T/ 0 T/ T ...
2 2

-1
Solución.
Del teorema de Parseval T /2 

 [ f (t )] dt   C
1 2 2
T n
T / 2 n  
y del ejemplo anterior
cn  1
n
[1  (1) n ]
sustituyendo

8  1 1 1 
c
2
 2 1     ...
  9 25 49
n
n   
16/11/2022 161
MSc. RAUL P. CASTRO VIDAL
UNMSM-FIEE
La serie numérica obtenida converge a
1 1 1
1    ...  1.2337
9 25 49
Por lo tanto,
T /2 
8
 [ f (t )] dt   c
2
1 2
 (1.2337)  1

T n 2
T / 2 n  

Como era de esperar.

16/11/2022 162
MSc. RAUL P. CASTRO VIDAL
UNMSM-FIEE
a) Sean c1 , c2   , con c1  c2 y la función:
c1 , x    ,0 )
f ( x)  
1. Calcúlese la serie de Fourier de f. c2 , x  0,  
2. Obténgase la identidad de Parseval en este caso y a partir de

1
ella calcule el valor de la serie:
 (2n  1)
n 1
2

3. ¿Converge la serie de Fourier de f puntualmente a f(0) en

x=0?

c2

c1

-π π
16/11/2022 MSc. RAUL P. CASTRO VIDAL 163
UNMSM-FIEE
 (c1  c2 )
1.
1 0 
a0    c1dx   c2 dx     c1  c2
   0  
c1  0  c2    c1  c2 
an    cos nxdx     cos nxdx    cos nxdx 
      0   0
c1  c2
 (senn  sen0)  0
n
c1  0  c2    c1  c2 
bn    sen(nx)dx     sen(nx)dx     sen(nx)dx 
     0   0


c1  c2
n
(cos n  cos 0)  c1  c2
n
( )
( 1)n  1 
n  2k  b2 k  0

 2(c2  c1 )
n  2k  1  b2 k 1  (2k  1)

c1  c2 2  (c2  c1 )
f ( x)    sen(2k  1)x
16/11/2022 2  k 1 (2k  1) MSc. RAUL P. CASTRO VIDAL
164
UNMSM-FIEE
 2.
c1  c2  1  2(c2  c1 )  1 sen(2k  1)x
f ( x) 
2
2  
 2  
k 1 (2k  1) 
 1 cos nx sen(nx) 
Como  , ,  es ortonormal en   ,   
 2   
4(c1  c2 )
( )
2
 c2  c1  
2

f (x ) dx 
1
  c2  c1  
2 2

2

 2 
 2 


k 1 (2k  1)
2



(c1  c2 )
2
 2 (c1  c2 ) 
4 2 1 


1

 2

 k 1 (2k  1) k 1 (2k  1)
2 2
2 8

c1  c2 2  (c2  c1 )
3.
c1  c2
No. Puesto que f (0)  c2 y   sen(2k  1)0 
2  k 1 (2k  1) 2
c1  c2
y en general c2  f es continua a trozos
2 y tiene derivadas laterales
16/11/2022 MSc. RAUL P. CASTRO VIDAL 165
UNMSM-FIEE
a) A partir de la serie de Fourier de la función f ( x)  x

definida en el intervalo   ,   : f ( x)     4 2 cos((2n  1)x )


determinar los valores de las series: 2 n 1  (2n  1)

 
1 1
1.  (2n  1)2 2.  (2n  1)4
n 1 n 1

1.
Particularizando para x  0, f(0)  0 :
 
4
0  cos((2n  1)0) 
2 n 1  (2n  1)2

 14 
0   
2  n 1 (2n  1)2


1   2
  2
n 1 (2 n  1)2
4 8

16/11/2022 MSc. RAUL P. CASTRO VIDAL 166
UNMSM-FIEE
2.
Aplicando la identidad de Parseval :

( )
 
2
1 a0
  an  bn
 
2
f ( x) dx 
2 2
 2 n 1
4
Sustituyendo f ( x)  x , a0   , an  , bn  0 :
 (2n  1)2

1  2 
16
 x dx 
2
 2 
  2 n 1  (2n  1)4

1 2 3  2 16  1
   2 
 3 2  n 1 (2n  1)4


1 4
 
n 1 (2n  1)
4
96

16/11/2022 MSc. RAUL P. CASTRO VIDAL 167

UNMSM-FIEE
16/11/2022 MSc. RAUL P. CASTRO VIDAL 168
UNMSM-FIEE