EVALUACIÓN MODULO III SyS 2022-II
EVALUACIÓN MODULO III SyS 2022-II
EVALUACIÓN MODULO III SyS 2022-II
SANCHEZ CARRION
% Señales originales
a = [2 5 7 9 11 15];
f = [30 40 50 60 60 60];
% Frecuencia de muestreo
fs = 300;
% Vector de tiempo
t = 0:T:1-T;
% Señal muestreada
s = zeros(1,length(t));
for i=1:length(a)
s = s + a(i)*sin(2*pi*f(i)*t);
end
% Gráficas
figure
subplot(2,1,1)
plot(t,s)
title('Señal muestreada en el dominio del tiempo')
xlabel('Tiempo (s)')
ylabel('Amplitud')
subplot(2,1,2)
stem(s)
title('Señal muestreada en el dominio de la secuencia')
xlabel('Número de muestra')
ylabel('Amplitud')
b) Si dispone de una señal senoidal de amplitud 10 y una frecuencia de
200 Hz y una frecuencia de muestreo de 4000 Hz.
Esto generará una figura con dos gráficas, una para la señal original y otra
para la nueva señal muestreada a 8000 Hz.
b-1) (4 puntos) De la anterior señal (amplitud 10 y una frecuencia de 200
Hz) se desea obtener una nueva frecuencia de muestreo de 1000 Hz.
Dibuje el diagrama de bloque de este nuevo sistema de compresión con
sus respectivos valores. Grafique y comente resultados.
% Filtro antialiasing
lpf_order = 20; % Orden del filtro pasa-bajos
lpf_cutoff = 0.4*Fs/2; % Frecuencia de corte del filtro en Hz
lpf_coeffs = fir1(lpf_order, lpf_cutoff/(Fs/2)); % Coeficientes del filtro
FIR
% Decimación
dec_factor = 4; % Factor de decimación
y = decimate(filter(lpf_coeffs, 1, x), dec_factor); % Señal decimada
b-2) (4 puntos) De la señal original (amplitud 10 y una frecuencia de 200
Hz) se desea obtener una nueva frecuencia de muestreo de 8000 Hz.
Dibuje el diagrama de bloques de este nuevo sistema de expansión con
sus respectivos valores. Grafique y comente resultados.
El resultado debería ser una gráfica de la señal original y la señal muestreada a la nueva
frecuencia de muestreo. Podemos observar que la señal muestreada a 8000 Hz tiene la
misma forma que la señal original, pero con una resolución temporal cuatro veces mayor
b-3) (4 puntos) De la señal original (amplitud 10 y una frecuencia de 300
Hz) se desea obtener una nueva frecuencia de muestreo de 6000 Hz.
Dibuje el diagrama de bloques de este nuevo sistema con sus respectivos
valores. Comente resultados.
% Interpolación
factor_interp = 20;
senal_interp = resample(senal_original, factor_interp, 1);
% Decimación
factor_decim = 5;
senal_final = downsample(senal_interp, factor_decim);
% Graficar señales
figure
subplot(3,1,1)
plot(tiempo, senal_original)
title('Señal original')
xlabel('Tiempo (s)')
ylabel('Amplitud')
tiempo_interp = 0:1/(factor_interp*6000):1-1/(factor_interp*6000);
subplot(3,1,2)
plot(tiempo_interp, senal_interp)
title('Señal interpolada')
xlabel('Tiempo (s)')
ylabel('Amplitud')
tiempo_final = 0:1/(factor_interp*factor_decim*6000):1-
1/(factor_interp*factor_decim*6000);
subplot(3,1,3)
plot(tiempo_final, senal_final)
title('Señal digitalizada')
xlabel('Tiempo (s)')
ylabel('Amplitud')
Podemos observar que la señal interpolada tiene una frecuencia de 6000
Hz y que la señal digitalizada tiene la frecuencia de muestreo deseada de
6000 Hz. La señal digitalizada sigue la forma de onda de la señal original,
aunque es importante mencionar que la frecuencia de muestreo debe ser
elegida cuidadosamente para evitar el efecto de aliasing en la señal
digitalizada.
% Graficar señales
figure
plot(tiempo, senal_original, 'b')
hold on
plot(tiempo, senal_cuantificada, 'r')
legend('Señal original', 'Señal cuantificada')
xlabel('Tiempo (s)')
ylabel('Amplitud')
NOTA
Los programas de MatLab que desarrolle, describa, especifique la función de
cada línea.