Teoria de Serie de Fourier

MATEMÁTICAS AVANZADAS
UNIDAD I
Series de Fourier
Semestre 2018-I
CONTENIDO
1. Funciones Periódicas
2. Serie trigonométrica de Fourier
3. Componente de directa, fundamental y armónicos
4. Ortogonalidad de las funciones seno y coseno
5. Cálculo de los coeficientes de la Serie de Fourier
6. Simetrías en señales periódicas
7. Fenómeno de Gibbs
8. Forma Compleja de las Series de Fourier
9. Espectros de frecuencia discreta
10. Potencia y Teorema de Parseval
11. De la serie a la Transformada de Fourier.
12. Obtención de la serie de Fourier usando FFT
13. Espectro de Frecuencia y medidores digitales
Preámbulo
El análisis de Fourier fue introducido en 1822 en la
“Théorie analyitique de la chaleur” para tratar la
solución de problemas de valores en la frontera en la
conducción del calor.
Más de siglo y medio después las aplicaciones de esta

teoría son muy bastas: Sistemas Lineales,
Comunicaciones, Física moderna, Electrónica,
Óptica y por supuesto, Redes Eléctricas entre
muchas otras.
Funciones Periódicas
Una Función Periódica f(t) cumple la siguiente

propiedad para todo valor de t.
f(t)=f(t+T)
A la constante mínima para la cual se cumple lo

anterior se le llama el periodo de la función
Repitiendo la propiedad se puede obtener:

f(t)=f(t+nT), donde n=0,1,  2, 3,...
Ejemplo: ¿Cuál es el período de la función f(t)  cos( 3t )  cos( 4t )?
Solución.- Si f(t) es periódica se debe cumplir:

f(t  T)  cos( t 3T )  cos( t 4T )  f(t)  cos( 3t )  cos( 4t )
Pero como se sabe cos(x+2kp)=cos(x) para cualquier entero k,
entonces para que se cumpla la igualdad se requiere que
T/3=2k1p, T/4=2k2p
Es decir,
T = 6k1p = 8k2p
Donde k1 y k2 son enteros,
El valor mínimo de T se obtiene con k1=4, k2=3, es
decir,T=24p
Gráfica de la función f(t)  cos( 3t )  cos( 4t )
3
T f(t)=cos(t/3)+cos(t/4)
2
1
f(t)
-1
-2
24p
-3
0 50 100 150 200
t
Podríamos pensar que cualquier suma de funciones seno
y coseno produce una función periódica.
Esto no es así, por ejemplo, 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= 2pm, w2T=2pn
De donde w1 m

w2 n
Es decir, la relación w1/ w2 debe ser un número racional.
Ejemplo: la función cos(3t)+cos(p+3)t no es periódica,
ya que w1  3 no es un número racional.
w2 3 p
f(t)=cos(3t)+cos((3+pi)t)
2
1
f(t)
-1
-2
0 5 10 15 20 25 30
t
ACTIVIDAD 1: 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(2pt)
3) f(t)= sen(t)+sen(t+p/2)
4) f(t)= sen(w1t)+cos(w2t)
5) f(t)= sen(2 t)
Serie Trigonométrica de Fourier
Algunas funciones periódicas f(t) de periodo T

pueden expresarse por la siguiente serie, llamada
f(t) = ½ a0 + a1cos(w0t)+a2cos(2w0t)+...
+ b1sen(w0t)+b2sen(2w0t)+...
Donde w0=2p/T.
Es decir,

f ( t )  12 a 0  [a n cos(nw0 t )  bnsen(nw0 t )]
n 1
Es posible escribir de una manera ligeramente
diferente la Serie de Fourier, si observamos que
el término ancos(nw0t)+bnsen(nw0t) se puede
escribir como
 an bn 
a 2n  b 2n  cos( nw t )  sen ( nw t ) 
 2
 2 0 2
 2 0 
 n
a b n a n b n 
Podemos encontrar una manera más compacta

para expresar estos coeficientes pensando en un
triángulo rectángulo:
an
 cos n
Cn  a 2n  b 2n a 2n  b 2n
bn
n bn
 senn
a 2n  b 2n
an
Con lo cual la expresión queda
Cn cos n cos(nw0 t )  senn sen (nw0 t )
 Cn cos( nw0 t  n )
Si además definimos C0=a0/2, la serie de Fourier

se puede escribir como

f ( t )  C 0   C n cos( nw0 t   n )
n 1
Así, Cn  a  b
2
n
2
n
1  b n 
y n  tan  
 an 
ACTIVIDAD 2:
Definir adecuadamente los coeficientes C0, Cn y
n, de manera que la serie de Fourier se pueda
escribir como

f ( t )  C0   C n sen (nw0 t  n )
n 1
Componentes y armónicas
Así, una función periódica f(t) se puede escribir como la
suma de componentes sinusoidales de diferentes
frecuencias wn=nw0.
A la componente sinusoidal de frecuencia nw0:
Cncos(nw0t+n) se le llama la enésima armónica de f(t).
A la primera armónica (n=1) se le llama la componente
fundamental y su periodo es el mismo que el de f(t)
A la frecuencia w0=2pf0=2p/T se le llama frecuencia angular
fundamental.
A la componente de frecuencia cero C0, se le

llama componente de corriente directa (cd) y
corresponde al valor promedio de f(t) en cada
periodo.
Los coeficientes Cn y los ángulos n son

respectiva-mente las amplitudes y los ángulos
de fase de las armónicas.
Ejemplo: La función f(t)  cos( t )  cos( t )
3 4
Como ya se mostró tiene un periodo T=24p, por lo tanto
su frecuencia fundamental es w0=1/12 rad/seg.
Componente fundamental es de la forma:
3
0*cos(t/12). f(t)=cos(t/3)+cos(t/4)
2
Tercer armónico:
1
cos(3t/12)=cos(t/4) 0
f(t)
Cuarto armónico: -1
Cos(4t/12)=cos(t/3) -2
24p
-3
0 50 100 150 200
t
Ejemplo: Como puede verse, la función anterior tiene
tantas partes positivas como negativas, por lo tanto su
componente de cd es cero, en cambio
f(t)  1  cos( 3t )  cos( 4t )
3
Tiene tantas partes
2
arriba como abajo
1
de 1 por lo tanto,
f(t)
0
su componente de
-1
cd es 1. f(t)=1+cos(t/3)+cos(t/4)
-2
24p
-3
0 50 100 150 200
t
ACTIVIDAD 3
: ¿Cuál es la componente fundamental, las
armónicas distintas de cero y la componente de
directa de
a) f(t) = sen2t
b) f(t) = cos2t ?
Justifícalo además mostrando la gráfica de las
funciones y marcando en ellas el periodo
fundamental y la componente de cd.
Ortogonalidad de senos y cosenos
Se dice que un conjunto de funciones 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
 f m (t)fn (t)dt  
a rn para m  n
Ejemplo: las funciones t y t2 son ortogonales en el
intervalo –1< t <1, ya que
1 1 4 1
t
 tt dt   t dt  0
2 3
1 1 4 1
Ejemplo: Las funciones sen t y cos t son ortogonales en

el intervalo –p/2< t <p/2, ya que
p 2 p
sen t
 sentcostdt  0
p 2 p
ACTIVIDAD 4:
Dar un ejemplo de un par de funciones que sean
ortogonales en el intervalo:
a) 0<t<1
b) 0<t<p
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,cosw0t, cos2w0t, cos3w0t,...,senw0t,sen2w0t,sen3w0t,...
(para cualquier valor de w0=2p/T).
Para verificar lo anterior podemos probar por pares:

1.- f(t)=1 Vs. cos(mw0t):
T/2 sen (mw0 t) T/2
2sen (mw0T/2) 2sen (mp)
 cos(mw0 t)dt    0
T / 2 mw0 T / 2 mw0 mw0
Ya que m es un entero.
2.- f(t)=1 Vs. sen(mw0t):
T/2  cos(mw0 t) T/2
 sen(mw0 t)dt  
T / 2 mw0 T / 2
1
 [cos(mw0T/2) - cos(mw0T/2)]  0
mw0
3.- cos(mw0t) Vs. cos(nw0t):
T/2  0 para m  n
 cos(mw0 t)cos(nw0 t)dt  
T / 2 T / 2 para m  n  0
4.- sen(mw0t) Vs. sen(nw0t):
T/2  0 para m  n
 sen(mw0 t)sen(nw0 t)dt  
T / 2 T / 2 para m  n  0
5.- sen(mw0t) Vs. cos(nw0t):
T/2
 sen(mw0 t)cos(nw0 t)dt  0 para cualquier m, n
T / 2
Para calcular las integrales de los casos 3, 4 y 5,
son útiles las siguientes identidades
trigonométricas:
cos A cos B = ½[cos(A+B)+cos(A-B)]

sen A sen B = ½[-cos(A+B)+cos(A-B)]
sen A cos B = ½[sen(A+B)+sen(A-B)]
Además:
sen2 = ½ (1-cos2)
cos2 = ½ (1+cos2)
Cálculo de los coeficientes de la Serie
Dada una función periódica f(t) ¿cómo se
obtiene su serie de Fourier?

n 1
Obviamente, el problema se resuelve si sabemos

como calcular los coeficientes a0,a1,a2,...,b1,b2,...
Esto se puede resolver considerando la

ortogonalidad de las funciones seno y coseno
comentada anteriormente.
Multiplicando ambos miembros por cos(nw0t) e

integrando de –T/2 a T/2, obtenemos:
T/2
an  2
T  f ( t ) cos(nw0 t )dt n  0,1,2,3,...
T / 2
Similarmente, multiplicando por sen(nw0t) e
integrando de –T/2 a T/2, obtenemos:
T/2
bn  T2  f ( t )sen (nw0 t )dt n  1,2,3,...
T / 2
Similarmente, integrando de –T/2 a T/2,
T/2
obtenemos: a  2  f ( t )dt
0 T
T / 2
El intervalo de integración no necesita ser
simétrico respecto al origen.
Como 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)
las fórmulas anteriores pueden calcularse en
cualquier intervalo que cumpla este requisito.
Ejemplo: Encontrar la Serie de Fourier para la

siguiente función de periodo T:
f(t)
1
t
... -T/ 0 T/ T ...
2 2
-1
Solución: La expresión para f(t) en –T/2<t<T/2 es

 1 para  T2  t  0
f (t)  
 1 para 0  t  T
2
T/2
Coeficientes an: a n  T2  f ( t ) cos(nw0 t )dt
T / 2
 0 T/2

 T    cos(nw0 t )dt   cos(nw0 t )dt 
2
 T / 2 0 
 1 0
1 T/2
 T2  sen (nw0 t )  sen (nw0 t ) 
 nw0 T / 2 nw0 0 

 0 para n  0
T/2
Coeficiente a0: a 0  T2  f ( t )dt
T / 2
 0 T/2

 T    dt   dt 
2
 T / 2 0 
 0 T/2
 T2  t t 
  T / 2 0 

0
T/2
Coeficientes bn: bn  T2  f ( t )sen (nw0 t )dt
T / 2
 0 T/2

 T    sen (nw0 t )dt   sen (nw0 t )dt 
2
 T / 2 0 
 1 0
1 T/2
 T2  cos(nw0 t )  cos(nw0 t ) 
 nw0 T / 2 nw0 0 

 (1  cos(np))  (cos(np)  1)

1
np

2
np
 
1  (1)n ) para n  0
Serie de Fourier: Finalmente la Serie de Fourier
queda como
f ( t )  sen (w0 t )  13 sen (3w0 t )  15 sen (5w0 t )  ...

4
p
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=p, es decir,
T=2:
Componentes de la Serie de Fourier
1.5
1
Componentes
0.5
-0.5
Suma
fundamental
-1 tercer armónico
quinto armónico
septimo armónico
-1.5
-1 -0.5 0 t 0.5 1
ACTIVIDAD 5: Encontrar la serie de Fourier

para la siguiente señal senoidal rectificada de
media onda de periodo 2p.
Senoidal rectificada de media onda
0.8
0.6
f(t)
0.4
0.2
-0.2
-6 -4 -2 0 2 4 6
t
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
2p p p 2p
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
2p p p 2p
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.
Ejemplo: ¿La función h(t)=f(1+t2) es par o
impar?, donde 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), por lo tanto h(t) es
función par, sin importar como sea f(t).
Ejemplo: De acuerdo al ejemplo anterior, todas
las siguientes funciones 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)
Como la función sen(nw0t) es una función impar
para todo n0 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
Por ejemplo, la señal cuadrada, ya analizada en
un ejemplo previo:
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:
4
p
Simetría de Media Onda
Una función periodica de periodo T se dice

simétrica de media onda, si cumple la propiedad
f ( t  12 T )   f ( t )
Es decir, si en su gráfica las partes negativas son
un reflejo de las positivas pero desplazadas
medio periodo:
f(t)
t
Simetría de Cuarto de Onda
Si una función tiene simetría de media onda y
además es función par o impar, se dice que tiene
simetría de cuarto de onda par o impar
Ejemplo: Función con simetría impar de cuarto

de onda:
f(t)
t
Ejemplo: Función con simetría par de cuarto de
onda:
f(t)
t
ACTIVIDAD 6: ¿Qué tipo de simetría tiene la

siguiente señal de voltaje producida por un triac
controlado por fase?
f(t)
t
Simetrías y Coeficientes de Fourier
Funciones
Simetría Coeficientes en la serie
T/2 T/2 Senos y
Ninguna an  2
T  f (t) cos(nw t)dt
T / 2
0 bn  2
T  f (t)sen(nw t)dt
T / 2
0 cosenos
T/2 únicamente
Par an  4
0
0
bn=0 cosenos
T/2
únicamente
Impar an=0 bn  4
0
0 senos
 0 n par  0 n par Senos y

media  T/2  T/2
an   4 bn   4 cosenos
T  T  f ( t )sen(nw0 t )dt n impar
f (t ) cos(nw0 t )dt n impar
onda
 0  impares
 0
Funciones
Simetría Coeficientes
en la serie
T/2 T/2
Senos y
Ninguna an  2
T  f (t) cos(nw0 t)dt
T / 2
bn  2
T  f (t)sen(nw0 t)dt
T / 2
cosenos
an=0 (n par)
T/4
Sólo
¼ de
onda par an  8
0
0
bn=0 cosenos
impares
(n impar )
bn=0 (n par)
¼ de T/4 Sólo
onda an=0 bn  8
0
0 senos
impar impares
(n impar )
Por ejemplo, la señal cuadrada, ya analizada en
un ejemplo previo:
f(t)
1
t
... -T/ 0 T/ T ...
2 2
-1
Es una función con simetría de ¼ de onda impar,
por ello su serie de Fourier sólo contiene
términos seno de frecuencia impar:
4
p
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, la sumatoria 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 anterior:
Fenómeno de Gibbs
1.5 Serie con 1 armónico
0.5
-0.5
-1
-1.5
-1 -0.5 0 0.5 1
Fenómeno de Gibbs
1.5 Serie con 3 armónicos
0.5
-0.5
-1
-1.5
-1 -0.5 0 0.5 1
Fenómeno de Gibbs
0.5
-0.5
-1
-1.5
-1 -0.5 0 0.5 1
Fenómeno de Gibbs
0.5
-0.5
-1
-1.5
-1 -0.5 0 0.5 1
Fenómeno de Gibbs
0.5
-0.5
-1
-1.5
-1 -0.5 0 0.5 1
Fenómeno de Gibbs
0.5
-0.5
-1
-1.5
-1 -0.5 0 0.5 1
Fenómeno de Gibbs
0.5
-0.5
-1
-1.5
-1 -0.5 0 0.5 1
Forma Compleja de la Serie de Fourier
Consideremos la serie de Fourier para una

función periodica f(t), con periodo T=2p/w0.

n 1
Es posible obtener una forma alternativa usando

las fórmulas de Euler:
cos( nw0 t )  12 (e jnw0 t  e  jnw0 t )
sen (nw0 t )  1
2j (e jnw0 t  e  jnw0 t )
Donde j  1
Sustituyendo

f ( t )  12 a 0   [a n 12 (e jnw0 t  e  jnw0 t )  b n 1
2j (e jnw0 t  e  jnw0 t )]
n 1
Y usando el hecho de que 1/j=-j


f ( t )  12 a 0   [ 12 (a n  jb n )e jnw0 t  12 (a n  jb n )e  jnw0 t ]
n 1
Y definiendo:
c0  12 a 0 , cn  12 (a n  jb n ), cn  12 (a n  jb n )
Lo cual es congruente con la fórmula para bn, ya
que b-n=-bn, ya que la función seno es impar.
La serie se puede escribir como


f ( t )  c 0   (c n e jnw0 t  c  n e  jnw0 t )
n 1
O bien,  
f (t )  c0   cn e jnw0 t
  cn e jnw0 t
n 1 n  1

Es decir,
f (t)  c e
n  
n
jnw0 t
A la expresión obtenida

f (t)  c e
n  
n
jnw0 t
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
 f ( t )e
 jnw0 t
cn  1
T dt
0
Para n=0, 1, 2, 3, ...

Los coeficientes cn son números complejos, y

también se pueden escribir en forma polar:
c n  c n e jn
Obviamente, cn  c*n  cn e jn
bn
Donde c n  a  b ,  n  arctan(  )
1 2 2
2 n n
an
Para todo n0,
Para n=0, c0 es un número real: c0  12 a 0
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):
an=0 para todo n
y b n  np [1  (1) ] para todo n
2 n
Podemos calcular los coeficientes cn de:
c n  [a n  jb n ]   j
1
2
1 2
2 np [1  (1) ] n
c n   j [1  (1) ]
1
np
n
Entonces la Serie Compleja de Fourier queda

f ( t )  p2 j(...  15 e  j5w0 t  13 e  j3w0 t  e  jw0 t
 e jw0 t  13 e j3w0 t  15 e j5w0 t  ...)
Solución 2. También podemos calcular los
coeficientes cn mediante la integral
T
 f ( t )e
 jnw0 t
cn  1
T dt
0
T/2 T
 T1 (  e  jnw0 t dt    e  jnw0 t
dt )
0 T/2
T/2 T
 jnw0 t  jnw0 t
 (1
T
1
 jnwo e  1
 jnwo e )
0 T/2
 jnw0T / 2  jnw0T  jnw0T / 2

 1
 jnwo T [(e  1)  (e e )]
Como w0T=2p y además e j  cos   jsen
cn  1
 jnwo T [( 1)  1)  (1  (1) )]
n n
 j 2
nwo T [1  (1) ]
n
  j [1  (1) ]
1
np
n
Lo cual coincide con el resultado ya obtenido.

ACTIVIDAD 7: Calcular los coeficientes cn

para la siguiente función de periodo 2p.
a) A partir de los coeficientes an,bn
b) Directamente de la integral
Senoidal rectificada de media onda
1
0.8
0.6
f(t)
0.4
0.2
0
-0.2
-6 -4 -2 0 t 2 4 6
Espectros de Frecuencia Discreta
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.
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.
Ejemplo. Para la función ya analizada:

f(t)
1
t
... -T/ 0 T/ T ...
2 2
-1
Se encontró que c n   j [1  (1) ]

1
np
n
Por lo tanto, cn  1
np
[1  (1) ]
n
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).
ACTIVIDAD 8. Dibujar el espectro de amplitud

para la función senoidal rectificada de ½ onda.
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=Th promedio
t
T
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
 [f (t)] dt
1 2
T
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.
El teorema de Parseval nos permite calcular la

integral de [f(t)]2 mediante los coeficientes com-
plejos cn de Fourier de la función periódica f(t):
T/2 
 [f (t)] dt   c
1 2 2
T n
T / 2 n  
O bien, en términos de los coeficientes an, bn:

T/2 
1
T  [f (t)] dt 
2 1
4 a 
2
0
1
2  (a
n 1
2
n b )
2
n
T / 2
Una consecuencia importante del teorema de Parseval
es el siguiente resultado:
El valor cuadrático medio de una función periódica f(t)

es igual a la suma de los valores cuadráticos medios de
sus armónicos, es decir,
2
T/2 
Cn
 [f (t )] dt  C  
1 2 2
T 0
T / 2 n 1 2
Donde Cn es la amplitud del armónico n-ésimo y C0 es
la componente de directa.
Para aclarar el resultado anterior es conveniente
encontrar la relación entre los coeficientes
complejos cn de la serie

f (t)  c e
n  
n
jnw0 t
Y los coeficientes reales Cn de la serie


f ( t )  C 0   C n cos( nw0 t   n )
n 1
Donde Cn es la amplitud del armónico n-ésimo y

C0 es la componente de directa.
Por un lado Cn  a 2n  b 2n ,
Mientras que c n  1
2 a 2n  b 2n
Entonces, cn  Cn Por lo tanto, c n  C
2 1 2
1
2 4 n
Además, para el armónico f n (t )  Cn cos(nw0 t  n )

Su valor rms es C n / 2 , por lo tanto su valor
cuadrático medio es C 2n / 2
Para la componente de directa C0, su valor rms
es C0, por lo tanto su valor cuadrático medio
será C02.
Ejemplo. Calcular el valor cuadrático medio de

la función f(t): f(t)
1
t
... -T/ 0 T/ T ...
2 2
Solución. -1
T/2 
 [f (t)] dt   c
2 2
Del teorema de Parseval 1
T n
T / 2 n  
y del ejemplo anterior c n  1

np
[1  (1) ]n

8  1 1 1 
sustituyendo n  2    
2
cn 1
 9 25 49 ...
  p
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.2337)  1
1 2 2
p
T n
T / 2 n  
Como era de esperarse.

ACTIVIDAD 9.
Calcular el valor cuadrático medio para la señal

senoidal rectificada de media onda de periodo
2p.
De la Serie a la Transformada de Fourier
La serie de Fourier nos permite obtener una

representación en el dominio de la frecuencia
para funciones periódicas f(t).
¿Es posible extender de alguna manera las series

de Fourier para obtener el dominio de la
frecuencia de funciones no periódicas?
Consideremos la siguiente función periodica de

periodo T
Tren de pulsos de amplitud 1, ancho p y periodo

T:
f(t)
1
... -T -T/ 0 T/ T ...

2 2
-p/ p/ t
2 2
0 T
2 t p
2
 p
f ( t )  1 2 t p
2
0 p
 t  T2
 2
Los coeficientes de la Serie Compleja de Fourier

en este caso resultan puramente reales:
sen ( nw p
0 2)
cn  ( T )
p
(nw0 2 )
p
El espectro de frecuencia correspondiente lo

obtenemos (en este caso) graficando cn contra
w=nw0.
Espectro del tren de pulsos para p=1, T=2

0.6
0.4
cn
0.2
-0.2
-60 -40 -20 0 20 40 60 w=nw
0
Si el periodo del tren de pulsos aumenta:
1.5
p=1, T=2
1
f(t)
0.5
0
-20 -10 0 t 10 20
1.5
p=1, T=5
1
f(t)
0.5
0
-20 -10 0 t 10 20
1.5
p=1, T=10
1
f(t)
0.5
0
-20 -10 0 t 10 20
1.5
p=1, T=20
1
f(t)
0.5
0
-20 -10 0 t 10 20
En el límite cuando T, la función deja de ser

periódica:
1.5
p=1, T=
1
f(t)
0.5
0
-20 -10 0 t 10 20
¿Qué pasa con los coeficientes de la serie de

Fourier?
0.6
p=1, T=2
0.4 cn
0.2
-0.2
w=nw0
-50 0 50
0.3
0.2 p=1, T=5

0.1
-0.1
-50 0 50
0.15
0.1
p=1, T=10
0.05
-0.05
-50 0 50
0.06
p=1, T=20
0.04
0.02
-0.02
-50 0 50
Si hace T muy grande (T): El espectro se

vuelve ¡continuo!
El razonamiento anterior nos lleva a reconsiderar
la expresión de una función f(t) no periódica en
el dominio de la frecuencia, no como una suma
de armónicos de frecuencia nw0, sino como una
función continua de la frecuencia w.

Así, la serie f (t)  c e
n  
n
jnw0 t
Al cambiar la variable discreta nw0 (cuando

T) por la variable continua w, se transforma
en una integral de la siguiente manera:
T/2
 f ( t )e
 jnw0 t
Como c n  1
T dt
T / 2
 1 T/2   jnw0 t
La serie queda f ( t )    T  f ( t )e
 jnw0 t
dt  e
n     T / 2 
 1 T/2
 
O bien, f ( t )    2p 
n   
f ( t ) e  jnw0 t
dt w
 0 e jnw0 t
T / 2 
cuando T, nw0w y w0dw y la sumatoria
se convierte en 
 
 jwt
   f (t )e
 jwt
1
f (t)  2p dt  e dw
  
Es decir,

 F(w)e
jwt
f (t)  1
2p dw Identidad
de Fourier

Donde

F(w)   f ( t )e  jwt Transformada
dt De Fourier

Estas expresiones nos permiten calcular la

expresión F(w) (dominio de la frecuencia) a
partir de f(t) (dominio del tiempo) y viceversa
Notación: A la función F(w) se le llama

transformada de Fourier de f(t) y se denota por
F, es decir

F[f ( t )]  F(w)   f ( t )e  jwt dt

En forma similar, a la expresión qu enos permite

obtener f(t) a partir de F(w) se le llama
transformada inversa de Fourier y se denota
por F –1 ,es decir

F 1[F(w)] f ( t )  
jwt
1
2p F( w) e dw

Ejemplo. Calcular F(w) para el pulso

rectangular f(t) siguiente f(t)
1
t
-p/ 0 p/
2 2
Solución. La expresión en el dominio del tiempo

de la función es 0 t p
2
 p
f ( t )  1 2 t p
2
0
2  t
p

 p/2
F(w)   f ( t )e  jwt dt   dt
e  jwt
 p / 2
Integrando  jwt

p/2
1
 jw e p / 2
 1jw (e jwp / 2  e jwp / 2 )

sen (wp / 2)
Usando la fórmula de Euler F(w)  p
wp / 2
Obsérvese que el resultado es igual al obtenido
para cn cuando T , pero multiplicado por T.
En forma Gráfica
F(w)
F(w) con p=1

1
0.5
-50 0 50 w
ACTIVIDAD 10. Calcular la Transformada de
Fourier de la función escalón unitario u(t):
u(t)
1
t
0
Graficar U(w)=F[u(t)]
¿Qué rango de frecuencias contiene U(w)?
¿Cuál es la frecuencia predominante?
La Transformada Rápida de Fourier
Cuando la función f(t) está dada por una lista de N
valores f(t1), f(t2), ...f(tN) se dice que está discretizada o
muestreada, entonces la integral que define la
Transformada de Fourier:

F(w)   f ( t )e  jwt dt

Se convierte en la sumatoria
N
F(n )   f ( t k )e
 j 2Npn ( k 1)
, para 1  n  N
k 1
(Donde k es la frecuencia discreta)

Llamada Transformada Discreta de Fourier
La Transformada Rápida de Fourier
La Transformada Discreta de Fourier (DFT)
requiere el cálculo de N funciones exponenciales
para obtener F(n), lo cual resulta un esfuerzo de
cálculo enorme para N grande.
Se han desarrollado métodos que permiten

ahorrar cálculos y evaluar de manera rápida la
Transformada discreta, a estos métodos se les
llama
Transformada Rápida de Fourier (FFT)

La FFT y la Serie de Fourier
Podemos hacer uso de la FFT para calcular los

coeficientes cn y c-n de la Serie compleja de
Fourier como sigue:
Ejemplo: Sea f(t) el tren de pulsos de ancho p y

periodo T.
f(t)
1
... -T -T/ 0 T/ T ...

2 2
-p/ p/ t
2 2
La versión muestreada f(k) de f(t) sólo puede

tomar un número finito de puntos. Tomemos por
ejemplo N=32 puntos cuidando que cubran el
intervalo de 0 a T (con p=1, T=2):
32 muestras de f(t), de 0 a T
1.5
1
f(k)
0.5
0
0 1 k 2
Para obtener estas 32 muestras usando Matlab se

puede hacer lo siguiente:
k=0:31
f=[(k<8)|(k>23)]
Plot(k,f,’o’)
Con los 32 puntos f(k) calculamos F(n) mediante

la FFT, por ejemplo, en Matlab:
F=fft(f)/N;
Con lo que obtenemos 32 valores complejos de
F(n). Estos valores son los coeficientes de la
serie compleja ordenados como sigue:
n 1 2 3 4 ... 16 17 18 19 ... 32
F(n) c0 c1 c2 c3 ... c15 c-16 c-15 c-14 ... c-1
Podemos graficar el espectro de amplitud

reordenando previamente F(n) como sigue
aux=F;
F(1:16)=aux(17:32);
F(17:32)=aux(1:16);
F(n) queda:
n 1 ... 13 14 15 16 17 18 19 ... 32
F(n) c-16 ... c-3 c-2 c-1 c0 c1 c2 c3 ... c15
Y para graficar el espectro de amplitud:
stem(abs(F))
Obteniéndose:
0.6
Espectro de Amplitud |F(n)|
|F(n) Para el tren de pulsos p=1,
|
T=2
0.4
0.2
0 n
0 10 20 30
Si deseamos una escala horizontal en unidades
de frecuencia (rad/seg):
w0=2*pi/T;
n=-16:15;
w=n*w0;
Stem(w,abs(F))
Espectro de Amplitud |F(n)|
0.6
para el tren de pulsos, p=1,T=2
|F(w)|
0.4
Obteniendo:
0.2
0
-50 0 50 w
También podemos obtener los coeficientes de la

forma trigonométrica, recordando que:
cn  (a n  jb n ), cn  (a n  jb n )
1
2
1
2
Podemos obtener
a 0  c0 , a n  2 Re( cn ), b  2 Im( cn )
Para el ejemplo se obtiene: a0=0.5, an=bn=0 (para
n par), además para n impar:
n 1 3 5 7 9 11 13 15
an 0.6346 -0.2060 0.1169 -0.0762 0.0513 -0.0334 0.0190 -0.0062

bn -0.0625 0.0625 -0.0625 0.0625 -0.0625 0.0625 -0.0625 0.0625
Como el tren de pulsos es una función par, se

esperaba que bn=0; (el resultado obtenido es
erróneo para bn, pero el error disminuye para N
grande): 1
a0
0.5
Coeficientes bn Coeficientes an
0
-0.5
0 10 20 30
ACTIVIDA 11: Usar el siguiente código para generar
128 puntos de una función periódica con frecuencia
fundamental w0=120p (60 hertz) y dos armónicos
impares en el intervalo [0,T]:
N=128;
w0=120*pi;
T=1/60;
t=0:T/(N-1):T;
f=sin(w0*t)+0.2*sin(3*w0*t)+0.1*sin(11*w0*t);
Usando una función periódica diferente a la subrayada:

a) Graficar la función.
b) Obtener y graficar el espectro de amplitud de la señal
usando la función FFT
Medidores Digitales
La FFT ha hecho posible el desarrollo de equipo

electrónico digital con la capacidad de cálculo de
espectros de frecuencia para señales del mundo
real, por ejemplo:
1) Osciloscopio digital Fuke 123 ($ 18,600.00 M.N.)

2) Osc. digital Tektronix THS720P ($3,796 dls)
3) Power Platform PP-4300
Medidores Digitales
El Fluke 123 scope meter

Medidores Digitales
Tektronix THS720P (osciloscopio digital)

Medidores Digitales
Analizador de potencia PP-4300
Es un equipo especializado en monitoreo de la

calidad de la energía: permite medición de 4
señales simultáneas (para sistemas trifásicos)
QUIZ 1
• 1.REALIZAR LAS ACTIVIDADES 1 A 11.
• 2 DESCRIBA EL PROCESO DE CONVOLUCION DIGITAL Y
ANALOGICO.
• 3 DESCRIBA LAS PROPIEDADES DE LA TRANSFORMADA DE
FOURIER Y FFT .
• 4. REPRESENTE EN MATLAB, LAS FUNCIONES TRIANGULAR,
EXPONENCIAL UN SOLO PULSO CUADRADO, Y UNA SEÑAL
ALEATORIA
• 5.DIGA COMO OBTIENE NUMEROS ALEATORIOS DE 1 A 1000 CON
MATLAB.

Teoria de Serie de Fourier

Cargado por

Copyright:

Formatos disponibles

Teoria de Serie de Fourier

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Teoria de Serie de Fourier

Cargado por

Copyright:

Formatos disponibles

MATEMÁTICAS AVANZADAS

Más de siglo y medio después las aplicaciones de esta

Una Función Periódica f(t) cumple la siguiente

A la constante mínima para la cual se cumple lo

Repitiendo la propiedad se puede obtener:

Solución.- Si f(t) es periódica se debe cumplir:

Esto no es así, por ejemplo, consideremos la función

ACTIVIDAD 1: Encontrar el periodo de las

Algunas funciones periódicas f(t) de periodo T

Podemos encontrar una manera más compacta

Si además definimos C0=a0/2, la serie de Fourier

A la componente de frecuencia cero C0, se le

Los coeficientes Cn y los ángulos n son

Se dice que un conjunto de funciones fk(t) son

Ejemplo: Las funciones sen t y cos t son ortogonales en

Para verificar lo anterior podemos probar por pares:

3.- cos(mw0t) Vs. cos(nw0t):

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

cos A cos B = ½[cos(A+B)+cos(A-B)]

Obviamente, el problema se resuelve si sabemos

Esto se puede resolver considerando la

Multiplicando ambos miembros por cos(nw0t) e

Ejemplo: Encontrar la Serie de Fourier para la

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

 (1  cos(np))  (cos(np)  1)

f ( t )  sen (w0 t )  13 sen (3w0 t )  15 sen (5w0 t )  ...

ACTIVIDAD 5: Encontrar la serie de Fourier

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

En forma similar, una función f(t) se dice

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

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

Una función periodica de periodo T se dice

Ejemplo: Función con simetría impar de cuarto

ACTIVIDAD 6: ¿Qué tipo de simetría tiene la

 0 n par  0 n par Senos y

1.5 Serie con 1 armónico

1.5 Serie con 3 armónicos

1.5 Serie con 5 armónicos

1.5 Serie con 7 armónicos

1.5 Serie con 13 armónicos

1.5 Serie con 50 armónicos

1.5 Serie con 100 armónicos

Consideremos la serie de Fourier para una

Es posible obtener una forma alternativa usando

Y usando el hecho de que 1/j=-j

La serie se puede escribir como

Se le llama forma compleja de la serie de

Para n=0, 1, 2, 3, ...

Los coeficientes cn son números complejos, y

Obviamente, cn  c*n  cn e jn

Entonces la Serie Compleja de Fourier queda

 jnw0T / 2  jnw0T  jnw0T / 2

Lo cual coincide con el resultado ya obtenido.

ACTIVIDAD 7: Calcular los coeficientes cn

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

Como n sólo toma valores enteros, la frecuencia angular

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

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

Ejemplo. Para la función ya analizada: