Rosario, Argentina, 22 al 26 de junio de 2015

La historia sin fin: ¡Mi estudio suena mal!

¿Cómo resolverlo?

Indio Gauvrona, María Andrea Farinaa-b

a
Audio Engineering Society - Sección Argentina, bFacultad de Bellas Artes - UNLP
E-mail: in_dio_ar@yahoo.com.ar - maria_afar@yahoo.com.ar

“Mi madre me enseñó que no hace falta ser

un genio para dedicarse a la ciencia. Lo im-
portante es ser feliz con la investigación y
si de paso puedes mejorar un poquito el
mundo con tu trabajo, todavía mejor”.
Hélène Langevin-Joliot

1 Resumen: “Por desgracia no se puede explicar lo que es Matrix.

Has de verla con tus propios ojos”
A partir de algunas anomalías de carácter acústico percibidas en un Control de es-
tudio de grabación y, bajo la premisa de que las modificaciones implicaran mínimos
cambios en el espacio, se trabajó para hallar una posible solución que contempló el de-
sarrollo de un software específico de predicción acústica.
El software permitió recrear y analizar el recinto. Luego de incorporar los datos
del proyecto se confrontaron los valores con los de mediciones realizadas en el lugar.
Finalmente se diseñaron, evaluaron y verificaron las soluciones.
Con todas estas herramientas y procedimientos se alcanzó el objetivo buscado. El
presente trabajo detalla las etapas en la búsqueda de la solución destacando la relevancia
de las disciplinas intervinientes para llegar al resultado.

2 Introducción: “En el principio de los tiempos…” o “Hace mucho

tiempo en una galaxia lejana, muy lejana…”
El trabajo tuvo tres etapas. En la primera se realizaron mediciones acústicas en el
lugar, un relevamiento perceptual y un registro gráfico del espacio. En la segunda etapa
se estudió la documentación disponible y se realizó el análisis de la obra a partir del
plano de arquitectura. Con toda esta información, durante la tercera etapa, se diseñó un
software de predicción acústica a medida para poder modelizar las condiciones iniciales

1
y, relacionando estos datos con los obtenidos a partir de las mediciones, elaborar una
solución.
Entre los problemas detectados encontramos que el espectro en el punto de audi-
ción sufría irregularidades ajenas a la respuesta de los monitores; dicho comportamiento
se manifestaba como atenuaciones y realces en determinadas bandas de frecuencia. Este
fenómeno introdujo cambios en la escucha objetiva que modificaban el resultado final
(mezcla, toma, etc.). Además, la permanencia temporal de la energía (tiempo de rever-
beración) desde los lados izquierdo y derecho en el sweet spot era diferente.

3 Primera etapa: A ver, a ver… ¿Qué tenemos para trabajar?

(o Relevamiento perceptual y medición acústica)
3.1 Visto y considerando que era de Viktor… (o Consideraciones)
El diagnóstico del estado acústico del espacio, en ese momento, requería indefec-
tiblemente realizar una mensura y reconocer el lugar con un oído objetivo.
Las mediciones y pruebas realizadas in situ permitieron prever que la resolución
de los problemas tendría cierto grado de dificultad y, a partir del relevamiento realizado
en el lugar, comenzó la diagramación del proyecto. Como situación a destacar, el Con-
trol está ubicado en el exterior (país: Krakozhia, propietario: Viktor Navorski), por lo
tanto se decidió trabajar con un modelo que simulara la situación para realizar gran can-
tidad de pruebas, correlacionar respuestas y llegar a un resultado satisfactorio.

3.2 Escuchando con el micrófono (o Mediciones acústicas)

En primera instancia se verificó la calibración de niveles del sistema electroacús-
tico en su totalidad. Esta etapa es esencial, no se puede verificar el comportamiento
acústico de un control con la cadena de monitoreo descalibrada. El proceso se realizó
con base en +4 dBu ajustados a -18 dBFS (norma ITU/EBU -Europea-). Este dato es
útil para los técnicos que utilizan instrumental de medición específico en grabación o
mezcla, ya sea externo (hardware) o interno (software); por ejemplo: vúmetros, PPM,
loudness meter, etc.
Las mediciones acústicas permitieron evaluar: coherencia, tiempo de reverbera-
ción, distribución espectral y respuesta al impulso; en diferentes posiciones, condiciones
y sistemas de monitores.
Una instancia previa al viaje consistió en realizar pruebas con el software de me-
dición1. El núcleo de la cuestión era buscar algún método que resultara eficiente, con-
fiable y conveniente para el trabajo. El que dio mejores resultados (se llevó todos los
créditos o levantó muchos porotos con un ¡falta envido!), fue el MLS filtrado con ruido
rosa por la mejor relación señal/ruido de las respuestas obtenidas, además de cierta can-
tidad de promedios de las muestras válidas.
Una de las primeras mediciones y/o procedimientos importantes para el relevo de
los datos y su posterior procesamiento es la coherencia. En sistemas ruidosos, con osci-
laciones internas o no linealidades, como la respuesta de un parlante en un recinto, apa-
recen componentes que capta el micrófono en régimen permanente que no proceden de
su respuesta lineal a la excitación de entrada. Estos componentes generan un error en el
cálculo de la función de transferencia.
La función de coherencia es muy eficiente y precisa. Esta función proporciona
una buena estimación de la salida de energía debida exclusivamente a la señal de
entrada. La función de coherencia es una relación, en el dominio de la frecuencia, adi-
mensional que tiene sólo la parte real y entrega valores en el rango de 0 a 1.

1
Se utilizó en esta ocasión el software: SATlive Versión 1.50.6 del Ing. (FH) Thomas Neumann.

2
 Los valores de la función de coherencia pueden ser menos de 1 en los siguientes
casos: mediciones contaminadas por ruido no relacionado, mediciones contaminadas
por alguna fuente de ruido externo, sistema no lineal, entradas adicionales presentes en
el sistema.
Los dos sistemas de monitoreo presentaban solos2, sin la ingerencia del recinto,
una coherencia excelente, tanto para el L como para el R de cada uno.
Ahora ya en el sweet spot, ambos sistemas y monitores (L y R) mostraban un de-
terioro (casi siempre debido a reflexiones específicas) en 100 Hz y 200 Hz. Si bien la
coherencia no es la herramienta adecuada para visualizar los problemas, da un panora-
ma previo de algo que puede estar en conflicto en esas zonas; tal cual ocurrió poste-
riormente.
Otra de las sorpresas, que posteriormente se pudo comprobar a partir del estudio
del material gráfico (planos), fue la diferencia de los valores de tiempo de reverberación
de un lado (L) con respecto al otro (R).

Figura 1. Gráfico del tiempo de reverberación en base al T30 y T20

del monitor L (arriba) y R (abajo) de uno de los juegos de monitores.

Los valores de tiempo de reverberación obtenidos se muestran en la tabla 1.

2
Este tipo de medición se realiza a corta distancia, siempre en campo lejano, para evitar la “contamina-
ción” por parte del recinto. Si bien el recinto sigue estando, la diferencia de energía entre esa distancia
cercana comparada con la de la posición de escucha (sweet spot), suele ser lo suficientemente grande
como para tomar en cuenta el ambiente y los resultados muestran claramente esta distinción entre una
posición y la otra.

3
 Tabla 1. Valores de los parámetros acústicos3 del recinto excitado con
monitor L (arriba) y R (abajo) de uno de los juegos de monitores.

L 62 Hz 125 Hz 250 Hz 500 Hz 1000 Hz 2000 Hz 4000 Hz 8000 Hz Full

RT60(T30) - 0.39 s 0.27 s 0.34 s 0.31 s 0.29 s 0.26 s 0.19 s 0.20 s
RT60(T20) 0.26 s 0.26 s 0.23 s 0.20 s 0.27 s 0.24 s 0.22 s 0.19 s 0.21 s
EDT - 0.39 s 0.27 s 0.34 s 0.31 s 0.29 s 0.26 s 0.19 s 0.20 s
C50 12.8 dB 18.0 dB 19.1 dB 16.6 dB 14.7 dB 13.3 dB 16.0 dB 16.5 dB 17.9 dB
D50 95 % 98 % 99 % 98 % 97 % 96 % 98 % 98 % 98 %
C80 24.0 dB 21.4 dB 24.1 dB 23.2 dB 21.0 dB 20.4 dB 23.6 dB 25.5 dB 26.1 dB
D80 100 % 99 % 100 % 100 % 99 % 99 % 100 % 100 % 100 %
Ts 10 ms 11 ms 9 ms 16 ms 15 ms 14 ms 10 ms 11 ms 7 ms
R 62 Hz 125 Hz 250 Hz 500 Hz 1000 Hz 2000 Hz 4000 Hz 8000 Hz Full
RT60(T30) - 0.23 s 0.33 s 0.36 s 0.33 s 0.42 s 0.26 s 0.18 s 0.17 s
RT60(T20) 0.24 s 0.21 s 0.23 s 0.19 s 0.25 s 0.27 s 0.20 s 0.18 s 0.18 s
EDT - 0.23 s 0.33 s 0.36 s 0.33 s 0.42 s 0.26 s 0.18 s 0.17 s
C50 18.0 dB 16.4 dB 16.5 dB 15.8 dB 16.5 dB 14.8 dB 17.0 dB 19.8 dB 16.4 dB
D50 98 % 98 % 98 % 97 % 98 % 97 % 98 % 99 % 98 %
C80 24.9 dB 26.2 dB 24.7 dB 23.0 dB 21.9 dB 20.4 dB 24.0 dB 28.5 dB 100.0 dB
D80 100 % 100 % 100 % 100 % 99 % 99 % 100 % 100 % 100 %
Ts 9 ms 10 ms 9 ms 15 ms 13 ms 11 ms 11 ms 7 ms 10 ms

4 Segunda etapa: estudio de la documentación gráfica

4.1 ¡A la perinola! ¡No a la simetría!
A partir del plano del lugar, suministrado por la arquitecta de Viktor, en formato
de AutoCAD, se realizó el análisis morfológico. En primera instancia se corroboró la
simetría. Con la planta y el corte del Control se pudo comprobar con respecto a un eje
central longitudinal que el mismo no poseía simetría axial. El muro (luego llamado
este) se encontraba desplazado unos centímetros hacia afuera y los muros SE y SO no
tenían la misma silueta. La materialización de las superficies internas del Control
tampoco eran simétricas, por ejemplo, la ventana que da a la sala se enfrenta a un ta-
bique y una puerta.

3
RT60 es el parámetro más usado históricamente en acústica de salas. Es el tiempo (teórico) que le toma
al sonido decaer 60 dB por debajo de su nivel inicial, una vez que la fuente cesó su emisión. De acuerdo
con la norma ISO 3382 se calcula como dos veces el tiempo que le toma al sonido decaer desde -5 dB a -
35 dB. Este valor lo etiquetamos RT60 (T30).
Si no es posible medir hasta -35 dB debido a ruido de fondo u otros problemas entonces la norma sugiere
el uso de tres veces el tiempo de decaimiento entre -5 dB y -25 dB. Este valor lo etiquetamos RT60 (T20).
El parámetro EDT -que significa Early Decay Time y fue definido por Jordan en 1970- es otro tiempo de
reverberación. Se calcula como seis veces el tiempo de decaimiento entre 0 dB y -10 dB. En situaciones
normales, el EDT es el tiempo de reverberación que uno siente o percibe en un recinto y el RT60 es un
valor más técnico utilizado para los cálculos.
C50 y D50 (Reichardt, 1975; Thiele, 1953) son parámetros cuyos valores relacionan la energía que llega
en los primeros 50 ms con la energía total de la señal (D50); o con la energía que llega después de los 50
ms iniciales (C50). Estos valores son de gran utilidad para la inteligibilidad de la palabra. Los valores
mayores a 0 dB (C50) o a 50 % (D50) indican una buena inteligibilidad.
El cálculo de los parámetros C80 y D80 es similar al del C50 y D50, su diferencia radica en que evalúan
la energía que llega en los primeros 80 ms. Estos valores se utilizan para juzgar la inteligibilidad musical
en un recinto. Un buen valor de estos parámetros se sitúa entre 0 dB y 3 dB (C80), que se corresponde
con un 50 % (D80).
El tiempo central (Ts), de acuerdo con Dietsch y Kraak (1986), indica el baricentro temporal de la res-
puesta al impulso de la sala. Este centro de tiempo define cierto “centro de gravedad” temporal de la
respuesta al impulso de la sala.

4
 Figura 2. Esquema del Control puesto con la orientación para trabajar.

La primera conclusión del análisis morfológico fue la asimetría.

La simetría en el monitoreo es uno de los pilares fundamentales para la escucha
objetiva en un control. Todo lo mal que “suene” de un lado debe ser igual del otro; por
supuesto que lo ideal sería que ambos lados tengan una buena respuesta.
Se suele comentar con excesiva ligereza que es muy conveniente “sacar de parale-
lismo las paredes” y de esa forma evitar los tan temidos modos de resonancia. Bien,
tenemos malas noticias, una vez que el sonido se halla confinado, los modos van a estar
allí. No se los puede eliminar, salvo que no tengamos paredes!! Piensen: lo que geomé-
tricamente no es paralelo, acústicamente sí lo es. Las longitudes de onda involucradas
en la zona de modos de resonancia comparadas con las pequeñas diferencias de los án-
gulos entre paredes, hacen que los valores de frecuencia de los modos cambien tan poco
como el ancho de banda de cada modo.
Sí, hay que reconocer que este tipo de “tratamiento” (es decir, hacer que todas las
paredes tengan un ángulo diferente) sirve para evitar el “flutter eco” -conocido como
“eco flotante”- que es el reflejo de la onda sonora contra dos superficies paralelas que se
mantiene hasta extinguirse. Si uno emite un sonido breve (como un chasquido de de-
dos), se percibe como una ráfaga de chasquidos. El mismo fenómeno, desde la perspec-
tiva de la óptica, se genera cuando se colocan dos espejos paralelos, en ese caso la ima-
gen se repite hasta el infinito (y más allá!!).
Dada la importancia de la simetría, en este caso, uno de los lados completos esta-
ba desplazado de su posición “correcta” más de un decímetro.
Uno de los primeros interrogantes que surgió era, si se podría con un adecuado di-
seño del tratamiento acústico interior y, estableciendo una ubicación para los monitores
y el oyente, mejorar esta condición inicial.

5
 4.2 Obteniendo las coordenadas del contorno
Para esta fase se utilizó un comando del AutoCAD llamado “list”. El mismo
muestra el listado de las coordenadas de los puntos del contorno del control que definen
el perímetro.

Figura 3. Esquema del contorno del control junto al archivo de salida

con el listado de las coordenadas de todos los puntos que definen di-
cho contorno; todo en AutoCAD R14.

5 Tercera etapa: desarrollo del software… o diciéndole a la máqui-

na que debería hacer por nosotros
5.1 Acerca del software a emplear
Un lenguaje de programación es una jerga o juego de expresiones formales dise-
ñadas para declarar procesos que pueden ser ejecutados por máquinas como las compu-
tadoras.
La decisión del tipo de software para implementar una solución depende de mu-
chos factores, pero, uno de los más importantes a la hora de llevar a cabo algún trabajo
de optimización o resolución de problemas es saber cómo funciona para poder aprove-
char sus características de manera rápida y eficiente.
Particularmente en este trabajo se utilizó MS Excel (2003) por su sencillez y rapi-
dez para implementar y “probar” algoritmos de forma eficaz, además, por su potente
aliado VBA (Visual Basic para Aplicaciones), que trae incorporado librerías para inter-
actuar con la API (Interfaz de Programación de Aplicaciones) de Windows y cualquier
otro programa de manera muy natural. La razón principal, asimismo, fue por la amplia
aplicación y conocimiento del mismo.
Si bien se utilizó la planilla para algunos cálculos, la mayor parte del software está
desarrollado en VBA y la planilla actuaba de interfaz para las entradas de datos y sali-
das de resultados.

6
 5.2 Cercando el rancho (o Perímetro, oyente y monitores)
La primera fase de esta etapa fue la incorporación de las coordenadas de los pun-
tos del contorno del Control, obtenidos partir de la documentación de arquitectura, para
definir el perímetro. Dicho contorno quedó entonces definido a partir de la ecuación de
las rectas que lo delimitan para de esta forma poder operar matemáticamente (establecer
ángulos, intersecciones, simetrías, etc.). Se ingresaron en el software los 11 pares orde-
nados correspondientes a los puntos del plano. Se trabajó con la Ecuación general o
implícita de la recta Ax + By + C = 0 para utilizar matemáticamente matrices. Además,
se usó la ecuación explícita de la recta y = mx + b para ir corroborando y graficando
resultados parciales de los procesos previos. De esta manera se definieron las rectas en
el plano que formaban las paredes o la silueta del control.
Se determinó el eje de simetría, línea por la que se mueve el oyente hacia el frente
y atrás, y las posiciones para ubicar cada uno de los monitores, además del oyente.
Contorno
X Y 10 m

0.369 m 3.888 m
0.568 m 2.833 m 9m
1.000 m 2.223 m
1.000 m 1.016 m
5.245 m 1.016 m 8m

5.245 m 2.249 m
6.199 m 3.594 m 7m
5.338 m 8.169 m
3.812 m 8.589 m
2.512 m 8.589 m 6m

1.185 m 8.223 m
0.369 m 3.888 m 5m

Mon L
2.170 m 7.490 m 4m

Mon R
4.155 m 7.490 m 3m
VERDADERO

Eje de simetría 2m

3.162 m 1.000 m
3.162 m 8.600 m 1m

Oyente 0m

3.162 m 5.300 m 0m 2m 4m 6m

Figura 4. Coordenadas de los puntos del contorno (extraídos por el co-

mando “list” del AutoCAD), posiciones de ambos monitores y oyente.
A la derecha, figura resultante del contorno del control junto a los moni-
tores y oyente (planteados en escala con un gráfico de dispersión).

5.3 La ecuación de ondas y yo caminando por el Palacio de Versalles

en la Galería de los Espejos4 (o Ecuación de onda esférica de las
imágenes espejo)
El siguiente paso fue calcular las posiciones de las fuentes fantasmas (FFn). Para
ello se obtuvieron las coordenadas de las imágenes espejo que tenía cada monitor en
relación al contorno y fue necesario utilizar ecuaciones de simetría con respecto a una

4
Obra del arquitecto Jules Hardouin-Mansart.

7
recta en el plano. De esta forma al tener dichas coordenadas y la del oyente, se podían
calcular las rectas que los unía (acústicamente fueron considerados como rayos sonoros)
y definir si cortaba o no a esa superficie del contorno (graficada como una línea, ya que
se trabajó todo en el plano por estar monitores y oyente a la misma altura).

⎛ ⎛ X ⎞ ⎞
Xp ⎜ m⎜ Yp + p − b ⎟ ⎟
⎜ ⎝⎜ ⎟
Yp + −b Y = 2⎜
m ⎠ + b⎟ − Y
m ⎟ p
X =2 − Xp ⎜
1
⎟
1 m+
m+ ⎜ m ⎟
m ⎝ ⎠

Espejito, espejito ¿quién es el más bonito?

Se estableció que las fuentes fantasmas primarias que intervienen, por la forma
que tiene el lugar, son cinco y esto tiene que ver con algunas de las posiciones relativas
del monitor y del oyente; y donde, además, algunas de las paredes del perímetro no in-
tervienen. Se trabajó con las ecuaciones de las rectas que definían las paredes y las tra-
yectorias.

FF3
FF2

FF4
FF1

FF5

Figura 5. Posición de cada imagen espejo (FFn) verificadas en Au-

toCAD.

8
 Posiciones fantasmas Distancias
0.001 m 7.898 m O 4.092 m
1.655 m 9.358 m NO 4.328 m
2.170 m 9.688 m N 4.499 m
11.208 m 1.079 m SE 9.085 m
8.536 m 8.688 m E 6.353 m
Oyente
3.162 m 5.300 m

Figura 6. Coordenadas de cada imagen espejo (FFn).

Como muestran las Figuras 5 y 6, se calcularon en el software los 5 pares ordena-

dos correspondientes a las fuentes fantasmas. Se establecieron las ecuaciones de las
rectas que pasaban por cada una de las coordenadas de las FFn y el oyente, se calculó la
distancia de cada una de esas rectas y los puntos de intersección entre la recta que le
corresponde a cada uno de los tabiques y la que le corresponde a cada una de las FF y el
oyente (Figura 7).

Ecuaciones rectas de paredes Distancia Mon-oyente

2.404 m
A B C
4.3356 -0.8161 1.572482 O
0.3657 -1.3271 10.479714 NO
0.0000 -1.3000 11.1657 N Pared O
-1.3455 0.9544 4.9105455 SE 1.185 m 8.223 m
-4.5744 -0.8610 31.451226 E 0.369 m 3.888 m
m = 5.31258424
Ecuaciones rectas de FF y Oy
b = 1.92682516
A B C
2.59824148 3.1611 -24.97034 O Pared NO
4.05753771 1.5069 -20.81785 NO 2.512 m 8.589 m
4.388 0.9923 -19.13536 N 1.185 m 8.223 m
-4.2206695 -8.0454 55.987504 SE m = 0.27556326
b = 7.89670243
3.3881968 -5.3736 17.765533 E
Coordenadas de intersección entre pared y FF-Oy Pared N
X Y 3.812 m 8.589 m
0.97357031 7.09899946 O OK 2.512 m 8.589 m
1.99385946 8.44613684 NO OK m= 0
2.41852728 8.589 N No cumple b= 8.589
6.25732425 3.67632472 SE No cumple
Pared SE
5.58981819 6.83061717 E OK 5.245 m 2.249 m
6.199 m 3.594 m
m = 1.40978625
b= -5.145165

Pared E
6.199 m 3.594 m
5.338 m 8.169 m
m= -5.312892
b = 36.5287174

Figura 7. Ecuaciones de las rectas de las paredes, las que pasan por
cada imagen espejo (FFn) y el oyente; y las coordenadas de la inter-
sección entre ambas para verificar si corta o no a dicha superficie.

Es muy importante verificar dicha intersección debido a que en determinadas po-

siciones de monitor y oyente la pared no refleja el sonido y, en consecuencia, no tendrá
influencia en la resultante final. Posteriormente cada superficie será tratada (o no) con
materiales absorbentes con coeficientes de absorción específicos para cada frecuencia.

9
 Lo último en agregarse fueron el piso y cielorraso respectivamente. Debido a que
todo el trabajo se realizó en el plano, como se mencionó anteriormente, la adición de
estas nuevas superficies se trataron independientemente; ya que los datos necesarios
para el cálculo eran las distancias a cada FFn y su intersección con la superficie en cues-
tión. Es muy lógico afirmar que tanto el piso como el cielorraso siempre van a estar
reflejando y debido a ello no se hicieron los sistemas de ecuaciones para dichas inter-
secciones. Cabe aclarar que al piso se lo trató de forma diferencial, ya que en la trayec-
toria de su imagen espejo se interponía la consola con su mobiliario.
Con todos estos datos reunidos, a saber:

• coordenada de la posición de los monitores

• coordenada de la posición del oyente
• coordenadas de las posiciones de cada FFn
• verificación de reflexiones por parte de cada FFn

se calculan las distancias del recorrido desde cada FFn hasta el oyente y con esos datos,
usando la solución general para una onda esférica simétrica, en el dominio de los núme-
ros complejos, se trabajó con la parte real y la imaginaria.
La suma de ambas partes (de cada una de las fuentes) nos da la resultante final
que posee la información de amplitud (asignando un valor de sensibilidad al monitor) y
fase (por tener la parte imaginaria). Este resultado, que se obtiene con el software, se
observa en trazo color rojo en la Figura 8. El cálculo completo, además, tomó en cuenta
la respuesta del monitor solo (a partir de una medición del mismo realizada muy cerca
con el objeto de obtener una respuesta sin interferencia del recinto, similar a medirlo en
campo libre; trazo color marrón) y, en base al resultado del software, obtuvimos la res-
puesta real del monitor en ese lugar específico; trazo color negro grueso.

100 dB Corr: 0.873163215

95 dB

90 dB

85 dB

80 dB

75 dB

70 dB

65 dB

60 dB

55 dB

50 dB
10 Hz 100 Hz 1000 Hz

Figura 8. Resultados parciales y finales juntos. Trazo marrón: moni-

tor en campo libre. Trazo rojo: resultado del cálculo mediante la ecua-
ción de onda con todas las superficies involucradas. Trazo azul: medi-
ción realizada en el sweet spot del estudio. Trazo negro: respuesta del
monitor en campo libre combinada con el resultado del cálculo; es de-
cir, resultado final para su evaluación.

10
 5.4 Funka o no funka (Volver al futuro o Verificación)
Para comprobar el funcionamiento del trabajo realizado hasta el momento, se bus-
có una posición cuyo resultado correspondiese con una de las mediciones realizadas en
el lugar. La idea era encontrar una posición donde esta aproximación teórica coincidiese
o se acercara a lo medido, teniendo en cuenta que existían diferencias entre el modelo
teórico (basado en ecuaciones) y la situación en obra.

Cachavacha barría… con el monitor

Una de las primeras acciones a tomar para esta corroboración era ir cambiando la
posición del monitor y con cada una de ellas, sacar una resultante.
El área sombreada en color azul de la Figura 9 determina la zona donde el monitor
se podía mover y donde podría llegar a ubicarse para rastrear la posición inicial.

10 m

9m

8m

7m

6m

5m

4m

3m

2m

1m

0m
0m 2m 4m 6m

Figura 9. Zona o área de consecutivas posiciones del monitor para

rastrear la posición que correspondiese con la medición en el lugar.

Esta región se utilizaría posteriormente para hallar la posición más adecuada, jun-
to con un tratamiento acústico específico, cuya respuesta sea similar a la del monitor en
campo libre. Debido a que se operaba en estéreo, se estableció como límite 8º de corri-
miento como máximo a partir del eje de escucha óptimo de 30º.
Se elaboró un algoritmo que recorría la zona y por cada posición del monitor se
obtenía un gráfico resultante. El paso mínimo que movía la posición del monitor era de
1 cm. A partir de indicar cada cuántos centímetros debía moverse el monitor, el mismo
se ubicaba en una posición y se exportaba, como resultado, un gráfico cuyo nombre
correspondía a las coordenadas de la posición (Figura 10).

11
 For X = Xmin to Xmax
For Y = Ymin to Ymax

no
X>RIS Sig X

si

no
X>RII Sig X

si

no
X<RDS Sig Y

si

no
X<RDI Sig Y

si
Área Gráfico Sig X

Figura 10. Zona y diagrama de flujo de las consecutivas posiciones

del monitor y exportar el gráfico.

Los resultados se exportaban a una carpeta o directorio establecido previamente desde el

soft y allí se iban amontonando… (ver Figura 11).

Figura 11. Ventana de: “Guardar como” (API de Windows) donde se

indicaba el directorio donde se exportarían los gráficos de resultados.
Ventana con los gráficos de resultados exportados en formato .gif con
sus nombres indicando las coordenadas de cada posición del monitor.
Dos gráficos elegidos al azar de resultados.

12
 Algo que vino del cielo... ¡mirando estrellas! (o Star Trek, Stargate o
Armageddon)
Los gráficos obtenidos se estudiaron con un método que se utiliza en astronomía5
que consiste en ir pasando las imágenes una a una, en forma rápida, y poder seleccionar
más fácilmente las de interés6.

¡Ufa! ¡cuántos gráficos! (o Los números y su historia7)

Tomando para la evaluación que el monitor se corriera de a 1 cm se obtenían co-
mo consecuencia más de 4300 gráficos. Como podemos observar en la figura 12 la va-
riación entre los resultados de las imágenes obtenidas con este valor era significativa.

Figura 12. No hay pregunta que por bien no venga: y un amigo nos
preguntó: ¿Tiene sentido cada 1 cm? Gráficos obtenidos cada 1 cm.

Y ahora… ¿quién podrá ayudarnos?

Debido a que la cantidad de gráficos obtenidos era muy grande implementamos la
utilización del coeficiente de correlación8, lo que nos permitió disminuir la cantidad de

5
Para medir la distancia a una estrella se observa la magnitud aparente que se relaciona con su brillo. Las
estrellas tipo “cefeidas” en particular, poseen un cambio de brillo en un determinado lapso de tiempo que
es constante y que permite determinar, midiendo dicho período de variación, a qué distancia se encuen-
tran. Con ese patrón los astrónomos pueden obtener la distancia al resto de las estrellas que componen
una galaxia lejana y así hallar su lugar en el universo.
El método utilizado para “cazar” cefeidas consistía en observar la fotografía de una parte del cielo, y
luego la misma imagen del cielo de una fecha posterior. Pasaban una y otra imagen en forma rápida en un
dispositivo que se llama: “comparador de destellos”, para determinar dónde estaba el cambio; y de ésta
manera las descubrían.
6
En el libro “Corazones solitarios en el cosmos”, de Dennis Overbye, se cuenta cómo usaban el compa-
rador de destellos. Ese “cómo lo usaban” quedó en forma de recuerdo latente en mi memoria y es la géne-
sis de donde surgió mi idea inicial de exportar los gráficos para verlos como una animación y poder deci-
dir entre ellos; aún muchos años después de haber leído el libro (Gauvron 1998).
7
Una muy linda obra del escritor y divulgador científico Isaac Asimov, donde hay interesantes compara-
ciones de números grandes y pequeños.

13
imágenes que cumplían cierto criterio (el de parecerse a la medición en el lugar). Para
ello se estableció un “umbral” a partir del cual se exportaba el gráfico, o no. Los pro-
blemas en el Control estaban en un rango de frecuencias específico y su verificación
sólo tenía sentido en una fracción de las frecuencias bajas, entonces, la correlación se
acotó a ese rango del espectro.

100 dB Corr: 0.873163215

Cov( X , Y )
ρ( x, y ) =
95 dB

σx·σy
90 dB

85 dB

80 dB

75 dB
Altura mon 1.150 m
Correlación
Altura techo 4.100 m 70 dB
M onitor sólo
65 dB

Límites p correleción Descartar por 60 dB

Lim inferior correlación
35.0 Hz VERDADERO
55 dB
Lim superior 400.0 Hz Corr mayor a 0.94
50 dB
10 Hz 100 Hz 1000 Hz

Figura 13. Sección del software donde se establece el umbral y rango

en el espectro del coeficiente de correlación y su descarte si no cum-
plía dicho criterio. Imagen de resultado con su correspondiente valor
del coeficiente de correlación.

Jugando a encontrar las diferencias

Como mencionamos anteriormente, el trazo azul correspondía a la medición en el
Control, que era inicialmente la curva que estábamos buscando y que manifestaba los
problemas detectados. El trazo negro es el resultado del software con el monitor. Lo que
buscábamos era la coincidencia entre los dos trazos, principalmente en los lugares que
habíamos detectado como significativos (meseta en 60 Hz y valles en 100 Hz y en 200
Hz). Al analizar los gráficos nos encontramos, en ocasiones, con valores altos de coefi-
ciente de correlación que no implicaban semejanzas de ese orden con la curva buscada,
por eso seguía siendo necesario mirar cada gráfico con el objeto de hallar el más coinci-
dente en las zonas de interés. Este comportamiento dio como finalidad el acotamiento
del espectro de evaluación.

Espejos siniestros 7…
Inicialmente trabajamos las superficies interiores del control como si las mismas
reflejaran toda la energía; pero, le podíamos incorporar el dato del coeficiente de absor-
ción alfa -entre 0 y 1- a cualquiera de los planos. En una primera instancia todas las su-
perficies se evaluaron con el alfa igual a 0. El siguiente paso consistió en comenzar a
probar sacando del cálculo alguna de las superficies (simplemente asignándole el valor
1 de alfa -la pared absorbe todo-) y como resultado las curvas se modificaron de una
manera significativa.
Considerando en el cálculo el piso, el cielorraso y las paredes (un total de 7 super-
ficies) armamos una combinatoria para que reflejen o absorban toda la energía con todas
las formas de combinación posible entre ellas. Esto daba un total de 127 combinatorias9.

8
En una distribución bidimensional puede ocurrir que las dos variables guarden algún tipo de relación
entre si. El coeficiente de correlación mide el grado de intensidad de esta posible relación entre las varia-
bles.
9
Resolver estas combinatorias fue tan sencillo como hacer un contador de 0 a 127 que pasara a binario de
7 dígitos el valor; cada dígito binario se trasladaba a cada una de los planos como su alfa (0 o 1).

14
Para cada una obteníamos un gráfico que nos indicaba además qué paredes o planos
estaban involucrados.

Ecuaciones rectas de paredes

A B C alfas
4.3356 -0.8161 1.572482 O 1.0
0.3657 -1.3271 10.479714 NO 1.0
0.0000 -1.3000 11.1657 N 1.0
-1.3455 0.9544 4.9105455 SE 1.0
-4.5744 -0.8610 31.451226 E 1.0
Ecuaciones rectas de FF y Oy 1.0 Te
A B C 0.5 Pi

Figura 14. Sección del software donde se asignan los alfas de cada
plano y parte de la combinatoria que asignaba automáticamente para
cada combinación los valores 0 y 1.

Como se muestra en la Figura 15, éste procedimiento nos permitió identificar una
posición donde la respuesta era bastante coincidente con la medición y los resultados;
además posibilitaron verificar que el software daba resultados compatibles con lo medi-
do (habiendo partido de la ecuación de onda, ecuaciones de las rectas, etc.) pudiendo
llegar al modelo.

Figura 15. Resultado hallado que corrobora el modelo con lo medido.

5.5 Disipando la energía acústica (o Tratamiento acústico interior)

La siguiente etapa consistió en trabajar con el modelo ya probado en la posición
coincidente con las mediciones, pero, ahora, incorporándole absorción a las superficies
interiores del control que necesitaban cierta regulación de sus reflexiones.
Las tablas que detallan los valores de los coeficientes de absorción (alfa) de los
materiales de construcción estudiados en la bibliografía, en general, tienen datos por
bandas de octava desde 125 Hz a 4 kHz. En este caso, necesitábamos valores que llega-
ran hasta los 40 Hz por el tipo de monitores que había en el Control y, además, la in-
formación la estábamos trabajando en 1/32 de octava. Todo el cálculo desarrollado está

15
realizado en 1/32 de octava porque la medición está hecha con esa resolución y, fue
necesario mantener la resolución con el objeto de poder comparar todo bajo iguales
condiciones.

Nuestra amiga, la Ingeniera, tenía guardado en el baúl…

Nuestra amiga en cuestión nos facilitó los datos correspondientes a una medición
de lana de vidrio de tres espesores diferentes con valores desde 100 Hz y hasta 5 kHz;
hechas en tercios de octava. Ver Figura 16.

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0
10 Hz 100 Hz 1000 Hz 10000 Hz

Figura 16. Valores de coeficientes de absorción (alfa) para lana de vi-

drio en 3 espesores diferentes. Las muestras están en resolución de
tercios de octava y van desde 100 Hz a 5 kHz.

Básicamente la solución de los problemas detectados en el Control se concluyó

con el agregado de absorción de lana de vidrio en determinados lugares específicos y la
ubicación más ajustada al cálculo de los monitores y el oyente. Una de las premisas
principales desde el inicio era que los cambios tenían que ser mínimos.

Un paseo por las nubes de Lagrange, Whittaker-Shannon, Hermite y

otros
A partir de estos valores originales, el primer paso fue aproximarlos a la resolu-
ción con la que estábamos trabajando (1/32 de octava), extender los datos hacia la parte
baja del espectro y ampliar el alcance para mayores espesores. Para ello se aplicaron
métodos de interpolación10 y extrapolación11. Después de probar varios métodos, empe-
zando por Lagrange, pasando por las polinómicas, revisando las ideas de Newton, ter-
minamos con trazadores cúbicos. Se interpoló, extrapoló y extendió por espesor.
Ahora conocíamos las ecuaciones que definían las curvas originales, lo que nos
permitiría trabajarla. Se obtuvieron todos los valores en 1/32 de octava y se extrapolaron
los datos hasta 40 Hz (Figura 17).

10
Obtención de nuevos puntos intermedios, partiendo de un conjunto discreto de valores.
11
Proceso de estimar más allá del intervalo de observación original, el valor de la variable en base a su
relación con otra variable.

16
 Para determinar el comportamiento del material con un mayor espesor también se
utilizó el trazador cúbico. Se extrapoló para obtener nuevos valores de lana de vidrio
que correspondían a un espesor mayor que el de las muestras de las cuáles teníamos los
datos de las mediciones (trazo rojo, figura 17).

1
0.9

0.8

0.7
0.6

0.5
0.4

0.3

0.2

0.1

0
10 Hz 50 mm 100 Hz 70 mm 100 mm 1000Int
Hz50 Int 70 10000
Int 100Hz
prom 50 prom 70 Espesor

Figura 17. Valores de coeficientes de absorción (alfa) como resultado

de las interpolaciones, extrapolaciones y extensión del espesor; todo
en resolución de 1/32 de octava.

Así se obtuvo el coeficiente de absorción final de esa lana de vidrio en varios es-
pesores a partir de la frecuencia de 40 Hz y con una resolución de 1/32 de octava. Estos
nuevos datos los incorporamos al programa.

1.00

0.80

0.60

0.40

0.20

0.00
10 Hz 100 Hz 1000 Hz

Figura 18. Curva final de coeficientes de absorción (alfa) usada en el

programa para los cálculos en 1/32 de octava.

17
 6 Y colorín, colorado te presento el resultado!!
El resultado final del cálculo se observa en la siguiente figura donde el trazo de
color negro es una copia muy aproximada del monitor solo (curva marrón). Si bien se
advierte que hay menos energía, no se observan las complicaciones con las que comen-
zamos el desarrollo del trabajo.
La curva negra (Figura 19) es el resultado al que arribamos con la implementación
en el modelo del tratamiento acústico interior en el control. Éste consistió en colocar
trampas de absorción con alfa y espesor específico en determinados lugares puntuales y
ubicar monitores y oyente en el espacio disponible.

Figura 19. Curva final como resultado del cálculo con el tratamiento
acústico completo.

Figura 20. Curva final, lugares de las trampas con alfa y espesor es-
pecífico, coordenadas de monitores y oyente en el control.

Este modelo se realizó en base únicamente a las reflexiones de primer orden. Para
no incorporar al cálculo las reflexiones de segundo orden se elaboraron dos hipótesis.
La primera se basaba en que tener una mayor cantidad de reflexiones provenientes des-
de distintos lugares iba a hacer que la respuesta fuera más plana. La segunda hipótesis
se relacionaba con la energía de esas reflexiones y, debido a los materiales absorbentes

18
que íbamos a incorporar en el proyecto podían resultar irrelevantes. Lo último que hici-
mos fue incorporar las reflexiones de segundo orden que dio validez a la primera hipó-
tesis. La gran cantidad de reflexiones, en especial si las paredes no tenían algún tipo de
tratamiento absorbente, hacía que la resultante fuera realmente más plana (Figura 21).

Figura 21. Comparación entre tomar en cuenta las reflexiones de pri-

mer orden (superior) y contemplar las de segundo orden (inferior); en
ambos casos con las superficies de alfa 0 (reflejo total).

Jugando a conocer otras disciplinas o Connections12

Para el desarrollo del presente trabajo se interactuó, estudió y trató con conceptos
y contenidos provenientes de diferentes disciplinas que nos interesa destacar: acústica,
álgebra, física, geometría analítica, métodos numéricos, astronomía, estadística, pro-
gramación, análisis matemático, geometría descriptiva y computación.

En el libro Música y Músicos de Robert Schumann figura entre las Reglas y

máximas para jóvenes músicos una que nos interesaría destacar en el final de este traba-
jo “El estudio no tiene fin”.

7 Paseando por la biblioteca real, no la de Babel13 (o bibliografía

consultada)
- Beckett, Royce (1983). “Numerical calculations and algorithms”. Robert E.
Krieger Publishing Company Inc. USA.
- Burden, Richard y Faires J. Douglas (1998). “Análisis numérico”. International
Thomson Editores. México.

12
Serie documental de la BBC creada, escrita y presentada por el locutor, historiador de la ciencia, autor
y productor británico de televisión James Burke. En Argentina se emitió con el nombre de: “Relaciones”.
13
La biblioteca de Babel es un cuento de Jorge Luis Borges. El relato es la especulación de un universo
compuesto de una biblioteca de todos los libros posibles, en la cual sus libros están arbitrariamente orde-
nados, o sin orden, y preexiste al hombre.

19
- Hildebrand, Francis Begnaud (1987). “Introduction to Numerical Analysis”.
Dover Publications. USA.
- Penick, Tom (2000). “Engineering acoustics”. Teicontrols. USA.
- Petrocelli, J. A. (1982). “Geometría analítica y métodos numéricos”. Apuntes
UTN, Regional Buenos Aires. Argentina.
- Sánchez Villegas, A. S. (1979). “Análisis numérico” -apuntes recopilados por A.
S. S. Villegas-. Universidad Politécnica de Madrid – Escuela Técnica Superior de
Ingenieros Industriales – Sección Publicaciones. Madrid.
- Scheid, Francis y Di Constanzo, Rosa Elena (1991). “Métodos numéricos”.
Schaum, Mc Graw-Hill. México.

Agradecimientos
A Adrián Rallap.
A la Ing. Nilda Vechiatti.
A Maximiliano Rispoli.
A Martín Wolmy.

20