Transformada Discreta de Fourier (DFT)
Seal Real x(n), n=0,1,N-1
Seal Discreta 1.5
Valores complejos
X(k ) x (n )e
n 0
N 1
2 kn N
k 0,..., N 1, (DFT)
0.5
X p (k )
x p (n), n 0,1,..., N 1
DFT
j ( k )
Fase
Procesamiento en Frecuencia
IDFT
-0.5
-1
8 n
10
12
14
16
18
X(k ) X(k ) e
Mdulo
Bloque de N=16 muestras
Regin de Inters
x p (n ) X p (k )e
k 0
Regin de Inters
N 1
2 kn N
n 0,....N 1 (IDFT)
Espectro de Fase
Espectro de Fase
Espectro de Mdulo
4.5 4 3.5 3
Espectro de Magnitud
1
2.5 2 1.5 1 0.5
-1
-2
0 -0.5 0 5 k 10 15
Centro de antisimetra
0 5 k 10 15
Centro de simetra
Bloque de N=16 frecuencias
Bloque de N=16 frecuencias
Prof. Dr. Guillermo Kemper
Transformada Discreta de Fourier (DFT)
Seal Real x(n), n=0,1,N-1
Seal Discreta 1.5
Valores complejos
X(k ) x (n )e
n 0
N 1
2 kn N
k 0,..., N 1, (DFT)
DSP
1 0.5
X p (k )
x p (n), n 0,1,..., N 1
DFT
Procesamiento en Frecuencia
IDFT
-0.5
-1
8 n
10
12
14
16
18
X( k ) X R (k ) j X I ( k )
Parte Re al Parte Im aginaria
Bloque de N=16 muestras
x p (n ) X p (k )e
k 0
N 1
2 kn N
n 0,....N 1 (IDFT)
Espectro de la parte Real
(simtrica)
8 6
Regin de Inters
Regin de Inters
Espectro de la parte Imaginaria
(antisimtrica)
Espectro de la Parte Real
8 6 4
Espectro de la Parte Imaginaria
4
2
2 0 -2 -4 -6 0 5 k 10 15
0 -2 -4 -6 -8 0
Centro de antisimetra
5 k 10 15
Centro de simetra
Bloque de N=16 frecuencias
Bloque de N=16 frecuencias
Prof. Dr. Guillermo Kemper
Transformada Discreta de Fourier (DFT): Algunos comentarios
La transformada de Fourier de un bloque de seal real discreta x(n) de N muestras resulta en una funcin discreta compleja X(k) de N valores (Tamao de la transformada: N muestras por periodo de espectro). El modulo de X(k) es un funcin discreta de N valores con centro de simetra ( funcin par). La fase de X(k) es un funcin discreta de N valores con centro de antisimtera (funcin impar). Si se desdobla X(k) en una funcin discreta correspondiente a la parte real y una funcin discreta correspondiente a la parte imaginaria, entonces la parte real presentar centro de simetra (funcin par) y la parte imaginaria centro de antisimetra (funcin impar). Si se procesa el espectro X(k) para obtener el espectro procesado Xp(k), este tiene que satisfacer las condiciones anteriores para que al aplicar la transformada inversa se obtenga una seal real discreta xp(n). De lo contrario la seal resultante xp(n) seria tambin compleja.
conjugado
Tambin se verifica que X(k)=X *(N-k) cuando x(n) es real.
El calculo de la DFT desgasta excesiva carga computacional por lo que su implementacin es hecha va la FFT (Fast Fourier Transform) y la IFFT (Inversa). Esto obliga a que el numero de componentes N sea potencia de 2.
Prof. Dr. Guillermo Kemper
Transformada Discreta de Fourier (DFT): Resolucin
N :Tamao de la Transformada (muestras por periodo) X(k) N=8
0
0 0
1 2
3 4 5
7 8
9 10 11 12 13
k (rad) f (Hz)
2 3 5 6 7 2 f 2f 3f fs/2 5f 6f 7f fs
Regin de inters
Numero de muestras en la regin de inters (incluyendo ) : N/2 + 1 Numero de muestras en la regin de inters (sin incluir ) : N/2
=2/N (rad)
f=fs/N (Hz)
Resolucin de la transformada
Prof. Dr. Guillermo Kemper
Transformada Discreta de Fourier (DFT): Carga Computacional
X(k ) x (n )e
n 0 N 1 j 2 kn N
k 0,..., N 1, (DFT)
DFT va Fuerza Bruta
No. de multiplicaciones complejas: N2 No. de adiciones complejas: N(N-1)
equivalente equivalente
No. de multiplicaciones reales: 4N2 No. de adiciones reales: 2N(N-1)
DFT va FFT (Fast Fourier Transform)
equivalente
No. de multiplicaciones complejas: (N/2)log2N No. de adiciones complejas: Nlog2N
equivalente
No. de multiplicaciones reales: 2Nlog2N No. de adiciones reales: 2Nlog2N
Para que funcione la FFT, el tamao de la transformada N tiene que ser una potencia de 2.
Prof. Dr. Guillermo Kemper
DFT (Fuerza Bruta) vs FFT
Prof. Dr. Guillermo Kemper
Proceso FFT (N=8)
Prof. Dr. Guillermo Kemper
Transformada Discreta de Fourier (DFT): Carga Computacional
rbol del Proceso FFT (decimacin en el tiempo) : N=8
1 Nodo
WN e W e
p N
2 N 2 p N
No. de Etapas de la FFT log 2 N Variacin del exp onente " p" en cada grupo de una det er min ada etapa: N p No. de la etapa 2
Etapa 1
Etapa 2
Etapa 3
Prof. Dr. Guillermo Kemper
