T.3 Transformada Dscreta de Fourier
T.3 Transformada Dscreta de Fourier
T.3 Transformada Dscreta de Fourier
DE FOURIER
Seales y sistemas digitales
ENC310
Ing Cristina Vilardell
Indice
Definicin e interpretacin
La DFT como transformacin lineal
Propiedades de la DFT
Enventanado de Secuencias
Algoritmos de Clculo de la DFT
Forma Directa
FFT
Otros Algoritmos
Filtrado lineal mediante la DFT
Filtrado de secuencias de larga duracin
Representacin frecuencial de una
seal
Esta representacin se hace a
travs de la transformada de
Fourier de la secuencia .
Porque????
La representacin frecuencial es
mas fcil de medir en laboratorio
Seales con diferentes
componentes frecuenciales
pueden compartir el mismo
soporte de transmisin sin
interferirse
Consiste en la descomposicion de x(n) en una
suma ponderada de exponenciales complejas
armonicas.
El conjunto de pesos y de fases de estas
exponenciales armonicas forman el espectro de
la seal.
Espectro: contenido frecuencial de una seal.
Si x(n) es periodica, esta suma sera una serie de
Fourier de tiempo discreto.
5. Si x(n) es una seal de energa, esta suma ser
una trasformada de Fourier de tiempo discreto.
Definicin e interpretacin
Antes de definir la DFT,
analizaremos primero la
Transformada de Fourier en
tiempo discreto (DTFT) y las
series discretas.
La DTFT describe el espectro
de seales discretas.
Deduciremos la DFT a partir
de la convolucin discreta
SERIES DE FOURIER PARA SEALES
DISCRETAS EN TIEMPO:
Las mismas motivaciones que nos condujeron al
desarrollo de las series y transformadas de Fourier
para seales contnuas, siguen siendo vlidas para
seales discretas en tiempo: Entre otras ventajas, si
una seal se puede descomponer en suma de
exponenciales complejas, cuando esta pase
por un sistema LTI, la salida
se calcular fcilmente
utilizando el principio de
superposicin.
SERIES DE FOURIER PARA SEALES
DISCRETAS EN TIEMPO:
Comencemos analizando el caso de seales discretas
peridicas. Sea x(n) una seal peridica con perodo N (5
en nuestro caso)
Si queremos representar x(n) en base a
funciones del tipo e
jnk(2/N)
lo primero que
debemos comprender es que solo existen
N distintas ( 0 < k <N-1).
Como se ve estas dos exponenciales son idnticas. Por lo tanto para
representar x(n) necesitaremos N exponenciales que podran ser aquellas que
se generan tomando 0 < k < N-1 1 < k < N A < k < A+N-1.
SERIES DE FOURIER PARA SEALES
DISCRETAS EN TIEMPO:
Esto se representar como k mod N.
Donde k mod N significa que la sumatoria
debe barrer , por ejemplo , desde n=0 hasta
N-1 o desde n=1 hasta N o cualquier otro
barrido de longitud N.
SERIES DE FOURIER PARA SEALES
DISCRETAS EN TIEMPO:
El valor de Ck se obtiene de la siguiente
forma:
Multiplicamos ambos lados de la ecuacin
anterior por y sumamos de acuerdo al
ndice n en un intervalo de N puntos:
Ahora bien la sumatoria
SERIES DE FOURIER PARA SEALES
DISCRETAS EN TIEMPO:
Se puede determinar
recordando que
Por lo tanto:
SERIES DE FOURIER PARA SEALES
DISCRETAS EN TIEMPO:
De esta forma se despeja Ck quedando igual a :
En conclusin : Cualquier secuencia discreta peridica x(n)
de perodo N, puede construrse en
base a N exponenciales de peso Ck tal y como sigue
SERIES DE FOURIER PARA SEALES
DISCRETAS . EJEMPLO:
Encuentre la serie de la seal x(n) mostrada en la figura:
Hagamos el cambio n = n + N1
SERIES DE FOURIER PARA SEALES
DISCRETAS .EJEMPLO:
DTFT Definicin e
interpretacin
Recordemos la definicin de convolucin discreta:
Se define la DTFT de una seal discreta x[n] como
Dualidad entre las series de
Fourier y la DTFT
Tenemos una seal peridica continua xp(t).
Mediante las series de Fourier transformamos
esa seal peridica continua en una funcin
aperidica y discreta (los coeficientes
espectrales X
S
[k]).
Dualidad entre las series de
Fourier y la DTFT
De una manera dual, podemos intercambiar tiempo
y frecuencia de forma donde S
F
=1/ts .
Ahora tenemos una seal aperidica discreta xs [k] y la
transformamos en una seal peridica continua
(Xp(f)) mediante la DTFT.
Series de Fourier y la DTFT
Todas las propiedades que se vieron para las
series de Fourier tienen su correspondientes
equivalencias en la DTFT.
Ejemplo : DTFT de la secuencia x[n]=[n]
Si tenemos una secuencia x[n]={1,0,3,-2}, a
partir de la anterior ecuacin y aplicando la
propiedad del desplazamiento,
Dualidad entre las series de
Fourier y la DTFT
En resumen : el comportamiento dual entre las
series de Fourier y la DTFT se manifiesta en lo
siguiente :
En las series de Fourier parto de una seal x(t),
temporal, continua y peridica (periodo T) y
obtengo los coeficientes X[k], que es una funcin
de la frecuencia, aperidica y discreta con una
distancia entre dos valores consecutivos de
f
0
=1/T.
Mas ejemplos DTFT
Ms ejemplos de la DTFT
DTFT (Discrete Time Fourier
Transform) versus DFT (Discrete
Fourier Transform )
La Transformada de Fourier (DTFT) es una funcin de
variable continua:
Imposible su uso mediante herramientas Informtica
Es preciso desarrollar una herramienta numrica: la DFT
Transformada Discreta de
Fourier
Si utilizamos la DTFT a la hora de realizar
operaciones tenemos los mismos problemas que
en las series de Fourier ya que seguimos
tratando con seales continuas o con series de
datos de longitud infinita.
La electrnica nos obliga a trabajar con un
nmero finito de datos discretos que adems
tienen una precisin finita.
DFT
Se muestrea el espectro de la seal x[n]:
X[k] recibe el nombre de Transformada Discreta
de Fourier (DFT)
Es necesario definir el valor de M
Repaso
Repaso
Ejemplos
RELACIONES DE LOS DSF
RELACION DSF
Los coeficientes del DSF de la extensin peridica
coinciden con las muestras del espectro de x[n]:
RELACION DSF
El muestreo X(e
j
) da como resultado la
extensin peridica de x[n] que en trminos de
seal finita
Se desea
Hay que evitar el solapamiento temporal
Transformada Discreta de Fourier
De lo que se trata es de conseguir discretizar
las variables continuas y de limitar el nmeros
de muestras en los dos dominios (temporal y
frecuencial).
Esto nos lleva a definir las series discretas de
Fourier y la Transformada Discreta de Fourier
(DFT).
Como lo hacemos??
Se muestrea el espectro de la seal x[n]:
X[k] recibe el nombre de Transformada Discreta de Fourier
(DFT)
Es necesario definir el valor de M
Muestreo de la respuesta en frecuencia
Relacion con el DSF
Los coeficientes del DSF de la extensin peridica coinciden con las
muestras del espectro de x[n]:
Relacin con el DSF
Formulacin
El muestreo X(ej) da como resultado la extensin peridica de x[n] que
en trminos de seal finita
Se desea
Hay que evitar el solapamiento temporal
La formulacin para la DFT:
En la prctica esto se implementa mediante un algoritmo FFT
Definicin e interpretacin
Solapamiento temporal
Definicin e interpretacin
Es
equivalente
a aadir
ceros
El comportamiento dual entre las series de
Fourier y la DTFT se manifiesta en lo
siguiente :
En las series de Fourier parto de una
seal
x(t), temporal, continua y peri-
dica (periodoT) y obtengo los coeficientes
X[k] , que es una funcin de la frecuencia,
aperidica y discreta con una distancia
entre dos valores
consecutivos de f0=1/T
En la DTFT parto de una seal discreta en
el tiempo x[n], con periodo de muestreo
ts=1/fs y aperidica y obtengo una funcin
X(f) , que es funcin continua de la
frecuencia y peridica con periodo
f
s
.
Todas las propiedades que se vieron para
las series de Fourier tienen su
correspondientes equivalencias en la
DTFT.
Ejemplo : DTFT de la secuencia
x[n]=
[n]
:
Si tenemos una secuencia
x[n]={1,0,3,-2}
, a partir de la anterior ecuacin
y aplicando la propiedad del
desplazamiento,
Se considera la base para las
exponenciales Complejas
Las frmulas se pueden expresar
Se definen los vectores x
N
y X
N
de N
puntos y la matriz cuadrada W
N
de
orden N.
Interpretacin matricial de la DFT
( ) G z
La frmula de la DFT:
Se puede rescribir
Interpretacin matricial de la DFT
De forma compacta:
Siendo W
N
matriz de transformacin lineal que es simtrica.
La frmula de inversin (IDFT):
Ahora bien, la IDFT se puede obtener mediante una forma alternativa si
volvemos a la expresin original.
DFT como transformacin lineal
La IDFT:
Rescrita en forma matricial:
DFT como transformacin lineal
De forma compacta, la IDFT:
Por lo que se deduce que:
Esto implica:
Donde I
N
es la matriz identidad de orden N. Se concluye que la
matriz de transformacin es ortogonal
Propiedades de la DFT
Desplazamiento Desplazamiento
circular
1
1
()
1
Uz
z
para
4) Serie geomtrica
() 0, 1, 2, 3, ... , . ==
k
fkak n
1
() ( )
k
Zafk Faz
(=
1
1
()
1
k
az
fka
az
Multiplicando y dividiendo
por a
()
z
Fz za
za
=
=
=
1 z
0
()
=
=
k
k
Fzkz
Para una secuencia geomtrica se tiene:
0
=
=
kk
k
z
az
za
Derivando con respecto a z:
22
0
()
()()
===
kk
k
ddzzaza
az
dzdzzazaza
1
2
0
( )
kk
k
a
kaz
za
=
=
( ) ( 2) y k u k = ( 1) u k ( ) u k
1
Z
++++
==
++++
Factorizando
1
01 1
12
1
.......... ()
()
()()..........()
()
mm
mm
n
n
i
i
bzbz bzbqz
Fz
zpzpzp
zp
=
++++
= =
H
Cuando todos los polos de en la ecuacin son diferentes
1
0 1 1
1 2
..........
()
( )( )...........( )
mm
mm
n
bzbz bzb
Fz
zpzp zp
++++
=
01
12
.........
n n
n
zzz
dddd
zpzpzp
=++++
=
i
i
i zp
zp
d Fz
z
1
0 1 1
1 2
...
( )( )... ( )
++++
=
mm
mm
n
i i i n
bzbz bzb
d
zpppppp
La secuencia resulta:
||
1
0 1122
() ()() .......
k k k
nn
fkZFzdkdpdpdp o
==++++
Con polos mltiples queda
1 2 1
1
0 1 1
2
..........
()
()()..........()
mm
mm
n n n
i n
bzbz bzb
Fz
zpzp zp
++++
=
La expansin de F(z), en este caso, tiene la forma:
1
11
2
0122
11 1
() .....
()()
=++++
n
nn
zzz
Fzdddd
zpzpzp
2
22
2
122
222
.....
()()
++++
n
nn
zzz
eee
zpzpzp
..... +
1
11
2
122
111
.....
()()
++++
n
nn
zzz
rrr
zpzpzp
Arreglar tamao en texto y frmulas
TABLA 8.II
PARES DE TRANSFORMADAS Z PARA RACES MLTIPLES
( ) F z ( ) f k
0 k >
z
z a
k
a
2
2
( )
z
z a
( 1)
k
k a +
3
3
( )
z
z a
(1)(2)
2!
k
kk
a
++
4
4
( )
z
z a
( 1)( 2)( 3)
3!
k
k k k
a
++ +
para
8.6 FUNCIN DE TRANSFERENCIA DE SISTEMAS DISCRETOS
El concepto de funcin de transferencia ; la cual se define
como la relacin de la transformada Z de la salida, , de un sistema
entre la transformada Z de su entrada,
( ) H z
( ) Y z
( ) U z
()
()
()
Yz
Hz
Uz
=
Arreglar tamao en texto y frmulas
La expresin general aplicable a la funcin de transferencia es:
1
011
1
011
... ()
()
() ...
++++
==
++++
mm
mm
nn
nn
bzbzbzb qz
Hz
pzazazaza
1
1 2
1 2
1
()
()()... ()
()()... ()
()
=
=
H
= =
H
m
j
j
m
n
n
i
i
zc
zczczc
zpzpzp
zp
Algunos sistemas tipicos:
1. Sistema en cascada
1 2
( ) ( ) y k u k =
1
( ) u k
2
( ) y k
1
( ) h k
2
( ) h k
En el domino de Z:
( ) u z
2
( ) y z
1 2
( ) ( ) ( ) = H z H z H z
2. Sistema inverso
1
( ) H Z
2
( ) H Z
( ) Y Z ( ) u z
() () YzUz =
1
2
1
()
()
Hz
Hz
=
1 2
() () ()1 Hz HzHz = =
La convolucin en este caso resulta:
00
()()()()()()
i i
ykhikiikik o
==
===
1
1 2
()
()
1 () ()
=
+
Hz
Yz
HzHz
8.7 ESTABILIDAD DE SISTEMAS DISCRETOS
Un sistema discreto es estable cuando produce una salida acotada al
aplicrsele una entrada acotada
Los sistemas discretos estables se caracterizan porque todos sus polos
se ubican en el plano complejo z , dentro de un crculo centrado en el
origen de radio unitario
8.7.1 POLOS DE H(z) Y RESPUESTA TRANSITORIA
La localizacin de los polos de H(z) en el plano z permite caracterizar
efectivamente las propiedades de la respuesta para un sistema discreto
lineal.
A.- Polo real en .
z a =
La respuesta caracterstica es de la forma
cos() u| +
k
Ark
Donde A y son constantes obtenidas de la expansin en fracciones
parciales y:
2 2
rab =+
1
tan
b
a
u
=
x x x x x x
(4)
Im z
(2)
(1)
(3)
(5)
(6)
Re z
Cambiar
dibujo
Casos:
1- . Sistema inestable.La respuesta a impulso es una
oscilacin creciente en magnitud.
2- . Sistema inestable.La respuesta es una oscilacin
parecida a un senoide con magnitud constante.
3-. Sistema estable. El resultado es una oscilacin parecida
a una senoide decreciente en magnitud.
2 2
1 ab +>
2 2
1 += ab
2 2
1 ab +<
x
Im z
Re z
(2)
(1)
(3)
(2)
x (3)
x (3)
(2)
x (1)
x (1)
Cambiar
dibujo
8.7.2. POLOS DOMINANTES
Son los que tienen una influencia de mayor importancia sobre la respuesta
transitoria.Son los polos que estn ms cerca del circulo unitario. Ej p1 y p2.
m z
3
p
4
p
Re z
x
x
1
p
x
x
2
p
8.8 RESPUESTA SENOIDAL PERMANENTE DE SISTEMAS
LINEALES (FILTROS DIGITALES)
Se asume que la entrada a un sistema es una seal senoidal
pura.
( ) u t ( ) u k
T
Sistema discreto
lineal H(z)
( ) Y k
()
1
()e = utsenT
()
1
()e = utsenkT
()
1 1
1
2
()
()()
ee
e
jTjT
senT
uz
zeze
1
0 1
1 2
...
()
()()... ()
+++
=
mm
m
n
bzbz b
Hz
zpzpzp
1
1 p <
Si consideramos que todos los polos son distintos
()()() YzHzHz =
11
012
12
() ...
uu
ee
ooo
=++++++
jj
njTjT
zzzzeze
Yzaab
zpzpznzeze
0
() ()
k
k
Hzhkz
=
=
Se tiene
1 1
0
( ) ()
e e
=
=
jT jTk
k
He hke
()()
1
11
0
()()cos
e
ee
=
= (
jT
k
HehkkTjsenkT
() ()
1 1
0 0
()cos () e e
= =
=
k k
hkkTjhksenkT
1.-
2.-
1 1
0
( ) ()
e e
=
=
jT jTk
k
He hke
() ()
1 1
0 0
()cos () e e
= =
= +
k k
hkkTjhksenkT
Arreglar tamao en texto y frmulas
Por ser complejas
111
()()()
eee
=Z
jTjTjT
HeHeHe
111
()()()
eee
=Z
jTjTjT
HeHeHe
1 1
1
()()()
e e
e
===
jT jT
MHHeHe
1
1
() ()
e
ue=Z
jT
He
1
1
() ( )
e
ue
=Z
jT
He
y
1 1
1
01 2
1 2
()
() ...
2
u u
ee
e
oooo
(
=+++++
(
j j
n jTjT
M zz z zeze
Yz
zpzpznjzeze
De ahi:
Antitransformando:
| |
01122 111
()() () () .... () () () oooooeeue ( =++++++
k k k
nn
ykkpppMsenT
||
111
()()() ykMsenkT eeue =+
Finalmente
1
()0
e
~
jT
He
Suprime la frecuencia
1
e
1
( )1
jwT
He>
Amplifica la frecuencia
1
e
8.8.1 PERIODICIDAD DE
( )
jw
He T
1
( ) Z
jwT
He
Factor de angulo fase
Una caracterstica particular en los sistemas discretos, es que los
factores de ganancia y ngulo son peridicos en relacin con la
frecuencia.
2
() j
jT
T
ee
t
e
e
+
=
Polar e
Forma
0
3
4 T
t
3 / 4 j
e
t
5 / 4 j
e
t
3 / 2 j
e
t
10
1 45
190
1135
1180
1 225
1 270
1360
3
2 T
t
5
4 T
t
2
T
t
T
t
j T
e
t
4 T
t
0 j
e
/4 j
e
t
/2 j
e
t
t
2 T
j
e
t
2 j
e
t
4T
t
e= 0
| | Re z
( )
0
o e e
o
+
=
= =
j T j T
c
Z e e
e =
2T
t
e =
T
t
e =
3
4T
t
e =
| | Im z
8.8.2 INTRODUCCIN A FILTROS DISCRETOS.
La caracterstica de ganancia de un filtro paso bajas ideal se
muestra abajo:
( )
j T
H e
e
1
1
e
2.- Filtro pasa altas:
( )
j T
H e
e
e
2
1
3.- Filtro pasa banda:
( )
j T
H e
e
1
2
(1)
( )
jT
jT
jT
e
He
e
e
e
e
o
o
(1)(cos )
(cos )
TjsenT
TjsenT
oee
eoe
+
=
+
222
(1)
()
cos2cos
jT
He
TTsenT
e
o
eoeoe
=
++
1
tan
(cos)
senT
T
T
e
ue
eo
=
+
2cos(3cos)(1cos)
c c c
TTT oeee =
Arreglar tamao en texto y frmulas