Facultad de Ingeniería

Ingeniería Electrónica - Ingeniería de Telecomunicaciones y Redes

Ingeniería Mecatrónica

2022-II

PROCESAMIENTO DIGITAL DE SEÑALES (EL222)

TRABAJO PARCIAL (TP)

Duración Total: 3 horas

Lectura: 10 minutos
Desarrollo: 2 horas 40 minutos
Envío: 10 minutos

Tema:
Señales y Sistemas Discretos — Respuesta Impulsiva Acústica

APELLIDOS NOMBRES CARRERA SECCIÓN GRUPO

EL76 1
EL76 1
EL76 1

Fecha: Martes 4 de octubre de 2022

Hora: 8:00 am – 11:00 am
Profesor: César D. Salvador

Objetivos
 Modelar la transmisión del sonido en un espacio acústico como un sistema lineal
 Calcular la respuesta impulsiva del espacio acústico a partir de mediciones
 Sintetizar la audición del espacio aplicando la respuesta impulsiva a una señal anecoica

Materiales
 Computador personal con tarjeta de sonido
 Sistema operativo Windows 10 (o superior)
 Matlab 2019 (o superior) con las siguientes toolboxes:
o Signal Processing Toolbox
o DSP System Toolbox
o Audio Toolbox

Instrucciones
Descargar los siguientes archivos:
 Plantilla de Word para el reporte: [Link] (este documento)
 Plantilla de script de Matlab: EL222PDSEL76G1TP.m
  Archivos de audio
o Señal del emisor: [Link]
o Señal del receptor: [Link]
o Señal de audio sin eco: [Link]

Seguir las instrucciones disponibles a lo largo de este documento.

Enviar como respuesta los siguientes archivos:

 Reporte en PDF: EL222PDSEL76G1TP_ApellidosDeQuienEnví[Link]
 Script de Matlab: EL222PDSEL76G1TP_ApellidosDeQuienEnvía.m
 Archivo de audio procesado con eco: AudioEco_ApellidosDeQuienEnví[Link]

El script de Matlab deberá ejecutarse sin errores y mostrar todos los resultados al presionar la tecla F5.
Las figuras deben mostrar valores, unidades, leyendas e indicaciones con tamaños de fuente legibles.
No serán considerados scripts que generen mensajes de error al ser ejecutados.

Cuando se requiera insertar imágenes en el reporte, estas deberán haber sido exportadas previamente
desde Matlab en formato PNG utilizando el comando PRINT. No serán calificadas las imágenes
insertadas mediante capturas de pantalla.

El audio procesado deberá ser creado con la misma frecuencia de muestreo asignada en los archivos
de audio originales. El audio procesado deberá estar correctamente escalado de manera que no hallan
saturaciones al reproducirlo. Archivos de audio defectuosos no serán considerados.

Bastará que un solo alumno de cada grupo suba los tres archivos. Tanto el reporte como el script
deberán tener escritos los datos completos de cada integrante del grupo.

Calificación
Todos los archivos enviados deben cumplir con las especificaciones mencionadas para acceder a la
siguiente calificación.

Experiencia Pregunta Puntos

1 4
1 2 7
3 5
2 1 4
Total 20
 Señales y Sistemas Discretos — Respuesta Impulsiva Acústica

Introducción
La respuesta impulsiva (RI) o función de transferencia (FT) de un espacio acústico caracteriza la
transmisión de sonido desde una posición en la cual se ubica un emisor hacia otra posición en la cual
se ubica un receptor. El proceso de transmisión acústica se modela mediante el siguiente sistema
lineal.

Emisor RI/FT Receptor

Una vez medida la transmisión acústica, la RI o FT se puede aplicar a cualquier otra señal de audio,
idealmente medida en condiciones anecoicas o sin eco, para sintetizar una nueva señal que incluya la
reverberación o eco del auditorio. Esta nueva señal de audio sintetizada se oirá como si hubiera sido
grabada en el auditorio.

Señal de audio anecoica RI/FT Señal de audio con eco

Experiencia 1: Cálculo de la Respuesta Impulsiva de un Espacio Acústico

Pregunta 1
En Matlab, utilizar la función audioread para leer el archivo de audio [Link]. Ingresar los valores
correspondientes en la siguiente tabla.

Señal del Emisor

Frecuencia de muestreo (Hz)
Número de muestras
Duración (s)
Número de canales

Utilizar las funciones plot y print para generar la siguiente gráfica.

Gráfica de la Señal del Emisor en el Dominio del Tiempo (s)

Utilizar nuevamente la función audioread para leer el archivo de audio [Link]. Ingresar los
valores correspondientes en la tabla a continuación.

Señal del Receptor

Frecuencia de muestreo (Hz)
Número de muestras
Duración (s)
Número de canales

Utilizar las funciones plot y print para generar la siguiente gráfica.

Gráfica de la señal del Receptor en el Dominio del Tiempo (s)

Escribir comentarios y observaciones acerca de las señales del emisor y del receptor. Describir las
señales indicando sus características y beneficios.

Observaciones sobre las Señales del Emisor y del Receptor

- El emisor es la señal emitida por algún altavoz de algún espacio repetido un número de
veces y la señal receptos es la señal capturada por algún micrófono dentro del espacio
en la que se encuentra el emisor `pero desde otro punto del espacio.
- Ambas señales tienen la misma frecuencia de muestreo.
- Comprar la duración de las señales.
- Ambas señales tienen canal uno porque al ejecutar la función audiowrite, el segundo
vector indica valor 1.

Pregunta 2
Utilizar las funciones xcorr y findpeaks para extraer el promedio de las tres mediciones contenidas en
la señal registrada por el receptor. En el dominio de la correlación, indicar los instantes en los cuales
ocurren las máximas correlaciones. Ingresar los valores en la siguiente tabla.

Picos hallados en la correlación entre las señales del emisor y del receptor
Pico Índice de muestra (lags) Instante de tiempo (s)
1
2
3

Utilizar la función tukeywin para aplicar una ventana en el dominio del tiempo a la señal obtenida al
promediar los tres registros del receptor.

Usando las funciones plot y print, generar la siguiente gráfica.

Gráfica de la Señal Receptor Promedio con Ventana en el Dominio del Tiempo (s)

Utilizar la función fft para calcular las transformaciones de Fourier de las señales emisor y receptor
promedio con ventana.

Usando las funciones plot y print, generar las gráficas a continuación con ejes horizontal y vertical
en escala logarítmica.

Gráfica de la Magnitud de la Transformación de Fourier de la Señal del Emisor, en el

Dominio de la Frecuencia (kHz)

Gráfica de la Magnitud de la Transformación de Fourier de la Señal Receptor Promedio

con Ventana, en el dominio de la Frecuencia (en kHz)

Escribir comentarios y observaciones acerca del propósito de calcular el promedio de las mediciones y
aplicar una ventana. Comentar también cada gráfica obtenida.
 Observaciones sobre el Procedimiento y sobre las Señales del Emisor y del Receptor

- El profesor dijo que posiblemente aquí no use la transformada de Foutier, sino un

análisis de espectrogramas.
- No importa si usan la función stft o espectogram. Si usan la función spectogram,
utilícenlo sobre la variable de receptor promedio.
- El receptorpromedio es la señal de correlación, en la que nos indica mediante la
función xcorr la señal repetida del emisor dentro de la señal receptor.
- Con respecto a la función STFT, pueden usar la misma entrada de la señal, pero
pueden agregarle la frecuencia de muestreo.
- Cuando el profesor dice que se tiene que usar el tamaño de ventana adecuado, se
refiere a que se debe usar el tiempo de la variable receptor promedio.
- Hay que recordar que, si el color de un spectograma cambia, significa que la
amplitud de la señal esta cambiando en el mismo orden. Si la señal cambia en el
ejey del espectrograma, significa que la frecuencia de la señal esta aumentando.
Es necesario recalcar que dependiendo del cambio de color de la señal del
espectorama y su cambios en el eje y, lo deben comprar con la grafica del recepto
promedio en el domingo del tiempo. Si sus observaciones coinciden, estarán
bien.
- Duracion = Numero de muestras/Fs tz = (0:1/Fs:Duracion)
- Pueden investigar más sobre el xcorr, pero la interpretación mas plausible es que
comprar la señal de entrada con respecto a la salida para ver cuantas veces e ha
repetido. Ahora, si quieren ser más técnicos, si encuentran alguna definición, aplíquenle
a la señal que están analizando. Ejemplo: se denomina covarianza, lo podrían usar
como “La función xcorr compara ambas señales mediante la correlación cruzada, la
cual consiste en obtener las covarianzas de ambas señales, para obtener la similitud.
La similitud, se le calculas sus picos. La función findspeak, lo que hace es calcular las
amplitudes máximas dentro de la ventana aplicada a la duración del emisor”
- Especifiquen a que se refieren cuando pongan “de la señal” o “de la gráfica”,
(emisor/espectograma)
- Si ven que en el espectograma hay un comportamiento lineal hacia arriba, sigan que es
un shirp.
- Deben poner por que eligen el valor de NFFT.

Pregunta 3
Utilizando el método de la división en frecuencia, calcular la función de transferencia del espacio
acústico.

Usando las funciones plot y print, gráficar la función de transferencia acústica con ambos ejes en
escala logarítmica.

Gráfica de la Magnitud de la Función de Transferencia Acústica en el Dominio de la

Frecuencia (kHz)

Usando la función ifft, calcular la respuesta impulsiva acústica. Centrar adecuadamente la respuesta
impulsiva utilizando la función circshift.

Usando las funciones plot y print, generar la siguiente gráfica.

Gráfica de la Respuesta Impulsiva Acústica en el Dominio del Tiempo (s)

Escribir comentarios y observaciones acerca de los resultados.

Observaciones sobre la Función de Transferencia y la Respuesta Impulsiva

- Pueden ver el video para ver el profesor, que va a pedir. Sin embargo, a partir de la
respuesta al impulso pueden saber si la señal es, causal, estable, invariante en el
tiempo. Entonces lo que van a hacer es decir si es causal, si es estable o si es
invariante en el tiempo, además de las observaciones del profesor.

Experiencia 2: Síntesis de Audio con Reverberación

Pregunta 1
En Matlab, utilizar la función audioread para leer el archivo [Link]. Ingresar los valores
correspondientes en la siguiente tabla.

Señal de Audio Anecoica (sin Eco)

Frecuencia de muestreo (Hz)
Número de muestras
Duración (s)
Número de canales

Utilizando la función fftfilt, hacer la convolución de la repuesta impulsiva acústica con la señal de
audio anecoico para sintetizar una nueva señal de audio con las características acústicas del espacio
en el cual fue realizada la medición.

Finalmente, usar la función audiowrite para guardar la señal de audio con eco en un archivo WAV.

Escribir comentarios y observaciones acerca de los resultados.

Observaciones sobre la Señal de Audio Sintetizada