2.
Considere la siguiente señal periódica
x ( n )= {… , 1 , 0 ,1 , 2 ,3 , 2 ,1 , 0 , 1, … }
a. Dibuje la señal x(n), el módulo y la fase de su espectro.
b. Utilizando los resultados del apartado (a), verifique la relación de Parseval.
Solución:
a. Se tiene la señal:
x ( n )= {… , 1 , 0 ,1 , 2 ,3 , 2 ,1 , 0 , 1, … }
Entonces se toman:
x ( n )= {3 , 2 ,1 , 0 , 1, 2 }
Con un periodo N=6.
Se sabe que:
N −1 2π
1 −j
N
kn
C k= ∑ x [n]e
N n=0
Se tendrán 6 coeficientes de tal modo que:
5
1 1
C 0= ∑ x [n]= ( x [ 0 ] + x [ 1 ] + x [ 2 ] + x [ 3 ] + x [ 4 ] + x [ 5 ] )
6 n=0 6
1 9
C 0= ( 3+2+1+0+ 1+ 2 )= =1.5
6 6
Entonces para el siguiente coeficiente:
5 π π π 2π 5π
1 −j n 1 −j −j −j −j
C 1= ∑ x [n]e 6 = (3+2 e 6 +e 3 + e 3 + 2e 6 )
6 n=0 6
Obteniendo así:
C 1=0.667
De igual modo para los siguientes coeficientes:
C 2=0
C 3=1.667
C 4=0
C 5=0.667
Utilizando el siguiente código en Matlab se obtuvo:
clc
clear all
close all
x=[3 2 1 0 1 2]
N=length(x);
[C,F]=Final(x)
modC=abs(C)
n=0:N-1;
subplot(311)
stem(n,x)
title('x[n]')
subplot(312)
stem(n,abs(C))
title('Espectro de amplitud')
subplot(313)
stem(n,F)
title('Espectro de fase')
Se pueden verificar que los valores de los coeficientes corresponden a los de la
gráfica.
b. Relación de Parseval
Se tiene que la relación de Parseval es:
N−1 N −1
∑ |x [ n ]| = N1 ∑ |X [ k ]|
2 2
n =0 n=0
Reemplazando:
5 5
∑|x [ n ]| = 16 ∑ |X [ k ]|
2 2
n=0 n=0
Con:
x [ n ] ={ 3 ,2 , 1 , 0 ,1 , 2 }
Y:
X [ k ]= {1.5 , 0.667,0,1 .667,0,0.667 }
Al resolver:
5
2
∑|x [ n ]| = [3 2+22 +1+0+ 1+ 22 ]=19
n=0
5
1
| X [ k ]| = 1 [ 1.52 +0.667 2+ 0+1.6672 +0+0.667 2 ]
2
∑
6 n=0 6
5
1
| X [ k ]| = 1 [ 0.98645 ]=0.1644
2
∑
6 n=0 6
Como se puede apreciar:
5 5
∑|x [ n ]| ≠ 1 ∑ |X [ k ]|
2 2
n=0 6 n=0
Al no cumplir la relación de Parseval, se demuestra que la Transformada de
Fourier de la señal no es unitaria.
4.- Considere la siguiente señal
x ( n )= {1 , 0 ,−1 , 2, 3 }
Cuya transformada de Fourier es X ( ω )=X R ( ω ) + j X I ( ω). Determine y dibuje la
señal y (n) cuya transformada de Fourier es
Y ( ω ) =X I ( ω ) + X R (ω) e j 2 ω
Solución:
Se muestra una función de simetría de la transformada de Fourier, el cual sus
funciones son complejas:
x ( n )=x R ( n )+ j x l ( n ) X ( w )=X R ( w )+ j X l ( w )
Del sistema que nos da en la señal anterior se tiene y en este caso analizaremos
la parte real y la parte imaginaria por separado:
x ( n )= {0,1 ,−1,2,3 }
x ( n ) + x (−n )
x R 1 ( n )= Analizamos el sistema anterior:
2
1 1
x R 1 ( n )={ , 0,1,2, 1,0 , }
2 2
En un sistema complejo:
x ( n )= {1,0 ,−1,2,3 }
x ( n )−x (−n )
x l 1 ( n )=
2
Obtenemos del sistema:
x l 1 ( n )= {12 , 0 ,−2 ,0 , 2 , 0 , 12 }
Analizándolos por las propiedades de la Transformada de Fourier para señales
discretas en el tiempo se tiene:
∞
X ( w )= ∑ x R 1 ( n )∗e− jwn
n=−∞
Para X R ( w ):
3
X R ( w )= ∑ x R 1 ( n )∗e− jwn
n=−3
Para j X l 1 ( w ) :
3
jX l 1 ( w )= ∑ x R 1 ( n )∗e− jwn
n =−3
El problema nos pide determinar y (n), si su transformada de Fourier esta dada por:
Y ( n )= X I ( w )+ X R ( w ) e j 2 w
Entonces para calcular y(n), en este caso aplicaremos la transformada inversa a
dicha función, el cual nos resultará:
y ( n )=F−1 {X I ( w ) }+ F−1 {X R ( w ) e j 2 w }
Como ya analizamos anteriormente el valor de X R ( w) y de X I ( w ) ; en este caso
solo nos quedaría analizar en la función para poder calcular y(n):
y ( n )=F−1 {X n ( w ) }
Entonces:
F−1 { X I ( w ) }= j x I 1 ( n ) F−1 { X R ( w ) }=x R 1 ( n+2 )
Por lo tanto y (n):
y ( n )=F−1 {X I ( w ) }+ F−1 {X R ( w ) e j 2 w } y ( n )= j x I 1 ( n ) + x R 1 ( n+2 )
Analizando el sistema con los resultados anteriores se tiene:
y ( n )= { 12 , 0 , 2−2 j , 2 , 2+2 j , 0 1−42 j , 0 2j }
Programa en Matlab:
clear all;
close all;
clc;
j=sqrt(-1);
n=-5:3;
y(n)=[1/2 0 (2-j)/2 2 (2+j)/2 0 (1-4*j)/2 0 j/2];
stem(n,y(n));
xlabel('t')
ylabel('y[n]')
title('VALORES DE LA SEÑAL')
grid on
5.- Determine y dibuje las transformadas de Fourier X 1 ( ω ) , X 2 ( ω ) , X 3 ( ω ) de las
siguientes señales
a. x 1 ( n )={ 1 ,1 , 1 ,1 , 1 }
b. x 2 ( n )={ 1 , 0 ,1 , 0 ,1 , 0 , 1 ,0 , 1 }
c. x 3 ( n )={ 1 , 0 ,0 , 1 , 0 , 0 ,1 , 0 , 0 ,1 , 0 , 0 ,1 }
d. ¿Existe alguna relación entre X 1 ( ω ) , X 2 ( ω ) , X 3 ( ω ) ? ¿Cuál es su significad
físico?
e. Demuestre que si
x k ( n )=
{ x ( nk ) n
si entero
k
0 en otro caso
Entonces X k (ω)= X (kω)
DESARROLLO:
a) x 1 ( n )={ 1 ,1 , 1 ,1 , 1 }
Aplicando la transformada de Fourier para la señal discreta:
∞
x 1 ( n )={ 1 ,1 , 1 ,1 , 1 } X 1 ( w )= ∑ x1 ( n ) e− jwn
n=−∞
X 1 ( w )=∑ x 1 ( n ) e− jwn
n
Entonces dicha sumatoria tomará todo el valor de n, los cuales nos detallan en
la parte inicial x 1 ( n )={ 1 ,1 , 1 ,1 , 1 }, tomando como valores el n de:
n=0 , ±1 , ±2
Reemplazando dichos valores en la sumatoria:
X 1 ( w )=∑ x 1 ( n ) e− jwn
n
n=0 , ±1 , ±2
Para n=−2 :
e− jwn=e− jw (−2 )=e 2 jw
Para n=−1:
e− jwn=e− jw (−1 )=e jw
Para n=0 :
e− jwn=e− jw (0 )=1
Para n=1:
e− jwn=e− jw (1 )=e− jw
Para n=2:
e− jwn=e− jw (2 )=e−2 jw
Por lo tanto, dicha sumatoria será:
X 1 ( w )=e 2 jw + e jw +1+e− jw +e−2 jw … … … ..(α)
Como sabemos la función trigonométrica coseno en su forma compleja es:
e jw +e− jw
cos ( w )=
2
2
Entonces reemplazando y multiplicando por ( ) a toda la ecuación (α ) :
2
2
X 1 ( w )= ( e ¿ ¿ 2 jw +e jw + 1+ e− jw + e−2 jw ) ¿
2
(e ¿ ¿ 2 jw +e jw +1+e− jw +e−2 jw )
X 1 ( w )=2∗ ¿
2
X 1 ( w )=2∗¿
X 1 ( w )=2 ¿
Por lo tanto tenemos que:
X 1 ( w )=2 cos w+ 2cos 2 w+1
c) x 3 ( n )={ 1 , 0 ,0 , 1 , 0 , 0 ,1 , 0 , 0 ,1 , 0 , 0 ,1 }
Aplicando las propiedades de la transformada de Fourier para señales discretas en
el tiempo se tiene:
x 3 ( n )={ 1 , 0 ,0 , 1 , 0 , 0 ,1 , 0 , 0 ,1 , 0 , 0 ,1 }
∞
X 3 ( w )= ∑ x3 ( n ) e− jwn
n=−∞
X 3 ( w )=∑ x 3 ( n ) e− jwn
n
Entonces dicha sumatoria tomará todos los valores de n, los cuales nos detallan en
la parte inicial x 3 ( n )={ 1 , 0 ,0 , 1 , 0 , 0 ,1 , 0 , 0 ,1 , 0 , 0 ,1 }; tomando como valores de n:
n=0 , ±1 , ±2 , ±3 , ± 4 ,± 5 , ±6
Reemplazando dichos valores en la sumatoria:
X 3 ( w )=∑ x 3 ( n ) e− jwn
n
Evaluando para cada valor de n:
Para n=−6 :
e− jwn=e− jw (−6 )=e 6 jw∗1=e 6 jw
Para n=−5 :
e− jwn=e− jw (−5 )=e 5 jw∗0=0
Para n=−4 :
e− jwn=e− jw (−4)=e 4 jw∗0=0
Para n=−3 :
e− jwn=e− jw (−3 )=e 3 jw∗1=e3 jw
Para n=−2 :
e− jwn=e− jw (−2 )=e 2 jw∗0=0
Para n=−1:
e− jwn=e− jw (−1 )=e jw∗0=0
Para n=0 :
e− jwn=e− jw (0 )∗1=1
Para n=1:
e− jwn=e− jw (1 )=e− jw∗0=0
Para n=2:
e− jwn=e− jw (2 )=e−2 jw∗0=0
Para n=3 :
e− jwn=e− jw (3 )=e−3 jw∗1=e−3 jw
Para n=4 :
e− jwn=e− jw (4 )=e−4 jw∗0=0
Para n=5 :
e− jwn=e− jw (5 )=e−5 jw∗0=0
Para n=6 :
e− jwn=e− jw (6 )=e−6 jw∗1=e−6 jw
Por lo tanto, dicha sumatoria será:
X 3 ( w )=e 6 jw +0+0+e 3 jw +0+0+1+0+0+e−3 jw +0+0+e−6 jw
X 3 ( w )=e 6 jw +e 3 jw +1+ e−3 jw +e−6 jw
Como sabemos la función trigonométrica coseno en su forma compleja es:
e jw +e− jw
cos ( w )=
2
2
Entonces reemplazando y multiplicando por ( 2 ) a toda la ecuación:
2
X 3 ( w )= ( e6 jw +e−6 jw +e 3 jw + e−3 jw +1)
2
(e 6 jw +e−6 jw ) (e 3 jw + e−3 jw ) 1
X 3 ( w )=2 [ 2
+
2
+ ¿
2 ]
X 3 ( w )=2 ¿
X 3 ( w )=2 cos 6 w+ 2cos 3 w+1
E) Demuestre que si
x k ( n )=
{ x ( nk ) n
si entero
k
0 en otro caso
Entonces X k (ω)= X (kω)
Desarrollo:
n
En este caso analizaremos para como valor entero y para los otros casos seria
k
0.
X k ( w )= ∑ x n ( n ) e− jwn
n
n a=
k
X k ( w )=∑ x n ( n ) e− jwnk
n
Por lo tanto, podemos decir que:
X k ( w )=∑ x n ( n ) e− jwnk = X (kw)
n
