CE725 PROCESAMIENTO AVANZADO DE SEÑALES DIGITALES

TAREA 1 (22 de enero de 2018)

1. Determine la señal periódica x(n) con N=8 si sus coeficientes de DTFS

ck= { ….0, 0.25, 0.5, 1,2, 1, 0.5, 0.25, 0 ….}

2. Considere el muestreo de la señal x(t) = (1/ t) sin( t) (a) Dibuja la FT de la señal muestreada
para los intervalos de muestreo =0.25 y =0.5 (b) Sea x(n) = x(n ). Dibuje la DTFT del
señal anterior para las tasas de muestreo dadas arriba.

3. Calcula la DFT de 8 puntos de la secuencia x (n) = { 2 -1 1 -1 2 3 4 3 }

4. Determina la convolución circular de 4 puntos de las secuencias x1(n) = {1, 2, 3, 1}y

x2(n) = {4, 3, 2, 2} usando (a) fórmula de dominio del tiempo (b) DFT

5. Considere la secuencia de longitud finita x[n] = {1, 1, 1, 1, 1, 1} y sea X(z) su transformada z.

2 k
j
Si muestreamos X (z) enzk e para4 k = 0, 1, 2, 3 obtenemos un conjunto de coeficientes DFT X-4(k).
Encuentra la secuencia, x4(n), la IDFT de X-4( k )

6. Encuentra los valores de X[k] y Y[k], dado que v[n] = x[n] + j y[n] y
V[k] = {3+j8, 2+j6, 1–j5, 4+j8, 5+j2, 3–j3, j5, - 3+j2}

7. Calcule la IDFT de 8 puntos de la siguiente secuencia utilizando los algoritmos DIT de radix 2.
Sigue exactamente los gráficos de flujo de señal y haz un seguimiento de todas las cantidades intermedias.
{ 2 1.707-j0.707 0.293-j0.707 1-j 0 1+j 0.293+j0.707 1.707-j0.707 } para
k= 0,1,2,…,7

8. Sea X(k) la DFT de N puntos de la secuencia de N puntos x(n). Calcula X(ej ) y X(k) si
(i) x(n) = x(N-1-n) y N es par
(ii) x(n)= x(N-1-n) y N es impar
(iii) x(n)= -x(N-1-n) y N es par
(iv) x(n) = -x(N-1-n) y N es impar

9. Dibuja el diagrama de flujo de señal para una FFT de 24 puntos. Optimízalo para el número mínimo de
multiplicaciones

10. Sea x(n) una señal arbitraria con DTFT X(ej ). Expresar la DTFT de
x1(n) = x(Dn) en términos de X(ej ) donde D es un entero -Decimación
x2(m) = x(m/U) m=0, U, 2U …
0 de lo contrario Interpolación
en términos de X(ej ).
Ejercicios de Matlab

Para los ejercicios de Matlab, necesitas hacer un informe con tus observaciones y
figuras y súbelo a Moodle. Puedes trabajar en grupos de dos.
Efecto del muestreo y cuantificación
1. Muestree la señal analógica x(t)=sin(10 t) a i) 1000Hz ii) 15 Hz y iii) 7.5 Hz. Grafica los discretos
señal de tiempo y observar el efecto de aliasing.

Solución de muestra

% Demostración de muestreo
t=0:0.001:2;
xa=sin(2*pi*5*t);
trazar(t, xa);
aguanta
%%
fs=15; ts=0:1/fs:2;
xs=sin(2*pi*5*ts);
grafica(ts,xs,'ro-');
%%
fs=7.5; ts=0:1/fs:2;
xs=sin(2*pi*5*ts);
trazar(ts,xs,'go-');

Modifica el código para muestrear la señal analógica x(t)=3cos(100 t) a i) 400Hz ii) 100 Hz y
iii) 60 Hz. Grafique la señal en el tiempo discreto y observe el efecto del aliasing.

Análisis en el dominio del tiempo

2. Considera la señal de tiempo discreto x(n) = cos (2πn/7).
(i) Obtén y grafica la señal muestreada (decimada) y(n) = x(2n) junto con x(n)
(ii) Obtenga y grafique la señal sobremuestreada (interpolada) z(n) = x(n/2) para n par y z(n) =
0de lo contrario junto con x(n)

3. En la marcación DTMF, un número se representa por un tono de doble frecuencia.

(i) Haz una búsqueda en la web y encuentra las frecuencias de cada dígito
(ii) Generar tonos DTMF correspondientes al número de teléfono 0824-2474040 mediante
muestreo la suma de sinusoidales a las frecuencias requeridas a Fs = 8192 Hz.
Concatena las señales poniendo 50 ceros entre cada señal (para representar
silencio) y escucha la señal.

4. Graba tu propia voz durante medio minuto y observa la forma de onda en el dominio del tiempo.
i) Encuentra el valor promedio y la energía de toda la señal.
ii) Divide la señal en segmentos de 20 ms de longitud y calcula la energía de cada segmento
y graficarlo.
5. Para los sistemas LTI descritos por las siguientes ecuaciones en diferencias, genera su impulso
respuesta y respuesta escalón unitario. Comente sobre las propiedades del sistema.

(i)y[n]-0.3695y[n-1]+0.1958y[n-2]=0.2066x[n]+0.4131x[n-1] +0.2066x[n-2]
(ii) y(n) = 0.77 y(n-1) + x(n) + x(n-1)
(iii) y(n) = 0.77 y(n-1) + 0.77 x(n) + x(n-1)

Análisis en el dominio de la frecuencia

Nota: Utilice freqz para todas las preguntas en esta sección sobre la Transformada de Fourier en Tiempo Discreto (DTFT).

No utilices fft

6. Para los sistemas LTI descritos por las siguientes ecuaciones en diferencias, genera su frecuencia
Comente sobre el tipo de respuesta.
y[n]-0.3695y[n-1]+0.1958y[n-2]=0.2066x[n]+0.4131x[n-1]+0.2066x[n-2]
Solución de muestra Respuesta en frecuencia
0

% ex6.m -20

% Respuesta en el dominio de la frecuencia de -40

ecuaciones de diferencias -60

% y[n]–0.3695y[n-1]+0.1958y[n-2] =
-80
0.2066x[n]+0.4131x[n-1]+0.2066x[n-2] 0 0.1 0.2 0.3 0.4 0.5 0.6
Frecuencia normalizada (Nyquist == 1)
0.7 0.8 0.9 1

0
b1 = [0.2066 0.4131 0.2066];
a1 = [1–0.3695 0.1958]; -50

freqz(b1,a1,64); -100

Respuesta de frecuencia -150

-200
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Frecuencia normalizada (Nyquist == 1)

7. Observe la relación entre la ubicación

de polos y ceros en el plano z, respuesta al impulso y respuesta de frecuencia para el
sistemas siguientes

iy(n) = 0.77y(n-1)+x(n)+x(n-1)
iiy(n) = 0.77y(n-1)+0.77x(n)-x(n-1)
iiiH(z) = 1-z -1/1+0.77z-1
ivH(z) = 1-z -1+z-2-z-3+z-4-z-5
vy(n) = x(n) + x(n-1) + x(n-2) + x(n-3) + x(n-4) + x(n-5)
viH(z)=3-3z -1

8. Sea x(n)=cos(0.2n ) + 0.5 cos(0.4n ); 0≤ n ≤ 99.

(i) Grafica el espectro de esta señal.
(ii) Generar una señal interpolada cero y(n)=x(n/2) y trazar su espectro. Comentar sobre
tus observaciones.
(iii) Genera la señal decimada d(n)=x(2n) y grafica el espectro. Comenta sobre tu
observaciones.
(iv) Genera la señal decimada d(n)=x(3n) y grafica el espectro. Comenta sobre tu resultado.
observaciones.
Transformada de Fourier Discreta
Nota: utiliza fft en esta sección

9. Calcula y traza la DFT de las siguientes secuencias y observa sus propiedades

a) Señal de impulso unitario; xi = {1 0 0 0 0 0 0 0}
b) Todos unos; x1 = {1 1 1 1 1 1 1 1}
c) Boxcar de tres puntos; xb = {1 1 1 0 0 0 0 0}
d) Caja simétrica; xbsy = {1 1 0 0 0 0 0 1}
e) Sinusoide; s[n] = A cos( 0n + Prueba con diferentes longitudes de secuencia.
f) Exponencial real; x[n] = (0.9)nu[n] 0 n 31. Compárelo con muestras de Y(ejw) =
(1-0.9 e-jw)-1 la DTFT de la exponencial infinitamente larga y explica la diferencia
en términos de ventana.

10. Compara las características de las siguientes funciones de ventana según su frecuencia.
gráficos de espectro. (usarwvtool)
a) Ventana rectangular w_rect[n] = 1 0 n M-1
b) Ventana de Hanning w_hann[n] = 0.5-0.5cos(2 n/M-1) 0 n M-1
c) Ventana de Hamming w_hamm[n] = 0.54-0.46cos(2 n/M-1) 0 n M-1
d) Ventana Blackmann w_blk[n] = 0.42-0.5cos(2 n/M-1)+0.08cos(4 n/M-1)
0 n M-1

11. Generar una señales[n] con tres componentes sinusoidales a 50, 120 y 240Hz corruptas
por AWGN. Traza el espectro e identifica los componentes de la señal.

Solución de muestra:
% identificación de sinusoides en ruido

fs=2000;
t = (0:199)/fs;
s = sen(2*pi*50.*t) + sen(2*pi*120.*t) + sen(2*pi*240.*t);

awgn = (0.5*randn(1,200)+.25); % N(0.25, 0.25)

sn = s + awgn;
subplot(211), plot(t,sn);
título(' Sinusoide con ruido'); cuadrícula;

Sn = fft(sn,200);
f = [Link];
sfmag = abs(Sn);
subplot(212), plot(f,sfmag(1:100));
título('Estimación espectral'); rejilla;

Repetir usando Simulink

12. Durante la grabación, una señal de ECG muestreada a 300Hz se contamina con un zumbido de 60Hz.
se proporcionan los latidos de la señal original y contaminada (600 muestras cada una) como ecgo
 y ECG. Carga estas señales y gráficas. Calcula la DFT de 600 puntos del
señal contaminada y anular los componentes DFT correspondientes a la señal de 60Hz.
Realice la IDFT para obtener la señal de ECG filtrada. Muestra las señales filtradas y originales.

13. a) Generar 600 muestras de la señal x(n) = cos(0.1n ) + cos(0.4n ) + cos(0.7n ). Trazar su
Magnitud de la DFT.

(b) Generar 200 muestras de cada una de las tres señales x1(n) = cos(0.1n ), x2(n)= cos(0.4n ),
y x3(n) = cos(0.7n ). Concaténalos {x1(n)2(n), x3(n)} para formar una muestra de 600
señal. Grafique su magnitud DFT.
(c) Obtener un gráfico de tiempo-frecuencia usando espectrograma

14. Deseamos calcular la FFT de radix 2 de las muestras de señal adquiridas de la señal x(t) =
Acos(2 F0t) + B Cos(2 (F0+ F)t) donde F0=100Hz a una frecuencia de muestreo Fs=480Hz.
(a) Si A=B=1, ¿cuál es el menor número de muestras de señal N?minrequerido para un
resolución de frecuencia de 2Hz si (i) ventana rectangular (ii) ventana de Hanning (iii)
Se utiliza la ventana Blackman. Traza la FFT en cada caso.
(b) Si A=1 y B=0.02, ¿cuál de las siguientes ventanas - Rectangular, Bartlett, Hanning,
Hamming, Blackman–¿puedes usar para obtener una resolución de frecuencia de? F=2Hz. ¿Cuál?
uno requerirá el número mínimo de muestras y por qué? Grafique la FFT en cada uno
caso.

15. Graficar (“gráfico de puntos”) tiempo vs longitud para la función fft para N=100 a 1000