Transformada de Fourier Discreta

Universidad Politcnica de Victoria
Procesamiento Digital de Seales

Transformada de Fourier Discreta
Yahir Hernndez Mier
Objetivos
Estudiar la definicin de la Transformada de Fourier Discreta
Estudiar las propiedades de la Transformada de Fourier Discreta
Aprender a interpretar los resultados espectrales de la DFT
Aplicar ventanas para reducir las fugas de la DFT
Aplicar la DFT a sinusoides reales y complejas
Transformada de Fourier continua
X ( f )= x (t) e
j 2 f t
Identidad de Euler
dt
Im
j
1 Re
Transformada de Fourier continua
X ( f )= x (t) e
j 2 f t
Identidad de Euler
Im
j
dt
e j =cos()+ j sen ()
sen
0 cos
1 Re
Transformada de Fourier discreta

N 1
X (m)= x (n) e
j 2n m
N
n=0
La DFT es un procedimiento matemtico utilizado para determinar el

contenido armnico o frecuencial de una seal discreta
La DFT es til para analizar secuencias discretas sin importar lo que
esta secuencia represente

N 1
X ( m)= x ( n)e
n=0
j 2m n
N

N 1
X ( m)= x ( n)e
n=0
j 2m n
N
N 1
= x ( n) e j n
n=0

Series geomtricas
2 4 8
s=1+ + + +
3 9 27
Cul es la suma total?

Series geomtricas
2 4 8
s=1+ + + +
3 9 27
Por qu valor se debe multiplicar el trmino anterior para obtener el trmino siguiente?

Series geomtricas
2 4 8
s=1+ + + +
3 9 27
r=
2
3
Razn comn

Series geomtricas
2 4 8
s=1+ + + +
3 9 27
r=
2
3
Razn comn
2
2 4 8 16
r s= s= + +
+ +
3
3 9 27 81

Series geomtricas
2 4 8
s=1+ + + +
3 9 27
r=
2
3
Razn comn
2
2 4 8 16
r s= s= + +
+ +
3
3 9 27 81
2
2 4 8
2 4 8 16
sr s=s s=1+ + + + + + + +
3
3 9 27
3 9 27 81

Series geomtricas
2 4 8
s=1+ + + +
3 9 27
r=
2
3
Razn comn
2
2 4 8 16
r s= s= + +
+ +
3
3 9 27 81
2
2 4 8
2 4 8 16
sr s=s s=1+ + + + + + + +
3
3 9 27
3 9 27 81
2
s s=1
3

Series geomtricas
2 4 8
s=1+ + + +
3 9 27
r=
2
3
Razn comn
2
2 4 8 16
r s= s= + +
+ +
3
3 9 27 81
2
2 4 8
2 4 8 16
sr s=s s=1+ + + + + + + +
3
3 9 27
3 9 27 81
2
s s=1
3
s=3

Series geomtricas
2 4 8
s=1+ + + +
3 9 27
r=
2
3
Razn comn
2
2 4 8 16
r s= s= + +
+ +
3
3 9 27 81
2
2 4 8
2 4 8 16
sr s=s s=1+ + + + + + + +
3
3 9 27
3 9 27 81
2
s s=1
3
s=3
Note que:
)
0
2 4 8
2
2
2
2
s=1+ + + + =
+
+
+
3 9 27
3
3
3
3
() () () ()
s=1r + 1r + 1r + 1r +

Series geomtricas
De manera general para r1 :
1
s=a+ ar + ar + ar + + ar
n1

Series geomtricas
1
n1
a: Es el primer trmino de la serie

r: Es la razn comn

Series geomtricas
1
n1
s= a r k
k =0
n1

r: Es la razn comn

Series geomtricas
1
n1
n1

r: Es la razn comn
s= a r k
k =0
s=ar + ar + ar + ar + + ar
n1
+ ar
n1
r =r
n1+ 1
=r

Series geomtricas
1
n1

r: Es la razn comn
n1
s= a r k
k =0
1
n1
sr s=a+ ar + ar + ar + + ar
+ ar
n1
r
1
n1
r =r
3
n1+ 1
=r
ar ar ar ar ar

Series geomtricas
1
n1

r: Es la razn comn
n1
s= a r k
k =0
1
n1
s (1r)=aar
+ ar
n1
r
1
n1
r =r
3
n1+ 1
=r
ar ar ar ar ar

Series geomtricas
1
n1

r: Es la razn comn
n1
s= a r k
k =0
1
n1
s (1r)=aar
n
1r
s=a
1r
+ ar
n1
r
1
n1
r =r
3
n1+ 1
=r
ar ar ar ar ar

Series geomtricas
1
n1

r: Es la razn comn
n1
s= a r k
k =0
1
n1
s (1r)=aar
n
1r
s=a
1r
+ ar
n1
r
1
r =r
3
n1+ 1
=r
ar ar ar ar ar
Si r< 1 y n :
n1
s=
a
1r

N 1
X ( m)= x ( n)e
n=0
N 1
s=
n=0
j 2m n
N
N 1
= x ( n) e j n
n=0
N 1
a
n
1r Si r< 1 y N : s=
a
r
=
a r =a
1r
1r
n=0
n

N 1
X ( m)= x ( n)e
n=0
N 1
s=
n=0
j2 m n
N
N 1
= x ( n) e
n=0
j n
N 1
= x (n) ( e
n=0
N 1
a
n
1r Si r< 1 y N : s=
a
r
=
a r =a
1r
1r
n=0
n
j n

N 1
X ( m)= x ( n)e
j2 m n
N
n=0
N 1
s=
n=0
N 1
= x ( n) e
n=0
j n
N 1
= x (n) ( e
j n
n=0
N 1
a
n
1r Si r< 1 y N : s=
a
r
=
a r =a
1r
1r
n=0
N
Transformada de Fourier discreta de funciones elementales

Funcin impulso unitario (n)
x (n t s )= 1,
0,
n=0
n0

N 1
s=
n=0
N 1
a
n
1r Si r< 1 y N : s=
a
r
=
a r =a
1r
1r
n=0
n
Funcin escaln unitario 1( n t s )
x (n t s )= 1 ,
0,
n=0,1,2,
n<0

N 1
s=
n=0
N 1
a
n
1r Si r< 1 y N : s=
a
r
=
a r =a
1r
1r
n=0
n
Funcin polinomial a
nts
x (n)= a ,
0,
n=0,1,2,
n<0

N 1
s=
n=0
N 1
a
n
1r Si r< 1 y N : s=
a
r
=
a r =a
1r
1r
n=0
n
Funcin exponencial e
a n t s
a n t s
x (n)= e
0,
n=0,1,2,
n<0

N 1
s=
n=0
N 1
a
n
1r Si r< 1 y N : s=
a
r
=
a r =a
1r
1r
n=0
N
Funcin senoidal
x (t)= sen(t ),
0,
Identidad de Euler
t0
t< 0
e j t e j t
sen ( t)=
2j
e j t =cos ( t)+ j sen ( t )

jt
e
=cos( t ) j sen ( t)
e j t + e j t
cos ( t)=
2

N 1
X (m)= x (n) e
n=0
j 2n m
N
Forma exponencial

N 1
X (m)= x (n) e
j 2n m
N
n=0
De la identidad de Euler:
e j =cos() j sen()
Forma exponencial

N 1
X (m)= x (n) e
j 2n m
N
Forma exponencial
n=0
e j =cos() j sen()
N 1
[ (
X (m)= x (n) cos

n=0
2nm
2nm
j sen
N
N
)]
Forma rectangular

N 1
X (m)= x (n) e
j 2n m
N
Forma exponencial
n=0
e j =cos() j sen()
N 1
[ (
X (m)= x (n) cos

n=0
2nm
2nm
j sen
N
N
)]
Forma rectangular
Donde:
X(m) es el m-simo componente de frecuencia, es decir, X(0), X(1), X(2), ...
m es el ndice de la salida de la DFT en el dominio de la frecuencia
m= 0, 1, 2, 3, ..., N-1
x(n) es la secuencia de muestras de entrada, x(0), x(1), x(2), x(3), ...
n es el ndice en el dominio del tiempo de las muestras de entrada, n=0, 1, 2, ..., N-1
j = 1
N es el nmero de muestras de la secuencia de entrada y el nmero de puntos de
frecuencia en la salida DFT
Significado de la Transformada de Fourier discreta

N 1
Ejemplo:
n=0
N =4
3
[ (
X (m)= x (n) cos

n=0
[ (
)]
X (m)= x (n) cos

2 nm
2 n m
j sen
4
4
2nm
2nm
j sen
N
N
)]
Cada trmino X(m) de la salida de la DFT es la suma del producto elemento por
elemento entre una seal de entrada y una sinusoide compleja cos () j sen()
m=0
200
200
210
210
j x (0)sen
+ x (1)cos
j x (1)sen
4
4
4
4
( )
( )
( )
( )
220
220
230
230
+ x (2)cos (
j
x
(2)sen
+
x(3)cos
j
x(3)sen
( 4 )
( 4 )
( 4 )
4 )
X (0)= x (0)cos
m=1

N 1
Ejemplo:
n=0
N =4
3
[ (
X (m)= x (n) cos

n=0
[ (
)]
X (m)= x (n) cos

2 nm
2 n m
j sen
4
4
2nm
2nm
j sen
N
N
)]
elemento entre una seal de entrada y una sinusoide compleja cos () j sen()
m=2
m=3

elemento entre una seal de entrada y una sinusoide compleja cos () j sen ()
La frecuencia exacta de cada sinusoide depende de f s y de N

Frecuencia a la cual la
seal original fue
muestreada
Nmero de muestras
que se toman en la
DFT

Ejemplo:
seal original fue
muestreada
f s =500 Hz (muestras/s)
DFT de 16 puntos en los datos muestreados
Nmero de muestras
que se toman en la
DFT

Ejemplo:
seal original fue
muestreada
La frecuencia fundamental de las sinusoides es:
f s 500
=
=31.25 Hz
N 16
Nmero de muestras
que se toman en la
DFT

Ejemplo:
seal original fue
muestreada
Nmero de muestras
que se toman en la
DFT
f s 500
=
=31.25 Hz
N 16
Las frecuencias de anlisis son mltiplos enteros de la frecuencia fundamental
Trmino DFT
Trmino de frecuencia
X(0)
X(1)
X(2)
X(3)
16o
0
1
2
3
31.25 = 0 Hz
31.25 = 31.25 Hz
31.25 = 62.5 Hz
31.25 = 93.75 Hz
...
...
...
X(15)
1o
2o
3o
4o
Frecuencia de anlisis
15 31.25 = 468.75 Hz

Ejemplo:
seal original fue
muestreada
Nmero de muestras
que se toman en la
DFT
f s 500
=
=31.25 Hz
N 16
Trmino DFT
X(0)
X(1)
X(2)
X(3)
16o
0
1
2
3
31.25 = 0 Hz
31.25 = 31.25 Hz
mfs
31.25 = 62.5 Hz f anlisis (m)= N
31.25 = 93.75 Hz
...
...
...
X(15)
1o
2o
3o
4o
15 31.25 = 468.75 Hz

Ejemplo:
seal original fue
muestreada
Nmero de muestras
que se toman en la
DFT
f s 500
=
=31.25 Hz
N 16
Trmino DFT
Magnitudes
de las
16o
0
1
2
3
31.25 = 0 Hz
31.25 = 31.25 Hz
mfs
31.25 = 62.5 Hz f anlisis (m)= N
31.25 = 93.75 Hz
...
X(15)
1o
2o
3o
4o
...
de la seal
X(0)
X(1)
X(2)
X(3)
...
componentes
15 31.25 = 468.75 Hz

Relacin trigonomtrica de un valor complejo X(m) individual de salida de la DFT
Eje imaginario j
Este punto representa el nmero

complejo
X (m)= X real (m)+ j X imag (m)
Eje real
X (m)= X real (m)+ j X imag (m)= X mag (m)a un ngulo X (m)

Eje imaginario j

complejo
Eje real

Eje imaginario j

complejo
X mag (m)
Eje real

X mag (m)=X (m)= X real (m) + X imag (m)
2
Magnitud X mag (m) de X (m)

Eje imaginario j

complejo
X mag (m)
Eje real

2
X (m)=tan
X imag (m)
X real (m)

ngulo de fase X (m) de X (m)

Eje imaginario j

complejo
X mag (m)
Eje real

2
X (m)=tan
X imag (m)
X real (m)
2

ngulo de fase X (m) de X (m)
X PS (m)= X mag (m) = X real (m) + X imag (m)
Espectro de potencia X PS (m) de X (m)

Ejemplo 1
Muestree y obtenga la DFT de 8 puntos de una seal continua conteniendo
componentes de frecuencia a 1 KHz y 2KHz, expresada como:
x in (t)=sen(2 1000 t )+ 0.5sen (2 2000 t+
3
)
4
3
=135
4

Ejemplo 1
3
3
)
=135
4
4
1
Muestreando a una frecuencia f s , se toman muestras cada
=t segundos:
fs s
3
x (n)=sen (2 1000 n t s )+ 0.5sen (2 2000 n t s + )
4
x in (t)=sen(2 1000 t )+ 0.5sen (2 2000 t+

Ejemplo 1
3
3
)
=135
4
4
1
=t segundos:
fs s
3
x (n)=sen (2 1000 n t s )+ 0.5sen (2 2000 n t s + )
4
Si f s=8000 muestras/s
m=0 m=1 m=2 m=7
Frecuencias de anlisis m f s : 0 KHz ,1 KHz , 2 KHz ,,7 KHz
N
x in (t)=sen(2 1000 t )+ 0.5sen (2 2000 t+

Ejemplo 1
3
3
x in (t )=sen(2 1000 t )+0.5 sen(2 2000 t + )
=135
4
4
1
=t segundos:
fs s
3
x (n)=sen (2 1000 n t s )+ 0.5sen (2 2000 n t s + )
4
m=0 m=1 m=2 m=7
N
3
Muestras x (n)=sen (2 1000 n t s )+ 0.5sen (2 2000 n t s + )
:
4 3
x (0)=sen (210000t s )+ 0.5 sen(22000nt s + )=0.3535534
4
x (1)=0.3535534
x (4)=0.3535534
x (7)=0.3535534
x (2)=0.6464466
x (5)=1.0606602
x (3)=1.0606602
x (6)=1.3535534

Ejemplo 1
3
:
)
4 3
x (0)=sen (210000t s )+ 0.5 sen(22000nt s + )=0.3535534
4
x (1)=0.3535534
x (4)=0.3535534
x (7)=0.3535534
x (2)=0.6464466
x (5)=1.0606602
x (3)=1.0606602
x (6)=1.3535534
Muestras x (n)=sen (2 1000 n t s )+ 0.5sen (2 2000 n t s +
sen(2 1000 t )
x in (t )
0.5sen (2 2000 t+
3
)
4

Ejemplo 1
7
[ (
X (m)= x (n) cos

n=0
2 nm
2 n m
j sen
8
8
X (1)=0 j 4.0=490
X (2)=1.4142+ j 1.4142=245
X (3)=0 j 0=00
X (4)=0 j 0=00
X (5)=0 j 0=00
X (6)=1.4142 j 1.4142=245
X (7)=0+ j 4.0=490
)]

Ejemplo 1
m=1
7
[ ( )
1
1
X (1)= x(n) cos 2 n j sen 2 n
8
8
n=0
f=
F
;
fs
)]
F= f f s
f s =8000 muestras/s ;
c 1 (t )=cos ( 2 1000 t )
1
F = 8000=1000 Hz
8
c 2 (t)=sen ( 2 1000 t )

Ejemplo 1
m=1
7
[ ( )
1
1
8
8
n=0
f=
F
;
fs
)]
F= f f s
c 1 (t )=cos ( 2 1000 t )
1
F = 8000=1000 Hz
8
c 2 (t)=sen ( 2 1000 t )
m=2
7
[ ( )
X (2)= x (n) cos 2

n=0
f=
F
;
fs
1
1
n j sen 2 n
4
4
)]
F= f f s
c 1 (t )=cos ( 2 2000 t )
1
F = 8000=2000 Hz
4
c 2 (t)=sen ( 2 2000 t )

Ejemplo 1
m=1
7
[ ( )
1
1
8
8
n=0
f=
F
;
fs
)]
F= f f s
c 1 (t )=cos ( 2 1000 t )
1
F = 8000=1000 Hz
8
c 2 (t)=sen ( 2 1000 t )
m=2
7
[ ( )
X (2)= x (n) cos 2

n=0
f=
F
;
fs
1
1
n j sen 2 n
4
4
)]
F= f f s
c 1 (t )=cos ( 2 2000 t )
1
F = 8000=2000 Hz
4
c 2 (t)=sen ( 2 2000 t )

Ejemplo 1
m=3
c 1 (t)=cos ( 2 3000 t )
m=7
c 2 (t)=sen ( 2 3000 t )
c 1 (t )=cos ( 2 2000 t )
c 2 (t)=sen ( 2 2000 t )

Ejemplo 1
X (4)=0 j 0=00
X (1)=0 j 4.0=490
X (5)=0 j 0=00
X (2)=1.4142+ j 1.4142=245
X (6)=1.4142 j 1.4142=245
X (3)=0 j 0=00
X (7)=0+ j 4.0=490

Ejemplo 1
X (4)=0 j 0=00
X (1)=0 j 4.0=490
X (5)=0 j 0=00
X (2)=1.4142+ j 1.4142=245
X (6)=1.4142 j 1.4142=245
X (3)=0 j 0=00
X (7)=0+ j 4.0=490

Ejemplo 1
X (4)=0 j 0=00
X (1)=0 j 4.0=490
X (5)=0 j 0=00
X (2)=1.4142+ j 1.4142=245
X (6)=1.4142 j 1.4142=245
X (3)=0 j 0=00
X (7)=0+ j 4.0=490
Propiedades de la Transformada de Fourier discreta

Prctica 4.1.1. Objetivo: Estudiar las grficas de la parte real, imaginaria, magnitud
y ngulo de fase de la DFT, mediante un programa de clculo numrico.
3
1. Tome 8 muestras de la funcin x (t)=sen(2 1000 t)+0.5 sen(2 2000 t + )
4
a una frecuencia de muestreo fs = 8000 Hz.
2. Calcule la DFT X(m) de 8 puntos sobre las muestras adquiridas x(n).
3. En una primera subgrfica, grafique la parte real de X(m).
4. En una segunda subgrfica, grafique la parte imaginaria de X(m).
5. En una tercera subgrfica, grafique la magnitud de X(m).
6. En una cuarta subgrfica, grafique el ngulo de fase de X(m)
7. Observe la simetra en las grficas, qu puede concluir?

Prctica 4.1.1 Realizar un programa en Scilab que calcule la DFT de 8 puntos
para la seal discreta del ejemplo 1:
x (n)=sen (2 1000 n t s )+ 0.5sen (2 2000 n t s+
x (0)=0.3535534
x (4)=0.3535534
x (1)=0.3535534
x (5)=1.0606602
x (2)=0.6464466
x (6)=1.3535534
x (3)=1.0606602
x (7)=0.3535534
3
)
4
f s =8000 muestras/s
N =8

x (n)=sen (2 1000 n t s )+ 0.5sen (2 2000 n t s+
x (0)=0.3535534
x (4)=0.3535534
x (1)=0.3535534
x (5)=1.0606602
x (2)=0.6464466
x (6)=1.3535534
x (3)=1.0606602
x (7)=0.3535534
DFT: N 1
2nm
2nm
X (m)= x (n) cos
j sen
N
N
n=0
[ (
3
)
4
f s =8000 muestras/s
N =8
)]
X (4)=0 j 0=00
X (1)=0 j 4.0=490
X (2)=1.4142+ j 1.4142=245
X (3)=0 j 0=00
X (5)=0 j 0=00
X (6)=1.4142 j 1.4142=245
X (7)=0+ j 4.0=490

DFT: N 1
2nm
2nm
X (m)= x (n) cos
j sen
N
N
n=0
[ (
)]

DFT: N 1
2nm
2nm
X (m)= x (n) cos
j sen
N
N
n=0
[ (
)]

Ejemplo 1
X (0)=0
X (4)=0 j 0=00
X (1)=0 j 4.0=490
X (5)=0 j 0=00
X (2)=1.4142+ j 1.4142=245
X (6)=1.4142 j 1.4142=245
X (3)=0 j 0=00
X (7)=0+ j 4.0=490

Ejemplo 1
3
x in (t )=sen (2 1000 t )+0.5 sen (2 2000 t + )
4
N =8
Relacin entre la amplitud de la seal y la
magnitud de su DFT
X (m)=A 0
N
2
2X (m)
A0 =
N

Simetra

Simetra

Simetra
La DFT acepta entrada compleja, pero la mayora de las seales fsicas son reales

Simetra
Cuando x (n) es real las salidas complejas de la DFT desde m=1 hasta m=
son redundantes con valores de salida de frecuencia para m>
N
2
N
1
2

Simetra
N
2
N
1
2
La m-sima salida (real) de la DFT tiene el mismo valor que la salida (N-m)-sima

Simetra
N
2
N
1
2
El ngulo de fase de la m-sima salida es el negativo del ngulo de la salida
(N-m)-sima

Simetra
N
2
N
1
2
(N-m)-sima
X (m)= X (m) a X (m) grados
X (m)=X ( N m) a X ( N m) grados

Simetra
N
2
N
1
2
(N-m)-sima
X (m)= X (m) a X (m) grados
X (m)=X ( N m) a X ( N m) grados
Cuando x(n) es real:
*
X (m)= X ( N m)
Recordar:
x=a+ j b
x =a j b ;
x=e
x =e

Simetra
Demostrar la propiedad de simetra de la DFT:
N 1
X ( N m)= x (n) e
n=0
j 2 n( N m)
N

Simetra
N 1
X ( N m)= x (n) e
j 2 n( N m)
N
n=0
N 1
X ( N m)= x (n) e
n=0
j 2 n N
N
j 2n m
N

Simetra
N 1
X ( N m)= x (n) e
j 2 n( N m)
N
n=0
N 1
X ( N m)= x (n) e
j 2 n N
N
j 2n m
N
n=0
N 1
X ( N m)= x (n) e j 2 n e
n=0
j 2n m
N

Simetra
N 1
X ( N m)= x (n) e
j 2 n( N m)
N
n=0
N 1
X ( N m)= x (n) e
j 2 n N
N
j 2n m
N
n=0
N 1
X ( N m)= x (n) e j 2 n e
n=0
e j =cos() j sen ()
j 2n m
N

Simetra
N 1
X ( N m)= x (n) e
j 2 n( N m)
N
n=0
N 1
X ( N m)= x (n) e
j 2 n N
N
j 2n m
N
n=0
N 1
X ( N m)= x (n) e j 2 n e
j 2n m
N
n=0
e j =cos() j sen ()
e j 2 n =cos (2 n) j sen(2 n)

Simetra
N 1
X ( N m)= x (n) e
j 2 n( N m)
N
n=0
N 1
X ( N m)= x (n) e
j 2 n N
N
j 2n m
N
n=0
N 1
X ( N m)= x (n) e j 2 n e
j 2n m
N
n=0
e j =cos() j sen ()
e j 2 n =1

Simetra
N 1
X ( N m)= x (n) e
j 2 n( N m)
N
n=0
N 1
X ( N m)= x (n) e
j 2 n N
N
j 2n m
N
n=0
N 1
X ( N m)= x (n) e j 2 n e
j 2n m
N
n=0
e j =cos() j sen ()
e j 2 n =1
N 1
X ( N m)= x (n)(1)e
n=0
j 2 nm
N

Simetra
N 1
X ( N m)= x (n) e
j 2 n( N m)
N
N 1
X ( N m)= x (n)(1)e
j 2 nm
N
n=0
n=0
N 1
X ( N m)= x (n) e
j 2 n N
N
j 2n m
N
n=0
N 1
X ( N m)= x (n) e j 2 n e
j 2n m
N
n=0
e j =cos() j sen ()
e j 2 n =1
N 1
X ( N m)= x (n) e
n=0
j 2 n m
N

Simetra
N 1
X ( N m)= x (n) e
j 2 n( N m)
N
N 1
X ( N m)= x (n)(1)e
j 2 nm
N
n=0
n=0
N 1
X ( N m)= x (n) e
j 2 n N
N
j 2n m
N
n=0
N 1
X ( N m)= x (n) e j 2 n e
j 2n m
N
n=0
e j =cos() j sen ()
e j 2 n =1
N 1
X ( N m)= x (n) e
j 2 n m
N
n=0
N 1
X ( m)= x ( n)e
n=0
j 2 n m
N

Simetra
N 1
X ( N m)= x (n) e
j 2 n( N m)
N
N 1
X ( N m)= x (n)(1)e
j 2 nm
N
n=0
n=0
N 1
X ( N m)= x (n) e
j 2 n N
N
j 2n m
N
n=0
N 1
X ( N m)= x (n) e j 2 n e
j 2n m
N
n=0
e j =cos() j sen ()
e j 2 n =1
N 1
X ( N m)= x (n) e
j 2 n m
N
n=0
N 1
X ( m)= x ( n)e
n=0
X (m)= X * ( N m)
j2 nm
N

Simetra
N 1
X ( N m)= x (n) e
j 2 n( N m)
N
N 1
X ( N m)= x (n)(1)e
j 2 nm
N
n=0
n=0
N 1
X ( N m)= x (n) e
j 2 n N
N
j 2n m
N
n=0
N 1
X ( N m)= x (n) e j 2 n e
j 2n m
N
n=0
e j =cos() j sen ()
e j 2 n =1
N 1
X ( N m)= x (n) e
j 2 n m
N
n=0
N 1
X ( m)= x ( n)e
j2 nm
N
n=0
X (m)= X * ( N m)
Complejo conjugado

Linealidad
La DFT de la suma de dos seales es igual a la suma de las DFT de cada seal
x suma (n)=x 1 (n)+ x 2 (n)

X suma (m)= X 1 (m)+ X 2 (m)
Demostracin:
N 1
X suma (m)= x suma (n) e

n=0
j 2n m
N

Linealidad
x suma (n)=x 1 (n)+ x 2 (n)

X suma (m)= X 1 (m)+ X 2 (m)
Demostracin:
N 1
X suma (m)= x suma (n)e
j 2 nm
N
n=0
N 1
X suma (m)= [ x 1 (n)+ x 2 (n) ] e

n=0
j 2n m
N

Linealidad
x suma (n)=x 1 (n)+ x 2 (n)

X suma (m)= X 1 (m)+ X 2 (m)
Demostracin:
N 1
j 2 nm
N
n=0
N 1
X suma (m)= [ x 1 (n)+ x 2 (n) ] e
j 2n m
N
n=0
N 1
X suma (m)= x 1 (n) e

n=0
j 2 n m
N
N 1
x2 (n) e
n=0
j 2n m
N

Linealidad
x suma (n)=x 1 (n)+ x 2 (n)

X suma (m)= X 1 (m)+ X 2 (m)
Demostracin:
N 1
j 2 nm
N
n=0
N 1
X suma (m)= [ x 1 (n)+ x 2 (n) ] e
j 2n m
N
n=0
N 1
X suma (m)= x 1 (n) e
j 2 n m
N
n=0
X suma (m)= X 1 (m)+ X 2 (m)
N 1
x2 (n) e
n=0
j 2n m
N

Periodicidad
Si la DFT es evaluada para todos los valores de n, y no solo para n=0, ..., N-1,
entonces la secuencia infinita resultante es una extensin peridica, con periodo
N, de la DFT.

Periodicidad
N, de la DFT.
Demostracin:
N 1
X ( m+ N )= x (n) e
n=0
j 2 n( m+ N )
N

Periodicidad
N, de la DFT.
Demostracin:
N 1
X ( m+ N )= x (n) e
j 2 n( m+ N )
N
n=0
N 1
X ( m+ N )= x (n) e
n=0
j 2 nm
N
j 2 n N
N

Periodicidad
N, de la DFT.
Demostracin:
N 1
X ( m+ N )= x (n) e
j 2 n( m+ N )
N
n=0
N 1
X ( m+ N )= x (n) e
j 2 nm
N
j 2 n N
N
n=0
N 1
X ( m+ N )= x (n) e
n=0
j 2 nm
N
e j 2 n

Periodicidad
N, de la DFT.
Demostracin:
N 1
X ( m+ N )= x (n) e
j 2 n( m+ N )
N
n=0
N 1
X ( m+ N )= x (n) e
j 2 nm
N
j 2 n N
N
n=0
N 1
X ( m+ N )= x (n) e
j 2 nm
N
e j 2 n
n=0
e j 2 n=cos(2 n) j sen ( 2 n)=1

Periodicidad
N, de la DFT.
Demostracin:
N 1
X ( m+ N )= x (n) e
j 2 n( m+ N )
N
X ( m+ N )= x (n) e
j 2 nm
N
j 2 n N
N
n=0
N 1
X ( m+ N )= x (n) e
j 2 nm
N
X ( m+ N )= x(n) e
n=0
n=0
N 1
N 1
e j 2 n
n=0
e j 2 n=cos(2 n) j sen ( 2 n)=1
j 2 nm
N

Periodicidad
N, de la DFT.
Demostracin:
N 1
X ( m+ N )= x (n) e
j 2 n( m+ N )
N
X ( m+ N )= x (n) e
j 2 nm
N
j 2 n N
N
n=0
N 1
X ( m+ N )= x (n) e
j 2 nm
N
X ( m+ N )= x(n) e
j 2 nm
N
n=0
n=0
N 1
N 1
N 1
X ( m)= x ( n)e
n=0
e j 2 n
n=0
e j 2 n=cos(2 n) j sen ( 2 n)=1
j 2nm
N

Periodicidad
N, de la DFT.
Demostracin:
N 1
X ( m+ N )= x (n) e
j 2 n( m+ N )
N
X ( m+ N )= x (n) e
j 2 nm
N
j 2 n N
N
n=0
N 1
X ( m+ N )= x (n) e
j 2 nm
N
N 1
X ( m)= x ( n)e
n=0
e j 2 n
n=0
j 2 n
X ( m+ N )= x(n) e
j 2 nm
N
n=0
n=0
N 1
N 1
=cos(2 n) j sen ( 2 n)=1
X ( m+ N )= X (m)
j 2nm
N
Mejora de la resolucin de la DFT

Rellenado con zeros:
F =3.0 Hz
16 muestras

F =3.0 Hz
16 muestras
16 ceros

F =3.0 Hz
16 muestras
16 ceros
48 ceros

F =3.0 Hz
16 muestras
16 ceros
48 ceros
112 ceros

Rellenado con zeros en Scilab:
-->xn = [2, 5, -6, 7, 2, 1]
xn =
2. 5. - 6. 7. 2. 1.
// Seal a rellenar
-->vc = zeros(1, 32)

vc =
column 1 to 18
0. 0. 0. 0. 0.
// Vector de ceros (ms elementos que xn)

0.
0.
0.
0.
0.
0.
0.
0.
0.
column 19 to 32
0. 0. 0. 0. 0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
-->vc( 1:length(xn) ) = xn
vc =
column 1 to 18
2. 5. - 6. 7. 2. 1.
column 19 to 32
0. 0. 0. 0. 0.
-->xnrellena = vc
0.
0.
0.
0.
0.
// Copiamos xn al inicio del vector de ceros

0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
// xnrellena toma el valor de vc
0.
0.
0.
0.

Prctica 4.1.2. Objetivo: Comprobar la propiedad de simetra mediante la
graficacin de la magnitud de la DFT en un programa de clculo numrico.
1
1. Considere la siguiente funcin discreta x (n)=3 sen (2 n) n=[0, 10]
10
2. Con la muestras obtenidas, calcule la DFT de 512 puntos.
3. Grafique la magnitud de la DFT y sus partes real e imaginaria (con plot, no con
sci_stem).
4. Escoja un valor de frecuencia m (por ejemplo m = 5). Muestre el resultado de
la siguiente diferencia: X(m) X(N-m+2). Recuerde que los ndices en los
programas de clculo numrico como Scilab y Matlab comienzan en 1.
5. Repita los pasos anteriores pero ahora con la funcin compleja
6. Qu puede concluir sobre la propiedad de simetra de la DFT?
x (n)=e
1
n
10

Prctica 4.1.3. Objetivo: Comprobar la propiedad de linealidad mediante la
1. Considere la siguientes funciones discretas y la suma de ellas:
1
x 1 (n)=3 cos (2 n)
7
2
x 2 ( n)=0.8 sen(2 n)
5
x s (n)= x 1 (n)+ x 2 ( n)
n=[0,10]
2. Calcule la DFT de 256 puntos sobre x1(n), x2(n) y xs(n).

3. En una primera subgrfica, grafique la magnitud de Xs(m) (con plot).
4. En una segunda subgrfica, grafique la magnitud de X2s(m) = X1(m) + X2(m).
5. Muestre la suma de la diferencia entre Xs(m) y X2s(m), mediante la instruccin
(en Scilab) disp( sum( abs(Xsm X2sm) ) ).
6. Con fundamento en las grficas y el resultado de la suma de la diferencia
entre las magnitudes de las DFT, qu puede concluir sobre la propiedad de
linealidad de la DFT?

Prctica 4.1.4. Objetivo: Comprobar la propiedad de periodicidad mediante la
1. Considere la siguiente funcin discreta x (n)=(0.9)
n=[0,10]
2. Con la muestras obtenidas, calcule la DFT de 512 puntos.

3. Con la muestras obtenidas, calcule la DFT de 512 puntos en dos ciclos con la
siguiente funcin:
4. Grafique en una primera subgrfica la magnitud de la DFT de un ciclo.

5. Grafique en una segunda subgrfica la magnitud de la DFT de dos ciclos.

6. Repita los pasos anteriores, pero para la siguiente funcin:
x (n)= 0.9 e
j n
3
7. En funcin de los resultados de las dos partes de la prctica, Qu puede

concluir sobre la propiedad de periodicidad de la DFT?
Teorema de desplazamiento para la DFT

Establece que un desplazamiento temporal de una secuencia peridica x(n) se
manifiesta como un desplazamiento de fase en los ngulos asociados con la DFT
N 1
X desp (m)= x (n+ k ) e

n=0
j 2 nm
N

N 1

n=0
q=n+ k
j 2 nm
N

N 1
j 2 nm
N
n=0
q=n+ k
n=qk

N 1
j 2 nm
N
n=0
q=n+ k
n=qk
N 1+ k
X desp (m)=
q=k
x (q)e
j 2 (qk ) m
N

N 1
j 2 nm
N
n=0
q=n+ k
n=qk
N 1+ k
X desp (m)=
x (q)e
j 2 (qk ) m
N
q=k
N 1
X desp (m)= x (q)e

q=0
j 2 q m
N
j 2 k m
N

N 1
j 2 nm
N
n=0
q=n+ k
n=qk
N 1+ k
X desp (m)=
x (q)e
j 2 (qk ) m
N
q=k
N 1
X desp (m)= x (q)e
j 2 q m
N
j 2 k m
N
q=0
X desp (m)=e
j 2 k m N 1
N
x (q)e
q=0
j 2 q m
N

N 1
j 2 nm
N
X desp (m)=e
n=0
q=n+ k
n=qk
N 1+ k
X desp (m)=
x (q)e
j 2 (qk ) m
N
q=k
N 1
X desp (m)= x (q)e
j 2 q m
N
j 2 k m
N
q=0
X desp (m)=e
j 2 k m N 1
N
x (q)e
q=0
j 2 q m
N
j 2 k m N 1
N
x (q)e
q=0
j 2 q m
N

N 1
j 2 nm
N
X desp (m)=e
n=0
q=n+ k
j 2 k m N 1
N
x (q)e
q=0
n=qk
N 1+ k
X desp (m)=
x (q)e
j 2 (qk ) m
N
q=k
N 1
X desp (m)= x (q)e
j 2 q m
N
j 2 k m
N
q=0
X desp (m)=e
j 2 k m N 1
N
x (q)e
q=0
j 2 q m
N
X desp (m)=e
j 2 k m
N
X (m)
j 2 q m
N

N 1
j 2 nm
N
X desp (m)=e
n=0
q=n+ k
j 2 k m N 1
N
x (q)e
j 2 q m
N
q=0
n=qk
N 1+ k
X desp (m)=
x (q)e
j 2 (qk ) m
N
X desp (m)=e
j 2 k m
N
X (m)
q=k
N 1
X desp (m)= x (q)e
j 2 q m
N
j 2 k m
N
q=0
X desp (m)=e
j 2 k m N 1
N
x (q)e
q=0
j 2 q m
N
Desplazamiento:
2 k m
360 k m
radianes o
grados
N
N

Ejemplo 2
Suponga que la seal del ejemplo 1 es muestreada con un adelanto en el
3
tiempo k =3:
x in (t)=sen(2 1000 t )+ 0.5sen (2 2000 t+ )
4
3
x (n+ 3)=sen (2 1000(n+ 3)t s )+ 0.5 sen(2 2000(n+ 3) t s + )
4

Ejemplo 2
Suponga que la seal del ejemplo 1 es muestreada con un adelanto en el
3
tiempo k =3:
x in (t)=sen(2 1000 t)+ 0.5sen (2 2000 t+ )
4
3
x (n+ 3)=sen (2 1000(n+ 3)t s )+ 0.5 sen(2 2000(n+ 3) t s + )
4
Muestras en el ejemplo 1
Muestras en el ejemplo 2

Ejemplo 2
3
Muestras x (n+ 3)=sen (2 1000(n+ 3)t s )+ 0.5 sen(2 2000(n+ 3) t s+ )
4
x (0)=1.0606602
x (1)=0.3535534
x (2)=1.0606602
x (3)=1.3535534
x (4)=0.3535534
x (5)=0.3535534
x (6)=0.3535534
x (7)=0.6464466

Ejemplo 2
Dominio de la frecuencia :

Ejemplo 2
X desp (m)=e
j 2 k m
N
m=1 N =8
X (m)
k =3

Ejemplo 2
X desp (m)=e
j 2 k m
N
m=1 N =8
X desp (1)=e
X (m)
k =3
j 2 (3)(1)
8
X (1)

Ejemplo 2
X desp (m)=e
j 2 k m
N
m=1 N =8
X desp (1)=e
X (m)
k =3
j 2 (3)(1)
8
X (1)
j
2
X (1)=490 =4 e

Ejemplo 2
X desp (m)=e
j 2 k m
N
m=1 N =8
X desp (1)=e
X (m)
k =3
j 2 (3)(1)
8
X (1)
j
2
X (1)=490 =4 e
X desp (1)=e
j 6
8
4e

Ejemplo 2
X desp (m)=e
j 2 k m
N
X (m)
m=1 N =8
X desp (1)=e
k =3
j 2 (3)(1)
8
X (1)
j
2
X (1)=490 =4 e
X desp (1)=e
X desp (1)=e
j 6
8
j 3
4
4e
4e

Ejemplo 2
X desp (m)=e
j 2 k m
N
X (m)
m=1 N =8
X desp (1)=e
k =3
j 2 (3)(1)
8
X (1)
j
2
X (1)=490 =4 e
X desp (1)=e
X desp (1)=e
j 6
8
j 3
4
X desp (1)=4 e
4e
4e
3 2
4 4

Ejemplo 2
X desp (m)=e
j 2 k m
N
X (m)
m=1 N =8
X desp (1)=e
k =3
j 2 (3)(1)
8
X (1)
j
2
X (1)=490 =4 e
X desp (1)=e
X desp (1)=e
j 6
8
4e
j 3
4
X desp (1)=4 e
X desp (1)=4 e
j
j
4e
3 2
4 4

Ejemplo 2
X desp (m)=e
j 2 k m
N
X (m)
m=1 N =8
X desp (1)=e
k =3
j 2 (3)(1)
8
X (1)
j
2
X (1)=490=4 e
X desp (1)=e
X desp (1)=e
j 6
8
4e
j 3
4
X desp (1)=4 e
X desp (1)=4 e
j
2
j
2
4e
3 2
4 4
j
4
X desp (1)=445

Ejemplo 2
X desp (1)=445
X (1)=490
Ejemplo 1

Prctica 4.1.5. Objetivo: Comprobar el teorema de desplazamiento mediante la
graficacin de la DFT de una seal desplazada k muestras.
1. Considere la funcin del ejemplo 1: x (t )=sen (2 1000 t )+0.5 sen ( 2 2000 t+
2. Muestree la seal a fs = 8000 Hz y obtenga 8 muestras n = [0,7], pero para n
adelantada tres muestras ( n+3 ).
3. Calcule la DFT X(m) de 8 puntos sobre las muestras adquiridas en n+3.
4. En una primera subgrfica, grafique la parte real de X(m).
5. En una segunda subgrfica, grafique la parte imaginaria de X(m).
6. En una tercera subgrfica, grafique la magnitud de X(m).
7. En una cuarta subgrfica, grafique el ngulo de fase de X(m).
8. Con fundamento en las grficas, qu puede concluir sobre el teorema de
desplazamiento de la DFT?
3
)
4
Transformada de Fourier Discreta Inversa (IDFT)

Definicin:
1
x (n)=
N
N 1
X (m)e
m=0
j 2m n
N

Definicin:
1
x (n)=
N
N 1
1
x (n)=
N
N 1
X (m)e
j 2m n
N
m=0
m=0
X (m) cos(
2mn
2mn
)+ j sen (
)
N
N

Definicin:
1
x (n)=
N
N 1
1
x (n)=
N
N 1
X (m)e
j 2m n
N
m=0
m=0
X (m) cos(
2mn
2mn
)+ j sen (
)
N
N
Una seal en tiempo discreto puede ser considerada como una suma de
varias frecuencias analticas sinusoidales

Definicin:
1
x (n)=
N
N 1
1
x (n)=
N
N 1
X (m)e
j 2m n
N
m=0
m=0
X (m) cos(
2mn
2mn
)+ j sen (
)
N
N
Las salidas X(m) de la DFT son un conjunto de N valores complejos que
indican la magnitud y la fase de cada frecuencia de anlisis involucrada en
la suma

Definicin:
1
x (n)=
N
N 1
1
x (n)=
N
N 1
X (m)e
j 2m n
N
m=0
m=0
X (m) cos(
2mn
2mn
)+ j sen (
)
N
N
Las salidas X(m) de la DFT son un conjunto de N valores complejos que
indican la magnitud y la fase de cada frecuencia de anlisis involucrada en
la suma
X (0)=0 j 0=00
X (4)=0 j 0=00
X (1)=0 j 4.0=490
X (5)=0 j 0=00
X (2)=1.4142+ j 1.4142=245
X (6)=1.4142 j 1.4142=245
X (3)=0 j 0=00
X (7)=0+ j 4.0=490

X (0)=0 j 0=00
X (4)=0 j 0=00
X (1)=0 j 4.0=490
X (5)=0 j 0=00
X (2)=1.4142+ j 1.4142=245
X (6)=1.4142 j 1.4142=245
X (3)=0 j 0=00
X (7)=0+ j 4.0=490
x (n)=
1
N
N 1
X (m) cos( 2 Nm n )+
m=0
j sen (
2mn
)
N

X (0)=0 j 0=00
X (4)=0 j 0=00
X (1)=0 j 4.0=490
X (5)=0 j 0=00
X (2)=1.4142+ j 1.4142=245
X (6)=1.4142 j 1.4142=245
X (3)=0 j 0=00
X (7)=0+ j 4.0=490
x (n)=
1
N
N 1
X (m) cos( 2 Nm n )+
m=0
j sen (
2mn
)
N
x (0)=0.3535534+ j 0.0
x (4)=0.3535534+ j 0.0
x (1)=0.3535534+ j 0.0
x (5)=1.0606602+ j 0.0
x (2)=0.6464466+ j 0.0
x (6)=1.3535534+ j 0.0
x (3)=1.0606602+ j 0.0
x (7)=0.3535534+ j 0.0
Transformada de Fourier discreta Inversa

Prctica 4.1.6. Objetivo: Comprender la relacin entre la DFT y la IDFT mediante
la programacin de esta ltima en un programa de clculo numrico.
1. Considere la funcin del ejemplo 1: x (t )=sen (2 1000 t )+0.5 sen ( 2 2000 t+
2. Muestree la seal a fs = 8000 Hz y obtenga 8 muestras n = [0,7].
3. Calcule la DFT de la seal con la funcin mi_DFT.sce.
4. Partiendo de la funcin mi_DFT.sce, programe la ecuacin de la IDFT:
5. Calcule la IDFT de 8 puntos sobre X(m) con la funcin mi_IDFT.sce.

6. Compare el resultado de la IDFT con las muestras x(n). qu puede concluir?
3
)
4
Fugas (leakage) en la DFT

3
3
=135
x in (t )=sen(2 1000 t )+0.5 sen(2 2000 t + )
4
4
1
=t segundos:
fs s
3
x (n)=sen (2 1000 n t s )+ 0.5sen (2 2000 n t s + )
4
m=0 m=1 m=2 m=7
N
Ejemplo 1
X (0)=0
X (1)=0 j 4.0=490
X (2)=1.4142+ j 1.4142=245
X (3)=0 j 0=00
X (4)=0 j 0=00
X (5)=0 j 0=00
X (6)=1.4142 j 1.4142=245
X (7)=0+ j 4.0=490

Las fugas hacen que los resultados de la DFT sean solo una aproximacin
al verdadero espectro de la seal antes del muestreo.
Se puede reducir mas no eliminar

Recordar que:
f anlisis (m)=
mfs
,
N
m=1, 2,, N 1

Recordar que:
mfs
, m=1, 2,, N 1
N
La DFT produce resultados correctos solo cuando la secuencia de datos de
f anlisis (m)=
entrada contiene energa precisamente en la frecuencias de anlisis

Considere una seal de frecuencia igual a 3 Hz, de la cual se extraen 64
muestras (n=0... 63) y es muestreada a f s =64 muestras/s . Calculando la DFT:

Frecuencias de anlisis m f s :
N
0 Hz ,1 Hz , 2 Hz , ,63 Hz

N
0 Hz ,1 Hz , 2 Hz , ,63 Hz

F =3.4 Hz
N
0 Hz ,1 Hz , 2 Hz , ,63 Hz

F =3.4 Hz Frecuencias de anlisis :
N
0 Hz ,1 Hz , 2 Hz , ,63 Hz
0 Hz ,1 Hz , 2 Hz , ,63 Hz

F =3.4 Hz Frecuencias de anlisis :
N
0 Hz ,1 Hz , 2 Hz , ,63 Hz
0 Hz ,1 Hz , 2 Hz , ,63 Hz

Qu causa ese comportamiento?
m=4 Output bin (Contenedor de salida)

Qu causa ese comportamiento?
m=4 Output bin (Contenedor de salida)
Las fugas hacen que cualquier seal cuya frecuencia no se encuentre en el

centro del contenedor (bin) se fuguen hacia los otros contenedores.

La respuesta en amplitud de una DFT de N-puntos sobre una funcin sinusoidal
en trminos del ndice del contenedor m puede ser aproximada por la funcin
sinc:
X (m)
N sen(( k m))
2
( k m)
Nx
)
sen( x)
sen( x)
2
DFT cos =
, o
, o
x
x
x
sen( )
sen ( )
N
2
sen (

N sen((k m))
Demostracin: X (m)
2
(k m)
Funcin real:
2nk
x r (n)=cos
N
Nx
)
sen( x)
sen ( x)
2
DFT cos =
, o
, o
x
x
x
sen( )
sen( )
N
2
sen(
Nx
)
sen(
x)
sen
(
x)
2
N sen((k m)) DFT =
,
o
,
o
Demostracin: X (m)
cos
x
x
x
2
(k m)
sen(
)
sen(
)
Funcin real:
N
2
2nk
Nmero de ciclos enteros
x r (n)=cos
en N muestras
N
sen(
Nx
)
sen(
x)
sen
(
x)
2
N sen((k m)) DFT =
,
o
,
o
Demostracin: X (m)
cos
x
x
x
2
(k m)
sen(
)
sen(
)
Funcin real:
N
2
2nk
x r (n)=cos
en N muestras
N
N 1
X r (m)= x r (n)e
n=0
2n m
N
sen(
Nx
)
sen(
x)
sen
(
x)
2
N sen((k m)) DFT =
,
o
,
o
Demostracin: X (m)
cos
x
x
x
2
(k m)
sen(
)
sen(
)
Funcin real:
N
2
2nk
x r (n)=cos
en N muestras
N
N 1
X r (m)= x r (n)e
n=0
sen(
2n m
N
N 1
=
n=0
2nk j
cos
e
N
2 nm
N
Nx
)
sen(
x)
sen
(
x)
2
N sen((k m)) DFT =
,
o
,
o
Demostracin: X (m)
cos
x
x
x
2
(k m)
sen(
)
sen(
)
Funcin real:
N
2
2nk
x r (n)=cos
en N muestras
N
N 1
X r (m)= x r (n)e
sen(
n=0
Identidad de Euler
2n m
N
N 1
=
n=0
2nk j
cos
e
N
e j t +e j t
cos( t )=
2
2 nm
N
Nx
)
sen(
x)
sen
(
x)
2
N sen((k m)) DFT =
,
o
,
o
Demostracin: X (m)
cos
x
x
x
2
(k m)
sen(
)
sen(
)
Funcin real:
N
2
2nk
x r (n)=cos
en N muestras
N
N 1
X r (m)= x r (n)e
2n m
N
n=0
Identidad de Euler
N 1
sen(
j
1
X r (m)= e
2 n=0
2 n k
N
N 1
=
n=0
2nk j
cos
e
N
e j t +e j t
cos( t )=
2
+e
2 nk
N
)e
2n m
N
2 nm
N
Nx
)
sen(
x)
sen
(
x)
2
N sen((k m)) DFT =
,
o
,
o
Demostracin: X (m)
cos
x
x
x
2
(k m)
sen(
)
sen(
)
Funcin real:
N
2
2nk
x r (n)=cos
en N muestras
N
N 1
X r (m)= x r (n)e
2n m
N
n=0
j
1
X r (m)= e
2 n=0
N 1
j
1
X r (m)= e
2 n=0
2nk j
cos
e
N
2 nm
N
e j t +e j t
cos( t )=
2
2 n k
N
2 n k
N
N 1
n=0
Identidad de Euler
N 1
sen(
+e
2 nk
N
2 nm
N
)e
1
+
2
2n m
N
N 1
e
n=0
2 nk
N
2n m
N
Nx
)
sen(
x)
sen
(
x)
2
N sen((k m)) DFT =
,
o
,
o
Demostracin: X (m)
cos
x
x
x
2
(k m)
sen(
)
sen(
)
Funcin real:
N
2
2nk
x r (n)=cos
en N muestras
N
N 1
X r (m)= x r (n)e
2n m
N
n=0
j
1
X r (m)= e
2 n=0
N 1
j
1
X r (m)= e
2 n=0
N 1
j
1
X r (m)= e
2 n=0
2nk j
cos
e
N
2 nm
N
e j t +e j t
cos( t )=
2
2 n k
N
2 n k
N
N 1
n=0
Identidad de Euler
N 1
sen(
+e
2 n (k m)
N
2 nk
N
2 nm
N
1
+
2
)e
1
+
2
N 1
2 nk
N
n=0
N 1
e
n=0
2n m
N
2 n(k + m)
N
2n m
N
Nx
)
sen(
x)
sen
(
x)
2
N sen((k m)) DFT =
,
o
,
o
Demostracin: X (m)
cos
x
x
x
2
(k m)
sen(
)
sen(
)
Funcin real:
N
2
2nk
x r (n)=cos
en N muestras
N
N 1
X r (m)= x r (n)e
2n m
N
n=0
j
1
X r (m)= e
2 n=0
N 1
j
1
X r (m)= e
2 n=0
N 1
j
1
X r (m)= e
2 n=0
2nk j
cos
e
N
2 nm
N
e j t +e j t
cos( t )=
2
2 n k
N
2 n k
N
N 1
n=0
Identidad de Euler
N 1
sen(
+e
2 n (k m)
N
2 nk
N
2 nm
N
1
+
2
)e
1
+
2
N 1
2 nk
N
n=0
N 1
e
n=0
2n m
N
2 n(k + m)
N
2n m
N
Nx
)
sen(
x)
sen
(
x)
2
N sen((k m)) DFT =
,
o
,
o
Demostracin: X (m)
cos
x
x
x
2
(k m)
sen(
)
sen(
)
Funcin real:
N
2
2nk
2 ( k m)
x r (n)=cos
q=
en N muestras
N
N
N 1
X r (m)= x r (n)e
2n m
N
n=0
j
1
X r (m)= e
2 n=0
N 1
j
1
X r (m)= e
2 n=0
N 1
j
1
X r (m)= e
2 n=0
2nk j
cos
e
N
2 nm
N
e j t +e j t
cos( t )=
2
2 n k
N
2 n k
N
N 1
n=0
Identidad de Euler
N 1
sen(
+e
2 n (k m)
N
2 nk
N
2 nm
N
1
+
2
)e
1
+
2
N 1
2 nk
N
n=0
N 1
e
n=0
2n m
N
2 n(k + m)
N
2n m
N
Nx
)
sen(
x)
sen
(
x)
2
N sen((k m)) DFT =
,
o
,
o
Demostracin: X (m)
cos
x
x
x
2
(k m)
sen(
)
sen(
)
Funcin real:
N
2
2nk
2 ( k m)
x r (n)=cos
q=
en N muestras
N
N
N 1
X r (m)= x r (n)e
2n m
N
n=0
j
1
X r (m)= e
2 n=0
N 1
j
1
X r (m)= e
2 n=0
N 1
j
1
X r (m)= e
2 n=0
2nk j
cos
e
N
2 nm
N
+e
2 n (k m)
N
2 nk
N
2 nm
N
1
+
2
)e
1
+
2
2 nk
N
n=0
N 1
e
n=0
2n m
N
N 1
2 n(k + m)
N
N 1
e
n=0
e j t +e j t
cos( t )=
2
2 n k
N
2 n k
N
N 1
n=0
Identidad de Euler
N 1
sen(
2n m
N
2 n(k m)
N
N 1
= e j n q
n=0
Nx
)
sen(
x)
sen
(
x)
2
N sen((k m)) DFT =
,
o
,
o
Demostracin: X (m)
cos
x
x
x
2
(k m)
sen(
)
sen(
)
Funcin real:
N
2
2nk
2 ( k m)
x r (n)=cos
q=
en N muestras
N
N
N 1
X r (m)= x r (n)e
2n m
N
n=0
j
1
X r (m)= e
2 n=0
N 1
j
1
X r (m)= e
2 n=0
N 1
j
1
X r (m)= e
2 n=0
2nk j
cos
e
N
2 nm
N
+e
2 n (k m)
N
2 nk
N
2 nm
N
1
+
2
)e
1
+
2
2 nk
N
n=0
N 1
e
n=0
2n m
N
N 1
2 n(k + m)
N
N 1
n=0
e j t +e j t
cos( t )=
2
2 n k
N
2 n k
N
N 1
n=0
Identidad de Euler
N 1
sen(
2n m
N
2 n(k m)
N
N 1
N 1
= e j n q = ( e j q )
n=0
n=0
Nx
)
sen(
x)
sen
(
x)
2
N sen((k m)) DFT =
,
o
,
o
Demostracin: X (m)
cos
x
x
x
2
(k m)
sen(
)
sen(
)
Funcin real:
N
2
2nk
2 ( k m)
x r (n)=cos
q=
en N muestras
N
N
N 1
X r (m)= x r (n)e
2n m
N
n=0
1
X r (m)= e
2 n=0
N 1
j
1
X r (m)= e
2 n=0
N 1
j
1
X r (m)= e
2 n=0
2nk j
cos
e
N
2 nm
N
+e
2 n (k m)
N
2 nk
j
N
2 nm
N
1
+
2
)e
1
+
2
2 n(k m)
N
2 nk
N
2 n(k + m)
N
s=
2n m
N
N 1
N 1
= e j n q = ( e j q )
n=0
n=0
n=0
N 1
n=0
N 1
2n m
j
N
N 1
N 1
n=0
e j t +e j t
cos( t )=
2
2 n k
j
N
2 n k
N
N 1
n=0
Identidad de Euler
N 1
sen(
1r N
a r =a
1r
n
n=0
Nx
)
sen(
x)
sen
(
x)
2
N sen((k m)) DFT =
,
o
,
o
Demostracin: X (m)
cos
x
x
x
2
(k m)
sen(
)
sen(
)
Funcin real:
N
2
2nk
2 ( k m)
x r (n)=cos
q=
en N muestras
N
N
N 1
X r (m)= x r (n)e
2n m
N
n=0
1
X r (m)= e
2 n=0
N 1
j
1
X r (m)= e
2 n=0
N 1
j
1
X r (m)= e
2 n=0
N 1
=
n=0
2nk j
cos
e
N
2 nm
N
2 n k
j
N
+e
2 nk
j
N
)e
N 1
2 n(k m)
N
n=0
e j t +e j t
cos( t )=
2
Identidad de Euler
N 1
sen(
s=
n=0
N 1
(e
n=0
2 n k
N
2 n (k m)
N
2 nm
N
1
+
2
1
+
2
N 1
2 nk
N
n=0
N 1
e
n=0
2 n(k + m)
N
2n m
N
N 1
= e j n q = ( e j q )
n=0
N 1
2n m
j
N
N 1
1r N
a r =a
1r
n
1e j q N
)=
1e j q
jq n
n=0
Nx
)
sen(
x)
sen
(
x)
2
N sen((k m)) DFT =
,
o
,
o
Demostracin: X (m)
cos
x
x
x
2
(k m)
sen(
)
sen(
)
Funcin real:
N
2
2nk
2 ( k m)
x r (n)=cos
q=
en N muestras
N
N
N 1
X r (m)= x r (n)e
2n m
N
n=0
j
1
X r (m)= e
2 n=0
N 1
j
1
X r (m)= e
2 n=0
N 1
j
1
X r (m)= e
2 n=0
2nk j
cos
e
N
2 nm
N
+e
2 nk
N
)e
N 1
2 n(k m)
N
n=0
e j t +e j t
cos( t )=
2
2 n k
N
2 n k
N
N 1
n=0
Identidad de Euler
N 1
sen(
s=
n=0
N 1
2 n (k m)
N
2 nm
N
1
+
2
1
+
2
N 1
2 nk
N
n=0
N 1
e
n=0
= e j n q = ( e j q )
n=0
1r N
a r =a
1r
n
jqN
( e ) = 1e
1e j q
n=0
N 1
n=0
N 1
2n m
N
N 1
2 n(k + m)
N
2n m
N
jq n
jq
N
2
jq
N
2
()
e
e
e
e
q
2
q
2
Nx
)
sen(
x)
sen
(
x)
2
N sen((k m)) DFT =
,
o
,
o
Demostracin: X (m)
cos
x
x
x
2
(k m)
sen(
)
sen(
)
Funcin real:
N
2
2nk
2 ( k m)
x r (n)=cos
q=
en N muestras
N
N
N 1
X r (m)= x r (n)e
2n m
N
n=0
j
1
X r (m)= e
2 n=0
N 1
j
1
X r (m)= e
2 n=0
N 1
X r (m)=
1
e
2 n=0
2nk j
cos
e
N
2 nm
N
+e
2 nk
N
)e
N 1
2 n(k m)
N
n=0
e j t +e j t
cos( t )=
2
2 n k
N
2 n k
N
N 1
n=0
Identidad de Euler
N 1
sen(
N 1
= e j n q = ( e j q )
n=0
N 1
s=
n=0
2n m
N
N 1
N 1
n=0
1r N
a r =a
1r
n
jqN
( e ) = 1e
1e j q
jq n
n=0
2 n (k m)
N
2 nm
N
1
2
1
+
2
N 1
n=0
2 nk
N
2n m
N
n=0
N 1
2 n(k + m)
N
jq
jq
e
e
N
2
jq N
e
N
2
j
j
e
e
q
2
jq
e
q
2
e
e
jq
N
2
jq
N
2
q
2
q
2
jq
N
2
jq
N
2
()
e
e
e
e
q
2
q
2

Demostracin: X (m)
e
N 1
jq n
(e ) =
n=0
jq
jq
e
e
N
2
jqN
e
N
2
j
j
e
e
q
2
jq
e
q
2
e
e
N sen((k m))
2
(k m)
jq
N
2
jq
N
2
q
2
q
2
Nx
)
sen( x)
sen ( x)
2
DFT cos =
, o
, o
x
x
x
sen( )
sen( )
N
2
sen(

Demostracin: X (m)
e
N 1
(e ) = e
n=0
jq n
jq
N
2
N
jq
2
e
e
jqN
e
e
q
2
q
j
2
jq
e
e
N sen((k m))
2
(k m)
jq
N
2
jq
N
2
q
2
q
2
e
=
jq
N
2
e
e
jqN
jq
e e j q e
e
jq
N
2
q
2
q
2
Nx
)
sen( x)
sen ( x)
2
DFT cos =
, o
, o
x
x
x
sen( )
sen( )
N
2
N
sen(
q
2

Demostracin: X (m)
e
N 1
(e ) = e
jq n
n=0
jq
N
2
N
jq
2
e
e
q
2
(e
jq
N
2
(e
jq
N
2
j
e
q
2
e
e
q
2
q
j
2
jq
e
e
jq
N
2
q
2
jqN
N sen((k m))
2
(k m)
jq
N
2
jq
N
2
q
2
q
2
e
=
jq
N
2
e
e
jqN
jq
e e j q e
e
jq
N
2
q
2
q
2
Nx
)
sen( x)
sen ( x)
2
DFT cos =
, o
, o
x
x
x
sen( )
sen( )
N
2
N
sen(
q
2

N sen((k m))
2
(k m)
Demostracin: X (m)
e
N 1
(e ) = e
jq n
n=0
jq
N
2
N
jq
2
e
e
q
2
(e
jq
N
2
(e
jq
N
2
j
e
q
2
e
e
q
2
q
j
2
jq
e
e
jq
N
2
q
2
jqN
) =e
jq
N
2
jq
N
2
q
2
q
2
q
2
jq
N
2
e
e
jqN
jq
N
2
e e j q e
e
e
jq
N
2
(e
(e
q
2
jq
N
2
q
2
jq
q
2
q
2
q
2
jq
N
2
Nx
)
sen( x)
sen ( x)
2
DFT cos =
, o
, o
x
x
x
sen( )
sen( )
N
2
N
sen(

N sen((k m))
2
(k m)
Demostracin: X (m)
e
N 1
(e ) = e
jq n
n=0
jq
N
2
N
jq
2
e
e
q
2
(e
jq
N
2
(e
jq
N
2
j
e
q
2
e
e
q
2
q
j
2
jq
e
e
jq
N
2
q
2
jqN
) =e
jq
N
2
jq
N
2
q
2
q
2
q
2
jq
N
2
e
e
jqN
jq
N
2
e e j q e
e
e
jq
N
2
(e
(e
q
2
jq
N
2
q
2
jq
q
2
q
2
q
2
jq
N
2
Nx
)
sen( x)
sen ( x)
2
DFT cos =
, o
, o
x
x
x
sen( )
sen( )
N
2
N
sen(

N sen((k m))
2
(k m)
Demostracin: X (m)
e
N 1
(e ) = e
jq n
n=0
jq
N
2
N
jq
2
e
e
q
2
(e
jq
N
2
(e
jq
N
2
j
e
q
2
j q jq
N
2
2
e
e
q
2
q
j
2
jq
e
e
jq
N
2
q
2
jqN
) =e
jq
N
2
jq
N
2
q
2
q
2
q
2
jq
N
2
e
e
jqN
jq
N
2
e e j q e
e
e
jq
N
2
(e
(e
q
2
jq
N
2
q
2
jq
q
2
q
2
q
2
jq
N
2
Nx
)
sen( x)
sen ( x)
2
DFT cos =
, o
, o
x
x
x
sen( )
sen( )
N
2
N
sen(

N sen((k m))
2
(k m)
Demostracin: X (m)
e
N 1
(e ) = e
jq n
n=0
jq
N
2
N
jq
2
e
e
q
2
(e
jq
N
2
(e
jq
N
2
j
e
q
2
j q jq
N
2
2
e
q
2
jq
e
e
jq
N
2
q
2
=e
e
e
q
j
2
jqN
) =e
jq
(1 N )
2
jq
N
2
jq
N
2
q
2
q
2
q
2
jq
N
2
e
e
jqN
jq
N
2
e e j q e
e
e
jq
N
2
(e
(e
q
2
jq
N
2
q
2
jq
q
2
q
2
q
2
jq
N
2
Nx
)
sen( x)
sen ( x)
2
DFT cos =
, o
, o
x
x
x
sen( )
sen( )
N
2
N
sen(

N sen((k m))
2
(k m)
Demostracin: X (m)
e
N 1
(e ) = e
jq n
n=0
jq
N
2
N
jq
2
e
e
q
2
(e
jq
N
2
(e
jq
N
2
j
e
q
2
j q jq
N
2
2
e
q
2
jq
N
2
q
2
jq
(1 N )
2
jq
N
2
q
2
q
2
q
2
) =e
N
2
jq
=e
e
e
q
j
2
jqN
jq
=e
jq
N
2
e
e
jqN
jq
N
2
e e j q e
e
e
jq
N
2
(e
(e
q
2
jq
[(N 1)]
2
jq
N
2
q
2
jq
q
2
q
2
q
2
jq
N
2
Nx
)
sen( x)
sen ( x)
2
DFT cos =
, o
, o
x
x
x
sen( )
sen( )
N
2
N
sen(

N sen((k m))
2
(k m)
Demostracin: X (m)
e
N 1
(e ) = e
jq n
n=0
jq
N
2
N
jq
2
e
e
q
2
(e
jq
N
2
(e
jq
N
2
j
e
q
2
j q jq
N
2
2
e
q
2
jq
N
2
q
2
jq
(1 N )
2
jq
N
2
q
2
q
2
q
2
) =e
N
2
jq
=e
e
e
q
j
2
jqN
jq
=e
jq
N
2
e
e
jqN
jq
N
2
e e j q e
e
e
jq
N
2
(e
(e
q
2
jq
[(N 1)]
2
N
2
jq
e
=e
q
2
jq
q
2
q
2
q
2
jq
j q (N 1)
2
N
2
Nx
)
sen( x)
sen ( x)
2
DFT cos =
, o
, o
x
x
x
sen( )
sen( )
N
2
N
sen(

N sen((k m))
2
(k m)
Demostracin: X (m)
e
N 1
(e ) = e
jq n
n=0
jq
N
2
e
e
q
2
N
2
(e
jq
=e
(e
jq
N
2
j
e
q
2
j q jq
N
2
2
j q (N 1)
2
N
jq
2
e
q
2
q
j
2
q
2
q
2
e
e
N
2
q
2
N
2
q
2
q
2
jq
N
2
e
e
jqN
jq
N
2
q
2
e
e
jq
N
2
j
(e
q
2
jq
[(N 1)]
2
N
2
jq
e
=e
q
2
jq
e e j q e
(e
=e
jq
jq
q
2
) =e
jq
(1 N )
2
N
2
jq
jq
N
2
e
N
2
=e
(e
(e
jq
jqN
jq
q
2
q
2
jq
j q (N 1)
2
N
2
Nx
)
sen( x)
sen ( x)
2
DFT cos =
, o
, o
x
x
x
sen( )
sen( )
N
2
N
sen(

N sen((k m))
2
(k m)
Demostracin: X (m)
e
N 1
(e ) = e
jq n
n=0
jq
N
2
e
e
q
2
N
2
(e
jq
=e
(e
jq
N
2
j
e
q
2
j q jq
N
2
2
j q (N 1)
2
N
jq
2
e
q
2
q
j
2
q
2
q
2
e
e
q
2
N
2
q
2
q
2
jq
N
2
e
e
jqN
jq
e
e
jq
N
2
j
(e
q
2
jq
[(N 1)]
2
N
2
jq
e
=e
q
2
N
2
q
2
jq
e e j q e
(e
=e
N
jq
2
jq
q
2
) =e
jq
(1 N )
2
N
jq
2
j
jq
N
2
e
N
2
=e
(e
(e
jq
jqN
jq
Nx
)
sen( x)
sen ( x)
2
DFT cos =
, o
, o
x
x
x
sen( )
sen( )
N
2
N
sen(
q
2
q
2
jq
N
2
j q (N 1)
2
Identidad de Euler
sen( t)=
jt
e
2j
jt
Nx
)
sen(
x)
sen
(
x)
2
N sen((k m)) DFT =
,
o
,
o
Demostracin: X (m)
cos
x
x
x
2
(k m)
sen( )
sen( )
N
2
N
N
N
N
jq
jq
jq
jq
2
N
e 2
e 2 e j q N e 2
jqN e
2
j
sen
q
e
j
q
(N
1)
N
N
N
jq
jq
jq
2
N 1
2
2
2
2
=e
n
e
e
=
q
( e j q ) = e j q
q
q
q
2
j
sen
j
j
j
n=0
2
e 2 jq e 2
e 2 e j q e 2
e
q
q
q
sen(
( )
()
e
=
q
2
N
2
(e
jq
=e
(e
jq
N
2
j
e
q
2
j q jq
N
2
2
j q (N 1)
2
2
N
2
q
2
=e
) =e
jq
(1 N )
2
N
jq
2
(e
(e
jq
q
2
e
e
q
2
q
2
e
e
jq
N
2
(e
=e
N
jq
2
(e
q
2
jq
[(N 1)]
2
N
2
jq
e
=e
q
2
jq
N
2
j q (N 1)
2
Identidad de Euler
sen( t)=
jt
e
2j
jt
Nx
)
sen(
x)
sen
(
x)
2
N sen((k m)) DFT =
,
o
,
o
Demostracin: X (m)
cos
x
x
x
2
(k m)
sen( )
sen( )
N
2
N
N
N
N
jq
jq
jq
jq
2
N
e 2
e 2 e j q N e 2
jqN e
2
j
sen
q
e
j
q
(N
1)
N
N
N
jq
jq
jq
2
N 1
2
2
2
2
=e
n
e
e
=
q
( e j q ) = e j q
q
q
q
2
j
sen
j
j
j
n=0
2
e 2 jq e 2
e 2 e j q e 2
e
2
(k
m)
q
q
q
sen(
( )
()
e
=
q
2
N
2
(e
jq
=e
(e
jq
N
2
j
e
q
2
j q jq
N
2
2
j q (N 1)
2
2
N
2
q
2
=e
) =e
jq
(1 N )
2
N
jq
2
(e
(e
jq
q
2
e
e
q
2
q
2
e
e
jq
N
2
(e
=e
N
jq
2
(e
q
2
jq
[(N 1)]
2
N
2
jq
e
=e
q
2
q=
jq
N
2
j q (N 1)
2
Identidad de Euler
sen( t)=
jt
e
2j
jt
Nx
)
sen(
x)
sen
(
x)
2
N sen((k m)) DFT =
,
o
,
o
Demostracin: X (m)
cos
x
x
x
2
(k m)
sen( )
sen( )
N
2
N
N
N
N
jq
jq
jq
jq
2
N
e 2
e 2 e j q N e 2
jqN e
2
j
sen
q
e
j
q
(N
1)
N
N
N
jq
jq
jq
2
N 1
2
2
2
2
=e
n
e
e
=
q
( e j q ) = e j q
q
q
q
2
j
sen
j
j
j
n=0
2
e 2 jq e 2
e 2 e j q e 2
e
2
(k
m)
q
q
q
sen(
( )
()
e
=
q
2
N
jq
2
=e
(e
jq
N
2
(e
q
j
2
j q jq
N
2
2
j q (N 1)
2
2
N
2
q
j
2
=e
) =e
jq
(1 N )
2
N
jq
2
(e
(e
jq
q
2
e
e
q
2
q
2
e
e
jq
N
2
(e
=e
N
jq
2
(e
q
j
2
jq
[(N 1)]
2
N
2
jq
e
=e
q
j
2
q=
jq
N
2
j q (N 1)
2
Identidad de Euler
sen( t)=
jt
e
2j
jt
=e
2 (k m)
( N 1)
N
2
2 (k m) N
N
2
(
)
2 (k m) 1
sen (
N
2)
sen
Nx
)
sen(
x)
sen
(
x)
2
N sen((k m)) DFT =
,
o
,
o
Demostracin: X (m)
cos
x
x
x
2
(k m)
sen( )
sen( )
N
2
N
N
N
N
jq
jq
jq
jq
2
N
e 2
e 2 e j q N e 2
jqN e
2
j
sen
q
e
j
q
(N
1)
N
N
N
jq
jq
jq
2
N 1
2
2
2
2
=e
n
e
e
=
q
( e j q ) = e j q
q
q
q
2
j
sen
j
j
j
n=0
2
e 2 jq e 2
e 2 e j q e 2
e
2
(k
m)
q
q
q
sen(
( )
()
e
=
q
2
N
jq
2
=e
(e
jq
N
2
(e
q
j
2
j q jq
N
2
2
j q (N 1)
2
2
N
2
q
j
2
=e
(e
(e
jq
) =e
jq
(1 N )
2
N
2
jq
q
2
e
e
q
2
N
2
q
2
e
e
jq
N
2
(e
=e
jq
(e
q
j
2
jq
[(N 1)]
2
N
2
jq
e
=e
q
j
2
q=
jq
N
2
j q (N 1)
2
sen( t)=
jt
e
2j
=e
=e
Identidad de Euler
e
jt
2 (k m)
( N 1)
N
2
(k m)
(N 1)
N
2 (k m) N
N
2
(
)
2 (k m) 1
sen (
N
2)
sen
sen (( k m))
(k m)
sen
N
Nx
)
sen(
x)
sen
(
x)
2
N sen((k m)) DFT =
,
o
,
o
Demostracin: X (m)
cos
x
x
x
2
(k m)
sen( )
sen( )
N
2
N
N
N
N
jq
jq
jq
jq
2
N
e 2
e 2 e j q N e 2
jqN e
2
j
sen
q
e
j
q
(N
1)
N
N
N
jq
jq
jq
2
N 1
2
2
2
2
=e
n
e
e
=
q
( e j q ) = e j q
q
q
q
2
j
sen
j
j
j
n=0
2
e 2 jq e 2
e 2 e j q e 2
e
2
(k
m)
q
q
q
sen(
( )
()
e
=
q
2
N
jq
2
=e
(e
jq
N
2
(e
q
j
2
j q jq
N
2
2
j q (N 1)
2
2
N
2
q
j
2
=e
(e
(e
jq
) =e
jq
(1 N )
2
N
2
jq
q
2
e
e
q
2
N
2
q
2
e
e
jq
N
2
(e
=e
jq
(e
q
j
2
jq
[(N 1)]
2
N
2
jq
e
=e
q
j
2
q=
jq
N
2
j q (N 1)
2
=e
=e
2 (k m)
( N 1)
N
2
(k m)
(N 1)
N
sen( t)=
jt
e
2j
jt
=e
sen (( k m))
(k m)
sen
N
Identidad de Euler
j (k m)
2 (k m) N
N
2
(
)
2 (k m) 1
sen (
N
2)
sen
(k m)
N
] sen (( k m))
sen
(k m)
N

Demostracin: X (m)
N 1
(e
[
) =e
j (k m)
jq n
(k m)
N
n=0
N 1
Xr
N sen((k m))
2
(k m)
] sen( (k m))
sen
j
1
X r (m)= e
2 n=0
[
(m)=e
2 n (k m)
N
j (k m)
(k m)
N
Nx
)
sen( x)
sen ( x)
2
DFT cos =
, o
, o
x
x
x
sen( )
sen( )
N
2
sen(
1
+
2
( k m)
N
N 1
q=
2 (k m)
N
2 n(k + m)
N
n=0
] 1 sen ((k m)) +e [
j (k m)
(k m)
sen
N
(k m)
N
] 1 sen ((k m))

2
sen
(k m)
N

Suponga que una sinusoide de frecuencia 8 KHz se muestrea a 32 KHz y se
obtienen 32 muestras


Suponga que la seal
muestreada es de
8.5 KHz

Suponga que la seal
muestreada es de
8.5 KHz

Suponga que la seal
muestreada es de
8.5 KHz
Suponga ahora que

la seal muestreada
es de 8.75 KHz

Suponga que la seal
muestreada es de
8.5 KHz
Suponga ahora que

la seal muestreada
es de 8.75 KHz

El espectro de frecuencia de una seal continua es peridico:
F =3.4 Hz

F =3.4 Hz

F =3.4 Hz

F =3.4 Hz
F =3.4 Hz
F =28.6 Hz
F =16.4 Hz
Ventanas
Las ventanas reducen las fugas en la DFT al minimizar la magnitud de los
lbulos laterales de la funcin sinc
N sen((k m))
X (m)
2
(k m)
Esto es lo que se busca reducir
Esto se logra al forzar a que la amplitud al inicio y al final de la seal de entrada

alcance suavemente un mismo valor de amplitud
Ventanas
Las ventanas reducen las fugas en la DFT al minimizar la magnitud de los
lbulos laterales de la funcin sinc
N sen((k m))
X (m)
2
(k m)
Esto es lo que se busca reducir
Esto se logra al forzar a que la amplitud al inicio y al final de la seal de entrada

alcance suavemente un mismo valor de amplitud
Ventanas
Rectangular
Triangular
Hanning
Hamming
Ventanas
Rectangular
w(n)=1 ;
n=0, 1, 2,, N 1
Triangular
n
; n=0, 1, 2, , N / 2
w(n)= N /2
n
2
; n= N /2, , N
N /2
Hanning
w(n)=0.50.5 cos(
2n
) ; n=0,1, 2,, N 1
N
Hamming
w(n)=0.540.46 cos(
2n
); n=0, 1, 2, , N 1
N
Ventanas
Rectangular
w(n)=1 ;
n=0, 1, 2,, N 1
Triangular
n=0, 1, 2, , N 1
Hanning
w(n)=0.50.5 cos(
2n
) ; n=0,1, 2,, N 1
N
Hamming
w(n)=0.540.46 cos(
2n
); n=0, 1, 2, , N 1
N
Ventanas: respuesta en frecuencia

N 1
W (m)= w (n)e
n=0
j2 nm
N

N 1
W (m)= w (n)e
j2 nm
N
n=0
W mag (m)=W (m)

Magnitud de la respuesta en frecuencia en
escala lineal
Rectangular
Hamming
Triangular
Hanning
fs
N
2 fs
N
3fs
N
4 fs
N
5fs
N

N 1
W (m)= w (n)e
j2 nm
N
n=0
W dB (m)=20log10
W mag (m)=W (m)

escala lineal
W (m)
W (0)

escala logartmica
Rectangular
Rectangular
Hanning
Hamming
Hamming
Triangular
Hanning
fs
N
2 fs
N
3fs
N
4 fs
N
5fs
N
fs
N
2 fs
N
3fs
N
4 fs
N
5fs
N
Triangular

Prctica 4.1.7. Objetivo: Observar la respuesta en frecuencia de distintas
funciones ventana
1. Calcule las funciones ventana Rectangular, Triangular, Hanning y Hamming.
Grafique su amplitud en subgrficas.
2. Calcule la DFT de 512 puntos de las funciones ventana.
3. Grafique la magnitud de las DFT de las funciones ventana
4. Grafique la magnitud de las DFT de las funciones ventana en escala logartmica
(en decibeles).
Ventanas: minimizacin de fugas en la DFT

N 1
X w (m)= w(n)x (n)e

n=0
j 2 nm
N

N 1
X w (m)= w(n)x (n)e
j 2 nm
N
n=0
Funcin ventana
x (n)=[0.50.5 cos (
Sinusoide de 3.4
2n
2 3.4 n
)][sen (
)]
64
64

N 1
X w (m)= w(n)x (n)e

n=0
Ventana Hanning
Funcin ventana
x (n)=[0.50.5 cos (
j 2 nm
N
Sinusoide de 3.4
2n
2 3.4 n
)][sen (
)]
64
64

N 1
X w (m)= w(n)x (n)e

n=0
Ventana Hanning
Funcin ventana
x (n)=[0.50.5 cos (
j 2 nm
N
Sinusoide de 3.4
2n
2 3.4 n
)][sen (
)]
64
64

N 1
X w (m)= w(n)x (n)e

n=0
j 2 nm
N
Ventana Hanning
Funcin ventana
Funcin ventana
x (n)=[0.50.5 cos (
Sinusoide de 3.4
2n
2 3.4 n
)][sen (
)]
64
64
x (n)=[0.50.5 cos (
Sinusoide de 3.4 y 7 ciclos
2n
2 3.4 n
2 7 n
)][sen (
)+ 0.1 sen(
)]
64
64
64

N 1
X w (m)= w(n)x (n)e

n=0
j 2 nm
N
Ventana Hanning
Funcin ventana
Funcin ventana
x (n)=[0.50.5 cos (
Sinusoide de 3.4
2n
2 3.4 n
)][sen (
)]
64
64
x (n)=[0.50.5 cos (
2n
2 3.4 n
2 7 n
)][sen (
)+ 0.1 sen(
)]
64
64
64

N 1
X w (m)= w(n)x (n)e

n=0
j 2 nm
N
Ventana Hanning
Funcin ventana
Funcin ventana
x (n)=[0.50.5 cos (
Sinusoide de 3.4
2n
2 3.4 n
)][sen (
)]
64
64
x (n)=[0.50.5 cos (
2n
2 3.4 n
2 7 n
)][sen (
)+ 0.1 sen(
)]
64
64
64
DFT de una funcin rectangular general

Nx
)
sen( x)
sen( x)
2
DFT rect =
, o
, o
x
x
x
sen( )
sen ( )
N
2
sen (
N
2
X (m)=
n=
x (n) e
j2 nm
N
x (n)
N
+1
2
N / 2+ 1
N /2
0

Nx
)
sen( x)
sen( x)
2
DFT rect =
, o
, o
x
x
x
sen( )
sen ( )
N
2
sen (
N
2
X (m)=
n=
x (n) e
j2 nm
N
x (n)
N
+1
2
K
N / 2+ 1
N /2
0
n=n 0
n=n0 n=n 0+ ( K 1)

Nx
)
sen( x)
sen( x)
2
DFT rect =
, o
, o
x
x
x
sen( )
sen ( )
N
2
sen (
N
2
X (m)=
n=
x (n) e
x (n)
N
+1
2
n 0+ ( K 1)
X (m)=
j2 nm
N
1e
j 2n m
N
N / 2+ 1
n=n0
N /2
0
n=n 0
n=n0 n=n 0+ ( K 1)

Nx
)
sen ( x )
sen ( x)
2
DFT rect =
, o
, o
x
x
x
sen ( )
sen ( )
N
2
sen (
N
2
X (m)=
n=
x (n) e
x (n)
N
+1
2
n 0+ ( K 1)
X (m)=
j2 nm
N
1e
j 2n m
N
N / 2+ 1
n=n0
0
n=n 0
2 m
q=
N
n 0+ ( K 1)
X (q)=
n=n 0
N /2
1e j q n
n=n 0+ ( K 1)

n 0+ ( K 1)
X (q)=
n=n 0
1e j q n

n 0+ ( K 1)
X (q)=
1e j q n
n=n 0
n 0+ ( K 1)
X (q)=
n=n 0
e j q n

n 0+ ( K 1)
X (q)=
1e j q n
n=n 0
n 0+ ( K 1)
X (q)=
e j q n
n=n 0
=e
j q (n 0)
+e
j q(n0+ 1)
+e
j q(n 0+ 2)
+ + e
j q (n0 + (K 1))

n 0+ ( K 1)
X (q)=
1e j q n
n=n 0
n 0+ ( K 1)
X (q)=
e j q n
n=n 0
=e
j q (n 0)
+e
j q(n0+ 1)
=e
j q (n 0)
+e
j q(n0 )
+e
j q(n 0+ 2)
j 1q
+ + e
+e
j q(n0 )
j q (n0 + (K 1))
j2q
+ + e
j q (n 0)
j (K 1)q

n 0+ ( K 1)
X (q)=
1e j q n
n=n 0
n 0+ ( K 1)
X (q)=
e j q n
n=n 0
=e
j q (n 0)
+e
j q(n0+ 1)
=e
j q (n 0)
+e
j q(n0 )
=e
j q (n 0)
j 0q
+e
+e
j q(n 0+ 2)
+ + e
j 1q
+e
j q(n0 )
j q(n0 )
j 1q
+e
j q (n0 + (K 1))
j2q
j q(n0 )
+ + e
j2q
j q (n 0)
+ + e
j (K 1)q
j q (n 0)
j(K 1)q

n 0+ ( K 1)
X (q)=
1e j q n
n=n 0
n 0+ ( K 1)
X (q)=
e j q n
n=n 0
=e
j q (n 0)
+e
j q(n0+ 1)
=e
j q (n 0)
+e
j q(n0 )
=e
j q (n 0)
=e
j q n0
[ e
j 0q
j0q
+e
+e
+e
j q(n 0+ 2)
+ + e
j 1q
+e
j q(n0 )
j 1q
+e
j2q
j q(n0 )
j 1q
+e
j q (n0 + (K 1))
j2q
+ + e
j q(n0 )
+ + e
j( K 1) q
j2q
j q (n 0)
+ + e
j (K 1)q
j q (n 0)
j(K 1)q

n 0+ ( K 1)
X (q)=
1e j q n
n=n 0
n 0+ ( K 1)
X (q)=
e j q n
n=n 0
=e
j q (n 0)
+e
j q(n0+ 1)
=e
j q (n 0)
+e
j q(n0 )
=e
j q (n 0)
=e
j q n0
X (q)=e
j q n0
[ e
j 0q
j0q
+e
+e
K 1
e j p q
p=0
+e
j q(n 0+ 2)
+ + e
j 1q
+e
j q(n0 )
j 1q
+e
j2q
j q(n0 )
j 1q
+e
j q (n0 + (K 1))
j2q
+ + e
j q(n0 )
+ + e
j( K 1) q
j2q
j q (n 0)
+ + e
j (K 1)q
j q (n 0)
j(K 1)q

n 0+ ( K 1)
X (q)=
1e j q n
n=n 0
n 0+ ( K 1)
X (q)=
e j q n
n=n 0
=e
j q (n 0)
+e
j q(n0+ 1)
=e
j q (n 0)
+e
j q(n0 )
=e
j q (n 0)
=e
j q n0
X (q)=e
j q n0
[ e
j 0q
j0q
+e
+e
K 1
e
p=0
j pq
+e
j q(n 0+ 2)
+ + e
j 1q
+e
j q(n0 )
j 1q
+e
j2q
j q(n0 )
j 1q
+e
+ + e
j( K 1) q
K 1
p=0
j2q
j q(n0 )
+ + e
j q (n0 + (K 1))
j2q
+ + e
K 1
j pq
= ( e j q )
p=0
j q (n 0)
j (K 1)q
j q (n 0)
j(K 1)q

n 0+ ( K 1)
X (q)=
1e j q n
n=n 0
n 0+ ( K 1)
X (q)=
e j q n
n=n 0
=e
j q (n 0)
+e
j q(n0+ 1)
=e
j q (n 0)
+e
j q(n0 )
=e
j q (n 0)
=e
j q n0
X (q)=e
j q n0
[ e
j 0q
j0q
+e
+e
K 1
e
p=0
j pq
+e
j q(n 0+ 2)
+ + e
j 1q
+e
j q(n0 )
j 1q
+e
j2q
j q(n0 )
j 1q
+e
+ + e
j( K 1) q
K 1
p=0
j2q
j q(n0 )
+ + e
j q (n0 + (K 1))
j2q
= ( e j q )
p=0
+ + e
K 1
j pq
j q (n 0)
j (K 1)q
j q (n 0)
j(K 1)q
Serie geomtrica
p

Series geomtricas
1
n1

r: Es la razn comn
n1
s= a r k
k =0
1
n1
s (1r)=aar
n
1r
s=a
1r
+ ar
n1
r
1
r =r
3
n1+ 1
=r
ar ar ar ar ar
Si r< 1 y n :
n1
s=
a
1r

K 1
jq p
(e )
p=0

K 1
jq p
(e )
p=0
n1
1r n
s= a r =a
1r
k =0
k

K 1
jq p
(e )
p=0
K 1
(e
p=0
1e j q K
)=
1e jq
jq p
n1
1r n
s= a r =a
1r
k =0
k

K 1
jq p
(e )
p=0
K 1
(e
p=0
1e j q K
)=
1e jq
jq p
Multiplicando por 1:
n1
1r n
s= a r =a
1r
k =0
k

K 1
n1
1r n
s= a r =a
1r
k =0
jq p
(e )
p=0
K 1
(e
p=0
1e j q K
)=
1e jq
jq p
e
jqK
jqK
2
jqK
2
1e
e
=
jq
jq
1e
e2
e
jq
2

K 1
n1
1r n
s= a r =a
1r
k =0
jq p
(e )
p=0
K 1
(e
1e j q K
)=
1e jq
jq p
p=0
jqK
2
jqK
2
jqK
1e
e
=
jq
jq
1e
e2
e
e
=
jqK
2
e
e
jqK
jq
2
jqK
2
jqK
2
jq
2
e e
e
jq
jq
2
jq
2

K 1
n1
1r n
s= a r =a
1r
k =0
jq p
(e )
p=0
K 1
(e
1e j q K
)=
1e jq
jq p
p=0
jqK
2
jqK
2
jqK
1e
e
=
jq
jq
1e
e2
e
e
=
jqK
2
e
e
jqK
jq
2
jqK
2
jqK
2
jq
2
e e
e
jq
jq
2
jq
2
jq
2
(e
jqK
2
jqK
2
(e
jq
2
jqK
2
jq
2

K 1
n1
1r n
s= a r =a
1r
k =0
jq p
(e )
p=0
K 1
(e
1e j q K
)=
1e jq
jq p
p=0
jqK
2
jqK
2
jqK
1e
e
=
jq
jq
1e
e2
e
e
=
jqK
2
e
e
jqK
jq
2
e e
e
jq
jq
2
jqK
2
jq
2
(e
jqK
2
jqK
2
jqK
2
jq
2
jq
2
jqK
2
(e
jq
2
jqK
2
(e
(e
jq
2
jqK
2
jq
2
jq
2
jqK
2
jq
2

K 1
n1
1r n
s= a r =a
1r
k =0
jq p
(e )
p=0
K 1
(e
1e j q K
)=
1e jq
jq p
p=0
jqK
2
jqK
2
jqK
1e
e
=
jq
jq
1e
e2
e
e
=
jqK
2
e
e
jqK
jqK
2
e e
e
jq
jq
2
jq
2
K 1
e
p=0
jqp
=e
(e
jqK
2
jqK
2
jqK
2
jq
2
jq
2
jq
2
jqK
2
jq
2
(e
jqK
2
(e
(e
jq
2
jqK
2
jq
2
jqK
2
jq
2
jqK
2
( e
(e
j q(K 1)
2
jq
2
jq
2
e
e
jqK
2
jq
2

K 1
e
p=0
jqp
=e
jqK
2
( e
(e
j q( K 1)
2
jq
2
e
e
jqK
2
jq
2

K 1
jqp
=e
p=0
e e
sen()=
2j
jqK
2
( e
(e
j q (K 1)
2
jq
2
e
e
jq K
2
jq
2

K 1
jqp
=e
p=0
e e
sen()=
2j
K 1
e
p=0
jqp
=e
jqK
2
( e
(e
j q( K 1)
2
jq
2
e
e
jqK
2
jq
2
jqK
2
( e
(e
j q( K 1)
2
jq
2
e
e
jqK
2
jq
2

K 1
jqp
=e
p=0
e e
sen()=
2j
K 1
jqp
=e
jq
2
e
e
jqK
2
jq
2
jqK
2
( e
(e
j q( K 1)
2
p=0
=e
jqK
2
( e
(e
j q( K 1)
2
jq
2
e
e
jqK
2
jq
2
qK
2
q
2 j sen
2
j q (K 1)
2
2 j sen
( )
()

K 1
jqp
=e
p=0
e e
sen()=
2j
K 1
jqp
=e
K 1
e
p=0
jqp
=e
jq
2
e
e
jqK
2
jq
2
jqK
2
( e
(e
j q( K 1)
2
p=0
=e
jqK
2
( e
(e
j q( K 1)
2
jq
2
e
e
jqK
2
jq
2
qK
2
q
2 j sen
2
( )
()
qK
sen (
2 )
q
sen ( )
2
j q (K 1)
2
j q( K 1)
2
2 j sen

K 1
e j q p=e
p=0
e e
sen()=
2j
K 1
jqp
=e
K 1
e
p=0
jqp
=e
jq
2
e
e
jqK
2
jq
2
jqK
2
( e
(e
j q( K 1)
2
jq
2
e
e
jqK
2
jq
2
qK
2
q
2 j sen
2
( )
()
qK
sen (
2 )
q
sen ( )
2
j q (K 1)
2
j q( K 1)
2
2 j sen
X (q)=e
j q n0
K 1
p=0
p=0
=e
jqK
2
( e
(e
j q( K 1)
2
e j p q

K 1
e j q p=e
p=0
e e
sen()=
2j
K 1
jqp
=e
K 1
e
p=0
jqp
=e
jq
2
e
e
jqK
2
jq
2
jqK
2
( e
(e
j q( K 1)
2
jq
2
e
e
jqK
2
jq
2
qK
2
q
2 j sen
2
( )
()
qK
sen (
2 )
q
sen ( )
2
j q (K 1)
2
j q( K 1)
2
2 j sen
X (q)=e
j q n0
K 1
p=0
p=0
=e
jqK
2
( e
(e
j q( K 1)
2
e j p q

K 1
e j q p=e
p=0
sen()=
K 1
jqp
e e
2j
=e
K 1
e
p=0
jqp
=e
jq
2
e
e
jqK
2
jq
2
X (q)=e
j q n0
( e
(e
j q( K 1)
2
jq
2
e
e
jqK
2
jq
2
qK
2
q
2 j sen
2
( )
()
qK
sen (
2 )
q
sen ( )
2
j q (K 1)
2
j q( K 1)
2
2 j sen
e j p q
qK
sen
j q (K 1)
2
jqn
2
X (q)=e e
q
sen
2
0
jqK
2
K 1
p=0
p=0
=e
jqK
2
( e
(e
j q( K 1)
2
( )
()

K 1
e j q p=e
p=0
sen()=
K 1
jqp
e e
2j
=e
K 1
e
p=0
jqp
=e
jq
2
e
e
jqK
2
jq
2
X (q)=e
j q n0
K 1
e j p q
p=0
qK
sen
j q (K 1)
2
jqn
2
X (q)=e e
q
sen
2
jqK
2
( e
(e
j q( K 1)
2
p=0
=e
jqK
2
( e
(e
j q( K 1)
2
jq
2
e
e
jqK
2
jq
2
qK
2
q
2 j sen
2
( )
()
qK
sen (
2 )
q
sen ( )
2
j q (K 1)
2
j q( K 1)
2
2 j sen
)
X (q)=e
j q n 0
K 1
2
( )
()
qK
sen (
2 )
)
q
sen ( )
2

K 1
e j q p=e
p=0
sen()=
K 1
jqp
e e
2j
=e
jqK
2
( e
(e
j q( K 1)
2
jq
2
e
e
jqK
2
jq
2
e
p=0
=e
e j p q
qK
sen
j q (K 1)
2
jqn
2
X (q)=e e
q
sen
2
0
jqK
2
( e
(e
j q( K 1)
2
jq
2
e
e
jqK
2
jq
2
qK
j q (K 1) 2 j sen
2
=e 2
q
2 j sen
2
jqp
X (q)=e
K 1
p=0
p=0
K 1
j q n0
( )
()
qK
sen (
2 )
q
sen ( )
2
j q( K 1)
2
)
X (q)=e
q=
2 m
N
j q n 0
K 1
2
( )
()
qK
sen (
2 )
)
q
sen ( )
2

K 1
e j q p=e
p=0
sen()=
K 1
jqp
e e
2j
=e
jqK
2
( e
(e
j q( K 1)
2
jq
2
e
e
jqK
2
jq
2
e
p=0
=e
e j p q
qK
sen
j q (K 1)
2
jqn
2
X (q)=e e
q
sen
2
0
jqK
2
( e
(e
j q( K 1)
2
jq
2
e
e
jqK
2
jq
2
qK
j q (K 1) 2 j sen
2
=e 2
q
2 j sen
2
jqp
X (q)=e
K 1
p=0
p=0
K 1
j q n0
( )
()
qK
sen (
2 )
q
sen ( )
2
j q( K 1)
2
)
X (q)=e
q=
j q n 0
K 1
2
2 m
N
( )
()
qK
sen (
2 )
)
q
sen ( )
2
2m K
sen (
2N )
)
2m
sen (
2N )
2m
K 1
n
(
)(
2
X (m)=e N
j

K 1
e j q p=e
p=0
sen()=
K 1
jqp
e e
2j
=e
jqK
2
( e
(e
j q( K 1)
2
jq
2
e
e
jqK
2
jq
2
e
p=0
=e
e j p q
qK
sen
j q (K 1)
2
jqn
2
X (q)=e e
q
sen
2
0
jqK
2
( e
(e
j q( K 1)
2
jq
2
e
e
jqK
2
jq
2
qK
j q (K 1) 2 j sen
2
=e 2
q
2 j sen
2
jqp
X (q)=e
K 1
p=0
p=0
K 1
j q n0
( )
()
qK
sen (
2 )
q
sen ( )
2
j q( K 1)
2
)
X (q)=e
q=
j q n 0
K 1
2
2 m
N
( )
()
qK
sen (
2 )
)
q
sen ( )
2
2m K
sen (
2N )
)
2m
sen (
2N )
mK
sen (
N )
)
m
sen (
N )
2m
K 1
n
(
)(
2
X (m)=e N
j
2m
K 1
n
(
)(
2
X (m)=e N
j

K 1
e j q p=e
p=0
sen()=
K 1
jqp
e e
2j
=e
jqK
2
( e
(e
j q( K 1)
2
jq
2
e
e
jqK
2
jq
2
e
p=0
=e
e j p q
qK
sen
j q (K 1)
2
jqn
2
X (q)=e e
q
sen
2
0
jqK
2
( e
(e
j q( K 1)
2
jq
2
e
e
jqK
2
jq
2
qK
j q (K 1) 2 j sen
2
=e 2
q
2 j sen
2
jqp
X (q)=e
K 1
p=0
p=0
K 1
j q n0
( )
()
qK
sen (
2 )
q
sen ( )
2
j q( K 1)
2
)
X (q)=e
q=
j q n 0
K 1
2
2 m
N
( )
()
qK
sen (
2 )
)
q
sen ( )
2
2m K
sen (
2N )
)
2m
sen (
2N )
mK
sen (
N )
)
m
sen (
N )
2m
K 1
n
(
)(
2
X (m)=e N
j
Forma General del

Kernel de Dirichlet
2m
K 1
n
(
)(
2
X (m)=e N
j
Tarea
Tarea 4.1. Objetivo: Obtener la DFT de una funcin rectangular simtrica y de
una funcin rectangular donde todos sus elementos son 1.
1. Basndose en la seccin 3.13.2 y 3.13.3 del libro Understanding Digital Signal
Processing de R. Lyons, obtenga la DFT de una funcin rectangular simtrica y de
una funcin rectangular donde todos sus muestras tienen amplitud 1.
2. Realice el desarrollo en su cuaderno. Escanee o fotografe el resultado,
fabrique un archivo en formato PDF y suba al sitio del curso.
DFT Inversa de una funcin rectangular general

X (m)
X (m) x (n)
K
N / 2+ 1
N /2
m
0
m=m0
m=m0 + ( K 1)

X (m)
X (m) x (n)
K
N / 2+ 1
1
x (n)=
N
N
2
m=
N
+1
2
X (m) e
j 2 mn
N
N /2
m
0
m=m0
m=m0 + ( K 1)

X (m)
X (m) x (n)
K
N / 2+ 1
1
x (n)=
N
N
2
m=
X (m) e
j 2 mn
N
N
+1
2
n K
sen
2n
K 1
j(
m
N
N )(
2 ) 1
x (n)=e

N
n
sen
N
0
( )
( )
N /2
m
0
m=m0
m=m0 + ( K 1)
DFT Inversa de una funcin rectangular simtrica

H (m)
K
N / 2+ 1
m=m0 =( K 1)/2
N /2
0
m
m=( K 1)/2

H (m)
n K
sen
2n
K 1
j(
m
N
N )(
2 ) 1
N / 2+ 1
x (n)=e

N
n
sen
N
0
( )
( )
IDFT de una funcin rectangular general
K
1
m=m0 =( K 1)/2
N /2
0
m
m=( K 1)/2

H (m)
n K
sen
2n
K 1
j(
m
N
N )(
2 ) 1
N / 2+ 1
x (n)=e

N
n
sen
N
( )
( )

Sustituyendo el valor de
m0 =
K 1
2
K
1
m=m0 =( K 1)/2
N /2
0
m
m=( K 1)/2

H (m)
n K
sen
2n
K 1
j(
m
N
N )(
2 ) 1
N / 2+ 1
x (n)=e

N
n
sen
N
( )
( )

m0 =
m=m0 =( K 1)/2
K 1
2
nK
sen
2 n K 1 K 1
j(
N
N )( 2
2 ) 1
h(n)=e

N
n
sen
N
( )
( )
N /2
0
m
m=( K 1)/2

H (m)
n K
sen
2n
K 1
j(
m
N
N )(
2 ) 1
N / 2+ 1
x (n)=e

N
n
sen
N
( )
( )

m0 =
( )
( )
( )
( )
m=m0 =( K 1)/2
K 1
2
nK
sen
2 n K 1 K 1
j(
N
N )( 2
2 ) 1
h(n)=e

N
n
sen
N
n K
sen
N
1
h(n)=
N
n
sen
N
K
N /2
0
m
m=( K 1)/2
Tarea
Tarea 4.2. Objetivo: Comprender el algoritmo de la transformada rpida de
Fourier (FFT).
1. Haga un resumen del captulo 4 del libro Understanding Digital Signal
Processing de R. Lyons. Ponga especial cuidado en incluir los diagramas y las
ecuaciones ms importantes.
2. Realice el resumen en su cuaderno. Escanee o fotografe el resultado,
fabrique un archivo en formato PDF y suba al sitio del curso.
Prctica integradora de la unidad

Prctica 4. Objetivo: Comprender la utilidad de la DFT en el anlisis de seales
mediante la obtencin del contenido frecuencial de una seal continua en el
tiempo formada por dos sinusoides de diferente frecuencia
1. Tome 201 puntos de una seal que est formada por la suma de dos sinusoides
de frecuencias 250 y 500 Hz, la cual es muestreada a 2000 Hz.
2. Calcule la DFT de 4096 puntos de la seal discreta.
3. Calcule el vector de frecuencias de anlisis correspondiente a la mitad de las
frecuencias de la DFT.
4. Grafique la seal continua en una primera subgrfica. Grafique la magnitud de
la mitad de las frecuencias de la DFT de la seal en otra subgrfica.
5. Aplique una funcin ventana sobre la seal y vuelva a calcular la DFT.
6. Grafique de nueva cuenta la DFT de la seal. Qu puede concluir?

Transformada de Fourier Discreta

Cargado por

Copyright:

Formatos disponibles

Transformada de Fourier Discreta

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

Transformada de Fourier Discreta

Cargado por

Copyright:

Formatos disponibles

Universidad Politcnica de Victoria

Procesamiento Digital de Seales

Yahir Hernndez Mier

Estudiar la definicin de la Transformada de Fourier Discreta

Estudiar las propiedades de la Transformada de Fourier Discreta

Aprender a interpretar los resultados espectrales de la DFT

Aplicar ventanas para reducir las fugas de la DFT

Aplicar la DFT a sinusoides reales y complejas

Transformada de Fourier continua

Transformada de Fourier continua

Transformada de Fourier discreta

La DFT es un procedimiento matemtico utilizado para determinar el

Transformada de Fourier discreta

Transformada de Fourier discreta

Transformada de Fourier discreta

Cul es la suma total?

Transformada de Fourier discreta

Cul es la suma total?

Transformada de Fourier discreta

Cul es la suma total?

Transformada de Fourier discreta

Cul es la suma total?

Transformada de Fourier discreta

Cul es la suma total?

Transformada de Fourier discreta

Cul es la suma total?

Transformada de Fourier discreta

Cul es la suma total?

Transformada de Fourier discreta

Cul es la suma total?

Transformada de Fourier discreta

Transformada de Fourier discreta

a: Es el primer trmino de la serie

Transformada de Fourier discreta

a: Es el primer trmino de la serie

Transformada de Fourier discreta

a: Es el primer trmino de la serie

Transformada de Fourier discreta

a: Es el primer trmino de la serie

Transformada de Fourier discreta

a: Es el primer trmino de la serie

Transformada de Fourier discreta

a: Es el primer trmino de la serie

Transformada de Fourier discreta

a: Es el primer trmino de la serie

Transformada de Fourier discreta

Transformada de Fourier discreta

Transformada de Fourier discreta

Transformada de Fourier discreta de funciones elementales

Transformada de Fourier discreta de funciones elementales

Funcin escaln unitario 1( n t s )

Transformada de Fourier discreta de funciones elementales

Transformada de Fourier discreta de funciones elementales

Transformada de Fourier discreta de funciones elementales

e j t =cos ( t)+ j sen ( t )

Transformada de Fourier discreta

Transformada de Fourier discreta

Transformada de Fourier discreta

X (m)= x (n) cos