Series de Fourier

Series de Fourier.
1
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. 2
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

Solucin.- Si f(t) es peridica se debe cumplir:

Pero como se sabe cos(x+2kt)=cos(x) para cualquier entero k,
entonces para que se cumpla la igualdad se requiere que
T/3=2k
1
t, T/4=2k
2
t
Es decir,
T = 6k
1
t = 8k
2
t
Donde k
1
y k
2
son enteros,
El valor mnimo de T se obtiene con k
1
=4, k
2
=3, es
decir,T=24t
)? cos( ) cos( f(t)
4
t
3
t
+ =
) cos( ) cos( T) f(t
4
T t
3
T t + +
+ = + ) cos( ) cos( f(t)
4
t
3
t
+ = =
Funciones Peridicas
Grfica de la funcin
0 50 100 150 200
-3
-2
-1
0
1
2
3
f(t)=cos(t/3)+cos(t/4)
t
f
(
t
)

24t
T
) cos( ) cos( f(t)
4
t
3
t
+ =
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(e
1
t)+cos(e
2
t).
Para que sea peridica se requiere encontrar dos enteros
m, n tales que
e
1
T= 2tm, e
2
T=2tn
De donde

Es decir, la relacin e
1
/ e
2
debe ser un nmero racional.
n
m
2
1
=
e
e
Funciones Peridicas
Ejemplo: la funcin cos(3t)+cos(t+3)t no es peridica,
ya que no es un nmero racional.
t +
=
e
e
3
3
2
1
0 5 10 15 20 25 30
-2
-1
0
1
2
f(t)=cos(3t)+cos((3+pi)t)
t
f
(
t
)

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)= sen
2
(2tt)
3) f(t)= sen(t)+sen(t+t/2)
4) f(t)= sen(e
1
t)+cos(e
2
t)
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) = a
0
+ a
1
cos(e
0
t)+a
2
cos(2e
0
t)+...
+ b
1
sen(e
0
t)+b
2
sen(2e
0
t)+...
Donde e
0
=2t/T.
Es decir,
] ) t n ( sen b ) t n cos( a [ a ) t ( f
1 n
0 n 0 n 0
2
1
e + e + =

=
Es posible escribir de una manera ligeramente
diferente la Serie de Fourier, si observamos que
el trmino a
n
cos(ne
0
t)+b
n
sen(ne
0
t) se puede
escribir como

Podemos encontrar una manera ms compacta
para expresar estos coeficientes pensando en un
tringulo rectngulo:
|
|
.
|
\
|
e
+
+ e
+
+ ) t n ( sen
b a
b
) t n cos(
b a
a
b a
0
2
n
2
n
n
0
2
n
2
n
n
2
n
2
n

Con lo cual la expresin queda
n
2
n
2
n
n
n
2
n
2
n
n
sen
b a
b
cos
b a
a
u =
+
u =
+
a
n
b
n
2
n
2
n n
b a C + =
u
n
| | ) t n ( sen sen ) t n cos( cos C
0 n 0 n n
e u + e u
| | ) t n cos( C
n 0 n
u e =
Si adems definimos C
0
=a
0
/2, la serie de Fourier
se puede escribir como

As,

y
| |
=
u e + =
1 n
n 0 n 0
) t n cos( C C ) t ( f
2
n
2
n n
b a C + =
|
|
.
|
\
|
= u

n
n
1
n
a
b
tan
Tarea:
Definir adecuadamente los coeficientes C
0
, C
n
y
u
n
, de manera que la serie de Fourier se pueda
escribir como

| |
=
u + e + =
1 n
n 0 n 0
) t n ( sen C C ) t ( f
Componentes y armnicas
As, una funcin peridica f(t) se puede escribir
como la suma de componentes sinusoidales de
diferentes frecuencias e
n
=ne
0
.

A la componente sinusoidal de frecuencia ne
0
:
C
n
cos(ne
0
t+u
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 e
0
=2tf
0
=2t/T se le llama frecuencia
angular fundamental.
A la componente de frecuencia cero C
0
, se le
llama componente de corriente directa (cd) y
corresponde al valor promedio de f(t) en cada
periodo.

Los coeficientes C
n
y los ngulos u
n
son
respectiva-mente las amplitudes y los ngulos
de fase de las armnicas.
Ejemplo: La funcin
Como ya se mostr tiene un periodo T=24t, por lo tanto
su frecuencia fundamental es e
0
=1/12 rad/seg.
Componente fundamental es de la forma:
0*cos(t/12).
Tercer armnico:
cos(3t/12)=cos(t/4)
Cuarto armnico:
Cos(4t/12)=cos(t/3)

0 50 100 150 200
-3
-2
-1
0
1
2
3
f(t)=cos(t/3)+cos(t/4)
t
f
(
t
)

24t
) cos( ) cos( f(t)
4
t
3
t
+ =
Ejemplo: Como puede verse, la funcin anterior tiene
tantas partes positivas como negativas, por lo tanto su
componente de cd es cero, en cambio

0 50 100 150 200
-3
-2
-1
0
1
2
3
f(t)=1+cos(t/3)+cos(t/4)
t
f
(
t
)

24t
) cos( ) cos( 1 f(t)
4
t
3
t
+ + =
Tiene tantas partes
arriba como abajo
de 1 por lo tanto,
su componente de
cd es 1.

Tarea: Cul es la componente fundamental, las
armnicas distintas de cero y la componente de
directa de
a) f(t) = sen
2
t
b) f(t) = cos
2
t ?
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 f
k
(t) son
ortogonales en el intervalo a<t<b si dos
funciones cualesquiera f
m
(t), f
n
(t) de dicho
conjunto cumplen

=
=
=
}
n m para r
n m para 0
dt (t) (t)f f
n
b
a
n m
Ejemplo: las funciones t y t
2
son ortogonales en el
intervalo 1< t <1, ya que

Ejemplo: Las funciones sen t y cos t son ortogonales en
el intervalo
t
/
2
< t <
t
/
2
, ya que
0
4
t
dt t dt tt
1
1 4
1
1
3
1
1
2
= =
}
=
}

0
2
t sen
sentcostdt
2
= =
}
t
t
t
t
Tarea:
Dar un ejemplo de un par de funciones que sean
ortogonales en el intervalo:
a) 0<t<1
b) 0<t<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,cose
0
t, cos2e
0
t, cos3e
0
t,...,sene
0
t,sen2e
0
t,sen3e
0
t,...
(para cualquier valor de e
0
=
2t
/
T
).

Para verificar lo anterior podemos probar por pares:
1.- f(t)=1 Vs. cos(me
0
t):

Ya que m es un entero.
0
m
) (m sen 2
m
T/2) (m sen 2
m
t) (m sen
t)dt cos(m
0 0
0
2 / T
2 / T
0
0
2 / T
2 / T
0
=
e
t
=
e
e
=
e
e
=
}
e
2.- f(t)=1 Vs. sen(me
0
t):

3.- cos(me
0
t) Vs. cos(ne
0
t):
0 T/2)] m ( cos - T/2) m [cos(
m
1
m
t) (m cos
t)dt sen(m
0 0
0
2 / T
2 / T
0
0
2 / T
2 / T
0
= e e
e
=
=
e
e
=
}
e
= =
=
=
}
e e
0 n m para 2 / T
n m para 0
t)dt t)cos(n cos(m
2 / T
2 / T
0 0
4.- sen(me
0
t) Vs. sen(ne
0
t):

5.- sen(me
0
t) Vs. cos(ne
0
t):
n , m cualquier para 0 t)dt t)cos(n sen(m
2 / T
2 / T
0 0
=
}
e e
= =
=
=
}
e e
0 n m para 2 / T
n m para 0
t)dt t)sen(n sen(m
2 / T
2 / T
0 0
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:
sen
2
u = (1-cos2u)
cos
2
u = (1+cos2u)
Clculo de los coeficientes de la Serie
Dada una funcin peridica f(t) cmo se
obtiene su serie de Fourier?

Obviamente, el problema se resuelve si sabemos
como calcular los coeficientes a
0
,a
1
,a
2
,...,b
1
,b
2
,...

Esto se puede resolver considerando la
ortogonalidad de las funciones seno y coseno
comentada anteriormente.
1 n
0 n 0 n 0
2
1
e + e + =

=
Multiplicando ambos miembros por cos(ne
0
t) e
integrando de T/2 a T/2, obtenemos:

Similarmente, multiplicando por sen(ne
0
t) e
integrando de T/2 a T/2, obtenemos:

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

,... 3 , 2 , 1 , 0 n dt ) t n cos( ) t ( f a
2 / T
2 / T
0
T
2
n
=
}
e =
,... 3 , 2 , 1 n dt ) t n ( sen ) t ( f b
2 / T
2 / T
0
T
2
n
=
}
e =
}
=
2 / T
2 / T
T
2
0
dt ) t ( f a
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 t
0
a t
0
+T, con t
0
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:

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

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

T
/
2
T . . .
-1
< <
< <
=
2
T
2
T
t 0 para 1
0 t para 1
) t ( f
Coeficientes a
n
:
}
e =
2 / T
2 / T
0
T
2
n
dt ) t n cos( ) t ( f a
(
}
e +
}
e =
2 / T
0
0
0
2 / T
0
T
2
dt ) t n cos( dt ) t n cos(
(
(
e
e
+ e
e
=
0
2 / T
0
0
2 / T
0
0
0
T
2
) t n ( sen
n
1
) t n ( sen
n
1
0 n para 0 = =
Coeficiente a
0
:
}
=
2 / T
2 / T
T
2
0
dt ) t ( f a
(
}
+
}
=
2 / T
0
0
2 / T
T
2
dt dt
(
(
+ =
0
2 / T
2 / T
0
T
2
t t
0 =
Coeficientes b
n
:
}
e =
2 / T
2 / T
0
T
2
n
dt ) t n ( sen ) t ( f b
(
}
e +
}
e =
2 / T
0
0
0
2 / T
0
T
2
dt ) t n ( sen dt ) t n ( sen
(
(
e
e
e
e
=
0
2 / T
0
0
2 / T
0
0
0
T
2
) t n cos(
n
1
) t n cos(
n
1
| | ) 1 ) n (cos( )) n cos( 1 (
n
1
t t
t
=
| | 0 n para ) ) 1 ( 1
n
2
n
=
t
=
Serie de Fourier: Finalmente la Serie de Fourier
queda como

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 e
0
=t, es decir,
T=2:
| | ... ) t 5 ( sen ) t 3 ( sen ) t ( sen
4
) t ( f
0
5
1
0
3
1
0
+ e + e + e
t
=
-1 -0.5 0 0.5 1
-1.5
-1
-0.5
0
0.5
1
1.5
Componentes de la Serie de Fourier
t
C
o
m
p
o
n
e
n
t
e
s

Suma
fundamental
tercer armnico
quinto armnico
septimo armnico
Tarea: Encontrar la serie de Fourier para la
siguiente seal senoidal rectificada de media
onda de periodo 2t.

-6 -4 -2 0 2 4 6
-0.2
0
0.2
0.4
0.6
0.8
1
Senoidal rectificada de media onda
t
f
(
t
)

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)

t 2t
f(t)
t
t 2t
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)

t 2t
f(t)
t
t 2t
Ejemplo: Las siguientes funciones son pares o
impares?
f(t) = t+1/t
g(t) = 1/(t
2
+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/(t
2
+1)=g(t), por lo
tanto g(t) es funcin par.
Ejemplo: La funcin h(t)=f(1+t
2
) es par o
impar?, donde f es una funcin arbitraria.
Solucin:
Sea g(t)= 1+t
2
, Entonces h(t)=f(g(t))
Por lo tanto h(-t) = f(g(-t)),
Pero g(-t)=1+(-t)
2
= 1+t
2
=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+t
2
)
h(t) = exp(1+t
2
)+5/ (1+t
2
)
h(t) = cos (2+t
2
)+1
h(t) = (10+t
2
)-(1+t
2
)1/2
etc...
Ya que todas tienen la forma f(1+t
2
)
Como la funcin sen(ne
0
t) es una funcin impar
para todo n=0 y la funcin cos(ne
0
t) 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 b
n
= 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:

Es una funcin impar, por ello su serie de
Fourier no contiene trminos coseno:
1
f(t)
t
. . . -T
/
2
0

T
/
2
T . . .
-1
| | ... ) t 5 ( sen ) t 3 ( sen ) t ( sen
4
) t ( f
0
5
1
0
3
1
0
+ e + e + e
t
=
Simetra de Media Onda
Una funcin periodica de periodo T se dice
simtrica de media onda, si cumple la propiedad

Es decir, si en su grfica las partes negativas son
un reflejo de las positivas pero desplazadas
medio periodo:
) t ( f ) T t ( f
2
1
= +

f(t)
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)
t
Ejemplo: Funcin con simetra par de cuarto de
onda:

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

f(t)
t
Simetras y Coeficientes de Fourier
Simetra Coeficientes
Funciones
en la serie
Ninguna
Senos y
cosenos
Par b
n
=0
nicamente
cosenos
Impar a
n
=0
nicamente
senos
media
onda
Senos y
cosenos
impares
}
e =
2 /
0
0
4
) cos( ) (
T
T n
dt t n t f a
}
e =
2 /
0
0
4
) ( ) (
T
T n
dt t n sen t f b
e
=
}
impar n dt t n t f
par n
a
T
T
n
2 /
0
0
4
) cos( ) (
0
e
=
}
impar n dt t n sen t f
par n
b
T
T
n
2 /
0
0
4
) ( ) (
0
}
e =
2 /
2 /
0
2
) cos( ) (
T
T
T n
dt t n t f a
}
e =
2 /
2 /
0
2
) ( ) (
T
T
T n
dt t n sen t f b
Simetra Coeficientes
Funciones
en la serie
Ninguna
Senos y
cosenos
de
onda par
a
n
=0 (n par)
b
n
=0
Slo
cosenos
impares
de
onda
impar
a
n
=0
b
n
=0 (n par)
Slo
senos
impares
}
e =
2 /
2 /
0
2
) cos( ) (
T
T
T n
dt t n t f a
}
e =
2 /
2 /
0
2
) ( ) (
T
T
T n
dt t n sen t f b
) (
) cos( ) (
4 /
0
0
8
impar n
dt t n t f a
T
T n
}
e =
) (
) ( ) (
4 /
0
0
8
impar n
dt t n sen t f b
T
T n
}
e =
Por ejemplo, la seal cuadrada, ya analizada en
un ejemplo previo:

Es una funcin con simetra de de onda impar,
por ello su serie de Fourier slo contiene
trminos seno de frecuencia impar:
1
f(t)
t
. . . -T
/
2
0

T
/
2
T . . .
-1
| | ... ) t 5 ( sen ) t 3 ( sen ) t ( sen
4
) t ( f
0
5
1
0
3
1
0
+ e + e + e
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 -0.5 0 0.5 1
-1.5
-1
-0.5
0
0.5
1
1.5
Serie con 1 armnico
Fenmeno de Gibbs
-1 -0.5 0 0.5 1
-1.5
-1
-0.5
0
0.5
1
1.5
Serie con 3 armnicos
Fenmeno de Gibbs
-1 -0.5 0 0.5 1
-1.5
-1
-0.5
0
0.5
1
1.5
Fenmeno de Gibbs
-1 -0.5 0 0.5 1
-1.5
-1
-0.5
0
0.5
1
1.5
Fenmeno de Gibbs
-1 -0.5 0 0.5 1
-1.5
-1
-0.5
0
0.5
1
1.5
Fenmeno de Gibbs
-1 -0.5 0 0.5 1
-1.5
-1
-0.5
0
0.5
1
1.5
Fenmeno de Gibbs
-1 -0.5 0 0.5 1
-1.5
-1
-0.5
0
0.5
1
1.5
Forma Compleja de la Serie de Fourier
Consideremos la serie de Fourier para una
funcin periodica f(t), con periodo T=2t/e
0
.

Es posible obtener una forma alternativa usando
las frmulas de Euler:

Donde
1 n
0 n 0 n 0
2
1
e + e + =

=
) e e ( ) t n ( sen
) e e ( ) t n cos(
t j n t j n
j 2
1
0
t j n t j n
2
1
0
0 0
0 0
e e
e e
= e
+ = e
1 j =
Sustituyendo

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

Y definiendo:

Lo cual es congruente con la frmula para b
n
, ya
que b
-n
=-b
n
, ya que la funcin seno es impar.
] ) e e ( b ) e e ( a [ a ) t ( f
1 n
t j n t j n
j 2
1
n
t j n t j n
2
1
n 0 2
1 0 0 0 0
=
e e e e
+ + + =
] e ) jb a ( e ) jb a ( [ a ) t ( f
1 n
t j n
n n 2
1
t j n
n n 2
1
0 2
1 0 0
=
e e
+ + + =
) jb a ( c ), jb a ( c , a c
n n 2
1
n n n 2
1
n 0 2
1
0
+ = = =

La serie se puede escribir como

O bien,

Es decir,

) e c e c ( c ) t ( f
1 n
t j n
n
t j n
n 0
0 0
=
e
e
+ + =

=
e
=
e
+ + =
1 n
t j n
n
1 n
t j n
n 0
0 0
e c e c c ) t ( f
=
e
=
n
t j n
n
0
e c ) t ( f
A la expresin obtenida

Se le llama forma compleja de la serie de
Fourier y sus coeficientes c
n
pueden obtenerse a
partir de los coeficientes a
n
, b
n
como ya se dijo, o
bien:

Para n=0, 1, 2, 3, ...
}
e
=
T
0
t j n
T
1
n
dt e ) t ( f c
0
=
e
=
n
t j n
n
0
e c ) t ( f
Los coeficientes c
n
son nmeros complejos, y
tambin se pueden escribir en forma polar:

Obviamente,

Donde ,
Para todo n=0,

Para n=0, c
0
es un nmero real:
n
j
n n
e c c
|
=
n
j
n
*
n n
e c c c
|
= =
2
n
2
n 2
1
n
b a c + =
)
a
b
arctan(
n
n
n
= |
0 2
1
0
a c =
Ejemplo. Encontrar la forma compleja de la
serie de Fourier para la funcin ya tratada:

Solucin 1. Como ya se calcularon los
coeficientes de la forma trigonomtrica (a
n
y b
n
):
a
n
=0 para todo n
y
1
f(t)
t
. . . -T
/
2
0

T
/
2
T . . .
-1
n todo para ] ) 1 ( 1 [ b
n
n
2
n
=
t
Podemos calcular los coeficientes c
n
de:

Entonces la Serie Compleja de Fourier queda
] ) 1 ( 1 [ j ] jb a [ c
n
n
2
2
1
n n 2
1
n
= =
t
] ) 1 ( 1 [ j c
n
n
1
n
=
t
...) e e e
e e e (... j ) t ( f
t 5 j
5
1
t 3 j
3
1
t j
t j t 3 j
3
1
t 5 j
5
1 2
0 0 0
0 0 0

+ + + =
e e e
e e e
t
Solucin 2. Tambin podemos calcular los
coeficientes c
n
mediante la integral
}
e
=
T
0
t j n
T
1
n
dt e ) t ( f c
0
) dt e dt e (
T
2 / T
t j n
2 / T
0
t j n
T
1 0 0
} }
e e
+ =
) e e (
2 / T
T
t j n
j n
1
0
2 / T
t j n
j n
1
T
1 0
o
0
o
e
e
e
e
=
)] e e ( ) 1 e [(
2 / T jn T jn 2 / T jn
T jn
1 0 0 0
o
e e e
e
=
Como e
0
T=2t y adems

Lo cual coincide con el resultado ya obtenido.
u u =
u
jsen cos e
j
)] ) 1 ( 1 ( ) 1 ) 1 [( c
n n
T jn
1
n
o
=
e
] ) 1 ( 1 [ j
n
T n
2
o
=
e
] ) 1 ( 1 [ j
n
n
1
=
t
Tarea: Calcular los coeficientes c
n
para la
siguiente funcin de periodo 2t.
a) A partir de los coeficientes a
n
,b
n
b) Directamente de la integral
-6 -4 -2 0 2 4 6
-0.2
0
0.2
0.4
0.6
0.8
1
Senoidal rectificada de media onda
t
f
(
t
)

Espectros de Frecuencia Discreta
A la grfica de la magnitud de los coeficientes c
n

contra la frecuencia angular e 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 c
n
contra e, se le llama el espectro
de fase de f(t).

Como n slo toma valores enteros, la frecuencia
angular e=ne
0
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
c
n
.

Por ello, los coeficientes c
n
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:

Se encontr que

Por lo tanto,
1
f(t)
t
. . . -T
/
2
0

T
/
2
T . . .
-1
] ) 1 ( 1 [ j c
n
n
1
n
=
t
] ) 1 ( 1 [ c
n
n
1
n
=
t
El espectro de amplitud se muestra a continuacin

Observacin: El eje horizontal es un eje de frecuencia,
(n=nmero de armnico = mltiplo de e
0
).
-30 -20 -10 0 10 20 30
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
Espectro de Amplitud de f(t)
n
|
C
n

|

Frecuencia negativa (?) Frecuencia
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)
1
f(t)
t
h=Altura
promedio
}
=
T
0
dt ) t ( f Area
T
Area=Th
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

Si f(t) es peridica, tambin lo ser [f(t)]
2
y el
promedio en un periodo ser el promedio en
cualquier otro periodo.
}
2 / T
2 / T
2
T
1
dt )] t ( f [
El teorema de Parseval nos permite calcular la
integral de [f(t)]
2
mediante los coeficientes com-
plejos c
n
de Fourier de la funcin peridica f(t):

O bien, en trminos de los coeficientes a
n
, b
n
:
}

=
=
n
2
n
2 / T
2 / T
2
T
1
c dt )] t ( f [
}

=
+ + =
1 n
2
n
2
n 2
1
2
0 4
1
2 / T
2 / T
2
T
1
) b a ( a dt )] t ( f [
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,

Donde C
n
es la amplitud del armnico n-simo y
C
0
es la componente de directa.
}

=
+ =
1 n
2
n
2
0
2 / T
2 / T
2
T
1
2
C
C dt )] t ( f [
Para aclarar el resultado anterior es conveniente
encontrar la relacin entre los coeficientes
complejos c
n
de la serie

Y los coeficientes reales C
n
de la serie

Donde C
n
es la amplitud del armnico n-simo y
C
0
es la componente de directa.
=
e
=
n
t j n
n
0
e c ) t ( f
| |
=
u e + =
1 n
n 0 n 0
) t n cos( C C ) t ( f
Por un lado

Mientras que

Entonces, Por lo tanto,

Adems, para el armnico
Su valor rms es , por lo tanto su valor
cuadrtico medio es

Para la componente de directa C
0
, su valor rms
es |C
0
|, por lo tanto su valor cuadrtico medio
ser |C
0
|
2
.
, b a C
2
n
2
n n
+ =
2
n
2
n 2
1
n
b a c + =
n 2
1
n
C c =
2
n 4
1
2
n
C c =
| | ) t n cos( C ) t ( f
n 0 n n
u e =
2 / C
n
2 / C
2
n
Ejemplo. Calcular el valor cuadrtico medio de
la funcin f(t):

Solucin.
Del teorema de Parseval

y del ejemplo anterior

sustituyendo
1
f(t)
t
. . . -T
/
2
0

T
/
2
T . . .
-1
}

=
=
n
2
n
2 / T
2 / T
2
T
1
c dt )] t ( f [
] ) 1 ( 1 [ c
n
n
1
n
=
t
(
+ + + +
t
=
=
...
49
1
25
1
9
1
1
8
c
2
n
2
n
La serie numrica obtenida converge a

Por lo tanto,

Como era de esperarse.
2337 . 1 ...
49
1
25
1
9
1
1 = + + + +
1 ) 2337 . 1 (
8
c dt )] t ( f [
2
n
2
n
2 / T
2 / T
2
T
1
=
t
= =

}

=
Tarea.

Calcular el valor cuadrtico medio para la seal
senoidal rectificada de media onda de periodo
2t.
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 -T
/
2
0

T
/
2
T . . .
p
-p
/
2

p
/
2
< <
< <
< <
=

2
T
2
p
2
p
2
p
2
p
2
T
t 0
t 1
t 0
) t ( f
Los coeficientes de la Serie Compleja de Fourier
en este caso resultan puramente reales:

El espectro de frecuencia correspondiente lo
obtenemos (en este caso) graficando c
n
contra
e=ne
0
.
) n (
) n ( sen
) ( c
2
p
0
2
p
0
T
p
n
e
e
=
Espectro del tren de pulsos para p=1, T=2
-60 -40 -20 0 20 40 60
-0.2
0
0.2
0.4
0.6
w=nw
0
c
n
Si el periodo del tren de pulsos aumenta:
-20 -10 0 10 20
0
0.5
1
1.5
p=1, T=2
t
f
(
t
)

t
-20 -10 0 10 20
0
0.5
1
1.5
p=1, T=5
f
(
t
)

-20 -10 0 10 20
0
0.5
1
1.5
p=1, T=10
t
f
(
t
)

-20 -10 0 10 20
0
0.5
1
1.5
p=1, T=20
t
f
(
t
)

En el lmite cuando T, la funcin deja de ser
peridica:

Qu pasa con los coeficientes de la serie de
Fourier?
-20 -10 0 10 20
0
0.5
1
1.5
p=1, T=
t
f
(
t
)

-50 0 50
-0.1
0
0.1
0.2
0.3
p=1, T=5
-50 0 50
-0.05
0
0.05
0.1
0.15
p=1, T=10
-50 0 50
-0.02
0
0.02
0.04
0.06
p=1, T=20
-50 0 50
-0.2
0
0.2
0.4
0.6
p=1, T=2
e=ne
0
c
n

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 ne
0
, sino como una
funcin continua de la frecuencia e.

As, la serie

Al cambiar la variable discreta ne
0
(cuando
T) por la variable continua e, se transforma
en una integral de la siguiente manera:
=
e
=
n
t j n
n
0
e c ) t ( f
Como

La serie queda

O bien,

cuando T, ne
0
e y e
0
de y la sumatoria
se convierte en
=
e
e
(
=
n
t j n
2 / T
2 / T
t j n
T
1 0 0
e dt e ) t ( f ) t ( f
}
e
=
2 / T
2 / T
t j n
T
1
n
dt e ) t ( f c
0
=
e
e
t
e
(
=
n
t j n
0
2 / T
2 / T
t j n
2
1 0 0
e dt e ) t ( f ) t ( f
} }

e
t
e
(
= d e dt e ) t ( f ) t ( f
t j t j
2
1
Es decir,

Donde

Estas expresiones nos permiten calcular la
expresin F(e) (dominio de la frecuencia) a
partir de f(t) (dominio del tiempo) y viceversa
}

e
t
e e = d e ) ( F ) t ( f
t j
2
1
}

e
= e dt e ) t ( f ) ( F
t j
Identidad
de Fourier
Transformada
De Fourier
Notacin: A la funcin F(e) se le llama
transformada de Fourier de f(t) y se denota por
F, es decir

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
}

e
t
e e = = e d e ) ( F ) t ( f )] ( F [
t j
2
1
1
F
}

e
= e = dt e ) t ( f ) ( F )] t ( f [
t j
F
Ejemplo. Calcular F(w) para el pulso
rectangular f(t) siguiente

Solucin. La expresin en el dominio del tiempo
de la funcin es
-p
/
2

0

p
/
2
1
f(t)
t
<
< <
<
=

t 0
t 1
t 0
) t ( f
2
p
2
p
2
p
2
p

Integrando

Usando la frmula de Euler

Obsrvese que el resultado es igual al obtenido
para cn cuando T , pero multiplicado por T.
} }

e

e
= = e
2 / p
2 / p
t j t j
dt e dt e ) t ( f ) ( F
2 / p
2 / p
t j
j
1
e
e
e
=
) e e (
2 / p j 2 / p j
j
1
e e
e
=
2 / p
) 2 / p ( sen
p ) ( F
e
e
= e
En forma Grfica
-50 0 50
0
0.5
1
F(w) con p=1
w
F
(
w
)
Tarea. Calcular la Transformada de Fourier de
la funcin escaln unitario u(t):

Graficar U(e)=F[u(t)]
Qu rango de frecuencias contiene U(e)?
Cul es la frecuencia predominante?
u(t)
0
1
t
La Transformada Rpida de Fourier
Cuando la funcin f(t) est dada por una lista de N
valores f(t
1
), f(t
2
), ...f(t
N
) se dice que est discretizada o
muestreada, entonces la integral que define la
Transformada de Fourier:

Se convierte en la sumatoria

(Donde k es la frecuencia discreta)
Llamada Transformada Discreta de Fourier

}

e
= e dt e ) t ( f ) ( F
t j
N n 1 para , e ) t ( f ) n ( F
N
1 k
) 1 k ( j
k
N
n 2
s s =
=

t
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 c
n
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.
1
f(t)
t
. . . -T -T
/
2
0

T
/
2
T . . .
p
-p
/
2

p
/
2
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):
0 1 2
0
0.5
1
1.5
32 muestras de f(t), de 0 a T
k
f
(
k
)

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) c
0
c
1
c
2
c
3
... c
15
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:

Y para graficar el espectro de amplitud:
stem(abs(F))
Obtenindose:
n 1 ... 13 14 15 16 17 18 19 ... 32
F(n) c
-16
...
c
-3
c
-2
c
-1
c
0
c
1
c
2
c
3
... c
15

Si deseamos una escala horizontal en unidades
de frecuencia (rad/seg):
0 10 20 30
0
0.2
0.4
0.6
Para el tren de pulsos p=1,
T=2
n
|
F
(
n
)
|

Espectro de Amplitud |F(n)|
w0=2*pi/T;
n=-16:15;
w=n*w0;
Stem(w,abs(F))

Obteniendo:
-50 0 50
0
0.2
0.4
0.6
para el tren de pulsos, p=1,T=2
w
|
F
(
w
)
|

Espectro de Amplitud |F(n)|
Tambin podemos obtener los coeficientes de la
forma trigonomtrica, recordando que:

Podemos obtener

Para el ejemplo se obtiene: a
0
=0.5, a
n
=b
n
=0 (para
n par), adems para n impar:
) jb a ( c ), jb a ( c
n n 2
1
n n n 2
1
n
+ = =

) c Im( 2 b ), c Re( 2 a , c a
n n n 0 0
= = =
n 1 3 5 7 9 11 13 15
a
n
0.6346 -0.2060 0.1169 -0.0762 0.0513 -0.0334 0.0190 -0.0062
b
n
-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 b
n
=0; (el resultado obtenido es
errneo para b
n
, pero el error disminuye para N
grande):
0 10 20 30
-0.5
0
0.5
1
Coeficientes b
n
Coeficientes a
n
a
0
Tarea: Usar el siguiente cdigo para generar 128
puntos de una funcin peridica con frecuencia
fundamental e
0
=120t (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

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

Series de Fourier.

También podría gustarte