Medición TR Mejía Rojas

Informe técnico medición de tiempo de reverberación
Mejía Soto Sebastián1, Rojas Sarria Camilo Andrés1.

Ingeniería de Sonido
Facultad de Ingenierías
Universidad de San Buenaventura
Medellín
Contextualización
Se realiza la medición para estimar el tiempo de reverberación del aula 312B, tercer piso de la
Universidad De San Buenaventura, de la ciudad de Medellín, sede San Benito, la medición se
realiza a las 01:00 P.M. el 21 de agosto de 2019. El proceso utilizado para el análisis de
resultados es la convolución de un audio sweep y el sweep inverso por medio del software
Matlab.
Objetivo General
Aplicar las directrices de la norma ISO 3382 -1/-2 para realizar una medición de tiempo de
reverberación para un recinto cerrado y comparar los resultados contra los estándares para un
aula de clase.
Objetivos específicos
Medir el tiempo de reverberación del aula 312B
Comparar el valor según estándares para aulas de clase con el valor obtenido en la medición
Conocer de forma práctica el sonómetro y la instrumentación requerida para su correcto uso

en las mediciones.
Utilizar softwares para calcular la respuesta al impulso de una sala mediante un barrido de
frecuencias.
Descripción de la edificación
El lugar elegido para llevar a cabo la medición fue el aula 312B de la Universidad San
Buenaventura sede San Benito, el recinto cuenta con un espacio abierto sobre la puerta,
televisor, tablero, dos ventanas y una puerta.
Documentación de referencia
Norma ISO 3382-1
Norma ISO 3382-2
Manual de CESVA SC310
Equipos empleados
● Sonómetro CESVA SC310

Sonómetro de clase 1, con ponderaciones frecuenciales A, B, C y Z, ponderaciones
temporales Fast,slow e impulse, cuenta con un micrófono ACO 7052E que posee una
sensibilidad de 35 mV/Pa, pre polarizado con un diámetro de 1/2”; Además cuenta
con los modo de medición SPL, Leq, SEL, Lden, Ltm3, Ltm5, LMax, LMin, LPeak
● Pistófono
● Flexómetro
● Trípode de aproximadamente 1.30 mts de altura al estar desplegado
● Micrófonos de medición dbx RTA-M (x3)
Es un micrófono de análisis de audio en tiempo real, tiene una respuesta en frecuencia
bastante plana, tiene un patrón polar omnidireccional, es de tipo condensador, por lo
tanto debe ser alimentado con phantom power (48V).
● Interfaces de audio M-Audio MobilePre USB (x2)
Es un preamplificador con una interfaz de audio integrada diseñada para grabaciones
y sesiones de muestreo móviles. Trabaja a 16 bits y con una frecuencia de muestreo
de 48 kHz.
● Cabina JBL PRX612M
Sistema de dos vías, capaz de generar hasta 134 dB SPL de pico, con un rango de
frecuencias de aproximadamente 60 Hz - 17.5 kHz.
● Software MATLAB
La calibración del sonómetro se verificó antes de la medición verificando con el pistófono

que la sensibilidad estuviera a 94 dB a 1kHz
Condiciones meteorológicas
Temperatura 26° C
Sensación térmica 29°C

Humedad 60%
Punto de Rocío 16°C
Índice UV Extremo
Visibilidad 9.2 Km
Velocidad del viento 10 Km/h
Presión 1019 hPa
Tabla 1: Condiciones meteorológicas primer día de medición.
Desarrollo
Tráfico denso por la hora (12:00 – 14:00). Se realiza el montaje del escenario para medir,
dejando el espacio totalmente vacío, para iniciar se realiza la medición de ruido de fondo
obteniendo un ruido de fondo promedio de 45 dB, esto mediante el sonómetro con un tiempo
de integración de 5s, se realiza la medición en 3 puntos.
Figura 1: Ruido de fondo punto A

Figura 2: Ruido de fondo punto B
Figura 3: Ruido de fondo punto
Debido las características y la ubicación del aula los resultados que se determinaron para el
nivel de ruido de fondo fueron considerablemente altos (tráfico denso), la fuente no permitió
subir el nivel de presión sonora 45 dB arriba del ruido de fondo con un correcto
funcionamiento, esto influyó en los resultados ya que mediante el software Matlab notamos
que los resultados fueren bajos para este recinto.
La medición se llevó a cabo mediante el método de respuesta impulsiva integrada que sugiere
la norma, que nos permite mediante un sine sweep de todo el rango de frecuencias requerido
(10 Hz a 20 kHz) representarlo en tiempo como si fuera un impulso.
Siguiendo la indicación de que se tomaran al menos tres puntos de medición por fuente con
los micrófonos dbx y mínimo dos posiciones de fuente realizamos en total 6 mediciones que
nos permitieron encontrar la respuesta al impulso de la sala, los puntos fueron determinados
aleatoriamente intentando simular la ubicación de los espectadores en el aula y a más de un
metro de cualquier superficie reflectante. Sin embargo, dentro de la medición surgió un
inconveniente pues uno de los audios no registro de manera correcta por lo que solo se
utilizaron 5 audios para el análisis de los resultados.
Figura 4: Resultados micrófono 1 fuente 1 en 1kHz
Figura 5: Resultados micrófono 2 fuente 1 1 kHz
Segundo punto de fuente
Figura 6: Resultados micrófono 1 fuente 2 1kHz

Frecuencia [Hz] 100 125 160 200 250 315 400 500 630 800 1000 1250 1600 2000 2500 3150 4000
Rt [s] F1P1 0.3 0.2 0.1 0.1 0.4 0.1 0.5 0.5 0.3 0.4 0.3 0.3 0.3 0.3 0.4 0.4 0.4
Rt [s] F1P2 0.3 0.2 0.1 0.1 0.4 0.1 0.5 0.5 0.3 0.4 0.3 0.3 0.3 0.3 0.4 0.4 0.4
Rt [s] F1P3 0.2 0.3 0.1 0.1 0.1 0.3 0.4 0.1 0.3 0.2 0.4 0.4 0.3 0.5 0.3 0.3 0.3
Rt [s] F2P1 0.2 0.2 0.2 0.4 0.3 0.2 0.4 0.6 0.4 0.4 0.6 0.3 0.3 0.4 0.3 0.4 0.3
Rt [s] F2P2 0.1 0.2 0.2 0.2 0.6 0.6 0.4 0.3 0.3 0.2 0.3 0.4 0.3 0.3 0.3 0.4 0.3
Promedio 0.22 0.22 0.14 0.18 0.36 0.26 0.44 0.4 0.32 0.32 0.38 0.34 0.3 0.36 0.34 0.38 0.34
Tabla 2: Resultados de la medición y rt promedio por bandas de octava

Para encontrar el tiempo de reverberación para las frecuencias medias utilizamos la siguiente
ecuación
Reemplazamos los valores obtenidos en la medición
(0.4 + 0.4)
𝑇𝑚𝑖𝑑 =
2
Y obtenemos el tiempo de reverberación para las frecuencias medias
Comparación del recinto
Según estándares el tiempo de reverberación medio de aulas de clase se debe encontrar entre
0,7 y 1,0 segundos, para el caso de los valores medidos los resultados obtenidos fue 0.4 un
caso de error puede deberse al correcto funcionamiento de la fuente y al nivel del ruido de
fondo.
Conclusiones
Como se evidencia en los datos obtenidos el tiempo de reverberación del recinto es muy bajo,
en principio sería un tiempo ideal ya que permite una escucha clara, sin las reflexiones tardías
que ocasionan difusión en la comunicación pero entrando al detalle histórico del aula se
evidencia un error dentro de la medición, en primera instancia un audio de la medición no
logro ser reproducido, esto por fallas del sistema y en segunda instancia el código utilizado
para realizar la convolución aún carece de precisión, sin embargo el procedimiento de medida
se realizó de forma correcta.
En los resultados vemos que existieron varios problemas en las frecuencias bajas, pero nos
apoyamos en que la norma contempla que para métodos ingenieriles el rango de frecuencia a
evaluar va desde los 125 Hz hasta los 4 KHz.
ANEXOS
Código.
clc; clear; close all;
fs = 44100;
%Loading the audio files

audio = audioread('F1P1.wav');
audio2 = audioread('F1P2.wav');
inv_sweep = audioread('inv_sweep_20-20kHz_10sec_44kHz.wav');
% Se debe cambiar la variable audio1 a audio2, después audio3 y así para cada
posición de fuente y mic
IR = conv(inv_sweep,audio#);
[aa,bb] = max(abs(IR)); % Se haya el máx de IR en abs para organizar los
límites de la grafica
IR = IR(bb-1000:bb+(0.2*fs)-1001);
IR_F1P1 = IR;
% filter IRs into one-third octave bands
IR_f = zeros(length(IR),10);
BandsPerOctave = 1;
N = 6; % Filter Order
F0 = 1000; % Center Frequency (Hz)
f = fdesign.octave(BandsPerOctave,'Class 1','N,F0',N,F0,fs);
f.BandsPerOctave = 3;
f.FilterOrder = 8;
F0 = validfrequencies(f);
Nfc = length(F0);
for n_idx = 1:Nfc,
f.F0 = F0(n_idx);
Hd(n_idx) = design(f,'butter');
IR_f(:,n_idx) = filter(Hd(n_idx),IR);
end
% Impulse integral
E = zeros(size(IR_f));
E_f = zeros(size(IR_f));
for n_idx = 1:Nfc

E(end:-1:1,n_idx) = cumsum(IR_f(end:-1:1,n_idx).^2); % Schroeder
integral:
E_f(:,n_idx) = 10.*log10(E(:,n_idx)./max(E(:,n_idx))); % energy decay
curve per band
end
t1 = (1:length(IR_f))'/fs;
T_60 = zeros(Nfc,1);
for n_idx = 1:Nfc-1;
i0_10 = find(E_f(:,n_idx)<-5 & E_f(:,n_idx)>-35); % Encuentra el rango
entre -5 a -35 dB medido en T30 o -5 a -25 en T20
p = polyfit(t1(i0_10),E_f(i0_10,n_idx),1); % linear curve fit
y_T(:,n_idx) = polyval(p,t1);
T_60(n_idx) = -60/p(1); % p(1) is gradient in
dB/sec
end
for n_plot = 22:-1:7

leg = [num2str(sprintf('%.0f',F0(n_plot))),' Hz'];
[a,b] = max(abs(IR));
xi = (b-2000)/fs;
xf = (b+(2*fs))/(fs*0.9);
IR_f(:,n_plot) = IR_f(:,n_plot)/max(abs(IR_f(:,n_plot)));
figure
subplot(1,2,1)
plot(t1,IR_f(:,n_plot),'k-')
title('Respuesta al Impulso F1-P1 Cámara sin material','FontSize',10)
xlabel('Tiempo (s)','FontSize',10); ylabel('Amplitud','FontSize',10)
set(gca, 'FontSize',10);
% xlim([xi xf]);
% ylim([-25 25]);
legend(leg);
grid;
subplot(1,2,2)
plot(t1,E_f(:,n_plot),'b-','LineWidth',1)
hold on;
plot(t1,y_T(:,n_plot),'r:','LineWidth',2);
title('Curva de Decaimiento F1-P1 Cámara sin material','FontSize',10)
% ylim([-70 30]);
% xlim([0 5]);
xlabel('Tiempo (s)','FontSize',10); ylabel('Amplitud (dB)','FontSize',10)
set(gca, 'FontSize',10);
legend(leg,'least squared fit');
grid;
end
Registro fotográfico:

Medición TR Mejía Rojas

Cargado por

Copyright:

Formatos disponibles

Medición TR Mejía Rojas

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Medición TR Mejía Rojas

Cargado por

Copyright:

Formatos disponibles

Informe técnico medición de tiempo de reverberación

Mejía Soto Sebastián1, Rojas Sarria Camilo Andrés1.

Medir el tiempo de reverberación del aula 312B

Conocer de forma práctica el sonómetro y la instrumentación requerida para su correcto uso

● Sonómetro CESVA SC310

La calibración del sonómetro se verificó antes de la medición verificando con el pistófono

Sensación térmica 29°C

Punto de Rocío 16°C

Velocidad del viento 10 Km/h

Presión 1019 hPa

Tabla 1: Condiciones meteorológicas primer día de medición.

Figura 1: Ruido de fondo punto A

Figura 3: Ruido de fondo punto

Figura 5: Resultados micrófono 2 fuente 1 1 kHz

Segundo punto de fuente

Figura 6: Resultados micrófono 1 fuente 2 1kHz

Figura 8: Resultados micrófono 2 fuente 2 1 kHz

Tabla 2: Resultados de la medición y rt promedio por bandas de octava

Reemplazamos los valores obtenidos en la medición

Comparación del recinto

%Loading the audio files

for n_idx = 1:Nfc

for n_plot = 22:-1:7

También podría gustarte