LABO1

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 7

Universidad Nacional de Ingeniería

Escuela Profesional de Ingeniería Mecatrónica

PROCESAMIENTO DIGITAL DE SEÑALES (MT417)


LABORATORIO 1: RECUANTIZACIÓN DE UNA SEÑAL

NOMBRE CODIGO
20191005G
Perez Cieza Piero Antonio
20204042H
Barzola Chancasanampa Bryan Cesar
20204017C
Brioso Rodriguez Janddy Kennedy

Quillatupa Vasquez Carlos 20200017i

Chinchay Vargas Ronald 20203006H

I.- Objetivos:

• Implementar el proceso de recuantización de señales mediante la técnica del factor de escala fijo.
• Comprobar los cambios en SNR producidos por el proceso de recuantización de señales.
• Detección de eventos sobre el espectro de módulo de una señal.

II.- Materiales

• Computador personal con tarjeta de sonido, SO Win7 (o mayor) y MATLAB 2015 (o mayor).
• Audífonos y micrófono

1
III.- Experiencia a desarrollar:

Experiencia 1
1. Capturar una frase vía micrófono. La frase a pronunciar durante el tiempo de captura será:
“Primer laboratorio calificado 2023 de procesamiento digital de señales”. Utilice los
parámetros de la siguiente tabla, en un rango de captura de 90%.
Frecuencia de muestreo 15000
Tiempo de captura 4 segundos
Número de canales 1
Bits/muestra 8
Formato ‘double’

2. Grabar la señal a utilizar en un archivo como se indica a continuación


𝑎𝑢𝑑𝑖𝑜𝑤𝑟𝑖𝑡𝑒(′𝑠𝑒𝑛𝑎𝑙. 𝑤𝑎𝑣,′ 𝑥, 𝐹𝑠)
Donde: senal, es el nombre del archivo con que se almacena la señal en el disco duro; x, es
el vector donde está la señal a almacenar; y Fs, es la frecuencia de muestreo con que se
almacena.

3. ¿Cuál es el valor máximo y mínimo de la señal capturada?


Valor máximo=0.33594 Valor mínimo=-0.39844

4. ¿Cuántas muestras tiene la señal capturada? Utilice la función 𝑡𝑎𝑚 = 𝑙𝑒𝑛𝑔𝑡ℎ(𝑥)


60000

5. Graficar la señal capturada en el tiempo, utilizando el comando “plot”. El eje de abscisas


debe estar en segundos.

2
6. Indique el factor de escala
𝒇𝒆 = 1 (POR TEORÍA MATLAB USA ESTE VALOR)

7. Los procesos de recuantización se deben realizar a partir de la señal capturada, es decir, no


se debe volver a capturar la señal para realizar la recuantización. Esto es muy importante
tomar en cuenta.

8. Para el proceso de recuantización y cálculo de la SNR considerar lo siguiente, utilice un rp=3,


y siga el siguiente procedimiento.

𝑥 ∙ (2(𝑟𝑝−1) − 1)
𝑦𝑝 = 𝑟𝑜𝑢𝑛𝑑 ( )
𝑓𝑒

𝑦𝑝 ∙ 𝑓𝑒
𝑦=
(2(𝑟𝑝−1) − 1)
e= y−x

𝐸𝑒 = 𝑠𝑢𝑚(𝑒.∗ 𝑒) + 10−20

Ex = sum( x. * x )

 Ex 
SNR = 10 * log 10 
 Ee 

9. Cuáles son los valores de:


3
Ex=352.7468 Ee=269.0662

10. Obtenga el intervalo de cuantización en el vector reproducible y en el vector de enteros.


Indique la operación para obtener lo requerido.
%Vector reproducible(y)doubles
rc_y=2*max(y)/((2^rp)-2)=0.1111
%Vector enteros(yp)
rc_yp=1; %teoria

rc_y=0.1111
rc_yp=1

11. Obtenga el factor de carga aplicando la formula teórica

3.6181

12. Obtenga el valor de la SNR


1.1760

13. Recuantizar la señal a rp=5 utilizando el procedimiento anterior.

14. Obtenga el intervalo de cuantización en el vector reproducible y en el vector de enteros.


Indique la operación para obtener lo requerido.

%Vector reproducible(y)doubles
rc_y=2*max(y)/((2^rp)-2)
%Vector enteros(yp)
rc_yp=1; %teoria

rc_y=0.0222
rc_yp=1

15. Cuáles son los valores de:


Ex= 352.7468 Ee=17.7532

16. Obtenga el factor de carga aplicando la formula teórica

4
3.6181

17. Obtenga el valor de la SNR


12.9819

18. Recuantizar la señal a rp=7 utilizando el procedimiento anterior

19. Obtenga el intervalo de cuantización en el vector reproducible y en el vector de enteros.


Indique la operación para obtener lo requerido.
%Vector reproducible(y)doubles
rc_y=2*max(y)/((2^rp)-2)=0.0053
%Vector enteros(yp)
rc_yp=1; %teoria

rc_y=0.0053
rc_yp=1

20. Cuáles son los valores de:


Ex=352.7468 Ee=1.3062

21. Obtenga el factor de carga aplicando la formula teórica


3.6181

22. Obtenga el valor de la SNR


24.3145

23. Realice en MATLAB la gráfica SNR (eje ordenadas) vs. rp (eje abscisas) con los valores de
SNR obtenidos para cada valor de rp. Utilizar el comando stem para graficar y en cada valor
de rp colocar el valor de SNR con “data tips”.

Gráfica SNR vs. Rp


Rojo:n=3;
Azul:n=5;
Verde:n=7;

5
24. ¿De acuerdo a los factores de carga obtenidos en cada proceso de recuantización? Pregunta:
¿Son iguales?, ¿Deben ser iguales o no? Justifique.
Sí, son iguales. Esto debido a que el factor de carga solo depende de la señal de
entrada mas no de la señal recuantizada por un proceso de recuantización.

Adjuntar el código al final de la guía


close all;clear all;clc
%% Parámetros de captura
Fs=15000;%frecuencia de muestreo
r=8;%bits/muestra
nc=1;%número de canales
tc=4;%tiempo de captura
%% Captura
laboratorio1=audiorecorder(Fs,r,nc);
record(laboratorio1,tc);
pause(tc+0.1);
x=getaudiodata(laboratorio1,'double');
%[x,Fs] = audioread("F:\CURSO OFIMATICA\POWER POINT\audio.wav",;
%%
audiowrite("F:\CURSO OFIMATICA\POWER POINT\AUDIOLABO1.wav", x,
Fs);
T=1/Fs; %Periodo de muestreo
t=0:T:tc-T;
figure, plot(t,x), grid, xlabel("Tiempo (s)"), ylabel("Amplitud")
tam=length(x);
%% artificio cuantizacion
rp=3;
fe=1;
yp=round(x*(2^(rp-1)-1)/fe);
y=yp*fe/(2^(rp-1)-1);
e=y-x;
Ee=sum(e.^2); 6
Ex=sum(x.^2);
SNR= 10*log10(ex/Ee); % Se calcula la relación señal-ruido
fcarga=max(x)/(std(x))
%% HALLAR INTERVALO DE CUANTIZACION
%Vector reproducible(y)doubles
rc_y=2*max(y)/((2^rp)-2);
%Vector enteros(yp)
rc_yp=1; %teoria
%%
t = [3];
x = [1.176];
hold on
stem(t, x,'r'),grid, xlabel("SNR"), ylabel("rp")
t = [5];
x = [12.9819];
hold on
stem(t, x,'b')
t = [7];
x = [24.3145];
hold on
stem(t, x,'g')

También podría gustarte

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy