Series de Fourier

Curso de Titulacin
Modelado y Anlisis de Sistemas Elctricos bajo

Condiciones de Operacin no Senoidales
Facultad de Ingeniera Elctrica

Universidad Michoacana
de San Nicols de Hidalgo
Febrero de 2003
Series de Fourier. 1
Series de Fourier
Contenido
1. Funciones Peridicas
2. Serie trigonomtrica de Fourier
3. Componente de directa, fundamental y armnicos
4. Ortogonalidad de las funciones seno y coseno
5. Clculo de los coeficientes de la Serie de Fourier
6. Simetras en seales peridicas
7. Fenmeno 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. Obtencin de la serie de Fourier usando FFT
13. Espectro de Frecuencia y medidores digitales
Prembulo
El anlisis de Fourier fue introducido en 1822 en la
Thorie analyitique de la chaleur para tratar la
solucin de problemas de valores en la frontera en la
conduccin del calor.
Ms de siglo y medio despus las aplicaciones de esta
teora son muy bastas: Sistemas Lineales,
Comunicaciones, Fsica moderna, Electrnica,
ptica y por supuesto, Redes Elctricas entre
muchas otras.
Funciones Peridicas
Una Funcin Peridica f(t) cumple la siguiente
propiedad para todo valor de t.
f(t)=f(t+T)
A la constante mnima para la cual se cumple lo
anterior se le llama el periodo de la funcin
Repitiendo la propiedad se puede obtener:
f(t)=f(t+nT), donde n=0,1, 2, 3,...
Funciones Peridicas
Ejemplo: Cul es el perodo de la funcin f(t) cos( 3t ) cos( 4t )?
Solucin.- Si f(t) es peridica se debe cumplir:
f(t T) cos( t 3T ) cos( t 4T ) f(t) cos( 3t ) cos( 4t )
Pero como se sabe cos(x+2k)=cos(x) para cualquier entero k,

entonces para que se cumpla la igualdad se requiere que
T/3=2k1, T/4=2k2
Es decir,
T = 6k1= 8k2
Donde k1 y k2 son enteros,
El valor mnimo de T se obtiene con k1=4, k2=3, es
decir,T=24
Funciones Peridicas
Grfica de la funcin
f(t) cos( 3t ) cos( 4t )
3
2
f(t)=cos(t/3)+cos(t/4)
f(t)
1
0
-1
-2
-3
24
0
50
100
150
200
Funciones Peridicas
Podramos pensar que cualquier suma de funciones seno
y coseno produce una funcin peridica.
Esto no es as, por ejemplo, consideremos la funcin
f(t) = cos(1t)+cos(2t).
Para que sea peridica se requiere encontrar dos enteros
m, n tales que
1T= 2m, 2T=2n
1 m
De donde
Es decir, la relacin 1/ 2 debe ser un nmero racional.

Funciones Peridicas
Ejemplo: la funcin cos(3t)+cos(+3)t no es peridica,
ya que 1 3 no es un nmero racional.
2
f(t)=cos(3t)+cos((3+pi)t)
f(t)
-1
-2
10
15
20
25
30
Funciones Peridicas
Tarea: Encontrar el periodo de las siguientes
funciones, si es que son peridicas:
1) f(t) = sen(nt), donde n es un entero.
2) f(t)= sen2(2t)
3) f(t)= sen(t)+sen(t+)
4) f(t)= sen(1t)+cos(2t)
5) f(t)= sen(2 t)
Serie Trigonomtrica de Fourier

Algunas funciones peridicas f(t) de periodo T
pueden expresarse por la siguiente serie, llamada
f(t) = a0 + a1cos(0t)+a2cos(20t)+...
+ b1sen(0t)+b2sen(20t)+...
Donde 0=2/T.
Es decir,
f ( t ) 12 a 0 [a n cos(n0 t ) b n sen (n0 t )]

n 1

Es posible escribir de una manera ligeramente
diferente la Serie de Fourier, si observamos que
el trmino ancos(n0t)+bnsen(n0t) se puede
escribir como
a 2n
an
b 2n
a 2n b 2n
cos(n0 t )
bn
a 2n b 2n
sen (n0 t )
Podemos encontrar una manera ms compacta

para expresar estos coeficientes pensando en un
tringulo rectngulo:

an
bn
C n a 2n b 2n
n
an
a 2n
b 2n
bn
a 2n
b 2n
cos n
senn
Con lo cual la expresin queda

C n cos n cos(n0 t ) sen n sen (n0 t )
C n cos(n0 t n )

Si adems definimos C0=a0/2, la serie de Fourier
se puede escribir como
f ( t ) C 0 C n cos(n0 t n )
n 1
As,
y
Cn a b
2
n
2
n
bn
n tan
an

Tarea:
Definir adecuadamente los coeficientes C 0, Cn y
n, de manera que la serie de Fourier se pueda
escribir como
f ( t ) C 0 C n sen (n0 t n )
n 1
Componentes y armnicas
As, una funcin peridica f(t) se puede escribir
como la suma de componentes sinusoidales de
diferentes frecuencias n=n0.
A la componente sinusoidal de frecuencia n0:
Cncos(n0t+n) se le llama la ensima armnica de
f(t).
A la primera armnica (n=1) se le llama la
componente fundamental y su periodo es el mismo
que el de f(t)
A la frecuencia 0=2f0=2/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 armnicas.
Ejemplo: La funcin f(t) cos( 3t ) cos( 4t )
Como ya se mostr tiene un periodo T=24, por lo tanto
su frecuencia fundamental es 0=1/12 rad/seg.
f(t)
Componente fundamental es de la forma:

3
0*cos(t/12).
f(t)=cos(t/3)+cos(t/4)
2
Tercer armnico:
1
cos(3t/12)=cos(t/4) 0
Cuarto armnico:
-1
Cos(4t/12)=cos(t/3) -2
-3
0
24
50
100
150
200
Ejemplo: Como puede verse, la funcin 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 )

2
1
f(t)
Tiene tantas partes

arriba como abajo
de 1 por lo tanto,
su componente de
cd es 1.
-1
-2
-3
0
f(t)=1+cos(t/3)+cos(t/4)
24
50
100
150
200
Tarea: Cul es la componente fundamental, las
armnicas distintas de cero y la componente de
directa de
a) f(t) = sen2t
b) f(t) = cos2t ?
Justifcalo adems mostrando la grfica 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)f n (t)dt
a
rn para m n

Ejemplo: las funciones t y t2 son ortogonales en el
intervalo 1< t <1, ya que
1
4 1
t
tt dt t dt
4
1
1
2
Ejemplo: Las funciones sen t y cos t son ortogonales en

el intervalo /2< t </2, ya que
sen t
sentcostdt
2

Tarea:
Dar un ejemplo de un par de funciones que sean
ortogonales en el intervalo:
a) 0<t<1
b) 0<t<

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,cos0t, cos20t, cos30t,...,sen0t,sen20t,sen30t,...
(para cualquier valor de 0=2/T).
Para verificar lo anterior podemos probar por pares:
1.- f(t)=1 Vs. cos(m0t):
sen (m0 t)
cos(m0 t)dt
m0
T / 2
T/2
T/2
T / 2
2sen (m0T/2) 2sen (m)
0
m0
m0
Ya que m es un entero.

2.- f(t)=1 Vs. sen(m0t):
T/2
cos(m0 t)
sen(m0 t)dt
m0
T / 2
T/2
T / 2
[cos(m0T/2) - cos(m0T/2)] 0
m0
3.- cos(m0t) Vs. cos(n0t):

para m n
T/2
T / 2
T / 2 para m n 0
cos(m0 t)cos(n0 t)dt

4.- sen(m0t) Vs. sen(n0t):
para m n
T/2
T / 2
T / 2 para m n 0
sen(m0 t)sen(n0 t)dt
5.- sen(m0t) Vs. cos(n0t):

T/2
sen(m0 t)cos(n0 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
trigonomtricas:
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)]
Adems:
sen2= (1-cos2)
cos2= (1+cos2)
Clculo de los coeficientes de la Serie

Dada una funcin peridica f(t) cmo se
obtiene su serie de Fourier?
f ( t ) 12 a 0 [a n cos(n0 t ) b n sen ( n0 t )]
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(n0t) e
integrando de T/2 a T/2, obtenemos:
an
2
T
T/2
f ( t ) cos(n0 t )dt
T / 2
n 0,1,2,3,...
Similarmente, multiplicando por sen(n0t) e

integrando de T/2 a T/2, obtenemos:
bn
2
T
T/2
f ( t )sen (n0 t )dt
T / 2
n 1,2,3,...
Similarmente, integrando de T/2 a T/2,

T/2
obtenemos:
a 0 T2 f ( t )dt
T / 2

El intervalo de integracin no necesita ser
simtrico respecto al origen.
Como la ortogonalidad de las funciones seno y
coseno no slo 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 frmulas anteriores pueden calcularse en
cualquier intervalo que cumpla este requisito.

Ejemplo: Encontrar la Serie de Fourier para la
siguiente funcin de periodo T:
1
...
-T
/2
f(t)
/2
T ...
-1
Solucin: La expresin para f(t) en T/2<t<T/2 es

1 para T2 t 0
f (t )
T
1
para
0
2

Coeficientes an:
T/2
a n T2 f ( t ) cos(n0 t )dt
T / 2
T/2
T / 2
cos(n0 t )dt cos(n0 t )dt
0
T/2
1
1
T
sen (n0 t )
sen (n0 t )
n0
n0
T / 2
0
0 para n 0

T/2
Coeficiente a0:
a 0 T2 f ( t )dt
T / 2
T/2
dt dt
T / 2
0
T t
t
T / 2
T/2

T/2
Coeficientes bn:
b n T2 f ( t )sen (n0 t )dt

T / 2
T/2
T / 2
sen ( n0 t )dt sen ( n0 t )dt
0
T/2
1
1
T
cos(n0 t )
cos(n0 t )
n0
n0
T / 2
0
1
(1 cos(n)) (cos(n) 1)
n
2
1 ( 1) n ) para n 0
n

Serie de Fourier: Finalmente la Serie de Fourier
queda como
4
f ( t ) sen (0 t ) 13 sen (30 t ) 15 sen (50 t ) ...
En la siguiente figura se muestran: la

componente fundamental y los armnicos 3, 5 y
7 as como la suma parcial de estos primeros
cuatro trminos de la serie para 0=, es decir,
T=2:

1.5
Componentes de la Serie de Fourier
Componentes
1
0.5
0
-0.5
Suma
fundamental
tercer armnico
quinto armnico
septimo armnico
-1
-1.5
-1
-0.5
0.5
1

Tarea: Encontrar la serie de Fourier para la
siguiente seal senoidal rectificada de media
onda de periodo 2.
Senoidal rectificada de media onda
1
f(t)
0.8
0.6
0.4
0.2
0
-0.2
-6
-4
-2
6
Funciones Pares e Impares

Una funcin (peridica o no) se dice funcin par
(o con simetra par) si su grfica es simtrica
respecto al eje vertical, es decir, la funcin f(t) es
par si
f(t) = f(-t)
f(t)

En forma similar, una funcin f(t) se dice
funcin impar o con simetra impar, si su grfica
es simtrica respecto al origen, es decir, si
cumple lo siguiente: -f(t) = f(-t)
f(t)

Ejemplo: Las siguientes funciones son pares o
impares?
f(t) = t+1/t
g(t) = 1/(t2+1),
Solucin:
Como f(-t) = -t-1/t = -f(t), por lo tanto f(t) es
funcin impar.
Como g(-t)=1/((-t)2+1) = 1/(t2+1)=g(t), por lo
tanto g(t) es funcin par.

Ejemplo: La funcin h(t)=f(1+t2) es par o
impar?, donde f es una funcin arbitraria.
Solucin:
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
funcin 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 funcin sen(n0t) es una funcin impar
para todo n0 y la funcin cos(n0t) es una
funcin par para todo n, es de esperar que:
Si f(t) es par, su serie de Fourier no contendr
trminos seno, por lo tanto bn= 0 para todo n
Si f(t) es impar, su serie de Fourier no
contendr trminos coseno, por lo tanto a n= 0
para todo n

Por ejemplo, la seal cuadrada, ya analizada en
un ejemplo previo:
1
...
-T
/2
f(t)
/2
T ...
-1
Es una funcin impar, por ello su serie de

Fourier no contiene trminos coseno:
4
f ( t ) sen (0 t ) 13 sen (30 t ) 15 sen (50 t ) ...
Simetra de Media Onda

Una funcin periodica de periodo T se dice
simtrica de media onda, si cumple la propiedad
f ( t 12 T ) f ( t )
Es decir, si en su grfica las partes negativas son
un reflejo de las positivas pero desplazadas
medio periodo:
f(t)
Simetra de Cuarto de Onda

Si una funcin tiene simetra de media onda y
adems es funcin par o impar, se dice que tiene
simetra de cuarto de onda par o impar
Ejemplo: Funcin con simetra impar de cuarto
de onda:
f(t)

Ejemplo: Funcin con simetra par de cuarto de
onda:
f(t)

Tarea: Qu tipo de simetra tiene la siguiente
seal de voltaje producida por un triac
controlado por fase?
f(t)
Simetras y Coeficientes de Fourier

Simetra
Funciones
en la serie
Coeficientes
T/2
an
Ninguna
2
T
T/2
f (t) cos(n t)dt
bn
2
T
T / 2
f (t)sen(n t)dt
0
T / 2
T/2
Par
an
4
T
Senos y
cosenos
nicamente
cosenos
bn=0
f (t) cos(n t)dt

0
T/2
bn
an=0
Impar
4
T
f (t)sen(n t)dt
0
nicamente
senos
media
onda
T/2
an
4
T
f (t) cos(n t)dt

0
n par
n impar
n par
T/2
bn
4
T
f (t)sen(n t)dt
0
n impar
Senos y
cosenos
impares

Simetra
Funciones
en la serie
Coeficientes
T/2
Ninguna
an
2
T
f (t) cos(n0 t)dt
T/2
bn
f (t)sen(n0 t)dt
2
T
T / 2
T / 2
an=0 (n par)
de
onda par
T/4
an
8
T
Slo
cosenos
impares
bn=0
f (t) cos(n t)dt

0
Senos y
cosenos
(n impar )
bn=0 (n par)
de
onda
impar
T/4
an=0
bn
8
T
f (t)sen(n t)dt
0
(n impar )
Slo
senos
impares

Por ejemplo, la seal cuadrada, ya analizada en
un ejemplo previo:
1
...
-T
/2
f(t)
/2
T ...
-1
Es una funcin con simetra de de onda impar,

por ello su serie de Fourier slo contiene
trminos seno de frecuencia impar:
4
f ( t ) sen (0 t ) 13 sen (30 t ) 15 sen (50 t ) ...
Fenmeno de Gibbs
Si la serie de Fourier para una funcin f(t) se
trunca para lograr una aproximacin en suma
finita de senos y cosenos, es natural pensar que a
medida que agreguemos ms armnicos, la
sumatoria se aproximar ms 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 armnicos.
Por ejemplo, consideremos el tren de pulsos
anterior:
Fenmeno de Gibbs
1.5
Serie con 1 armnico
1
0.5
0
-0.5
-1
-1.5
-1
-0.5
0.5
1
Fenmeno de Gibbs
1.5
Serie con 3 armnicos
1
0.5
0
-0.5
-1
-1.5
-1
-0.5
0.5
1
Fenmeno de Gibbs
1.5
1
0.5
0
-0.5
-1
-1.5
-1
-0.5
0.5
1
Fenmeno de Gibbs
1.5
1
0.5
0
-0.5
-1
-1.5
-1
-0.5
0.5
1
Fenmeno de Gibbs
1.5
1
0.5
0
-0.5
-1
-1.5
-1
-0.5
0.5
1
Fenmeno de Gibbs
1.5
1
0.5
0
-0.5
-1
-1.5
-1
-0.5
0.5
1
Fenmeno de Gibbs
1.5
1
0.5
0
-0.5
-1
-1.5
-1
-0.5
0.5
1
Forma Compleja de la Serie de Fourier

Consideremos la serie de Fourier para una
funcin periodica f(t), con periodo T=2/0.
f (t)
1
a
2 0
[a n cos(n0 t ) b n sen (n0 t )]

n 1
Es posible obtener una forma alternativa usando

las frmulas de Euler:
cos(n0 t ) 12 (e jn0 t e jn0 t )
sen (n0 t )
Donde
1
2j
(e jn0 t e jn0 t )
j 1

Sustituyendo
f ( t ) 12 a 0 [a n 12 (e jn0 t e jn0 t ) b n
n 1
1
2j
(e jn0 t e jn0 t )]
Y usando el hecho de que 1/j=-j
f ( t ) 12 a 0 [ 12 (a n jb n )e jn0 t 12 (a n jb n )e jn0 t ]
n 1
Y definiendo:
c0 12 a 0 , c n 12 (a n jb n ), c n 12 (a n jb n )
Lo cual es congruente con la frmula para b n, ya
que b-n=-bn, ya que la funcin seno es impar.

La serie se puede escribir como
f ( t ) c 0 (c n e jn0 t c n e jn0 t )
n 1
O bien,
f (t ) c0 cn e
jn0 t
n 1
Es decir,
f (t)
cn e
jn0 t
n 1
c e
jn0 t

A la expresin obtenida
f (t)
c e
jn0 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
cn
1
T
f ( t )e
jn0 t
dt
Para n=0, 1, 2, 3, ...

0

Los coeficientes cn son nmeros complejos, y
tambin se pueden escribir en forma polar:
c n c n e j n
c n c*n c n e jn
Obviamente,
cn
1
2
Donde
Para todo n0,
a b
2
n
2
n
bn
n arctan( )
an
Para n=0, c0 es un nmero real:
c 0 12 a 0

Ejemplo. Encontrar la forma compleja de la
serie de Fourier para la funcin ya tratada:
1
...
-T
/2
f(t)
/2
T ...
-1
Solucin 1. Como ya se calcularon los

coeficientes de la forma trigonomtrica (a n y bn):
an=0 para todo n
n
2
b n n [1 (1) ] para todo n
y

Podemos calcular los coeficientes cn de:
c n [a n jb n ] j
1
2
1 2
2 n
[1 (1) ]
n
c n j [1 (1) ]
n
1
n
Entonces la Serie Compleja de Fourier queda

f ( t ) 2 j(... 15 e j50 t 13 e j30 t e j0 t
e
j0 t
e
1
3
j3 0 t
e
1
5
j5 0 t
...)

Solucin 2. Tambin podemos calcular los
coeficientes cn mediante la integral
T
cn
1
T
f ( t )e
jn0 t
dt
T/2
T1 ( e jn0 t dt
0
T1 ( jn1o e jn0 t
T/2
jn0 t
e
dt )
T/2
jn1o e jn0 t
1
jno T
[(e
jn0 T / 2
1) (e
)
T/2
jn0 T
jn0 T / 2
)]

Como 0T=2 y adems e
cn
1
jno T
cos jsen
[(1) 1) (1 (1) )]
2
no T
[1 (1) ]
n
j [1 (1) ]
1
n
Lo cual coincide con el resultado ya obtenido.

Tarea: Calcular los coeficientes cn para la
siguiente funcin de periodo 2.
a) A partir de los coeficientes an,bn
b) Directamente de la integral
Senoidal rectificada de media onda
1
f(t)
0.8
0.6
0.4
0.2
0
-0.2
-6
-4
-2
6
Espectros de Frecuencia Discreta

A la grfica de la magnitud de los coeficientes c n
contra la frecuencia angular de la componente
correspondiente se le llama el espectro de
amplitud de f(t).
A la grfica del ngulo de fase n de los
coeficientes cn contra , se le llama el espectro de
fase de f(t).
Como n slo toma valores enteros, la frecuencia
angular =n0 es una variable discreta y los
espectros mencionados son grficas discretas.

Dada una funcin peridica f(t), le corresponde
una y slo 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 funcin en el dominio del
tiempo.

Ejemplo. Para la funcin ya analizada:
1
...
-T
/2
f(t)
/2
T ...
-1
Se encontr que c n j [1 (1) ]

n
1
n
Por lo tanto,
cn
1
n
[1 (1) ]
n

El espectro de amplitud se muestra a continuacin
0.7
Espectro de Amplitud de f(t)
0.6
Cn
0.5
0.4
0.3
0.2
0.1
0
-30
-20
-10
10
20
30
Frecuencia negativa
Frecuencia
(?) El eje horizontal es un eje de frecuencia,
Observacin:
(n=nmero de armnico = mltiplo de 0).

Tarea. Dibujar el espectro de amplitud para la
funcin senoidal rectificada de onda.
Potencia y Teorema de Parseval

El promedio o valor medio de una seal
cualquiera f(t) en un periodo dado (T) se puede
calcular como la altura de un rectngulo que
tenga la misma rea que el rea bajo la curva de
f(t)
T
f(t)
Area f ( t )dt
0
1
Area=Th
T
h=Altura
promedio
t

De acuerdo a lo anterior, si la funcin peridica
f(t) representa una seal de voltaje o corriente, la
potencia promedio entregada a una carga
resistiva de 1 ohm en un periodo est dada por
T/2
1
T
2
[
f
(
t
)]
dt
T / 2
Si f(t) es peridica, tambin 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 complejos cn de Fourier de la funcin peridica f(t):
T/2
1
T
[f (t )] dt c
2
T / 2
2
n
O bien, en trminos de los coeficientes a n, bn:

1
T
T/2
T / 2
n 1
2
2
2
1 2
1
[
f
(
t
)]
dt
(
a
b
n
n)
4 0
2

Una consecuencia importante del teorema de
Parseval es el siguiente resultado:
El valor cuadrtico medio de una funcin
peridica f(t) es igual a la suma de los valores
cuadrticos medios de sus armnicos, es decir,
T/2
1
T
[f (t )] dt C
2
T / 2
2
0
n 1
Cn
2
Donde Cn es la amplitud del armnico n-simo y

C0 es la componente de directa.

Para aclarar el resultado anterior es conveniente
encontrar la relacin entre los coeficientes
complejos cn de la serie
f (t)
jn0 t
c
e
n
Y los coeficientes reales Cn de la serie
f ( t ) C 0 C n cos(n0 t n )
n 1
Donde Cn es la amplitud del armnico n-simo y

C0 es la componente de directa.

2
2
C
b
Por un lado n
n
n,
Mientras que c n
1
2
a 2n b 2n
C nlo tanto,
Entonces, c n Por
1
2
cn C
1
4
2
n
Adems, para el armnico f n ( t ) C n cos(n0 t n )

C n / 2 , por lo tanto su valor
Su valor rms es
cuadrtico medio es C 2n / 2
Para la componente de directa C 0, su valor rms es
C0 , por lo tanto su valor cuadrtico medio ser
C0 2.

Ejemplo. Calcular el valor cuadrtico medio de
la funcin f(t):
f(t)
1
...
-T
/2
/2
-1
Solucin.
Del teorema de Parseval
T/2
1
T
y del ejemplo anterior c n

c
sustituyendo
n
2
n
8
2
[f (t )] dt c
2
T / 2
T ...
1
n
2
n
[1 (1) n ]
1 1
1
...
9 25 49

La serie numrica obtenida converge a
1 1
1
1
... 1.2337
9 25 49
Por lo tanto,
T/2
1
T
[f (t )] dt c
2
T / 2
2
n
8
2 (1.2337) 1
Como era de esperarse.

Tarea.
Calcular el valor cuadrtico medio para la seal
senoidal rectificada de media onda de periodo
2.
De la Serie a la Transformada de Fourier

La serie de Fourier nos permite obtener una
representacin en el dominio de la frecuencia
para funciones peridicas f(t).
Es posible extender de alguna manera las series
de Fourier para obtener el dominio de la
frecuencia de funciones no peridicas?
Consideremos la siguiente funcin periodica de
periodo T

Tren de pulsos de amplitud 1, ancho p y periodo
T:
1
f(t)
...
-T
-T
/2
0
-p
f (t) 1
0
/2
p
/2
T ...
p
2
p
2
/2
T
2
p
2
p
2
t
t
T
2

Los coeficientes de la Serie Compleja de Fourier
en este caso resultan puramente reales:
p
sen
(
n
p
0 2)
cn ( T )
p
(n0 2 )
El espectro de frecuencia correspondiente lo

obtenemos (en este caso) graficando c n contra
=n0.

Espectro del tren de pulsos para p=1, T=2
0.6
cn
0.4
0.2
0
-0.2
-60
-40
-20
20
40
60 w=nw0

Si el periodo del tren de pulsos aumenta:
1.5
p=1, T=2
f(t)
1
0.5
0
-20
-10
1.5
10
20
10
20
p=1, T=5
f(t)
1
0.5
0
-20
-10
1.5
p=1, T=10
f(t)
0.5
0
-20
-10
1.5
10
20
10
20
p=1, T=20
f(t)
0.5
0
-20
-10

En el lmite cuando T, la funcin deja de ser
peridica:
1.5
p=1, T=
f(t)
0.5
0
-20
-10
10
20
Qu pasa con los coeficientes de la serie de

Fourier?

0.6
cn
p=1, T=2
0.4
0.2
0
-0.2
-50
=n0
50
0.3
p=1, T=5
0.2
0.1
0
-0.1
-50
0.15
50
p=1, T=10
0.1
0.05
0
-0.05
-50
0.06
50
p=1, T=20
0.04
0.02
0
-0.02
-50
50

Si hace T muy grande (T): El espectro se
vuelve continuo!

El razonamiento anterior nos lleva a reconsiderar
la expresin de una funcin f(t) no peridica en
el dominio de la frecuencia, no como una suma
de armnicos de frecuencia n0, sino como una
funcin continua de la frecuencia .
As, la serie
f (t)
c e
jn0 t
Al cambiar la variable discreta n0 (cuando

T) por la variable continua , se transforma
en una integral de la siguiente manera:

T/2
Como c n
1
T
jn0 t
f
(
t
)
e
dt
T / 2
La serie queda f ( t )
O bien, f ( t )
T/2
1
2
T/2
f ( t )e
1
T
f ( t )e
T / 2
jn0 t
T / 2
jn0 t
dt e jn0 t
dt 0 e jn0 t
cuando T, n0 y 0d y la sumatoria

se convierte en
j t
jt
1
f (t)
f (t )e
dt e d

Es decir,
f (t)
1
2
F()e
jt
Donde
F() f ( t )e
j t
dt
Identidad
de Fourier
Transformada
De Fourier
Estas expresiones nos permiten calcular la

expresin F() (dominio de la frecuencia) a
partir de f(t) (dominio del tiempo) y viceversa

Notacin: A la funcin F() se le llama
transformada de Fourier de f(t) y se denota por
F, es decir
F[f ( t )] F() f ( t )e jt dt
En forma similar, a la expresin 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()] f ( t )
1
2
j t
F
(
)
e
d

Ejemplo. Calcular F(w)
rectangular f(t) siguiente f(t)
para
el
pulso
t
/2
-p
/2
Solucin. La expresin en el dominio del tiempo

de la funcin es
0
t p
f (t) 1
0
p
2
t
p
2 t
p
2
p/2
p / 2
F() f ( t )e jt dt
Integrando
1
j
1
j
j t
jt
e
dt
p/2
p / 2
( e j p / 2 e j p / 2 )
sen (p / 2)
Usando la frmula de Euler F() p
p / 2
Obsrvese que el resultado es igual al obtenido

para cn cuando T , pero multiplicado por T.

En forma Grfica
F(w)
F(w) con p=1
0.5
0
-50
50

Tarea. Calcular la Transformada de Fourier de la
funcin escaln unitario u(t):
u(t)
1
t
0
Graficar U()=F[u(t)]
Qu rango de frecuencias contiene U()?
Cul es la frecuencia predominante?
La Transformada Rpida de Fourier

Cuando la funcin 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() f ( t )e jt dt
Se convierte en la sumatoria
N
F(n ) f ( t k )e
j 2Nn ( k 1)
, para 1 n N
k 1
(Donde k es la frecuencia discreta)

Llamada Transformada Discreta de Fourier
La Transformada Rpida de Fourier

La Transformada Discreta de Fourier (DFT)
requiere el clculo de N funciones exponenciales
para obtener F(n), lo cual resulta un esfuerzo de
clculo enorme para N grande.
Se han desarrollado mtodos que permiten
ahorrar clculos y evaluar de manera rpida la
Transformada discreta, a estos mtodos se les
llama
Transformada Rpida 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
/2
0
-p
/2
/2
p
/2
T
T ...
t

La versin muestreada f(k) de f(t) slo puede
tomar un nmero 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
f(k)
1
0.5
0
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
...
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
F(n) c-16
...
13
14
15
16
17
18
19
...
32
...
c-3 c-2 c-1 c0
c1
c2
c3
... c15
Y para graficar el espectro de amplitud:

stem(abs(F))
Obtenindose:
|
F(n)|
0.6
Espectro de Amplitud |F(n)|

Para el tren de pulsos p=1,
T=2
0.4
0.2
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
|F(w)|
para el tren de pulsos, p=1,T=2
0.4
Obteniendo:
0.2
-50
50

Tambin podemos obtener los coeficientes de la
forma trigonomtrica, recordando que:
c n (a n jb n ), c n (a n jb n )
1
2
1
2
Podemos obtener
a 0 c 0 , a n 2 Re(c n ), b 2 Im(c n )
Para el ejemplo se obtiene: a 0=0.5, an=bn=0 (para

n par), adems para n impar:
n
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 funcin par, se
esperaba que bn=0; (el resultado obtenido es
errneo para bn, pero el error disminuye para N
grande): 1
a0
0.5
0
-0.5
Coeficientes bn Coeficientes an
10
20
30

Tarea: Usar el siguiente cdigo para generar 128 puntos
de una funcin peridica con frecuencia fundamental
0=120 (60 hertz) y dos armnicos 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 funcin peridica diferente a la subrayada:

a) Graficar la funcin.
b) Obtener y graficar el espectro de amplitud de la seal
usando la funcin FFT
Medidores Digitales
La FFT ha hecho posible el desarrollo de equipo
electrnico digital con la capacidad de clculo de
espectros de frecuencia para seales 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 energa: permite medicin de 4
seales simultneas (para sistemas trifsicos)

Series de Fourier

Cargado por

Copyright:

Formatos disponibles

Series 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

Series de Fourier

Cargado por

Copyright:

Formatos disponibles

Curso de Titulacin

Modelado y Anlisis de Sistemas Elctricos bajo

Facultad de Ingeniera Elctrica

Pero como se sabe cos(x+2k)=cos(x) para cualquier entero k,

f(t) cos( 3t ) cos( 4t )

Es decir, la relacin 1/ 2 debe ser un nmero racional.

Serie Trigonomtrica de Fourier

f ( t ) 12 a 0 [a n cos(n0 t ) b n sen (n0 t )]

Serie Trigonomtrica de Fourier

Podemos encontrar una manera ms compacta

Serie Trigonomtrica de Fourier

Con lo cual la expresin queda

Serie Trigonomtrica de Fourier

Serie Trigonomtrica de Fourier

Componente fundamental es de la forma:

f(t) 1 cos( 3t ) cos( 4t )

Tiene tantas partes

Ortogonalidad de senos y cosenos

Ortogonalidad de senos y cosenos

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

Ortogonalidad de senos y cosenos

Ortogonalidad de senos y cosenos

2sen (m0T/2) 2sen (m)

Ortogonalidad de senos y cosenos

3.- cos(m0t) Vs. cos(n0t):

cos(m0 t)cos(n0 t)dt

Ortogonalidad de senos y cosenos

sen(m0 t)sen(n0 t)dt

5.- sen(m0t) Vs. cos(n0t):

sen(m0 t)cos(n0 t)dt 0 para cualquier m, n

Ortogonalidad de senos y cosenos

Clculo de los coeficientes de la Serie

Obviamente, el problema se resuelve si sabemos

Clculo de los coeficientes de la Serie

Similarmente, multiplicando por sen(n0t) e

f ( t )sen (n0 t )dt

Similarmente, integrando de T/2 a T/2,

Clculo de los coeficientes de la Serie

Clculo de los coeficientes de la Serie

Solucin: La expresin para f(t) en T/2<t<T/2 es

Clculo de los coeficientes de la Serie

cos(n0 t )dt cos(n0 t )dt

Clculo de los coeficientes de la Serie

Clculo de los coeficientes de la Serie

b n T2 f ( t )sen (n0 t )dt

sen ( n0 t )dt sen ( n0 t )dt

Clculo de los coeficientes de la Serie

En la siguiente figura se muestran: la

Clculo de los coeficientes de la Serie

Componentes de la Serie de Fourier

Clculo de los coeficientes de la Serie

Funciones Pares e Impares

Funciones Pares e Impares

Funciones Pares e Impares

Funciones Pares e Impares

Funciones Pares e Impares

Funciones Pares e Impares

Funciones Pares e Impares

Es una funcin impar, por ello su serie de