Universidad Nacional Mayor de San Marcos

Facultad de Ingeniería Electrónica y Eléctrica

Escuela Profesional de Ingeniería de Telecomunicaciones
“Año del Bicentenario del Perú: 200 años de Independencia”

TRABAJO FINAL
“PREGUNTA N°01”

• Alumno: Cornelio Cosisilla, Luis Angel

• Código: 18190327

• Profesor: Flavio Carrillo

• Horario: Jueves 2:00-4:00 p.m.

• Curso: Laboratorio de Procesamiento

Digital de Señales en TR

2021
 UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
FACULTAD DE INGENIERIA ELECTRICA Y ELECTRONICA

I. ACTIVIDAD

La figura nos muestra un modelo de convertir la frecuencia de muestreo de la

secuencia x[n] por un factor que está relacionado con los valores de M y D. La
finalidad del filtro pasa bajo H(z) es para filtrar las componentes no deseadas de
proceso del up sampling
Considere la señal x[n] contenida en el archivo texto anexo, la misma que es
aplicada al sistema discreto mostrada en la figura. La secuencia x[n] ha sido
obtenida con una frecuencia de muestreo Fm = 44100 Hz. El objetivo de este
sistema discreto es que la salida y[n] se reproduzca con una frecuencia de muestreo
de 29400 Hz. Analiza la secuencia de entrada y de acuerdo con la información que
obtengas, diseña los valores D, M y H(z) tal que la frecuencia de muestreo de la
secuencia de salida al nuevo valor deseado.
(a) Determina y prueba la expresión de H(z) que corresponde a un filtro digital pasa
bajo, desarrollando parta ello dos modelos, uno de tipo FIR y otro de tipo IIR, donde
el rizado de la banda de paso δp = 0.01 dB, y el rizado de la banda de atenuación δp
= -50 dB para ambos modelos. Las frecuencias de corte de la banda de paso y de
atenuación lo deberás determinar de acuerdo con el análisis del espectro y el
objetivo del sistema discreto.
(b) Graficar en el tiempo x[n], r[n], v[n] e y[n] así como en frecuencia el módulo y
fase de la respuesta en frecuencia de X(ejΩ), R(ejΩ), V(ejΩ) y Y(ejΩ). La escala de
frecuencia deberá estar expresado en Hz y el tiempo en segundos.
(c) Entregar en Octave el código de todos los procesos de cálculo que corrobore los
resultados y especificaciones propuesto.
(d) ¿Qué diferencias y similitudes observas entre las señales x[n] y v[n] con respecto
al tiempo, así como en frecuencia?
 UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
FACULTAD DE INGENIERIA ELECTRICA Y ELECTRONICA

II. DESARROLLO
Para comenzar con este ejercicio lo primero que debemos hacer es cargar el archivo
proporcionado por el profesor a Matlab para poder darle el respectivo
procesamiento, para ello mediante el comando audioread obtenemos el vector de
audio de el archivo Halcon_0602.wav y lo almacenamos en un vector x, asimismo
obtenemos con este mismo comando la frecuencia de muestreo Fm, que es igual a
44.1 kHz, ahora bien el problema del ejercicio nos pide que la señal de salida tenga
una frecuencia de reproducción (muestreo) de 29.4 kHz aplicando técnicas de
upsampling y downsampling. Como bien sabemos la técnica del upsampling
depende de un factor N que agrega N-1 ceros entre muestras y se utiliza
principalmente para aumentar N veces la frecuencia de muestreo, por otro lado, el
downsampling de factor M quita M-1 muestras entre dos muestras y se utiliza para
reducir M veces la frecuencia de muestreo, con esto en mente hallamos la razón
entre ambos valores de frecuencia para saber en qué factor k debemos modificar la
muestra, entonces tenemos:
29.4 𝑘𝐻𝑧 2
𝑘= =
44.1 𝐾𝐻𝑧 3
Como podemos ver la frecuencia original varia en 2/3 respecto a la frecuencia de
salida, esto quiere decir que si aumentamos la frecuencia original de 44.1 kHz por
un factor de 2 y luego la reducimos en un factor de 3 obtendremos la frecuencia
deseada de 29. kHz, entonces podemos concluir que en el bloque M se hará un
upsampling de factor 2 mientras que el bloque D un downsampling de factor 3.
Ahora bien, también se nos pide diseñar 2 filtros, pero primero avanzaremos bloque
por bloque ya que para diseñar el filtro debemos conocer el espectro de la señal de
salida del upsampling, es por ello por lo que tenemos la primera parte del código
donde se realizara el ploteo en el tiempo y frecuencia de la señal de entrada y la que
sale luego de realizar el upsampling de factor 2.
 UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
FACULTAD DE INGENIERIA ELECTRICA Y ELECTRONICA

Fig. 01 Código para graficar la señal original y con upsampling de factor 2 en el tiempo y frecuencia
 UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
FACULTAD DE INGENIERIA ELECTRICA Y ELECTRONICA

Fig. 02 Señal original y con upsampling de factor 2 en el tiempo

Fig. 03 Señal original y con upsampling de factor 2 en la frecuencia

Podemos apreciar que en cuanto a la figura 2 en donde se muestran ambas señales

en el dominio del tiempo su duración en segundos se duplica, cabe mencionar que
ambas se han muestreado a la misma frecuencia de muestreo que es 44.1 kHz, el
hecho de que suceda esto es que como ya se mencionó, la técnica del upsampling lo
que hace es aumentar la frecuencia de muestreo para poder recuperar la señal
original es por ello que si queremos ver graficas iguales deberíamos graficar la
 UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
FACULTAD DE INGENIERIA ELECTRICA Y ELECTRONICA

segunda señal con una frecuencia de muestreo de 88.2 kHz, por otro lado en la
figura 3 vemos que en el espectro ha aparecido la componente espejo de la señal,
esto se debe a que al realizar el upsampling la frecuencia pasa a reducirse en un
factor M a la frecuencia original y como la componente espejo ahora tiene menos
frecuencia que la frecuencia de muestreo puede aparecer en el espectro debido a
dicha compresión, ahora lo que debemos realizar es un filtrado que nos permita
seleccionar solo la primera componente que es la perteneciente a la señal original,
entonces establecemos que el filtro será pasa bajo con una frecuencia corte 6 kHz.
III. DISEÑO DEL FILTRO
Según lo planteado por el ejercicio se nos pide diseñar dos filtros uno FIR y otro IIR
ambos a una misma frecuencia de corte, con rizado en la banda de paso de δp = 0.01
dB y rizado de la banda de atenuación δp = 50 dB. Tomando esto en cuenta se
utilizo lo visto en clase para diseñar el filtro FIR usando la herramienta fdatool, para
ello solo escribimos en la ventana de comandos fdatool y se nos abrirá una ventana
igual a la figura 4, luego en ella seleccionamos el tipo de filtro que se desea y el
método, en este caso pasa bajo(Lowpass) y FIR, asimismo llenamos los otros
parámetros para la frecuencia de muestreo en 44.1 kHz, frecuencia de paso de 6
kHz, de atenuación de 6.5 kHz y finalmente los rizados para la banda de paso y
atenuación, finalmente pulsamos el botón Design Filter y el filtro se diseñará
automáticamente, ahora lo que debemos hacer para usarlo en nuestro ejercicio es
exportar los valores de los coeficientes del filtro desde el mismo lugar

Fig. 04 Herramienta fdatool para generar coeficientes de filtro FIR

 UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
FACULTAD DE INGENIERIA ELECTRICA Y ELECTRONICA

Por otro lado, para realizar el diseño del filtro IIR se utilizo un filtro elíptico con las
mismas características que el filtro anterior de orden 6, cabe resaltar para esta parte
que los coeficientes obtenidos para el filtro FIR fueron almacenados en un
documento de texto para que no se pierdan al cerrar la aplicación, luego de esto se
calculo la respuesta en frecuencia para ambos filtros con 100 puntos a una
frecuencia de muestreo de 44.1 kHz

Fig. 05 Código para el diseño del filtro IIR y respuesta en frecuencia

Fig. 06 Respuesta en frecuencia para ambos filtros

 UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
FACULTAD DE INGENIERIA ELECTRICA Y ELECTRONICA

IV. FILTRADO DE LA SEÑAL

Luego de diseñado ambos filtros se procede a realizar el filtrado de la señal usando
el comando filter, luego se procede a graficarlos en el tiempo y frecuencia usando el
código

Fig. 07 Código para el filtrado y grafica en el tiempo y frecuencia usando ambos filtros
 UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
FACULTAD DE INGENIERIA ELECTRICA Y ELECTRONICA

Fig. 08 Señal original y con filtrada en el dominio del tiempo

Fig. 09 Señal original y con filtrada en el dominio de la frecuencia

En cuanto a los resultados obtenidos podemos ver que en el tiempo no hay ninguna
variación, la señal sigue teniendo una duración del doble de la original lo que
significa que sigue manteniendo su frecuencia de muestreo al doble, es decir, a 88.2
kHz por otro lado el espectro en frecuencia si ha cambiado, ya no observamos la
frecuencia espejo de hace un momento en ambas salidas por lo que se puede notar
que en ambos casos el filtrado ha sido el correcto.
 UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
FACULTAD DE INGENIERIA ELECTRICA Y ELECTRONICA

V. DOWNSAMPLING
Y finalmente luego del filtrado procedemos a realizar el downsampling de factor 3 y
realizamos las graficas respectivas para representarlas en el tiempo y frecuencia,
para realizar el downsampling usaremos el comando downsample y se realizará este
proceso para la salida de cada filtro (FIR e IIR)

Fig. 10 Código para graficar la señal original y con downsampling de factor 3 en el tiempo y frecuencia
 UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
FACULTAD DE INGENIERIA ELECTRICA Y ELECTRONICA

Fig. 11 Señales filtradas con downsampling de factor 3 en el tiempo

Fig. 12 Señales filtradas con downsampling de factor 3 en la frecuencia

Lo que podemos observar en este último proceso en la figura 11 es que la duración

ha vuelto a cambiar, esto se debe nuevamente a que la frecuencia de muestreo se ha
vuelto a modificar en este caso ya obtuvimos lo deseado que es una frecuencia de
reproducción de 29.4 kHz lo cual se comprobará mas adelante, por otro lado en el
espectro en frecuencia de ambas señales con upsampling se observa que se ha
modificado nuevamente su frecuencia pero esta vez se ha expandido por lo que
 UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
FACULTAD DE INGENIERIA ELECTRICA Y ELECTRONICA

supone que la frecuencia de muestreo se ha reducido desde la anterior vez, asimismo

se observa el mismo espectro que la señal original solo que comprimida un factor
k=2/3 veces.
VI. PRUEBAS DE REPRODUCCIÓN
Finalmente, en el último paso pasamos a reproducir las señales obtenidas en la
salida del sistema para ello usaremos el comando sound que nos permite reproducir
datos almacenados en vectores o matrices, para usarlo debemos colocar en sus
parámetros primero el vector a reproducir y seguidamente la frecuencia de
reproducción, además también se adicionó pause que nos permite realizar una pausa
entre comandos para que puedan escucharse por separado, las pruebas de sonido
podrán ser comprobadas ejecutando el código.

Fig. 12 Prueba de reproducción