Manual ARTA
Manual ARTA
Manual ARTA
Programa para Medicin de Respuesta de Impulsos y Anlisis en Tiempo Real del Espectro y de la Respuesta en Frecuencia
Manual de usuario
Versin 1.2.0
Ivo Mateljan Laboratorio de Electroacstica Facultad de Ingeniera Elctrica R. Boskovica bb. 21000 Split, Croacia Traduccin de: Andrea Gutirrez Prez andreag@agptraducciones.com
1 Introduccin
ARTA es un programa para la medicin de la respuesta de impulsos, el anlisis de espectros en tiempo real y la medicin en tiempo real de la respuesta en frecuencia. Es una herramienta para realizar mediciones acsticas y pruebas "punto por punto" de la calidad del sonido en los sistemas de comunicacin. ARTA dispone de funciones de los siguientes sistemas de medicin: 1. Sistema de medicin de respuesta de impulsos con generadores de seales: ruido peridico blanco, ruido peridico rosa, MLS, ondas sinusoidales de barrido lineal y logartmico. 2. Analizador de Fourier de doble canal con generadores de seales: ruido blanco, ruido rosa, ruido peridico blanco y ruido peridico rosa. 3. Analizador de Fourier de canal sencillo con generadores de seales: ruido peridico blanco y ruido peridico rosa. 4. Analizador de espectros, de banda de octavas y THD con generadores de seales: onda sinusoidal, dos rfagas sinusoidales, multitono, ruido blanco, ruido rosa, ruido peridico blanco y ruido peridico rosa. Nota: El modo 2 y 3 tambin se pueden usar para estimar la respuesta de impulsos.
ARTA es un poderoso analizador de: 1. 2. 3. 4. 5. 6. 7. 8. Respuesta de frecuencias limitada, Respuesta de frecuencia suavizada (en bandas de 1/n octavas) Respuesta de etapa, Envolvente de la respuesta de impulsos (curva - ETC) Curva de deterioro espectral acumulado, Deterioro de la energa en los entornos reverberantes, Parmetros acsticos de la sala, Medidas de inteligibilidad del habla: MTF, STI, RASTI, %AL.
1.1 Requisitos
Los requisitos para usar ARTA son:
Sistemas operativos: Windows 98/Me/2000/XP Procesador Pentium, Pentium III o superior, con una frecuencia de reloj de 400 MHz o
ARTA se ha usado con resultados satisfactorios con las siguientes tarjetas de sonido:
RME Fireface 800, RME DIGI96, RME HDSP Terratec EWX 24/96 Firewire FW X24 M-audio Audiophile 2496, USB Transit, Delta 44, YAMAHA GO 46 Digigram VxPocket 440 una tarjeta PCMCIA para porttil Echo Layla 24 Echo Indigo TASCAM US-122 - USB audio ESI Quatafire 610, U24 USB and Waveterminal, Soundblaster X-Fi Soundblaster Live 24 y Audigy y Extigy-USB (pero slo a la frecuencia de muestreo de 48kHz), Tarjetas Turtle Beach Pinnacle y Fuji, ARTA se puede usar con un rendimiento ligeramente menor con las siguientes tarjetas:
Soundblaster MP3+ USB (nota: no instale un controlador USB, use el controlador por defecto
de Windows XP),
Tarjetas de sonido y audio incorporado con cdecs AC97 (problema con ruido alto en el modo
FTT). La instalacin de este software es simple: Siga el programa de configuracin de ARTA o simplemente copie los archivos "ARTA.exe" y "ARTA.chm" en una carpeta y haga un acceso directo a "ARTA.exe". Todos los datos del registro se guardarn automticamente la primera vez que se ejecute el programa. Los archivos con la extensin ".PIR" se registran para que se abran con el programa ARTA. Contienen los datos de la respuesta de impulso peridica (PIR). Los resultados de otros tipos de mediciones (respuesta en frecuencia y espectro) se pueden guardar como archivos de texto con formato ASCII, o como archivo de superposicin. ARTA puede exportar e importar archivos de diferentes formatos (wav, mlssa, tim y txt). ARTA no enva los grficos a la impresora, pero se pueden copiar todos los grficos en el Portapapeles de Windows y pegar en otras aplicaciones de Windows.
La configuracin general para medicin de las pruebas de sistema es la que aparece en la Ilustracin 1.1. El canal de salida de seal izquierdo de la tarjeta de sonido se usa como salida del generador de seales. La entrada de lnea izquierda se usa para medir el voltaje de salida del dispositivo a prueba y la entrada de lnea derecha de la tarjeta de sonido se usa para medir el voltaje de entrada del dispositivo a prueba. En una configuracin de canal sencillo slo se registra un voltaje de salida del dispositivo a prueba. En una configuracin de canal semi-doble la entrada de la lnea derecha se usa para medir el voltaje de la salida de la lnea derecha. En la modalidad de bucle, la salida de lnea izquierda est conectada a la entrada de lnea izquierda y la salida de lnea derecha se conecta a la entrada de la lnea derecha.
Ilustracin 1.1 Configuracin de medicin general para pruebas de respuesta del sistema (D.U.T. = dispositivo a prueba)
Las configuraciones para realizar mediciones acsticas se muestran en las Ilustraciones 1.2, 1.3, 1.4 y 1.5.
Para proteger la entrada de la tarjeta de sonido frente a los voltajes elevados que se generan en el amplificador de potencia, se recomienda usar el circuito de sonda de voltaje, como aparece en la Ilustracin 1.3. Los valores de las resistencias R1 y R2 se han de elegir para una atenuacin arbitraria 5
ARTA User Manual (es decir que R1=8200 y R2=910 ohmios nos da una sonda con una atenuacin de -20,7 dB (0,0923) si la tarjeta de sonido tiene una impedancia normal de entrada de - 10k).
Ilustracin 1.3 Sonda de voltaje con proteccin frente a sobrecarga del canal de entrada de la tarjeta de sonido.
ARTA User Manual ARTA est diseado para realizar pruebas "punto por punto" de la calidad del sonido en los sistemas de comunicacin. La Ilustracin 1.7 muestra la configuracin que se usa para probar esos sistemas. Se puede obtener una interfaz para telfonos mviles usando la entrada/salida de unos cascos auriculares. La interfaz para la lnea de telfono normal (POTS) se muestra en la Ilustracin 1.8. ARTA puede medir la respuesta de impulso y en frecuencia, las distorsiones de seales sinusoides, sinusoides dobles, multitonos, estimar retardos, ecos e ndices de transmisin del habla. Se aplica una tcnica especial de medicin, con una excitacin de ruido interrumpida, para evitar el comportamiento variable en el tiempo de estos sistemas (control de ganancia automtica, reduccin de ruido, activacin de la voz).
Ilustracin 1.8 Interfaz entre la salida/entrada de la tarjeta de sonido y la lnea de telfono normal (POTS)
Ilustracin 1.9 Ventana de respuesta de impulsos Usando el men Mode [Modo], podr cambiar entre las tres ventanas del dominio de la frecuencia para realizar anlisis en tiempo real.
Ventana de mediciones de respuesta en frecuencia de canal doble Ventana de mediciones de respuesta en frecuencia de canal sencillo Ventana del analizador de espectros
El modo de medicin se puede escoger tambin haciendo clic sobre los siguientes iconos de la barra de herramientas: FR2 - Ventana de mediciones de respuesta en frecuencia de canal doble FR1 - Ventana de mediciones de respuesta en frecuencia de canal sencillo SPA - Ventana del analizador de espectros IMP - Ventana de respuesta de impulsos Ahora haga clic sobre estos mens o iconos de la barra de herramientas para ver cmo funcionan las ventanas de medicin.
Ilustracin 1.10 Ventana de respuesta en frecuencia en canal doble - FR2 (la ventana respuesta en frecuencia en canal sencillo - FR1 es igual)
ARTA User Manual La ventana de la respuesta de impulsos es la ms importante para el anlisis de la respuesta de un sistema. Se describir con ms detalle una vez que veamos cmo se analiza el espectro y la respuesta de frecuencia del sistema.
Ilustracin 1.12 Cuadro de dilogo para la configuracin de los dispositivos de audio El cuadro de dilogo Configuracin de dispositivos de audio tiene los siguientes controles:
En la seccin Sound Card [Tarjeta de sonido]: Input device [Dispositivo de entrada] elija la tarjeta de sonido que funciona como dispositivo de entrada. Output device [Dispositivo de salida] elija la tarjeta de sonido que funciona como dispositivo de salida. Generalmente habr de escoger la misma tarjeta como dispositivo de entrada y salida (no ser necesario con los modos de Analizador Fourier de Doble Canal o de Analizador de Espectros). Wave format [Formato de la onda] elija el formato de onda Windows de 16 bits, 24 bits o 32 bits. ARTA utiliza los controladores de sonido de Windows habituales (WMME para Windows 98 y WDM formato extensible para Windows XP). Le recomendamos el modo WMME y 16 bits. As que dejaremos la casilla de verificacin 'Extensible' sin activar. Por supuesto podr experimentar con modos diferentes, pero tenga en cuenta que el modo de 32 bits no funciona con muchas tarjetas de sonido en formato WDM Extensible. Se recomienda usar los modos de 24 o 32 bits slo si tiene una tarjeta de sonido de muy alta calidad (muchas tarjetas de sonido se anuncian como de 24 bits, pero su resolucin real es inferior a 16 bits). En la seccin I/O amplifier interface [Interfaz del amplificador de Entrada/salida]:
10
Antes de continuar con el programa ARTA, asegure los siguientes valores por defecto.
Nota importante: Por favor, enmudezca los canales de seal y micrfono en la mezcladora de salida de la tarjeta de sonido. De no hacerlo as puede obtener una retroalimentacin positiva durante las mediciones. Si usa una tarjeta de sonido profesional, desactive el control directo o de latencia cero de las entradas de lnea.
Ilustracin 1.13 Configuracin tpica de la salida de la mezcladora de la tarjeta de sonido (la encontrar en el cuadro de dilogo del panel de control 'Sound and Audio Devices' [Dispositivos de sonido y audio])
11
1.5 Calibracin
El comando de men Setup > Calibrate devices [Configuracin > Calibracin de dispositivos] abre el cuadro de dilogo 'Soundcard and Microphone Calibration' [Calibracin de la tarjeta de sonido y del micrfono] que aparece en la ilustracin 1.14.
Ilustracin 1.14 Cuadro de dilogo para calibrar la tarjeta de sonido y el micrfono Hay tres secciones que nos conducen a la calibracin de (a) canal izquierdo de salida de la tarjeta de sonido, (b) canales izquierdo y derecho de entrada a la tarjeta de sonido y (c) calibracin del micrfono.
12
ARTA User Manual Por ejemplo, el archivo "MB550-B.mic" (que aparece en la Ilustracin 1.15) tiene el siguiente contenido:
Micrfono mb550 frec(Hz) Magn(dB) 48.280 0.34 48.936 0.28 49.601 0.21 .....
Ilustracin 1.15 Respuesta en frecuencia tpica para un micrfono electret. La casilla de verificacin "Show spline interpolated data" [Mostrar datos interpolados de estra] nos permite ver la curva de la respuesta en frecuencia interpolada que se usar en la compensacin de la respuesta en frecuencia. El botn 'Copy' [Copiar] copia la imagen del grfico actual al portapapeles. El cuadro de lista desplegable 'Range (dB)' [Rango (dB)] establece el rango dinmico de la magnitud del grfico (10-100dB). El botn 'Use frequency response compensation' [Usar compensacin de la respuesta en frecuencia] activa/desactiva la compensacin de la respuesta en frecuencia.
o o
ARTA User Manual Para conseguir la copia de la imagen del grfico que hay dentro de la ventana, el usuario ha de pulsar al mismo tiempo las teclas Ctrl+C o activar el comando de men 'Edit->Copy' [Edicin > Copiar], o pulsar el botn correspondiente de 'Copy' [Copiar]. En la barra de herramientas de la ventana principal, el botn 'Copy' [Copiar] aparece con el icono . Este comando abre el cuadro de dilogo 'Copy to Clipboard with Extended Information' [Copiar al portapapeles con informacin ampliada] que se indica en la Ilustracin 1.16. Aqu el usuario ha de configurar tres opciones de copiado:
15
Ilustracin 1.17 Copia del grfico, que aparece en la Ilustracin 1.15. La informacin extendida muestra el nombre del archivo, la fecha, el tiempo y el texto definido por el usuario
16
2 Analizador de espectros
El analizador de espectros de ARTA se presenta como analizador de espectros en base FFT a tiempo real. Un generador incorporado proporciona las siguientes seales: onda sinusoidal, dos rfagas sinusoidales, multitono, ruido blanco, ruido rosa, ruido peridico blanco (PN blanco), ruido peridico rosa (PN rosa) y ruido hablado peridico (PN habla). La forma de trabajo con el analizador de espectros se explicar a travs del procedimiento de pruebas de la tarjeta de sonido. La forma ms sencilla de probar la calidad de la tarjeta de sonido es en el modo analizador de espectros. 1. Haga una conexin en bucle para probar la tarjeta de sonido. 2. Haga clic sobre el elemento del men Mode > Spectrum Analyzer [Modo > Analizador de espectros] o haga clic sobre el icono de la barra de herramientas SPA. 3. Haga clic sobre el elemento del men Generator > Setup [Generador > Configuracin] o haga . Obtendr un grfico como el que aparece en la clic sobre el icono de la barra de herramientas ilustracin 2.1.
Ilustracin 2.1 Cuadro de dilogo para la configuracin del generador de seal El cuadro de dilogo tiene los siguientes controles:
La seccin Signal generator type [Tipo de generador de seales] presenta dos cuadros desplegables. El cuadro desplegable Spectrum mode [Modo de espectro] establece el tipo de generador de corriente para el anlisis de espectros: onda sinusoidal, dos rfagas sinusoidales (para pruebas de intermodulacin), multitono, ruido blanco continuo, ruido rosa continuo y ruido peridico con espectro blanco, rosa y hablado (PN blanco, PN rosa, PN habla). El cuadro desplegable FR mode [Modo de respuesta de frecuencia] establece el tipo de generador de corriente para las mediciones de la respuesta en frecuencia: ruido blanco continuo, ruido rosa continuo y ruido peridico con espectro blanco, rosa y hablado (PN blanco, PN rosa, PN habla). Nota 1: El PN (ruido peridico) es una seal peridica similar a un ruido con un nivel de espectro controlado y una fase aleatoria. El ruido peridico y el multitono pertenecen a la clase de seales multisinusoidales (que se explicarn ms adelante).
17
Ahora elija: Spectrum mode generator type [Tipo de generador de espectros]: Sinusoide Frequency [Frecuencia]: 1000Hz Peak level [Nivel pico]: -3dB Dither level [Nivel de ruido aleatorio]: 16 bit. 4. Usando la barra del dilogo, elegimos:
(frecuencia de muestreo) (numero se muestras en el cuadro de anlisis de FFT) (ventana de seal para evitar prdidas en el anlisis de FFT) (compensacin de la seal)
Los mismos parmetros se pueden configurar en un cuadro de dilogo 'Spectrum Analysis Setup' [Configuracin de anlisis de espectros] que aparece en la Ilustracin 2.2. (se obtiene haciendo clic sobre el elemento de men Setup > Measurement [Configuracin > medicin]). Usando este cuadro de dilogo se puede configurar (1) el canal de entrada preferido, (2) los parmetros de compensacin y (3) la resolucin FFT.
5. Elija: Input channel [Canal de entrada]: Izquierdo. 6. Prepare la mezcladora de sonido de Windows: Active el canal de entrada de lnea Enmudezca el canal de entrada de lnea de la mezcladora de salida. Ajuste el volumen de la lnea de salida a la mxima sensibilidad de salida. Ajuste el volumen de la lnea de entrada cerca de la sensibilidad de entrada mnima. 7. Con el comando de men Setup > Spectrum Scaling [Configuracin > Escala de espectros] ( ), o haciendo clic con el botn derecho del ratn en la zona del nombre del grfico, se abre el cuadro de dilogo 'Spectrum Scaling' [Escala de espectros] (que se muestra en la ilustracin 2.3). Use este cuadro de dilogo para configurar (1) la escala de magnitud, (2) la ponderacin de potencia y (3) las medidas de distorsin.
Ilustracin 2.3. Escala de espectros Seccin Scaling [Graduacin]: Magnitude Scaling [Graduacin de magnitud]: dBFS (dB en la escala completa), dBV o SPL (nivel de presin sonora), PSD (modo de densidad espectral de potencia en dB V/Hz). Voltaje units [Unidades de voltaje]: dBW o dBu. 19
ARTA User Manual dB re 20uPa o dB re 1Pa (vlido slo si se conecta y activa un micrfono).
Seccin Power [Potencia]: Power Weighting combo box [Ponderacin de potencia] Elija: Sin filtro, filtro A, B o C para estimacin de la potencia ponderada de la seal. Show RMS level [Mostrar nivel RMS] active esta casilla para que se muestre el nivel de potencia en la parte inferior del grfico. Seccin Distortion [Distorsin]: THD = distorsiones armnicas totales (en pruebas de respuesta sinusoidal) THD+N = distorsiones armnicas totales + ruido (en pruebas de respuesta sinusoidal) IMD = distorsiones de intermodulacin (en pruebas de respuesta de dos rfagas sinusoidales) Multitone TD+N [Multitono TD + N] active esta casilla para que se muestre la distorsin total + ruido (TD+N) en pruebas de respuesta de multitono. Normalize with full power [Normalizar con plena potencia] active esta casilla para obtener una THD normalizada con potencia de seal que incluya los armnicos superiores. Low cut-off (Hz) combo box [Cuadro combinado de corte bajo (Hz)] establece un corte bajo de frecuencia en las mediciones de THD+N. 2nd and 3rd order IMD [IMD de 2 y 3er orden] active esta casilla para mostrar la distorsin de intermodulacin de 2 y 3er orden, tal como se define en las normativas SMPTE, DIN, CCIF y IEC. Frequency weighting [Ponderacin de frecuencia] active esta casilla para usar la ponderacin de frecuencia (A,B,C) en las mediciones THD+N y TD+N.
8. Active las siguientes casillas: THD, THD+M y Mostrar nivel RMS. 9. Comience la grabacin haciendo clic sobre el icono de la barra de herramientas (o a travs del men Recorder > Run [Grabadora > Iniciar]. Obtendr una respuesta como la que aparece en la ilustracin 2.4. Esta ilustracin se puede obtener con los comandos copiar/pegar (men Edit > Copy [Edicin > Copiar]). Suba lentamente el volumen del canal de entrada de lnea (con el mezclador de la tarjeta de sonido) hasta que llegue a un nivel de pico cercano a -3dB FS.
Ilustracin 2.4 Espectro de un generador sinusoidal de 1 kHz para la tarjeta de sonido Terratec EWX 24/96 en modo bucle. Ventana de seal: Kaiser, tamao FFT: 16384, Fs: 48000Hz. La parte inferior de la Ilustracin 2.4 muestra el valor del espectro en la posicin del cursor (frecuencia y magnitud), el nivel de RMS y las distorsiones. El cursor aparece como una lnea fina que se puede mover pulsando el botn izquierdo del ratn o las teclas de desplazamiento izquierda y derecha en el teclado. 20
ARTA User Manual Si obtiene un THD+N menor del 0,1%, tiene una tarjeta de sonido que podr usar. Si obtiene un THD+N menor del 0,01%, tiene una tarjeta de sonido buena. Nota: Durante la medicin podr usar la barra de control para cambiar el tipo de compensacin, borrar la compensacin, cambiar la frecuencia de muestreo, cambiar el tipo de seal de excitacin y el tamao de FFT. Podr cambiar cualquier parmetro del grfico (rango dinmico, gama de frecuencias y eje) en el cuadro de dilogo 'Spectrum graph setup' [Configuracin del grfico del espectro] (llegar a este cuadro haciendo clic sobre el elemento de men Setup > Graph setup [Configuracin > Configuracin del grfico] o haciendo clic sobre el botn derecho del ratn en la zona del grfico). La forma ms sencilla de ajustar los mrgenes del grfico es usando la barra de control derecha. Las funciones de los botones de la barra se explican en la Ilustracin 2.5.
Ilustracin 2.5 Barra de control para configurar los mrgenes del grfico (tambin se usan para las ventanas de las respuestas en frecuencia) Nota: Las teclas rpidas para cambiar el margen superior de la magnitud del grfico son las teclas "Arriba" y "Abajo" y la rueda de desplazamiento del ratn (el grfico se mover arriba y abajo).
X k = wn xn e j 2kn / N
n =0
N 1
se calcula usando el algoritmo FFT. Este clculo nos da los componentes espectrales en la forma de valores complejos en las frecuencias discontinuas. fk = k f , donde f es una resolucin espectral DFT f = fs / N. Para seales reales, hay N/2 componentes espectrales de potencia de una sola cara Gk:
G0 =| X 0 / N |2
dc component
Gk = 2 | X k / N |2 , k = 1,2,..N / 2 1
4. El espectro de magnitud se expresar en uno de los siguientes modos de graduacin:
Modo de graduacin Nivel pico (referencia a la escala completa) Nivel RMS (espectro de potencia) Densidad espectral de potencia Nivel 10 log(2Gk) 10 log( Gk (sensibilidad_entrada / ganancia_preamp)2 ) 10 log( Gk (sensibilidad_entrada / ganancia_preamp)2 / f) Unidades dBFS dBV (o dBu) dBV/Hz (o dBu/Hz)
Nota: Si se activa la casilla Use Microphone [Usar micrfono] del cuadro de dilogo 'Audio device setup' [Configuracin de dispositivos de audio], los niveles de RMS o PSD se elevan en 20log10( 210-5 Pa) sensibilidad_micrfono(mV/Pa) ). 5. El grfico del espectro muestra los niveles de las magnitudes espectrales como grficos de lneas. Nota: El espectro DFT se define como un conjunto discontinuo de frecuencias, de forma que sera ms apropiado mostrar el espectro como un grfico de barras discontinuo. Sin embargo, cuando tratamos con una cantidad grande de componentes espectrales, como es el caso en ARTA, el grfico de lneas nos proporciona una mejor impresin visual de las magnitudes espectrales. La parte inferior de la Ilustracin 2.4 muestra: RMS Nivel RMS de una seal de entrada definida como 10 log10 (suma de todos los componentes del espectro de potencia DFT) Si la ponderacin de potencia, en el cuadro de dilogo 'Spectrum Scaling' [Graduacin de espectro], est definida como filtro A, B o C, se ponderar cada uno de los componentes espectrales, antes de realizar la suma de espectros, con una respuesta de magnitud de los filtros A, B o C (consulte la definicin de estos filtros en la seccin 2.4). THD distorsin armnica totalse define como un porcentaje de la raz cuadrada del ratio de la suma de potencia de los armnicos superiores (H2, H3,..) respecto de la potencia del armnico fundamental de la seal (H1).
2 2 H 2 + H 32 + .. + H n (%) THD = 100 H 12
22
En el denominador se usa la potencia de la seal distorsionada completa. Esta definicin se encuentra ms cerca del valor de las distorsiones armnicas que se miden con instrumentaciones analgicas en sistemas de ruido bajo. En un sistema con un ruido elevado, una medicin de mayor calidad, equivalente a la instrumentacin analgica, es la de THD*N. THD+N distorsin armnica total ms ruidose define como un porcentaje de la raz cuadrada del ratio de la suma de potencia de los armnicos superiores y la potencia del ruido (N) que tambin incluye distorsin y potencia de ruido.
Con la instrumentacin analgica, se obtiene Potencia Armnica + Potencia del Ruido aplicando un filtro de muesca a la frecuencia fundamental. El valor RMS de la seal y la seal con el armnico fundamental en muesca se miden en algunas bandas de frecuencia predefinidas, generalmente desde un corte de frecuencia bajo (10, 20 0 100 Hz) hasta un corte de frecuencia alto (22, 30 u 80kHz). ARTA no usa la limitacin de alta frecuencia. Se realiza de forma automtica con el filtro antisolapamiento de un conversor AD de entrada. El corte de frecuencia bajo puede ser definido por el usuario. Nota: Si no hay seal en la entrada de la tarjeta, el RMS mostrar la relacin Seal/ruido en el canal de entrada.
YkM
1 M
M i =1
X ki
Exponential averaging [Compensacin exponencial] se suele emplear para hacer el seguimiento de espectros con variacin lenta. Pone relevancia en los acontecimientos recientes, suaviza las variaciones en las altas frecuencias y pone en evidencia las tendencias a largo plazo. ARTA usa un filtro de suavizado que simula un filtro analgico de primer orden y paso bajo con una constante temporal T: 23
M 2 k
= a Y
M 1 2 k
+ (1 a ) X k
La constante a ser igual a T fs /( T fs +1). En ARTA, T est predefinido como 10 veces el tiempo de adquisicin de una secuencia sencilla. Peak hold [Mantenimiento de pico] realmente no se trata de compensacin, sino que simplemente YkM son iguales a los valores mximos de los componentes espectrales,
24
Kaiser
xk I 0 ( x) = k k = 0 2 k!
Tabla 2.1 Definicin de las ventanas de seal empleadas en ARTA Para aplicar la ventana de seal, se recomienda seguir estas normas: Para seales no peridicas continuas (ruido) use la ventana Hanning Para medir las distorsiones armnicas y de intermodulacin use la ventana Kaiser o Blackman4 Para calibrar con una seal sinusoidal, use la ventana de techo plano Para el ruido peridico, los multitonos y otras seales peridicas dentro de la ventana de adquisicin, use la ventana Uniforme Ejercicio: Cambie la ventana de seal y repita las mediciones. Los resultados normales se muestran en la Ilustracin 2.6.
Ilustracin 2.6 Espectro de una seal sinusoidal de 1kHz con las ventanas siguientes aplicadas:
a) Uniforme b) Hanning c) Techo plano
a) Uniforme, b) Hanning, c) Techo Plano, d) Blackman3, e) Blackman4, f) Kaiser (Magnitudes mostradas de 20dB/div).
d) Blackman3
e) Blackman4
f) Kaiser
25
El usuario cambia el color del fondo de "Negro" a "Blanco" haciendo clic sobre el comando del men Edit > B/W background color [Edicin > Color del fondo B/N] o sobre el icono de la barra de herramientas . 26
ARTA User Manual o El usuario indica un color de primer plano arbitrario para cada elemento del grfico usando el cuadro de dilogo 'Edit->Colors and grid style' [Edicin > Colores y estilo de cuadrcula]. As se abre el cuadro de dilogo 'Color Setup' [Configuracin del color sonido] que aparece en la Ilustracin 2.8. Haciendo clic con el botn izquierdo del ratn sobre un rectngulo de color se abre el cuadro de dilogo estndar de Windows 'Color' que se muestra en la Ilustracin 2.9. Nota 1: Si se activa la casilla de verificacin 'All overlays with same color' [Todas las superposiciones con el mismo color], todas las superposiciones se dibujarn con el mismo color. Nota 2: Si se activa la casilla de verificacin 'Dotted graph grid' [Cuadrcula con lnea de puntos], la cuadrcula del grfico se dibujar con lneas de puntos.
Ilustracin 2.8 Cuadros de dilogo para la configuracin de los colores del grfico a) del fondo en negro y b) del fondo en blanco
27
Tabla 2.2 Caractersticas de las ventanas de seal ARTA En los analizadores de bandas de octava el espectro de potencia se mide en una frecuencia fk en una banda de frecuencia que tiene un ancho de banda relativo constante. En un filtro de 1/n octavas, el ancho de banda relativo es igual a 28
Por ejemplo, el filtro de 1/3 octava tiene un ancho de banda que es el 23% de la frecuencia central. En ARTA, las frecuencias centrales se determinan con la expresin:
1
f k = 2 n f k 1 ,
f k 1 = 2
1 n
fk
con 1 kHz como valor de referencia. Esta frmula nos da unos valores cercanos a las frecuencias estndar de la tabla 2. En ARTA, la energa de la banda de octavas se estima aadiendo las potencias espectrales de depsitos de DFT en una banda de octavas 1/n, pero slo si al menos hay tres componentes espectrales DFT que contribuyen a dicha banda. Este requisito quiere decir que el ancho de banda de una banda de 1/n octava ha de ser mayor que el doble del ancho de banda de resolucin de DFT, lo que nos da que la frecuencia ms baja de una banda de 1/n octava es:
f lowest 2
fs N
1 (2
1 2n
1 2n
Por ejemplo, para una frecuencia de muestreo fs = 48000Hz y un nmero de muestras N=16384, la frecuencia menor de los espectros de DFT es igual a 2,93 Hz, la banda de 1/3 octava menor es 25 Hz y la banda de 1/12 octava menor es 100Hz.
16.0 31.5 63.0 125 250 500 1000 2000 4000 8000 16000 20.0 40.0 80.0 160 315 630 1250 2500 5000 10000 20000 25.0 50.0 100 200 400 800 1600 3150 6300 12500 25000
Tabla 2.3. ISO 266 Frecuencias centrales preferidas de bandas de 1/1 octava y 1/3 octava. (La primera columna muestra las frecuencias de banda de 1/1 octava)
29
ARTA User Manual Ilustracin 2.10 Suma de potencias espectrales en un paso de banda de a) filtro de pared de ladrillos, b) filtro de paso de banda Butterworth de 6 polos. En ARTA, se puede calcular power Pk [Potencia Pk] en la banda f1k fk f2k, de dos formas: 1. Potencia de salida de un filtro de paso de banda de pared de ladrillos Primero, se asume que cada componentes DFT consigue una densidad espectral de potencia constante Gn/ f en la regin de frecuencia nf-f/2 f nf +f /2 (de esta forma obtenemos una densidad espectral continua a trozos). A continuacin, se obtiene la potencia de una banda como la integral de la funcin de densidad espectral continua, de f1k a f2k Este proceso se muestra en la Ilustracin 2.10a). 2. Potencia de salida de un filtro de paso de banda Butterworth de 6 polos - Primero, el espectro de potencia se pondera con una magnitud cuadrada de la respuesta del filtro de paso de banda. A continuacin, se calcula Pk como la suma de los componentes espectrales de potencia entre frecuencias en las que la respuesta del filtro es de -20dB. Este proceso se muestra en la Ilustracin 2.10b). El primero mtodo es el que se prefiere en ARTA. La densidad espectral de potencia de la banda kava es Pk /( f2k - f1k).
30
1 RMS = T
T0 +T
x (t )
dt
T0
ARTA usa una constante de integracin T igual a la duracin de un bloque FFT (los ejemplos se muestran en la Tabla 2.4).
FFT - longitud T (fs = 48000Hz) segundos T (fs = 44.100Hz) segundos 4096 0.085 0.093 8192 0.170 0.185 16384 0.341 0.371 32768 0.682 0.743 65636 1.365 1.486 131072 2.730 2.972
Ilustracin 2.4 Duracin del bloque de FFT (para frecuencias de muestreo de 48000 Hz a 44100 Hz)
31
Ilustracin 2.12 Registro de tiempo de la ltima seal capturada La curva muestra la seal correctamente graduada, el valor RMS y el factor de cresta (ratio de valor pico y valor RMS). La lnea amarilla indica la posicin del cursor, y la lnea roja indica la posicin del marcador. Puede cambiar el marcador y la posicin del cursor pulsando y arrastrando los botones izquierdo y derecho del ratn. Haga doble clic sobre el botn derecho del ratn para activar o desactivar el marcador. La etiqueta 'Cursor:' [Cursor:] indica el informe de la magnitud de la seal en la posicin del cursor (tiempo en ms o posicin de la muestra con llaves). La etiqueta 'Gate:' [Limitacin] indica el informe para la diferencia en tiempo (y muestras) entre el cursor y el marcador. Los botones del cristal derecho sirven como comandos para desplazarse por la curva de la seal, para ampliar la imagen, para cambiar la ganancia y la compensacin vertical.
Comandos de zoom: Up [Arriba] aumenta la ampliacin Down [Abajo] reduce la ampliacin All [Todos] muestra todos los puntos Max [Mx] muestra los puntos entre el cursor y el marcador, pero si se desactiva el marcador, ampla el grfico a un nmero de puntos de la pantalla. Comandos de Gain [Ganancia]: Up [Arriba] incrementa el factor de ganancia
32
El cuadro combinado 'Channel' [Canal] muestra el canal que se est empleando actualmente (izquierdo o derecho). Se pueden usar las siguientes teclas de acceso directo:
Teclas arriba y abajo - para cambiar la ganancia Teclas Ctrl+arriba y Ctrl+abajo - para cambiar la compensacin vertical Teclas izquierda y Ctrl+izquierda para desplazar la curva hacia la izquierda. Teclas derecha y Ctrl+derecha para desplazar la curva hacia la derecha Teclas Shift+izquierda y shift+derecha para mover el cursor hacia la izquierda y la derecha Tecla Repag y Avpag - para cambiar el factor de ampliacin
Arrastrando el ratn a la zona de las etiquetas se desplaza la curva horizontal y verticalmente. Haciendo doble clic sobre el botn izquierdo del botn sobre la zona del eje de tiempo cambia la etiqueta de posicin tiempo/muestra.
33
Nota: Los analizadores clsicos de sonido en tiempo real usan la excitacin del espectro rosa para un anlisis de la respuesta del altavoz en banda de octava o en octava suavizada. En un caso ideal (despus de la compensacin) la excitacin del espectro rosa da una respuesta plana en el modo de Espectro de Potencia (dBFS o dBV RMS). Si usamos la excitacin del espectro blanco, el anlisis en banda de octava o en suavizado de octava nos da la respuesta plana en el modo PSD (modo de densidad espectral de potencia). Es importante estudiar las caractersticas de las seales de espectro blanco y rosa, puesto que se usarn para estimar la respuesta en frecuencia y la respuesta de impulso. Ejercicio: Primero, configure: Generador: Graduacin Tamao FFT: Fs (Hz): Ventana: Eje de frecuencia: Compensacin:
Ruido blanco PSD 32768 48000Hz Hanning Suavizado de octava, 1/3 oct, 20Hz-20000Hz Ninguno
Ilustracin 2.13 Densidad espectral con suavizado de octava del generador de ruido blanco (graduacin PSD) Tenga en cuenta que la densidad espectral de una breve secuencia "blanca" no es plana. La ondulacin es muy alta (10 dB). Si repetimos las mediciones con la modalidad de compensacin en Lineal, despus de 20 compensaciones obtenemos el espectro que se muestra en la Ilustracin 2.14. La ondulacin se reduce a 1 dB. Usando 100 compensaciones, la ondulacin se puede reducir a 0,2dB. 34
Ilustracin 2.14 Densidad espectral con suavizado de octava del ruido blanco despus de 20 compensaciones (graduacin PSD) Ahora, si cambiamos a: Generador: Ventana: Ruido peridico blanco (PN Blanco). Uniforme
Ilustracin 2.15 Densidad espectral con suavizado de octava del ruido blanco peridico Es obvio que la seal de ruido peridico tiene una densidad espectral perfectamente plana sin necesidad de compensacin ni de ventanas de seal. Nota: El mismo resultado se puede obtener con una excitacin de ruido rosa (PN rosa) y una graduacin a dBFS o dBV.
35
En ARTA, una seal multisinusoidal de banda ancha se genera usando el DFT inverso:
gn =
donde Ak' =
1 N
A e e
k =0 ' k j
k
N 1
j 2kn / N
0, k = 0 , Ak , k > 0
N k = k | aleatorio [0,2 ] , N / 2 = 0
ARTA introduce tres tipos de ruido peridico: 1. PN Blanco tiene Ak = const, y un espectro blanco 2. PN Rosa - tiene Ak =2A2k y sus magnitudes espectrales se reducen 3dB/oct (despus de cierta frecuencia de corte vase la Ilustracin 2.16). En ARTA, la frecuencia de corte de baja frecuencia variable se puede cambiar en el cuadro de dilogo 'Signal Generator Setup' [Configuracin del generador de seales] El ruido rosa se suele usar en las modalidades de banda de octava o suavizado de 1/n octava con graduacin dBV (vase la Ilustracin 2.15), puesto que proporciona una respuesta de frecuencia lineal en bandas con un ancho de banda relativo constante. 3. PN Habla tiene magnitudes espectrales determinadas por la norma IEC 60268-16 sobre medicin del ndice de Transmisin del Habla (vase la Ilustracin 2.18). Las amplitudes de estas seales tienen una distribucin normal y un factor de cresta de 12-13dB. Si se sigue procesando, el factor de cresta se puede reducir a valores inferiores a 6dB. En ARTA la generacin de ruido peridico se introduce con un factor de cresta inferior a 10dB.
Ilustracin 2.16 Densidad espectral con suavizado de octava del ruido peridico rosa (graduacin PSD)
36
Ilustracin 2.17 Densidad espectral con suavizado de octava del ruido peridico rosa (graduacin dBV)
Ilustracin 2.18 Densidad espectral con suavizado de 1/3 octava del ruido peridico del "habla"
37
Ilustracin 2.19 3. En el cuadro de dilogo Signal generador setup [Configuracin del generador de seal], establezca dos frecuencias sinusoidales en la seccin de Usuario en 11.000 Hz y 12.000 Hz, y un ratio de amplitud de 1:1, tal como se muestra en la Ilustracin 2.20.
Ilustracin 2.21 Distorsiones de intermodulacin de la tarjeta de sonido Terratec EWX 24/96 (medida con dos rfagas sinusoidales: 11kHz y 12kHz) 38
ARTA User Manual El ejemplo de la Ilustracin 2.21 muestra el espectro de intermodulacin para la excitacin de 11kHz y 12 kHz. Observe el componente dominante de intermodulacin a 1 kHz (es la diferencia entre 11 kHz y 12 kHz).
En ARTA, slo se usan los componentes mayores, para 1 < m,n < 8.
Estos factores muestran distorsiones de intermodulacin dominantes cuando f2 >> f1, es decir, para las pruebas de altavoces usaremos f2 = 8.5 f1. Para las mediciones de la distorsin del amplificador, cuando f2 f1, la normativa IEC 60268 3 define dos factores:
39
En esta expresin, las amplitudes de las bandas laterales se suman (en rms) y se expresan como porcentaje del nivel de frecuencia superior. Este factor de intermodulacin est muy cerca del valor de la distorsin de intermodulacin que se puede medir con instrumentacin analgica SMPTE.
40
ARTA User Manual Si f2 / f1 > 7 ARTA usa el mtodo DIN (SMPTE) y genera el informe sobre la distorsin de modulacin: IMDDIN, MD2 y MD3. Si 2 < f2 / f1 < 7 ARTA usa el mtodo de Potencia y genera el informe IMD La relacin de amplitudes, si difieren de las recomendaciones en las normativas, deberan ser incluidos en un informe adicional para el usuario.
ITU_T O.81
Rango de la dcada inferior Diez seales sinusoidales, espaciadas 1/3 octava en una dcada inferior al rango de muestreo. Factor de cresta 10 2dB Rango de la dcada superior Diez seales sinusoidales, espaciadas 1/3 octava en una dcada superior al rango de muestreo. Factor de cresta 10 2dB
Todos los componentes multitonos se generan de forma que cada frecuencia sinusoidal coincide con una frecuencia de depsitos de DFT. Y por esto, a la hora de analizar la respuesta del multitono no se necesita aplicar una ventana de seal. La Ilustracin 2.22 muestra el espectro del multitono del "habla" a travs de un sistema GSM. Observe un alto porcentaje de distorsiones. Nota: El trmino TD+N (distorsin total + ruido) es bastante nuevo. An no est aceptado en ninguna norma. Significa: TD+N es el porcentaje de la raz cuadrada de la relacin entre la potencia de todos los depsitos de ruido y la potencia del multitono.
41
Ilustracin 2.22 Espectro del multitono del "habla" que pasa a travs de dos telfonos mviles de un sistema GSM. Nota: Las pruebas de distorsin con seales multitonos parece que son la nica medida significativa de las distorsiones en sistemas codificados.
Ilustracin 2.23 Barra de estado de las pantallas ARTA El texto de la izquierda muestra el contexto actual del men. El cuadro de texto de la derecha muestra la modalidad de medicin. Cuadros de texto izquierdo y derecho: Muestran los niveles mximos a escala completa de los canales de entrada izquierdo y derecho durante la medicin. Estos "medidores de pico" tienen un rango dinmico de -100 a 0dB. Si se sobrecargan los canales de entrada, los niveles de pico muestran el texto "OVRL".
42
Ilustracin 2.24 Espectro despus de 47 compensaciones (negro) y superposicin que mantiene los valores del espectro despus de 2 compensaciones. La superposicin se puede almacenar en un disco (comando del men Overlay > Save [Superposicin > Guardar]) o cargar de un disco (comando del men Overlay > Load [Superposicin > Cargar]). Las superposiciones se almacenan en archivos de formato binario con la extensin OVS. La superposicin se elimina con el comando del men Overlay > Delete [Superposicin > Borrar]. En la ventana de respuesta en frecuencia se usan las superposiciones casi de la misma forma, pero se almacenan en archivos con extensin .OVF. Nota: Los valores se la curva de espectro actual se han de guardar en archivos de texto (comando del men File > Export ASCII [Archivo > Exportar ASCII]) El archivo almacenado contiene lneas de texto con valores de frecuencia en Hz y magnitud en dB.
43
Ilustracin 3.1. Sistema LTI La relacin entre la entrada y la salida de un sistema LTI, en el dominio de la frecuencia, se puede expresar como:
Y ( f ) = X ( f )H ( f )
en la que la funcin completa H(f) se denomina una respuesta en frecuencia:
H( f ) =
Y( f ) = H ( f ) e j ( f ) X( f )
|H(f)| se denomina una respuesta de magnitud, y (f) una respuesta de fase. La respuesta en frecuencia muestra cmo el sistema cambia el espectro de la magnitud y de la fase de una seal de entrada. La transformada de Fourier inversa de la respuesta de frecuencia se denomina respuesta de impulso. Se escribe como h(t). El producto X(f) H(f) tiene una pareja Fourier en el dominio temporal definida por la circunvolucin x(t)h(t). Esta circunvolucin es igual a la seal de salida y(t):
La funcin h(t) se llama la respuesta de impulso del sistema, puesto que es una respuesta del sistema a una excitacin de la funcin- de un impulso. Esto es obvio, puesto que analizando la circunvolucin (t)h(t), se obtiene:
h(t ) = h( ) (t )d
La respuesta en frecuencia del sistema se suele calcular usando el espectro cruzado de entrada-salida y el espectro automtico de entrada. Al volver a escribir la expresin de la funcin de transferencia de la siguiente forma: 44
H( f ) =
Y ( f ) Y ( f ) X ( f ) S xy ( f ) = = X ( f ) X ( f ) X * ( f ) S xx ( f )
*
obtenemos la respuesta de frecuencia dividiendo el espectro cruzado de entrada-salida con el espectro automtico de entrada (el asterisco indica un valor conjugado complejo) Esta ecuacin se suele denominar el estimador H1. Los pares de transformadas de Fourier del espectro cruzado Sxy(f) y del espectro automtico de entrada Sxx(f) son la correlacin cruzada Rxy(t) y la correlacin automtica (Rxx(t)), es decir.
Rxy (t ) S xy ( f ) Rxx (t ) S xx ( f )
Si la entrada del sistema tiene un espectro blanco ( Sxx(f)=1), entonces Rxx(t)=(t), la respuesta de impulso es igual a la correlacin cruzada de entrada-salida;
Ilustracin 3.2 Diagrama de bloques para el sistema de medicin Nota: En las mediciones acsticas obviamos la influencia del ruido del generador y el ruido del canal de entrada x, puesto que son mucho menores que el ruido y las distorsiones en el canal de salida y. En un sistema de canal doble la entrada al dispositivo de prueba es adquirida por el ordenador como secuencia discreta x. En un sistema de canal sencillo no se mide la seal en la entrada del sistema, y se considera la seal conocida g como excitacin del sistema. En las siguientes secciones trataremos los sistemas de medicin de canal doble y canal sencillo.
H e ( ) =
Y ( f )X ( f )
*
X ( f )X ( f )
* i =1 i i
i =1 N
S xy ( f ) S xx ( f )
( H1 estimator )
Donde He(f) denota la respuesta en frecuencia estimada. Los corchetes <> denotan el valor compensador. El estimador H1 proporciona una previsin polarizada de la funcin de transferencia real H(f), que depende del ruido, las distorsiones y el retardo entre el canal de entrada y salida. Cuando slo el ruido aporta a la polarizacin, el efecto de la compensacin se puede expresar con la ecuacin.
He ( f ) H ( f )+
n N s ( f )A( f )X * ( f ) n X * ( f )X ( f )
* 1 N s ( f )G ( f ) D* ( f ) H ( f )+ , 2 n G ( f )G * ( f ) D( f )
Observe que el trmino seal se aade de forma coherente, mientras que la parte aleatoria del ruido no se suma la potencia. La conclusin es que las compensaciones reducen el nivel de ruido proporcionalmente a una raz cuadrada del nmero de compensaciones, mejorando as la medicin S/N por 10log(n). Si existen distorsiones no lineales, entonces parte del ruido del sistema es coherente con una seal generada, y una mejor medida de la proporcionalidad del ruido+distorsin y un nmero de compensaciones es 1 / n, donde es la funcin de coherencia entrada-salida, definida como;
=
2
| S xy ( f ) |2 S xx ( f ) S yy ( f )
La funcin de coherencia es una medida de la proporcin de la potencia en y que se debe a las operaciones lineales de la seal x. Al estimar la funcin de transferencia, la funcin de coherencia es una comprobacin til de la calidad de los datos empleados. El valor de coherencia mximo es 1. En ARTA se puede ver la coherencia, as que ser posible comprobar la coherencia asociada con las mediciones de "canal doble". En la prctica, necesitamos que 2 est cerca de 1 para garantizar una buena estimacin, pero hemos de tener en cuenta que la coherencia tiene sentido slo si el nmero de compensaciones es mayor de 1.
Ilustracin 3.3 Ilustracin de estimacin sin correlacin en un analizador de Fourier clsico Los principales problemas en el analizador de Fourier clsico con la excitacin de ruido continuo son: La seal de excitacin no tiene un espectro constante. Esto genera una polarizacin del ruido selectivo de frecuencia. Es elevado en las frecuencias en las que el espectro del generador tiene muescas. Esta polarizacin de la resolucin se puede reducir enormemente incrementando el nmero de ciclos de compensaciones. Se recomienda hacer al menos 8 compensaciones de espectro y controlar la funcin de coherencia. 46
ARTA User Manual En un sistema con un retardo grande entre la entrada y la salida (vase la Ilustracin 3.2), es decir, al medir la respuesta del altavoz en la sala, o la respuesta de sistemas de comunicacin con un retardo grande, habr una correlacin baja entre las seales de entrada y salida medidas. En ARTA es posible retardar la adquisicin del canal de entrada, de forma que se elimina este tipo de error. Pero, si medimos la respuesta de frecuencia en un entorno altamente reverberante, no ser posible compensar todos los retardos posibles.
He ( f ) =
1 N
Y ( f )X ( f ) Xi ( f )Xi * ( f ) i =1 i i
N *
Este tipo de compensacin se denomina la compensacin asncrona en el dominio de la frecuencia. Tericamente tiene la misma potencia de reduccin de ruido y distorsiones que el estimador H1, pero se prefiere usar el estimador H1 puesto que nos permite controlar la funcin de coherencia. En ARTA nos referimos a ambos mtodos como la compensacin en el dominio de la frecuencia.
Ilustracin 3.4 Ilustracin de la generacin y adquisicin de seal durante la compensacin asincrnica en el dominio de la frecuencia. Para su aplicacin correcta, se han de cumplir dos condiciones: El inicio de la adquisicin ha de ser posterior al ciclo de pre-compensacin que es necesario para alcanzar la respuesta en estado estacionario. Despus de cada bloque adquirido, se ha de detener la generacin de seal, y se ha de generar la nueva secuencia de ruido peridico. La longitud del bloque FFT ha de ser igual a la longitud de la secuencia generada de ruido peridico. Esto garantiza que las seales generadas y adquiridas estn siempre correlacionadas, as que no habr polarizacin debida al retardo entrada/salida.
Este mtodo de excitacin con el ruido peridico interrumpido es la mejor forma de realizar mediciones de la respuesta de frecuencia en sistemas de comunicacin que se activan con la voz y que tienen un procesado de seal que vara con el tiempo (control de ganancia y reduccin de ruido automtico). El ruido interrumpido mantiene el canal de comunicacin en estado "activo", mientras que se hacen las mediciones en un intervalo pequeo de tiempo para garantizar la estacionalidad del sistema. Si la excitacin se realiza con una secuencia peridica nica, N veces repetidas (vase la Ilustracin 3.5), el estimador puede tener tambin la forma:
47
y (t ) = yi (t ),
i =1
x(t ) = xi (t ), H e ( ) =
i =1
Y ( f )X
(f ) X ( f )X ( f )
* *
Este tipo de compensacin se denomina la compensacin sincrnica en el dominio temporal. Este estimador reduce el ruido aleatorio del sistema, pero no puede reducir las distorsiones no lineales y el ruido estacionario que es peridico dentro del periodo de excitacin.
Ilustracin 3.5 Ilustracin de la generacin y adquisicin de seal durante el proceso de compensacin sincrnica en el dominio temporal. La calidad de la estimacin de la respuesta de frecuencia depende del ruido, tal como se indica a continuacin: 1. La influencia del ruido es elevada en frecuencias cercanas a fs/2, en las que echamos de menos la energa que excitacin, que se filtra con un filtro de suavizado D/A. El problema se soluciona en ARTA con el filtrado digital de la respuesta de impulso estimada en la frecuencia inferior a la frecuencia de corte del filtro antisolapamiento del convertidor. 2. La influencia del ruido es tambin alta en frecuencias extremadamente bajas, puesto que el sistema est acoplado en CA. Esto implica que se han de usar tarjetas de sonido con una frecuencia de corte muy baja. En las mediciones acsticas, el periodo del multisinusoide ha de ser mayor al tiempo de reverberacin - T60. El siguiente razonamiento demuestra este requisito. La respuesta acstica de la sala tiene el ancho de banda de picos de resonancia igual a 2.2/T60.. Si escogemos que la diferencia de frecuencia entre dos componentes multisinusoidales sea inferior a la mitad de este valor, para permitir que se aumenten todas las resonancias de la sala, podremos llegar a la conclusin de que el periodo del ruido peridico haya de ser igual o mayor al tiempo de reverberacin. Adems, tambin se concluye que la longitud del ciclo de pre-compensacin ha de ser mayor o igual al tiempo de reverberacin.
48
Ilustracin 3.6 Diagrama del sistema de medicin de canal sencillo (D conversor digital a analgico, A - conversor de analgico a digital, g - seal discreta generada) La respuesta en frecuencia estimada He es la relacin entre el espectro cruzado YG* y el espectro automtico GG*:
He ( f )
Y ( f )G * ( f ) A( f )N ( f )G * ( f ) = H ( f )A( f )D( f ) + G ( f )G * ( f ) G ( f )G * ( f )
H e ( f ) = H ( f )A( f )D( f )
Vemos que la respuesta en frecuencia estimada est siempre polarizada con funciones de transferencia de los conversores A/D y D/A. No es nada malo, como se podra pensar en un principio, puesto que de esta forma obtenemos una buena relacin Seal/Ruido en frecuencias extremadamente bajas y altas (los conversores A/D y D/A eliminan los espectros en las bajas frecuencias y cerca de fs/2). Esto es especialmente importante si usamos tarjetas de sonido baratas que pueden tener una frecuencia de corte inferior alta. Tericamente, la tcnica de compensaciones de un sistema de canal doble se puede aplicar a un sistema de canal sencillo, pero debido a las limitaciones del controlador de sonido de Windows, que no mantiene la sincronicidad entre la reproduccin y grabacin de la seal, slo se puede aplicar la compensacin sincrnica en el dominio temporal. La respuesta de fase no se puede estimar con precisin absoluta. Para finalizar: El sistema de canal sencillo se recomienda como sistema de medicin cuando usamos tarjetas de sonido de baja calidad. Es apropiado para medir la magnitud de la respuesta de frecuencia, pero no proporciona una respuesta de fase absoluta y precisa. Si deseamos usar el canal del micrfono para la tarjeta de sonido, slo podremos aplicar un sistema de canal sencillo, puesto que la entrada del micrfono de todas las tarjetas de sonido es monofnica.
49
Para hacer una medicin en tiempo real de la respuesta en frecuencia, sern necesarios los siguientes pasos:
50
ARTA User Manual 1. Introduzca la ventana de respuesta en frecuencia de canal sencillo o doble, y conecte el sistema de medicin tal como se indica en el captulo 1. 2. Active el cuadro de dilogo Audio device setup [Configuracin del dispositivo de audio] haciendo clic en el men Setup > Audio devices [Configuracin > Dispositivos de audio] o . Use este cuadro de dilogo para haciendo clic sobre el icono de la barra de herramientas configurar los parmetros de entrada y salida de la tarjeta de sonido, la ganancia de la sonda de voltaje, la ganancia del amplificador de potencia (slo para la modalidad de 1 canal) y la sensibilidad del micrfono (si usa el micrfono). 3. Active el cuadro de dilogo Signal generador setup [Configuracin del generador de seal] haciendo clic en el men Generator > Configure [Generador > Configurar] o haciendo clic sobre el icono de la barra de herramientas . Use este cuadro de dilogo para elegir el tipo de seal de excitacin (ruido blanco, ruido rosa, PN blanco, PN rosa). En la modalidad de canal simple, slo se permite la excitacin con ruido peridico (PN blanco o PN rosa). Si usa la excitacin de PN rosa, establezca la frecuencia de corte inferior. Adems, determine el volumen de salida de la tarjeta de sonido. Nota: La eleccin del tipo de generador se suele hacer usando la barra de control. 4. Active el cuadro de dilogo Frequency Response Measurement setup [Configuracin de medicin de la respuesta en frecuencia] haciendo clic en el men Setup > Measurement [Configuracin > Medicin] o haciendo clic sobre el icono de la barra de herramientas . Este cuadro de dilogo se muestra en la Ilustracin 4.3.
51
5. Active el cuadro de dilogo Frequency response graph setup [Configuracin del grfico de la respuesta en frecuencia] haciendo clic en el men Setup > Graph setup [Configuracin > Configuracin del grfico] o haciendo clic con el ratn en la zona del grfico. Este cuadro de dilogo se muestra en la Ilustracin 4,4. Use este cuadro de dilogo para ajustar (1) el rango dinmico que se muestra, (2) el rango de frecuencias que se muestra y (3) la resolucin del eje de frecuencia.
52
6. Finalmente, podr iniciar la medicin haciendo clic sobre el comando de men Recorder > Run [Grabadora > Iniciar] o haciendo clic sobre el icono de la barra de herramientas . Durante las mediciones podr usar la barra de control para cambiar el tipo de compensacin, borrar la compensacin, cambiar la frecuencia de muestreo, y cambiar el tipo de seal de excitacin y el tamao de FFT. Tambin podr cambiar cualquier parmetro del grfico (rango dinmico, rango de frecuencias y eje) activando el cuadro de dilogo Graph setup [Configuracin del grfico] (haga clic sobre el men Setup > Graph setup [Configuracin > Configuracin del grfico] o haciendo clic con el botn derecho del ratn en la zona del grfico). Las teclas de acceso directo para cambiar el margen superior de la magnitud del grfico son las teclas Arriba y Abajo del teclado y la rueda de desplazamiento del ratn. Con estas teclas se mueve el grfico hacia arriba y hacia abajo. 7. Podr detener las mediciones en cualquier momento haciendo clic sobre Grabadora > Fin o sobre el icono de la barra de herramientas Error! Nombre de archivo no especificado.. La duracin de la medicin depende del tipo de compensacin. Si no se usa compensacin, las mediciones se repiten hasta que el usuario detiene la grabacin. Si se usa compensacin, la medicin se detiene cuando se alcanza el nmero mximo de ciclos de compensacin o cuando el usuario detiene la grabacin.
53
Ilustracin 4.5 Respuesta en frecuencia de Terratec EWX 24/96 en modo bucle Ahora desconecte el canal de entrada izquierdo y mida la respuesta. Lo que mide ahora es un suelo de ruido del canal de entrada. Haga una medicin sin la compensacin y guarde la curva medida como superposicin (comando de men File > Set as overlay [Archivo > Establecer superposicin]). A continuacin repita la medicin con una compensacin = 100. En ambos casos, configure el eje de Frecuencia en la Frequency Response Graph Setup [Configuracin del grfico de respuesta de frecuencia] a Octave-smoothed [Suavizado de octava] con una resolucin 1/3 octave' [1/3 octava]. Obtendr la curva del suelo de ruido que aparece en la Ilustracin 4.6. Observe que el nivel absoluto del suelo de ruido nos da realmente el mximo rango dinmico que se puede obtener con una tarjeta de sonido. Lo que vemos es que la compensacin reduce un suelo de ruido y aumenta un rango dinmico, pero no de igual forma en todas las frecuencias.
54
ARTA User Manual Ilustracin 4.6 Rango dinmico en la modalidad FR1 para la tarjeta de sonido EWX 24/96, sin compensacin (curva superior) y con compensacin (curva inferior). Ahora conecte el canal derecho en el modo bucle y repita la medicin en la modalidad FR2. Compruebe la modalidad Frequency domain averaging [Compensacin en el dominio de la frecuencia] en el cuadro de dilogo Frequency Response Measurement Setup [Configuracin de la medicin de la respuesta en frecuencia]. Obtendr una respuesta en frecuencia casi perfectamente plana y un rango dinmico tal como se muestra en la Ilustracin 4.7.
Ilustracin 4.7 Rango dinmico en la modalidad FR2 para la tarjeta de sonido EWX 24/96, sin compensacin (curva superior) y con compensacin (curva inferior). Vemos que la compensacin en el dominio de la frecuencia, en la modalidad FR2, nos da un rango dinmico ligeramente ms grande que en la modalidad FR1. Esta ganancia en el rango dinmico ser ms pronunciada al medir sistemas que exhiben distorsiones no lineales. La compensacin permite mediciones con un rango dinmico grande, pero slo si usamos la tarjeta de sonido de alta calidad. La Ilustracin 4.8 muestra la comparacin del suelo de ruido para dos tarjetas de sonido (Audio integrado Intel y sistema de sonido de muy alta calidad RME Fireface 800). Durante las mediciones ambas tarjetas estaban en modo de 16 bits, pero la diferencia en el rango dinmico es de 30-50 dB.
Ilustracin 4.8 Rango dinmico en modo FR1 sin compensacin para sonido integrado Intel (curva superior) y RME Fireface 800 (curva inferior).
55
ARTA User Manual Nota: Los anteriores ejemplos muestran curvas con un suavizado 1/3 octava. El suavizado de la curva de la respuesta de frecuencia se hace de forma diferente cuando se realiza el suavizado de la magnitud del espectro. ARTA usa un planteamiento de primero interpolar y compensar la respuesta en el eje de la frecuencia logartmica, a continuacin hace el suavizado mediante la circunvolucin de esa respuesta con la respuesta del filtro Butterworth de seis polos de paso de banda.
Aqu tratamos con dos tipos de superposiciones. El primer tipo se puede fijar a partir de la curva actual usando el comando Overlay > Set overlay' [Superposicin > Establecer superposicin] o cargndolo de un disco con el comando 'Overlay->Load overlay' [Superposicin > Cargar superposicin]. Se puede mostrar slo una superposicin de este tipo. Este tipo de superposicin se puede guardar en el disco usando el comando Overlay > Save overlay' [Superposicin > Guardar superposicin]. Los nombres de los archivos de superposicin tienen la extensin '.ovf'. El segundo tipo de superposiciones son respuesta de filtros objetivo definidos por el usuario. Se pueden generar usando el cuadro de dilogo Target Filter response [Respuesta de filtro objetivo] que aparece en la Ilustracin 4.9. Se activa con el comando 'Overlay > Generate target response [Superposicin > Generar respuesta objetivo].
56
Ilustracin 4.9 El dilogo para generar una superposicin con la respuesta de filtro objetivo del filtro de cruce del altavoz La seleccin de respuestas filtro objetivo se hace a partir del conjunto de las respuestas de filtro de cruce del altavoz. El usuario elige el tipo del filtro de cruce ptimo del cuadro desplegable de la
Un tercer tipo de superposiciones son respuestas arbitrarias definidas por el usuario provenientes de un archivo ASCII (con el comando de men Overlay > Load target response [Superposicin > Cargar respuesta objetivo]). El contenido debera tener el formato .frd.
* se puede usar un asterisco al inicio para designar un comentario opcional * magnitud de frecuencia (dB) 1000 32 1020 33 1040 34 .... ...
57
Ilustracin 4.10 Respuesta en frecuencia del sistema GSM. (fs = 16000Hz, compensacin lineal, tamao del bloque FFT 4096 muestras o 235 ms, generador: PN rosa con frecuencia de corte 200Hz). Haciendo clic sobre el icono de la barra de herramientas IMP (o el men Mode > Impulse response [Modalidad > Respuesta de Impulso]) abre el cuadro de dilogo que aparece en la Ilustracin 4.11. Una confirmacin a la pregunta: Convertir a respuesta de impulso, con el botn 'Yes' [S] se abre la ventana de respuesta de impulso (Ilustracin 4.12) que contiene la vista temporal de la respuesta de impulso.
58
ARTA User Manual Ilustracin 4.11 Cuadro de dilogo para obtener la respuesta de impulso a partir de la respuesta en frecuencia Nota: La conversin a la respuesta en impulso se puede hacer usando el comando de men 'File > Save as PIR' [Archivo > Guardar como PIR]. Este comando transforma la respuesta en frecuencia actual en una respuesta en impulso que la archiva como un archivo .PIR.
Ilustracin 4.12 Ventana de respuesta de impulsos En la ventana de respuesta de impulso podremos modificar la vista de una respuesta de impulso de la misma forma y con los mismos controles que en la ventana de registro de tiempo.
La lnea amarilla indica la posicin del cursor, y la lnea roja indica la posicin del marcador. El Cursor: La etiqueta 'Cursor:' [Cursor:] indica el informe de la magnitud de la seal en la posicin del cursor (tiempo en ms o posicin de la muestra con llaves). La etiqueta 'Gate:' [Limitacin] indica el informe para la diferencia en tiempo (y muestras) entre el cursor y el marcador. De forma alternativa, se puede establecer, en el comando de men View > Gate Time (m @ 344m/s) [Ver > Tiempo de limitacin (m @ 344m/s)] para mostrar la distancia equivalente de la propagacin de la onda sonora. Puede cambiar la posicin del marcador y del cursor pulsando y arrastrando las teclas izquierda y derecha del ratn. Haga doble clic sobre el botn derecho del ratn para activar y desactivar el marcador. Lo mismo se puede hacer con los botones Set [Fijar] y Del [Borrar] del panel derecho. Otros botones sirven para hacer zoom sobre el grfico, para cambiar la ganancia [gain] y la desviacin [offset] vertical. El comando Zoom All [Ampliar todos] muestra todos los puntos de la seal. El comando Zoom Max [Mx Ampliacin] muestra los puntos entre el cursor y el marcador, pero si se desactiva el marcador, ampla el grfico a un nmero de puntos de la pantalla. Se pueden usar las siguientes teclas de acceso directo: Teclas arriba y abajo - para cambiar la ganancia Teclas Ctrl+arriba y Ctrl+abajo - para cambiar la compensacin vertical Teclas izquierda y Ctrl+izquierda para desplazar la curva hacia la izquierda. Teclas derecha y Ctrl+derecha para desplazar la curva hacia la derecha Teclas Shift+izquierda y shift+derecha para mover el cursor hacia la izquierda y la derecha
59
En este ejemplo, el marcador (lnea roja) est en la posicin del mximo. Es decir, a 207,7 ms desde el cursor (lnea amarilla). La etiqueta Limitacin muestra esta diferencia (realmente nos da la cantidad de retardo del sistema). Nota: El cursor est colocado en una posicin de referencia para mediciones de canal doble. Es la posicin de retardo cero. En ARTA, la posicin referente es el ndice de muestra 300. En la modalidad de medicin de canal sencillo, la posicin de retardo cero no est fija ARTA siempre pone en la posicin de referencia la primera muestra cuya amplitud sea 20dB inferior al mximo de la respuesta de impulso.
En la configuracin de mediciones, un parmetro el retardo del sistema - es bastante arbitrario. El retardo est siempre presente en las mediciones acsticas (debido a la propagacin del sonido desde el altavoz al micrfono). Tambin est siempre presente en todos los sistemas procesados digitalmente, especialmente, es bastante grande en un sistema GSM. Si repetimos las mediciones, pero con la casilla de edicin Propagation Delay [Retardo de propagacin] en la ventana FR2 configurado a 204,7 ms, obtenemos la respuesta de impulso que se muestra en la Ilustracin 4.13.
Ilustracin 4.13 Respuesta de impulso medida con una compensacin del retardo de propagacin Ahora, el mximo de la respuesta de impulso estar cercana a la posicin referente de "retardo cero". Quiere decir que las mediciones se realizan con unas seales de entrada y salida mximamente correlacionadas. Es extremadamente importante para estimar el retardo adecuadamente al medir con excitacin con ruido continuo. Nota: En ARTA se pueden guardar la respuesta de impulso como archivo de formato *.PIR usando el comando de men File > Save [Archivo > Guardar] o File > Save as [Archivo > Guardar como]. La descripcin de las operaciones del archivo se incluye en el siguiente captulo. 60
Ilustracin 4.14 Cuadro de dilogo para ver e introducir la informacin del archivo El formato del archivo binario .PIR es el siguiente:
char filesignature[4] unsigned int version; int infosize; al final del archivo int reserved1; int reserved2 float fskHz; int samplerate; int length; int inputdevice; float devicesens; int numchannels; int avgtype; frecuencia) int numavg; mediciones int bfiltered; 2 int genused; float peakleft; derecho float peakleft; derecho float reserved1; float reserved2; // cuatro caracteres de firma: 'P','I','R','\0' // versin de formato de archivo // longitud del texto definido por el usuario // // // // // // // // // 0; 0; frecuencia de muestreo en kHz frecuencia de muestreo en Hz longitud pir 0 - voltprobe, 1-mic, 2-acc V/V o V/Pa (para entrada de micrfono) medicin en canal sencillo o doble tipo de compensacin (0-tiempo o 1-
// nmero de compensaciones empleadas en las // filtro antisolapamiento forzado en el canal // tipo de generador // valor pico (ref. 1.0) en el canal de entrada // valor pico (ref. 1.0) en el canal de entrada // 0 // 0
61
ARTA puede exportar datos de la respuesta de impulso en tres formatos: 1) Archivo .wav de 16 bits de Microsoft (mximo de PIR de 32000) 2) Archivo ASCII MLSSA
// Todas las lneas comienzan con espacio inicial, en el siguiente orden: Cero (0). Intervalo de muestreo en milisegundos. Nmero total de puntos de datos para seguir. Puntos de datos. Cadena de nombre.
3) Archivo ASCII Archivo de texto con lneas de comentario y lneas con dos columnas de datos (tiempo en segundos y amplitud en voltios)
// Algn comentarios // tiempo (segundos) amplitud (V) 0.01 0.234 0.02 0.45 0.03 -0.98
Archivo .wav de 16 bits de Microsoft (mximo de .wav establecido en un valor 1.0 en el archivo .pir) Archivo .tim MLSSA binario Archivo ASCII MLSSA tal como se describe antes Archivo ASCII tal como se describe antes (tambin acepta formatos de tiempo ASCII CLIO)
2. 3. 4.
62
Cada hoja sirve para configurar y hacer la medicin de una respuesta de impulsos. Se debatir ms adelante cundo usar cada tipo concreto de excitacin.
63
Ilustracin 5.1 Cuadro de dilogo para la medicin de la respuesta de impulsos usando la excitacin de ruido peridico
Los controles que se usan para la medicin de la respuesta de impulso con la excitacin del ruido peridico son: Sequence length [Longitud de la secuencia] elija el nmero de muestras en un periodo del ruido peridico Sampling Rate [Frecuencia de muestreo] elija la frecuencia de muestreo, Noise Spectrum [Espectro de ruido] elija la forma espectral del ruido (Rosa, Blanco o Habla), Output Volume [Volumen de salida] elija el volumen de salida en dB, Pink cutoff-Hz [Corte rosa-Hz] introduzca la frecuencia de corte del ruido peridico rosa, Preferred input channel [Canal de entrada preferido] elija el canal de entrada de la tarjeta de sonido que se emplea para medir la salida de un dispositivo a prueba, Dual channel mode [Modo de canal doble] elija el modo de canal sencillo o doble, Invert Phase of input channel [Invertir la fase del canal de entrada] activar esta casilla para cambiar la polaridad de la seal de entrada. Number of averages [Nmero de compensaciones] introduzca el nmero de compensaciones Frequency domain 2ch averaging [Compensacin de canal doble en el dominio de frecuencia] active esta casilla para usar la compensacin en el dominio de la frecuencia Filter dual channel response [Filtrado de la respuesta de canal doble] Active esta casilla para usar el filtro antisolapamiento de la respuesta de impulso (se elimina el ruido cerca de fs/2 en el modo de canal doble), Generate [Generar] iniciar o detener el generador, Record [Grabacin] iniciar o detener la medicin (grabacin y generacin de seal), Default [Defecto] seleccionar la configuracin por defecto OK cierra el cuadro de dilogo y vuelve a la ventana de respuesta de impulsos, que contendr la nueva respuesta de impulsos peridica grabada. Esc - cierra el cuadro de dilogo y vuelve a la ventana de la respuesta de impulsos, sin modificar su contenido.
64
g (t ) = sin(2 (t ))
La frecuencia instantnea de esta seal se define como:
f (t ) =
d (t ) dt
En ARTA se usan dos tipos de ondas sinusoidales de barrido; una con una dependencia de tiempofrecuencia logartmica y otra lineal. Una seal sinusoidal de barrido lineal se define con la siguiente funcin de fase:
(t ) = f1t + ( f 2 f1 )t 2 / 2T
En la que T denota la duracin total del barrido (en segundos), f1 es la frecuencia de inicio y f2 es la frecuencia de parada. Se define una onda sinusoidal de barrido logartmico, conforme a Farina, con la siguiente funcin de fase:
t ln( f 2 / f1 ) f1T T (t ) = (e 1) ln( f 2 / f1 )
En ambos casos, el factor de cresta es 3dB, un valor mucho menor que el factor de cresta del ruido. Quiere decir que la medicin con una onda sinusoidal de barrido nos da una relacin seal/ruido elevada. Cmo usa ARTA la onda sinusoidal de barrido en una medicin de respuesta de impulso? ARTA trata las ondas sinusoidales de barrido como seales no peridicas y las usa como seales de excitacin en un analizador Fourier con un estimador H1. La idea bsica se muestra en la Ilustracin 5,2. Primero, se genera la secuencia de la onda de barrido de longitud N. En el mismo momento, ARTA inicia la adquisicin del bloque de 2N muestras para un anlisis de FFT. Ser importante duplicar la longitud de la secuencia adquirida en una medicin acstica puesto que garantiza que se recogen todos los reflejos de la sala. Un requisito adicional en las mediciones acsticas es que la duracin de la secuencia generada sea mayor que el tiempo de reverberacin de la sala.
Ilustracin 5.2 Principios de la generacin de las ondas sinusoidales de barrido y la adquisicin de seal en ARTA
65
Ilustracin 5.3 Cuadro de dilogo para la medicin de la respuesta de impulsos usando una onda sinusoidal de barrido
Los siguientes controles se usan para la medicin de la respuesta de impulsos usando la excitacin de la onda sinusoidal de barrido: Sequence length [Longitud de la secuencia] elija el nmero de muestras en un periodo de la onda sinusoidal de barrido Sampling Rate [Frecuencia de muestreo] elija la frecuencia de muestreo, Output Volume [Volumen de salida] elija el volumen de salida en dB, Log-frequency sweep [Barrido de frecuencia logartmica] elija la dependencia tiempo-frecuencia logartmica o lineal Generate voice activation [Generar la activacin de voz] seleccione esta casilla para activar una seal especial de tipo voz Ser necesario en un sistema activado por la voz. Preferred input channel [Canal de entrada preferido] elija el canal de entrada de la tarjeta de sonido que se emplea para medir la salida de un dispositivo a prueba, Dual channel measurement mode [Modo de medicin de canal doble] seleccione esta casilla si desea el modo de medicin de canal doble, Invert Phase of input channel [Invertir la fase del canal de entrada] activar esta casilla para cambiar la polaridad de la seal de entrada. Number of averages [Nmero de compensaciones] introduzca el nmero de compensaciones, aunque no se suele usar, Record [Grabacin] iniciar o detener la medicin (grabacin y generacin de seal), Filter dual channel response [Filtrado de la respuesta de canal doble] Active esta casilla para usar el filtro antisolapamiento de la respuesta de impulso (para eliminar el ruido cerca de fs/2 en las mediciones en modo de canal doble), Center impulse response [Centrado de la respuesta de impulsos] seleccione esta casilla para activar la grabacin de la parte de la respuesta de impulso no causal (inicio anterior). Esta parte de la respuesta de impulso (en la estimacin con un barrido logartmico) contiene las imgenes de respuesta de impulsos inducidas por distorsin. Default [Defecto] seleccionar la configuracin por defecto OK cierra el cuadro de dilogo y vuelve a la ventana de respuesta de impulsos, que contendr la nueva respuesta de impulsos peridica medida. Esc - cierra el cuadro de dilogo y vuelve a la ventana de la respuesta de impulsos, sin modificar su contenido.
66
ARTA User Manual Notas: La onda sinusoidal de barrido es una seal de excitacin ptima para la medicin rpida de una respuesta de impulsos acstica sin compensacin. Nos proporciona una estimacin mejor que otras seales de excitacin en entornos acsticos que varan con el tiempo y sistemas ligeramente no lineales. La onda sinusoidal de barrido no una seal de excitacin buena si el entorno genera un gran nivel de ruido coloreado o impulsivo. Tambin nos da una mala estimacin en el sistema que tiene un control de ganancia automtica sensible a la frecuencia o la supresin automtica del ruido. En aquellos casos en los que la excitacin con ruido peridico nos da una estimacin mejor.
b( x ) = x 4 + x 3 + 1
Estos son polinomios de cuarto orden y la secuencia MLS generada tiene una longitud (periodo) N = 24 1 = 15.
Estados de registros de cambio ---------------------a3 a2 a1 a0 ---------------------1 0 0 0 1 1 0 0 1 1 1 0 1 1 1 1 0 1 1 1 1 0 1 1 0 1 0 1 1 0 1 0 1 1 0 1 0 1 1 0 0 0 1 1 1 0 0 1 0 1 0 0 0 0 1 0 0 0 0 1 1 0 0 0
Cada secuencia MLS tiene las siguientes caractersticas: 1. La funcin de autocorrelacin es constante Rkxx = 1 2. El valor CC es igual a 1/N 3. El factor de cresta es 1 (0dB) 67
k =0 k 0
ARTA User Manual Para N mayores, un valor CC (1/N) se acerca a cero. Entonces, la autocorrelacin es igual a 1 para k=1, sino es igual a cero. El espectro de potencia Sn y la autocorrelacin Rk son parejas Fourier:
S nxx = Rkxx e j 2 nk / N = 1
k =0
N 1
Este espectro de potencia es una constante, que significa que la secuencia MLS tiene un espectro blanco. Cuando la excitacin del sistema tiene un espectro blanco la correlacin-cruzada de una seal de salida con una seal de entrada es proporcional a la respuesta de impulsos del sistema ( hk Rkxy ). Una simple generacin de hardware y una computacin de correlacin rpida fueron la razn principal para la popularidad de una instrumentacin de base MLS (la correlacin con una secuencia MLS se puede hacer con la transformada Hadamard que es un algoritmo ms rpido que un FFT). La segunda razn para la popularidad de MLS es una propiedad terica MLS que tiene el factor de cresta ms bajo posible. En realidad, cuando se genera el MLS con una tarjeta de sonido, esto no es verdad, puesto que la seal MLS cambia en la salida del filtro "antisolapamiento" conversor D/A y al pasar a travs de cualquier otro filtro. Un valor de factor de cresta de 6dB a 9dB es comn en la salida de las tarjetas de sonido de PC. El mayor problema con la seal MLS es que algunas de las subsecuencias de MLS estn correlacionadas y pueden generar graves distorsiones al medir la respuesta de sistemas no lineales. Y es por eso que las ondas sinusoidales de barrido y el ruido peridico son seales mejores para medir la respuesta en frecuencia de sistemas que presentan una ligera no linealidad.
Ilustraciones 5.5 Cuadro de dilogo para la medicin de la respuesta de impulsos usando MLS
68
Al usar una seal MLS, slo se permite la configuracin de canal sencillo o de canal semi-doble. En la configuracin del canal semi-doble el segundo canal se usa como referencia de tiempo, pero la estimacin de la respuesta de impulsos se hace como en el sistema de canal sencillo (de esta forma se puede estimar un retardo de entrada/salida de precisin de una muestra, pero la funcin de transferencia estimada est polarizada por las funciones de transferencia de los conversores A/D y D/A). Algunas recomendaciones generales:
El sistema con base MLS es inferior a los sistemas accionados por ondas sinusoidales de barrido o por ruido peridico cuando se implemente con una tarjeta de sonido de PC normal. El sistema con base de ondas sinusoidales de barrido nos proporcionan la mejor estimacin en un entorno de bajo ruido. El sistema de medicin con una excitacin de ruido rosa peridico nos da la estimacin ms robusta, y se puede considerar como sistema de uso general.
69
Ilustracin 5.6 Cuadro de dilogo para la graduacin arbitraria de una respuesta de impulsos En el cuadro de edicin el usuario introduce una constante del punto flotante o una expresin aritmtica compuesta por: Nmeros enteros, Nmeros de coma flotante, Operadores por orden de prioridad: Exponenciacin (^), multiplicacin (*, /), suma (+, -) Corchetes ( ) para agrupado. Ejemplo: La expresin vlida para introducir la escala equivalente de 0,7 dB es la expresin 10^ (0.7/20). La graduacin se suele usar para escalar la respuesta de campo cercano de un altavoz para obtener la estimacin de la respuesta de campo lejano. Ejemplo: La medicin de los altavoces en la respuesta de la sala no nos puede dar una estimacin de la respuesta de bajos en campo abierto. Para obtenerlo medimos la respuesta de bajos en el campo cercano de los altavoces pNF y escalar la respuesta con un factor que es igual a la relacin del radio de la membrana con el doble de la distancia normal de medicin r. a p (r ) = p NF 2r Las respuestas obtenidas de esta forma son una buena aproximacin de la respuesta de campo abierto del altavoz montado en una caja infinita (2), para frecuencias por debajo de 200Hz.
70
Estos anlisis se pueden activar en el submen Analysis [Anlisis] o haciendo clic con el botn izquierdo del ratn sobre uno de los iconos de la barra de herramientas "verde" que se muestra en la Ilustracin 6.1.
Ilustracin 6.1 Ventana de la respuesta de impulso iconos de la barra de herramientas Las siguientes secciones describen los mtodos de anlisis de la respuesta de impulsos.
ARTA User Manual o activar el cuadro de dilogo 'Impulse Response Analysis Setup' [Configuracin del anlisis de la respuesta de impulsos]. Las partes limitadas y no limitadas de la respuesta de impulsos se dibujan con diferentes colores.
Ilustracin 6.2 Componentes que determinen la parte "limitada" de la respuesta de impulsos. El cursor, el marcador y la longitud del bloque de FFT.
72
Ilustracin 6.4 Cuadro de dilogo para la configuracin del anlisis de la respuesta de impulsos
Controles de la seccin Gated time window [Ventana de limitacin temporal]: Maximum length (FFT length) [Longitud mxma (longitud FFT)] - elija la longitud (tamao) del bloque FFT. Warn if FFT length < gate [Aviso si la longitud FFT < limitacin] active esta casilla si desea que se le avise cuando la longitud FFT es menor que la limitacin Window type [Tipo de ventana] seleccione el tipo de ventana: Uniforme, Hann12%, Hann25% o Hann 50%. Seccin ETC window [Ventana ETC]: Type [Tipo] elija la ventana para la estimacin del ETC (envolvente de la respuesta de impulsos: Uniforme, Half-Hann, Habla, Causal Seccin de Delay for phase estimation [Retardo para la estimacin de fase]: Time (ms) [Tiempo (ms)] establecer el intervalo de tiempo desde la posicin del cursor hasta la posicin de referencia para estimacin de la fase (normalmente es una distancia desde la posicin del cursor hasta la posicin del mximo de la respuesta de impulsos).
Antes de realizar la estimacin de la respuesta de frecuencia limitada a partir de una respuesta de impulsos limitada en tiempo, ha de establecer el tipo de ventana que se le aplicar a la respuesta limitada. Podr aplicar las siguientes ventanas: Uniforme, Hann12%, Hann25% o Hann 50%. El porcentaje que aparece detrs del nombre Hann quiere decir que se aplica la mitad de la ventana Hanning en el porcentaje de la parte final de la ventana para atenuar levemente la respuesta de impulsos. En la Ilustracin 6.5 se muestra un ejemplo de la ventana Hann12%.
Ilustracin 6.5 Ejemplo de una respuesta de impulsos de un altavoz en una sala pequea (adems, tambin se muestra la ventana Hann12% - determina la parte de la respuesta de impulsos que se usar para estimar la respuesta de frecuencia limitada)
73
Ilustracin 6.6 Ventana de respuesta de frecuencia suavizada (La barra horizontal, dibujada en color de cursor, denota la regin de la frecuencia en la que el requisito de ancho de bandatiempo no se satisface). La Ilustracin 6.6 muestra la Ventana de respuesta de frecuencia suavizado. La barra
derecha contiene controles para la configuracin del margen del grfico (los mismos controles que estn en las ventanas de respuesta de espectro y de frecuencia). Se puede usar un cuadro desplegable adicional Smoothing [Suavizado] para establecer la resolucin de suavizado. La barra inferior tiene varios botones:
Botn Mag elija el grfico de magnitud Botn M+P elija el grfico de magnitud + fase Botn Ph elija el grfico de fase Botn Gd elija el grfico de retardo de grupo Botn Record para repetir la medicin, Botn Overlay para gestionar las superposiciones, Botn Copy para copiar el grfico al portapapeles, Botn B/W para cambiar el color de fondo.
La misma presentacin se encuentra en la ventana de respuesta de frecuencia DFT (pero sin los mens y controles para Suavizado y Superposiciones). La ilustracin 6.7 muestra ejemplos de las respuesta de frecuencia no suavizados y con suavizado de 1/3 octava. La comparacin de las dos curvas muestra que la respuesta suavizada nos da una tendencia general y realmente la mejor perspectiva de la respuesta de frecuencia. Las grandes ondulaciones en la respuesta de frecuencia no suavizada (DFT) son la consecuencia de los reflejos de la sala. 74
a)
b)
Ilustracin 6.7 Ejemplos de a) respuesta de frecuencia no suavizados y b) con suavizado de 1/3 octava. ARTA ofrece suavizado de 1/1, 1/2, 1/3, 1/6, 1/12 y 1/24 octava. Los filtros de suavizado se describen en el captulo 2.3.
Ilustracin 6.8 Cuadro de dilogo para la configuracin del grfico. La configuracin de los mrgenes del grfico y de la presentacin se podr ajustar con la barra de control derecha y con el dilogo de configuracin del grfico, que aparece en la Ilustracin 6.8. (comando de men View > Setup [Ver > Configuracin]. Los controles son los mismos que en las ventanas de respuesta de frecuencia y de espectro, adems de eso, el usuario puede ajustar el rango para el grfico de fase y de retardo de grupo, y usar las dos casillas de verificacin adicionales.
Thick plot lines [Lneas gruesas] permite dibujar curvas con un trazo grueso,
75
La aplicacin de la limitacin, de duracin T, influye en la resolucin de frecuencia. Si se activa la casilla Time-Bandwidth [Ancho de banda tiempo] (o men View > Time-Bandwidth requirement [Ver > Requisitos ancho de banda tiempo]), ARTa slo mostrar los componentes de la respuesta de frecuencia para las frecuencias en las que fk > 1/T. Si no se activa esta opcin, entonces ARTA mostrar todos los depsitos de frecuencia DFT, y la barra horizontal de la parte inferior del grfico indicar la regin de la frecuencia en la que el requisito de ancho de banda-tiempo no se ha satisfecho (Ilustracin 6.6). El Apndice nos da una descripcin completa de los mens en la ventana de Respuesta de frecuencia suavizada y la ventana Respuesta de frecuencia DFT.
Una definicin sencilla de fase mnima es: Una caracterstica de fase del sistema para la que se puede realizar un sistema equivalente con las mismas caractersticas de magnitud y un cambio de fase mnimo (sobre todas las frecuencias). La diferencia entre la fase y las caractersticas de la fase mnima se suele llamar exceso de fase. Matemticamente, la fase mnima se puede estimar a partir de la magnitud de la respuesta de frecuencia usando la transformada de Hilbert. ARTA, as como otros programas similares, usa el DFT para calcular la transformada de Hilbert. Introduce la periodicidad en la estimacin de la fase mnima y nos da un resultado que se acerca a la fase mnima verdadera slo en las frecuencias inferiores a fs/4. Para definir el retardo de grupo, primero analizamos la respuesta de frecuencia H (j) de la siguiente forma:
H ( j ) = A( )e j ( ) = e ( )+ j ( )
Donde: A()=H(j) es la respuesta de magnitud, () es la respuesta de fase y () es el logaritmo de la respuesta de magnitud. Si asumimos que un sistema ideal tiene una respuesta de frecuencia H(j) = const., la distorsin lineal de la respuesta de frecuencia se puede expresar como:
76
Esta expresin muestra que hay dos tipos de distorsiones lineales debido a los cambios de magnitud y fase. La parte real de esta expresin se define como un retardo de grupo Td:
Td = F{t h(t )} d ( ) = Re d ( ) F{h(t )}
Lo que muestra esta expresin es que el retardo de grupo se puede estimar a partir de la respuesta de impulso del sistema. Para obtener una estimacin aceptable del retardo de grupo, se han de cumplir dos condiciones: 1. La respuesta de impulsos ha de deteriorarse completamente dentro de la regin limitada, 2. El retardo de grupo ha de ser menor de la mitad del tamao de FFT. Si no se cumple esta condicin, obtenemos el valor negativo del retardo de grupo. A veces es til analizar el exceso de retardo de grupo. Es el retardo de grupo obtenido del exceso de fase.
El comando de men Generar respuesta objetivo sirve para generar superposiciones con las caractersticas ptimas de las respuestas de filtros de cruce, y el comando de men Cargar respuesta objetivo sirve para cargar el objetivo arbitrario que se encuentra en un archivo ASCII (formato .frd).
77
Ilustracin 6.9 Cuadros de dilogo 'FR Overlay Manager' [Gestor de superposiciones de respuestas en frecuencia] y 'Overlay colors' [Colores de la superposicin] Otras funciones avanzadas con superposiciones se podrn hacer usando el cuadro de dilogo 'FR Overlay Manager' [Gestor de superposiciones de respuestas en frecuencia] (que se muestra en la Ilustracin 6.9). Se activa con el comando 'Overlay > Manage Overlays' [Superposicin > Gestor de superposicin]. Algunos botones del cuadro de dilogo duplican los comandos de men:
Add [Aadir] establece la curva actual como superposicin Add above crs [Aadir encima del cursor] guarda la parte de la curva actual por debajo del cursor como superposicin Add below crs [Aadir debajo del cursor] guarda la parte de la curva actual por encima del cursor como superposicin Delete all [Borrar todo] Borra todas las superposiciones
Otros botones permiten realizar funciones avanzadas en los elementos del men desplegable seleccionados:
Replace sel [Sustituir la seleccin] sustituye la superposicin seleccionada por la curva activa actual Delete sel [Borrar la seleccin] Borra todas las superposiciones seleccionadas Color cambia el color de los elementos seleccionados usando el cuadro de dilogo 'Overlay Colors' [Colores de superposicin] que se muestra en la Ilustracin 6.10.
Haga clic con el ratn sobre el elemento de la lista desplegable y suceder lo siguiente:
Un solo clic - selecciona los elementos Un solo clic sobre la casilla de verificacin hace que sea visible o invisible la superposicin Doble clic - permite editar el nombre de la superposicin
Todos los elementos de la lista se pueden hacer visibles pulsando sobre el botn 'Check All' [Activar todo].
78
Los elementos 'Corte por debajo del cursor', 'Corte por encima del cursor' y 'Nivel de escala' se usan normalmente para combinar dos grficos, uno para la respuesta de frecuencia alta y otro para la respuesta de frecuencia baja. Lo mismo se puede hacer con el comando Merge' [Fusionar]. En este caso, la curva resultante se puede exportar en formato de archivo ASCII (Magnitud + Fase). El elemento de men 'Nivel de escala' abre un cuadro de dilogo (Ilustracin 6,10) en el que el usuario introduce un nivel arbitrario (en dB) para escalar la respuesta de magnitud. Esta operacin no cambia la respuesta de impulso. Slo cambia la respuesta en frecuencia que se est mostrando actualmente.
Ilustracin 6.10 Cuadro de dilogo 'Magnitud de Escala' Los elementos de men Subtract overlay [Resta superposicin] y Subtract from overlay [Resta de la superposicin] se puede usar para obtener la diferencia entre dos respuesta (es decir, la calibracin de la respuesta de micrfonos con la otra calibrada). Despus de estas operaciones la curva de
ARTA User Manual (esfrica, cuadrada o rectangular), el ancho y la altura de la caja del altavoz. Estos valores se usan para definir la funcin de transferencia de graduacin W(f) que se usa para estimar la respuesta en campo abierto para la respuesta del altavoz que se mide montado en una caja infinita (en el campo cercano).
Ilustracin 6.11 Cuadro de dilogo 'Difraccin de caja en bajas frecuencias' ARTA usa la siguiente expresin para la funcin de transferencia de graduacin de difraccin en bajas frecuencias.
W( f ) =
1 + j f / f0 2 + j f / f0
siendo f0 = 42.7 / d para una caja esfrica de dimetro d, o f0 = 34.16 / d para una caja cuadrada con un ancho d. Estos valores se obtienen aplicando numricamente la funcin de transferencia W(f) con una funcin de transferencia de una caja de altavoz esfrica. Esta funcin de transferencia se llama tambin ecualizador 2/4 puesto que da la diferencia entre la respuesta de un altavoz de baja frecuencia a media esfera (2) y la respuesta en el espacio completo (4). Para una caja rectangular, con un ancho de caja w y alto h, ARTA usa como aproximacin una caja cuadrada equivalente de ancho d = w (h/w)1/3. La Ilustracin 6,12 muestra un ejemplo de la respuesta de altavoz de campo cercano medida (curva superior) y la respuesta en campo abierto estimada (curva inferior). A unas frecuencias muy bajas, la diferencia de nivel es de 6dB.
Ilustracin 6.12 Respuesta de altavoz en campo cercano (curva superior) y la respuesta ecualizada 2/4 (curva inferior). 80
81
ventana haciendo clic sobre Analysis > Step Response [Anlisis > Respuesta de etapa] o sobre el icono de la barra de herramientas .
La Ilustracin 6.13 muestra un ejemplo de la respuesta de etapa en un sistema GSM. La oscilacin transitoria con el periodo 3.5 ms es habitual para el filtro de entrada del canal de voz a 300 Hz.
Ilustracin 6.13 Ejemplo de la respuesta de etapa en el sistema GSM Los controles del panel derecho son los mismos que en las ventanas de la respuesta de impulsos o de registro de tiempo. Los comandos de men permiten almacenar la respuesta de etapa en un archivo de texto (File > Export ASCII [Archivo > Exportar a ASCII]), guardar el mapa de bits del grfico al portapapeles (Edit > Copy [Edicin > Copiar]) y cambiar el color del fondo del grfico (Edit > B/W background color [Edicin > Color de fondo B/N]). Podr determinar o borrar el marcador con un clic sencillo o doble del botn derecho del ratn.
82
t d = x ( t ) t .
x ( )
La transformada de Hilbert es una circunvolucin de x(t) y 1/(t). Una forma ms simple de obtenerlo es con la multiplicacin en el dominio de la frecuencia:
~(t ) = F 1{ j sgn( ) X ( j )} = F 1{ j sgn( ) F{x (t )}} , x Donde: 1 za > 0 1 sgn( ) = 0 za = 0 , F = j sgn( ) - 1 za < 0 t
En ARTA, obtenemos la envolvente de la respuesta de impulso en una ventana independiente haciendo clic sobre el elemento de men Analysis > ETC [Anlisis > ETC] o haciendo clic sobre el icono de la barra de herramientas . La Ilustracin 6.14 muestra el ETC de la respuesta del altavoz. El cursor se coloca en la posicin del primer sonido directo y el marcador se encuentra en la posicin del primer reflejo fuerte.
Ilustracin 6.14 Envolvente de la respuesta de impulso (ETC) de una respuesta de altavoz La parte derecha de la ventana tiene los siguientes controles:
dB range [Gama dB] - establece la gama de magnitud
83
Notas: La mayora de los sistemas de medicin usan la notacin ETC (Curva Energa Tiempo), pero preferimos el nombre de envolvente de respuesta de impulso. El nombre ETC lo puso Richard Heyser, que observ que la ortogonal de los componentes de la seal analstica es anloga al intercambio de la energa potencial y cintica en las ondas acsticas. Su conclusin no tiene justificacin terica puesto que los componentes analsticos generan la funcin no causal mientras que la energa ha de ser una funcin causal. Para obtener una curva ETC ms cercana a la funcin causal, ARTA podr aplicar las siguientes ventanas a los datos del dominio de la frecuencia: Ventana Half-Hann de uso general (ligera supresin de frecuencia extremadamente bajas y altas), Ventana de habla se usa para mejorar la gama hablada (suprime las frecuencias altas y bajas), Ventana causal se usa si la respuesta de impulsos se obtiene con la excitacin MLS. Uniforme se usa sin ventana
C (t , f ) = h( )u ( t )e j2f d
Donde h(t) es la funcin de respuesta de impulso y u(t) es la funcin de etapa de la unidad Realmente C(t, f) es una transformada de Fourier de la parte de la respuesta de impulse definida a partir del tiempo =t hasta el infinito, tal como se muestra en la Ilustracin 6.15. Una aplicacin repetida de la Transformada de Fourier, cada vez para una parte de la respuesta de impulso que est por delante en el tiempo de un intervalo dt, obtendremos la funcin frecuencia-tiempo como en la Ilustracin 6.16. Para estimar la Transformada de Fourier, ARTA usa el DFT y sustituye la funcin de etapa de la unidad por una funcin de ventana de longitud finita. Para evitar el corte abrupto de la respuesta de impulso con la ventana rectangular, ARTA aplica una ventana que sube y baja gradualmente (el usuario elige el tiempo de subida de 0,02 a 1ms).
84
Ilustracin 6.16 Deterioro espectral acumulado en el sistema GSM ARTA muestra el Deterioro Espectral Acumulado en una ventana independiente, pero primero se han de establecer algunos parmetros para la transformada tiempo-frecuencia. El comando de men 'Analysis-Cumulative spectral decay [Anlisis-Deterioro espectral acumulado] abre el cuadro de dilogo Configuracin del deterioro espectral acumulado (Ilustracin 6.17). El Deterioro espectral acumulado aparece despus de pulsar el botn OK.
85
ARTA User Manual Ilustracin 6.17 Cuadro de dilogo de configuracin del Deterioro espectral acumulado
El espectro de deterioro acumulado puede dibujarse como un grfico en cascada o como un sonograma. Las curvas del grfico en cascada se dibujan con distintos colores dependiendo de la magnitud. El sonograma es un grfico de una funcin cuya magnitud est expresada por la cantidad y orden predefinido de colores, lo que se denomina paleta (Ilustracin 6.19). El uso de una paleta con una pequea cantidad de colores discretos mejora las curvas en el plano de color (Ilustracin 6.20). La Ilustracin 6.18 muestra un grfico de cascada coloreada de CSD obtenida de la respuesta de impulsos de los auriculares AKG 240.
La ventana CSD (Ilustracin 6.18) dispone de diversos controles para realizar la manipulacin del grfico. Caja combinada Mode [Modo] elija entre tipo de grfico de cascada o sonograma. Caja combinada Palette [Paleta] elija entre varias paletas de colores (Azabache, Gris, Cobre y Fro). Colored Waterfall [Cascada coloreada] active esta casilla para elegir entre un grfico de cascada coloreado (o son un solo color). Contoured [Contorneado] active esta casilla para elegir entre un cambio de color gradual o contorneado (vea las Ilustraciones 6.19 y 6.20). Control de efecto Rango (dB) - cambia el rango dinmico del rango de 5 a 70dB. Botn Copy [Copiar] copia el grfico actual al portapapeles Botn B/W [Blanco y negro] cambia el color de fondo del grfico entre blanco y negro.
87
Smbolo de cantidad Tiempo de reverberacin T60, T30, T20 (s) Tiempo de deterioro temprano EDT (s) Tiempo central TS (s) claridad C80 (dB) definicin D50 (%) Fuerza G (dB) Fraccin de energa lateral LF, LEF (%) Correlacin cruzada interaural IACC
Aspecto subjetivo Reverberacin, volumen, implicacin Reverberacin, claridad claridad claridad Definicin del habla Nivel relativo Impresin espacial Impresin espacial
Tabla 7.1 Resumen de los parmetros acsticos generales de una sala (ISO 3382) El parmetro ms importante de la sala es el tiempo de reverberacin - T 60 Se define como el intervalo de tiempo necesario para que la energa sonora se deteriore 60dB una vez que se ha detenido la excitacin. Para obtener el tiempo de reverberacin, necesitamos medir o estimar la curva de deterioro de energa una vez que se haya desconectado la fuente de sonido. La curva de deterioro de energa es irregular y ruidosa r(t) que solemos aproximar con el deterioro lineal, tal como se muestra en la Ilustracin 7.1.
Ilustracin 7.1 Curva de deterioro de energa para una fuente de sonido que se desconecta a t=0. El tiempo de reverberacin se determina con la pendiente del deterioro lineal, de la siguiente forma:
T60 = 60
ISO 3382 define:
dt dr
88
ARTA User Manual T20 es el tiempo de reverberacin determinado a partir de la pendiente media de la curva de deterioro de energa obtenida a partir de parte de la curva de deterioro entre -5dB y -25dB.
T30 es el tiempo de reverberacin determinado a partir de la pendiente media de la curva de deterioro de energa obtenida a partir de parte de la curva de deterioro entre -5dB y -35dB.
r (t ) h 2 ( )d
t
Tenga en cuenta que en esta expresin, el denominador representa la energa total. En ARTA, obtenemos la curva de deterioro de energa haciendo clic sobre el elemento de men Analysis > Energy Decay' [Anlisis > Deterioro de energa] o haciendo clic sobre el icono de la barra de herramientas . La Ilustracin 7.2 muestra la ventana de deterioro de energa con dos curvas. La curva superior es la curva de deterioro de energa, y la curva inferior gris es la envolvente de la respuesta de impulsos. ISO 3383 recomienda el anlisis del deterioro de la energa en seis bandas 1/1 octava (125 Hz, 250Hz, 500Hz, 1kHz, 2kHz y 4kHz). ARTA permite medir los parmetros acsticos en dos bandas adicionales: 62,5Hz y 8kHz. Para obtener esto, ARTA filtra la respuesta de impulsos con filtros de bandas de octava de clase I IEC. La misma norma recomienda que durante las mediciones: La posicin del micrfono ha de estar al menos a 1 m de las superficies reflectantes y no demasiado cerca de ninguna fuente; la distancia mnima desde la fuente se puede calcular con la ecuacin:
d min = 2
V [m] cT
Donde: V es el volumen de la sala, c la velocidad del sonido, T es una estimacin del tiempo de reverberacin esperado. La fuente de sonido debera ser tan omni-direccional como sea posible. El micrfono ha de ser omni-direccional.
Tambin es importante que el ciclo pre-compensacin de las mediciones con un ruido peridico o con MLS sea mayor al tiempo de reverberacin.
89
Ilustracin 7.1 Ventana de deterioro de energa La curva superior es la curva de deterioro de energa, y la curva inferior es la envolvente de la respuesta de impulsos.
Los controles de la derecha se usan para configurar el grfico y estimar los parmetros acsticos de la sala, de la siguiente forma:
Filtering [Filtrado] elija la frecuencia de una respuesta filtrada en bandas de octava o la respuesta de banda ancha. dB range [Gama dB] introduce la gama de magnitud T60 inicia los clculos de los parmetros acsticos, el valor de los cuales se indica en la parte inferior de la ventana. Noise Tail [Cola de ruido] presenta dos controles. En el primer control se elige el porcentaje de la curva de deterioro que se estima como ruido. El otro control se usa para elegir el mtodo de reduccin de ruido: Trunc quiere decir que no se incluir el porcentaje de la cola de la curva en la estimacin del deterioro. Sub quiere decir que el nivel medio de potencia de ruido de cola se resta de la curva de deterioro. Log [Registro] se obtiene un informe de los parmetros acsticos de la sala que se han calculado. Zoom - establece el factor de zoom horizontal a Mx o Todos. Scroll [Desplazamiento] mueve el grfico de izquierda o derecha
ARTA necesita cierta interaccin por parte del usuario para estimar el deterioro de energa y los parmetros acsticos adecuadamente. El usuario ha de estimar el porcentaje de la cola de la curva ETC en la que el ruido es el principal contribuyente, y el mtodo para reducir ese ruido. El mtodo ms simple para reducir el ruido es truncar la respuesta de impulso. ISO 3382 sugiere definir un punto de truncamiento como el punto en el que el nivel de la seal + ruido es 5dB superior al nivel medio de ruido de la cola de la respuesta de impulsos. En ARTA el punto de truncamiento se determina como el porcentaje de la cola de la respuesta de impulsos que est detrs de dicho punto de truncamiento. Hay dos mtodos para reducir el ruido: 90
ARTA User Manual Mtodo Trunc trunca (elimina) la cola del ruido de la estimacin del deterioro de energa
Mtodo Sub resta el valor medio de potencia de ruido de cola con la integracin de
Schroeder. La Ilustracin 7.3 muestra el caso en el que se trunca el 40% de la cola de ruido de la estimacin del deterioro de energa. Un criterio que se puede usar para validar la estimacin es: en una sala sencilla no acoplada, el deterioro de energa logartmico debera tener un deterioro lineal.
Ilustracin 7.3 Ventana de deterioro de energa La curva superior es la curva de deterioro de energa con el 40% de la cola de ruido eliminada (truncada) de la estimacin del deterioro, y la curva inferior gris es la envolvente de la respuesta de impulsos.
Los parmetros acsticos se estiman a partir de la curva de deterioro de energa. El usuario ha de elegir una parte ms o menos lineal de la curva de deterioro y fijar el cursor y el marcador en esa parte de la curva. Segn ISO 3382, el cursor debera estar cerca del nivel -5dB, mientras que el marcador debera estar detrs de la posicin en la que el nivel es 30dB por debajo del nivel del cursor (o al menos 20dB). Si el nivel de deterioro en la posicin del marcador es inferior a 25dB, ARTA estima T20, y si el nivel de deterioro en esa posicin es inferior a -35dB, ARTA tambin estima T30. Haga clic sobre el botn 'T60' para obtener la estimacin del tiempo de reverberacin y otros parmetros acsticos. Para obtener el tiempo de reverberacin la curva de deterioro se aproxima con una ecuacin lineal y = ax + b y se estima la pendiente de la curva mediante la regresin lineal. La curva a que mejor se ajusta y los valores b polarizados se evalan con el mtodo de los mnimos cuadrados. La calidad de la estimacin ajustada a la lnea que se produce con la regresin lineal se describe con el coeficiente de 91
ARTA User Manual correlacin r. El coeficiente de correlacin tiene el rango [-1,1] en la que una alta correlacin produce unos valores cercanos a la unidad.
El mismo procedimiento se debera repetir para cada respuesta filtrada en bandas de octava y para la respuesta en banda ancha. ARTA memoriza los valores de los parmetros acsticos estimados, de forma que cuando se hace clic sobre el botn 'Log' aparecer el informe de los valores de los parmetros acsticos en todas las bandas de octavas.
EDT Tiempo de deterioro temprano (EDT) se define como el intervalo de tiempo necesario para que la energa sonora se deteriore 10dB una vez que se ha detenido la excitacin. Para permitir una comparacin directa con el tiempo de reverberacin, el resultado se multiplica por 6 El deterioro exponencial idneo en un campo difuso, el valor de EDT esperado sera igual a T60.
92
Claridad La claridad o "ndice temprano a tarde" C80 (C50) se define como la relacin logartmica de la energa de la respuesta de impulso antes del tiempo te, y la energa despus de te. El valor te = 50ms se usa para expresar la claridad del habla, mientras que te = 80ms ser ms apropiado para la msica.
80 ms
C80 = 10 log
p p
50 ms 2
p p
(t )dt dB (t )dt
80 ms
50 ms
El nombre original en alemn para el trmino claridad es "klarheitsmass". Unos altos valores de claridad indican una gran cantidad de energa temprana, que corresponde con una sensacin subjetiva de claridad. Por el contrario, unos valores de claridad bajos indican que el sonido no es claro, y es excesivamente reverberante.
Definicin La Definicin D50 o "relacin de energa sonora temprana a total" es una medida de la definicin del habla. Tambin se conoce por su nombre alemn, Deutlichkeit. Se define como:
50 ms
D50 = 100
(t )dt (%)
h
0
(t )dt
TS-Tiempo central El tiempo central TS se corresponde con el centro de gravedad del cuadrado de la respuesta de impulso:
TS =
t h
0
(t )dt
h
0
(t )dt
Los lmites superiores de integracin se consideran los puntos de truncamiento, o el final de la respuesta de impulsos, dependiendo de la opcin de tratamiento de ruido especfica. El subndice S en el nombre TS es la abreviatura del nombre alemn Schwerpunktzeit. El valor de TS se expresa en milisegundos. Un TS bajo sugiere una sensacin de claridad, mientras que un valor TS alto sugiere un sonido reverberante. TS est muy correlacionado con EDT (rara vez contiene alguna informacin adicional si lo comparamos con EDT).
93
m( F ) =
h (t )e
2 0 2 0
j 2 F t
dt
h (t )dt
En la que F es la frecuencia de modulacin. A partir de esta definicin, MTF es proporcional a la transformada de Fourier del cuadrado de la respuesta de impulsos (el denominador, que representa la energa total, normaliza la expresin). Puesto 94
ARTA User Manual que MTF tambin es dependiente del ruido, durante una medicin de la respuesta de impulso no debemos usar la compensacin.
ARTA muestra el MTF para la respuesta de impulso filtrada en bandas de octavas en una ventana independiente (Ilustracin 8.2) que se activa haciendo clic sobre el elemento de men Analysis > Modulation Transfer Function' [Anlisis > Funcin de Transferencia de la Modulacin] o haciendo clic sobre el icono de la barra de herramientas .
m X i = 10 log( i ) , 1 mi
2. Los valores de Xi estn limitados a 15 95
si Xi >15 entonces Xi =15, si Xi >-15 entonces Xi =-15, El mtodo STI indica que la relacin sonido /ruido en la gama de frecuencias de -15dB a 15dB es linealmente dependiente del ndice de inteligibilidad en la gama de 0 a 1. Es por esto que la relacin sonido/ruido se convierte en el ndice de transmisin TI.
TI i =
X i + 15 30
3. El valor medio de TIi para cada banda de octava (ndice de transmisin de octava OTI) se define con la siguiente ecuacin:
OTI n =
1 14 TI i , n = 1,2,..7 14 i =1
En la que los factores ponderados k y k estn experimentalmente determinados para la voz masculina o femenina (consulte la Tabla 8.1). Frecuencia (Hz) masculina masculina femenina femenina
125 0.085 0.085 0 0 250 0.127 0.078 0.117 0.099 500 0.230 0.065 0.223 0.066 1000 0.233 0.011 0.216 0.062 2000 0.309 0.047 0.328 0.025 4000 0.224 0.095 0.250 0.076 8000 0.173 0 0.194 0
Ilustracin 8.1 Factores ponderados para el habla masculina y femenina IEC 60268-16.
Los valores STI se encuentran en la gama de 0 a 1. En la Tabla 8.2 se incluyen los correspondientes ndices subjetivos.
STI 0,0 < STI < 0,3 0,3 < STI < 0,45 0,45 < STI < 0,6 0,6 < STI < 0,75 0,75 < STI < 1,00 ndice subjetivo equivalente Malo Pobre Suficiente Bueno Excelente
96
3. 4.
5. 6. 7.
ARTA User Manual Ajustar la ganancia del amplificador de potencia para que obtenga una SPL =68dBA. Controlar el espectro. Ha de tener los niveles de banda de octava que se determinan en la Tabla 8.3 (si fuera necesario usar el ecualizador de frecuencia). Colocar ARTA en el Modo de Respuesta de Impulso con una seal de excitacin del tipo Habla PN. El periodo del ruido peridico ha de ser mayor al periodo de la frecuencia de modulacin ms baja. (es decir, mayor de 1/0,63Hz=1.58 segundos). Quiere decir que el nmero de muestras en un periodo (longitud FFT) ha de ser mayor de 1,58Fs. La Tabla 8.4 nos indica la longitud FFT adecuada. Hacer la medicin de la respuesta de impulsos sin la compensacin, en una posicin de micrfono arbitraria. En la ventana de respuesta de impulsos active el comando de men 'Analysis-STI' o haga clic sobre el icono de la barra de herramientas. Finalmente, ARTA generar un informe con los valores MTF, los valores OTI y los valores STI, como en la Ilustracin 8.3.
Banda de octavas (Hz) Ruido del habla masculina (dB) Ruido del habla femenina (dB) 125 9.7 250 9.7 14.4 500 6.0 7.2 1000 0.0 0.0 2000 -6.0 -6.7 4000 -12.0 -7.6 8000 -18.0 -8.9
Tabla 8.3 Niveles de bandas de octava del ruido del habla (referidos a un nivel de 1kHz).
Fs (Hz) Longitud FFT 48000 131072 44100 131072 32000 65536 22020 65536
97
98
Ilustracin 8.4 %ALcons como funcin de la relacin sonido/ruido y el tiempo de reverberacin (datos experimentales para la posicin del oyente en el campo de reverberacin difuso).
El %ALcons se puede estimar a partir del tiempo de reverberacin medido T60 y de la relacin sonido/ruido usando los datos experimentales de Peutz de la Ilustracin 8.4 si el oyente est colocado en el campo de reverberacin difuso). Si no es as, si el oyente est en el campo del hablante directo, %ALcons se puede estimar usando la expresin: % ALcons 200d 2T60 = , VD
2
VD T60
En la que V es el volumen de la sala, D es la directividad del hablante, d es la distancia del hablante al oyente. Esta frmula se suele usar en clculos arquitectnicos, pero no tiene en cuenta la influencia del ruido. Tanto los ndices de inteligibilidad del habla, STI y %ALcons, son parmetros acsticos de la sala tiles. La comparacin con los ndices de inteligibilidad subjetivos muestra que STI nos da una estimacin dentro de 5.6% del ndice subjetivo, un resultado mejor que el 10% para %ALcons.
99
9 Herramientas
Esta seccin describe las herramientas que ARTA aade para la presentacin de resultados de mediciones.
La primera definicin representa un grupo de respuestas en frecuencia para una distancia de medicin constante r. La segunda definicin representa el grupo de respuestas en frecuencia normalizadas con respuesta en frecuencias en el grado cero. Puede ser mayor de uno. La tercera definicin representa el grupo de respuestas en frecuencia normalizadas con respuesta en frecuencia en un ngulo en el que la respuesta en frecuencia es mxima. Siempre es igual o inferior a uno. Estas tres definiciones se pueden usar en ARTA para mostrar los grficos de patrones de directividad. Normalmente, la magnitud de respuesta en el patrn de directividad se muestra en la escala de nivel de dB. Adems de la representacin grfica de las caractersticas direccionales, se definen tres parmetros de directividad [44]: 1) El factor de directividad Q(f) es la relacin entre la intensidad en un eje determinado de un radiador de sonido a una distancia r concreta respecto de la intensidad que se producira en la misma posicin por fuente de punto si estuviera radiando la misma potencia acstica total que el radiador. 2) El ndice de directividad DI(f) se expresa en dB como resultado de la expresin DI = 10log(Q). 3) El ngulo de ancho del haz de un patrn de directividad se define como el ngulo entre dos puntos a cada lado del eje principal (normalmente en el grado cero) en el que el nivel de presin de sonido es inferior a 6dB desde su valor al grado cero.
100
ARTA User Manual La construccin del grfico de cascada es simple, slo necesitamos medir las respuestas en frecuencia de los altavoces en ngulos diferentes, y de forma opcional normalizar la respuesta con una respuesta medida en el eje del grado cero.
El segundo tipo de grfico de patrn de directividad que se usa en ARTA es un sonograma. El sonograma grfico bidimensional que muestra planos de color. En cada plano, el nivel de la magnitud de respuesta se presenta como un valor del color de la paleta de colores predefinida. La Ilustracin 9.3 muestra el sonograma con colores que cambian gradualmente, mientras que la Ilustracin 9.4 muestra el sonograma contorneado con colores variables en tramos predefinidos. Normalmente, el sonograma contorneado muestra los patrones de directividad ms claramente que un sonograma gradual.
101
El tercer tipo de patrn de directividad en ARTA es un diagrama polar. Muestra el patrn de directividad para una frecuencia sencilla en el sistema coordinado polar (Ilustracin 9.5). El fondo del grfico muestra tambin el factor de directividad Q, el ndice de directividad DI y el ngulo del ancho del haz para la frecuencia de corriente. Se asume que el altavoz radia a un espacio sin lmites (condiciones de campo abierto). Se pueden mostrar varios patrones de directividad polar en el mismo grfico, como curvas de superposicin (como en la Ilustracin 9.6).
102
Para la radiacin en la media esfera suele ser ms apropiado usar un diagrama semi-polar, como en la Ilustracin 9.7.
103
104
La ventana del Patrn de directividad dispone de un men y varios controles. Los comandos del men son:
File [Archivo] Create directivity pattern file [Crear un archivo de patrn de directividad] abre un cuadro de dilogo para crear un archivo de patrn de directividad a partir de los archivos .pir. Save directivity pattern file [Guardar el archivo de patrn de directividad] guardar el archivo de patrn de directividad (.dpf). Load directivity pattern file [Cargar el archivo de patrn de directividad] carga el archivo de patrn de directividad (.dpf). Export ASCII 1/3 octave data [Exportar datos ASCII de 1/3 octava] exporta en valores de archivo ASCII en frecuencias estndar de 1/3 octava. Export ASCII & 1/1 octave data [Exportar datos ASCII y 1/1 octava] exporta en valores de archivo ASCII en frecuencias estndar de 1/1 octava. Export Directivity Index and Angle (-6dB) [Exporta el ndice y ngulo de directividad (6dB)] exporta en un archive DI, Q y el ngulo en ASCII (-6dB) File Info [Informacin de archive] proporciona informacin y texto definido por el usuario originaria del archive .dpf actual. Edit [Edicin] Copy [Copiar] copia el mapa de bits del grfico y el texto definido por el usuario al portapapeles B/W background color [Color de fondo B/N] determina el color de fondo en blanco o en negro
Los controles inferiores de la ventana son: Cuadro combinado Mode [Modo] Selecciona entre el tipo de grfico cascada1, cascada2, sonograma, diagrama polar, o semi-polar. Caja combinada Palette [Paleta] elija entre varias paletas de colores (Azabache, Gris, Cobre y Fro). 105
ARTA User Manual Colored Waterfall [Cascada coloreada] active esta casilla para elegir entre un grfico de cascada coloreado (o son un solo color). Contoured [Contorneado] active esta casilla para elegir entre un cambio de color gradual o contorneado (vea las Ilustraciones 6.19 y 6.20). Control de efecto Rango (dB) - cambia el rango dinmico del rango de 5 a 70dB. Botn Copy [Copiar] copia el grfico actual al portapapeles. Botn B/W [Blanco y negro] cambia el color de fondo del grfico entre blanco y negro.
Los controles del lado izquierdo de la ventana son: Control de efecto Frecuencia cambia la frecuencia del patrn polar actual a la siguiente banda de octava estndar de 1/3 Botn Superposicin abre el cuadro de dilogo del gestor de superposiciones. 0 => 0dB active esta casilla para usar el patrn de directividad con una magnitud normalizada con el valor de la magnitud en el grado cero, si no se hace as, el diagrama polar muestra el patrn con una magnitud normalizada con un valor de magnitud mximo. Contoured [Contorneado] active esta casilla para elegir entre un cambio de color gradual o contorneado (vea las Ilustraciones 6.19 y 6.20). Show DI [Mostrar DI] active esta casilla para ver el valor de los parmetros de directividad en el fondo del diagrama polar. Thick line [Trazo grueso] active esta casilla para usar una lnea de trazo grueso. El gestor de superposiciones del diagrama polar aparece en la Ilustracin 9.9. Es similar al gestor de superposiciones de una ventana FR, aunque hay una diferencia significativa. Hay una ventana emergente que se puede usar simultneamente con los comandos de las ventanas del patrn de directividad.
Los botones de la ventana del gestor de superposiciones activan los siguientes comandos:
Add [Aadir] establece la curva actual como superposicin Replace sel [Sustituir la seleccin] sustituye la superposicin seleccionada por la curva actual Delete sel [Borrar la seleccin] Borra todas las superposiciones seleccionadas Delete all [Borrar todo] Borra todas las superposiciones Color cambia el color de los elementos seleccionados usando el cuadro de dilogo 'Overlay Colors' [Colores de superposicin] que se muestra en la Ilustracin. 6.10. Update [Actualizar] actualiza las etiquetas del diagrama polar.
Haga clic con el ratn sobre el elemento de la lista desplegable y suceder lo siguiente:
Un solo clic - selecciona los elementos Un solo clic sobre la casilla de verificacin hace que sea visible o invisible la superposicin Doble clic - permite editar el nombre de la superposicin
106
El procedimiento para crear un archivo DPF comienza haciendo clic en el men 'File > Create directivity pattern file [Archivo > Crear un archivo de patrn de directividad]. Abre el cuadro de dilogo Directivity data definition' [Definicin de datos de directividad] (Ilustracin 9.11). Primero, hemos de pulsar sobre el botn 'Load Files [Cargar archivos] para obtener el directorio y nombre principal de los archivos PIR. As se abre el cuadro de dilogo estndar para abrir archivos, que se muestra en la Ilustracin 9.10. Despus de seleccionar uno de los archivos, hemos de pulsar sobre el botn OK. A continuacin se abrir el cuadro de dilogo 'Directivity data definition' [Definicin de datos de directividad] tal como se indica en la Ilustracin 9.11.
107
Antes de pulsar sobre el botn 'OK que confirma el proceso de generar datos del patrn de directividad, hemos de configurar los parmetros de creacin de datos usando los siguientes controles: Seccin Magnitude [Magnitud]: Cuadro combinado Suavizado elija el suavizado de la respuesta en frecuencia de 1/1 octava a 1/12 octava. Normalize to zero degree [Normalizar al grado cero] active esta casilla para hacer que las respuestas del patrn de directividad se normalicen con una respuesta en el grado cero. Symmetrical for neg. angles [Simtrico para ngulos negativos] - active esta casilla si slo dispone de archivos PIR definidos para ngulos positivos y desee que el patrn sea simtrico para los ngulos negativos. Seccin Frecuency range [Rango de frecuencias]: Low (Hz) [Bajo] introduzca el lmite de las bajas frecuencias para la estimacin de la respuesta en frecuencia High (Hz) [Alto] introduzca el lmite de las altas frecuencias para la estimacin de la respuesta en frecuencia Seccin FFT analysis [Anlisis FFT]: Ungated from estimated PIR maximum position [Sin limitacin desde la posicin PIR mxima estimada] active esta casilla para estimar la respuesta en frecuencia a partir de la respuesta limitada que se determina con el requisito de baja frecuencia y tiempo-ancho de banda y el punto de inicio del archive PIR que representa la primera muestra cuya magnitud es 20 dB inferior a la magnitud PIR mxima. Ungated from PIR start position [Sin limitacin desde la posicin de inicio PIR] seleccione esta casilla para la estimacin de la respuesta en frecuencia desde la primera muestra en el archive PIR y con una tamao FFT determinado por el requisito de baja frecuencia y tiempo-ancho de banda. 108
ARTA User Manual Gated [Limitada] seleccione esta casilla para obtener una estimacin de respuesta en frecuencia a partir de la posicin de inicio definida por el usuario en el archive PIR y la limitacin de Longitud de tiempo definida por el usuario. El Inicio y la Longitud se han de introducir en milisegundos.
La seccin de User info [Informacin de usuario] contiene un cuadro de edicin en el que el usuario introduce un texto arbitrario que se archivar como archivo DPF.
109
Bibliografa
[1] Mateljan, I., Signal Selection for the Room Acoustics Measurement, Proc. 1999 IEEE Workshop on Applications of Signal Processing to Audio and Acoustics, New Paltz, New York, 1999, ISBN - 0-78035612-8 Mateljan, I., Ugrinovi K. The Comparison of Room Impulse Response Measuring Systems, Proceedings of the First Congress of Alps Adria Acoustics Association, Portoroz, Slovenia, 2003, ISBN 961-6238-73-6
[2]
[3] Mateljan I., Audio Quality Measurements in Communication Systems, Proceedings of the Second Congress of Alps Adria Acoustics Association, Opatija, 2005, ISBN 953-95097-0-X. [4] Mateljan I., Models for the Estimation of the Loudspeaker In-Room Response, Int. Journal for Engineering Modelling, vol. 6., No.1-4, 1993, ISSN 1330-1365 [5] [6] [7] [8] [9] Rife, D.D., Vanderkooy, J., Transfer Function Measurement with Maximum-Length Sequences, J. Audio Eng. Soc., vol. 37, June 1989. Rife, D. D., Modulation Transfer Function Measurement with Maximum-Length Sequences, J. Audio Eng. Soc., vol. 40, October 1992. Farina, A., Simultaneous measurement of impulse response and distortion with a swept sine technique, 108 AES Convention, Paris, 2000. Pinelton, R., Schoukens, J., Measurement and Modeling of Linear Systems in the Presence of Non-Linear Distortions , Mechanical Systems and Signal Processing, 16(5), 2002. Schoukens, J., Pinelton, R., Ven der Ouderaa, E., and Renneboog, J., Survey of Excitation Signals for FFT Based Signal Analysers, IEEE Trans. Instrumentation and Measurement, vol. 37, September 1988.
[10] Ven der Ouderaa, E., Schoukens, J., and Renneboog, J., Peak Factor Minimisation of Input and Output Signals of Linear Systems, IEEE Trans. Instrumentation and Measurement, vol. 37, June 1988. [11] Dunn, C., and Hawksford, M.O., Distortion Immunity of MLS-Derived Impulse Response Measurement, J. Audio Eng. Soc., vol. 41, May 1993. [12] Vanderkooy, J., Aspects of MLS Measuring Systems, J. Audio Eng. Soc., vol. 42, April 1993. [13] MacWilliams, F. J., and Sloane, N. J. A., Pseudo Random Sequences and Arrays, Proc. IEEE, vol. 64, December 1976. [14] Bendat, J. S., and Piersol, A. G., Engineering applications of Correlation and Spectral Analysis, Wiley, New York, 1980. [15] Pierce, A. D., Acoustics An Introduction to Its Physical Principles and Applications, McGraw-Hill, NewYork, 1981. [16] Tan, Moore, Zacharov, The Effect of Nonlinear Distortion on Percived Quality of Music and Speech Signals, JAES, vol. 5, November, 2003 [17] IEC-268-16, Sound system equipment: the objective rating of speech intelligibility in auditoria by the RASTI method, International Electronical Committee, Geneva, IEC-268-16, 1988. [18] ITU-T Recommendation P.501, Test signals for use in telephonometry, 1996. [19] ISO -3382, Acoustics measurement of the reveberation time of rooms with reference to other acoustical parameters. 1997. [20] ISO Publication 266, Acoustics Preferred frequencies for measurements, 1975. [21] IEC 1260, Electroacoustics Octave-band and fractional octave-band filters, 1995. [22] IEC 61672-1:2002, Electroacoustics Sound level meters - Part 1 Specifications, 2002. [23] Keele, D.B., Low-Frequency Loudspeaker Assessment by Nearfield Sound-Preassure Measurement, JAES, Vol. 22, IV, 1974. [24] Heyser, R.C., Loudspeaker Phase Characteristics and Time Delay Distortion: Part 1, JAES, I, 1969. [25] Heyser, R.C., Loudspeaker Phase Characteristics and Time Delay Distortion: Part 2, JAES, IV, 1969. [26] Heyser, R.C., Determination of Loudspeaker Signal Arrival Times: Part I, JAES, X, 1971.
110
111
File [Archivo] New [Nuevo] - crea un nuevo archivo llamado Untitled.pit Open [Abrir] Abre el archivo Save [Guardar] Save As [Guardar como] guarda el archivo con un nombre nuevo Info [Informacin] muestra la informacin sobre el archivo .pir actual Export [Exportar] guarda los datos de la respuesta de impulso en los siguientes formatos: ASCII file [Archivo ASCII] guarda los datos de la respuesta de impulso en formato ASCII. MLSSA ASCII file [Archivo MLSSA ASCII] guarda los datos de la respuesta de impulso en formato MLSSA ASCII. WAV file [Archivo WAV] guarda los datos de la respuesta de impulso en formato .wav. Import [Importar] importa los datos de la respuesta de impulso de los siguientes formatos: WAV file [Archivo WAV] importa los datos de la respuesta de impulso en formato .wav. MLSSA .TIM file [Archivo MLSSA .TIM] importa los datos de la respuesta de impulso en formato MLSSA .TIM. MLSSA ASCII file [Archivo MLSSA ASCII] importa los datos de la respuesta de impulso en formato MLSSA ASCII. ASCII file [Archivo ASCII] importa los datos de la respuesta de impulso (tiempo-amplitud) en formato ASCII. Load and sum [Cargar y sumar] carga el archivo .pir y lo suma a la respuesta de impulso actual Recent File [Archivos recientes] abre un archivo que se haya abierto recientemente Exit [Salir] Sale del programa Edit [Edicin] Copy [Copiar] copia el mapa de bits del grfico y el texto definido por el usuario al portapapeles Colors and grid styles [Estilos de colores y cuadrcula] activa el cuadro de dilogo de configuracin del color B/W background color [Color de fondo B/N] determina el color de fondo en blanco o en negro Set Marker [Establecer marcador] - establece un marcador en la posicin del cursor Delete Market [Borrar marcador] Borra el marcador Invert [Invertir] - cambia la polaridad de la respuesta de impulsos Rotate at cursor [Girar en el cursor] gira la respuesta de impulso peridica (el cursor se convierte en la primera muestra) Scale [Escalar] multiplica la respuesta de impulsos por una constante arbitraria View [Ver] Toolbar [Barra de herramientas] muestra u oculta la barra de herramientas Status Bar [Barra de estado] muestra u oculta la barra de estado Record [Registro] Impulse response [Respuesta de impulsos] - activa el cuadro de dilogo de medicin de la respuesta de impulso Analysis [Anlisis] Smoothed Frequency response [Respuesta de frecuencia suavizada] - muestra la respuesta de frecuencia suavizada en bandas de 1/n octavas DFT frequency response [Respuesta de frecuencia DFT] muestra los componentes DFT de la respuesta de frecuencia. Step response [Respuesta de etapa] muestra la respuesta de etapa ETC Impulse Response Envelope [ETC - Envolvente de la respuesta de impulso] - muestra la envolvente de la respuesta de impulso (ETC) Energy decay [Deterioro de energa] muestra la curva de deterioro de energa de Schroeder Cumulative spectral decay [Deterioro espectral acumulado] muestra el cuadro de dilogo del deterioro espectral acumulado Modulation transfer function [Funcin de transferencia de modulacin] Muestra el MTF de octava
112
El texto de la izquierda muestra el contexto del men. El cuadro de texto de la derecha muestra la modalidad de medicin. Cuadros de texto izquierdo y derecho: Muestran los niveles mximos a escala completa de los canales de entrada izquierdo y derecho durante la ltima medicin.
113
FFT - establece el tamao de la ventana de anlisis FFT Window [Ventana] establece el tipo de ventana de seal Delay for phase estimation (ms) [Retardo para la estimacin de fase (ms)] establece el intervalo de tiempo, desde la posicin del cursor, de referencia para la estimacin de fase. Get [Obtener] - obtiene el intervalo de tiempo entre el cursor y el marcador Zero [Cero] establece un retardo cero Inv - cambia la polaridad de la respuesta de impulsos Acceso rpido con el ratn Puede cambiar la posicin del marcador y del cursor pulsando y arrastrando los botones izquierdo y derecho del ratn. Haga doble clic sobre el botn derecho del ratn para activar o desactivar el marcador. Arrastrando el ratn a la zona de las etiquetas se desplaza la curva horizontal o verticalmente. Haciendo doble clic sobre el botn izquierdo del botn sobre el rea del eje de tiempo cambia la etiqueta de posicin tiempo/muestra. Acceso rpido con el teclado Teclas arriba y abajo - para cambiar la ganancia Teclas Ctrl+arriba y Ctrl+abajo - para cambiar la compensacin vertical Teclas izquierda y Ctrl+izquierda para desplazar la curva hacia la izquierda. Teclas derecha y Ctrl+derecha para desplazar la curva hacia la derecha Teclas Shift+izquierda y shift+derecha para mover el cursor hacia la izquierda y la derecha Tecla Repag y Avpag - para cambiar el factor de ampliacin Tecla Ctrl + S guardan el archivo Tecla Ctrl + N crean un nuevo archivo Tecla Ctrl + O abren el archivo Tecla Ctrl + C - copia un grfico de mapa de bits al portapapeles (con las opciones definidas por el usuario) Tecla Ctrl + P copia toda la ventana como mapa de bits al portapapeles Tecla Ctrl + B cambia el color del fondo
114
Men de la ventana del analizador de espectros File [Archivo] Export ASCII [Exportar a ASCII] exporta los valores del grfico a un archivo de formato ASCII. Exit [Salir] Sale del programa Overlay [Superposicin] Set overlay [Establecer superposicin] establecer la respuesta en frecuencia actual como superposicin Delete overlay [Borrar superposicin] Eliminar la superposicin Load overlay [Cargar superposicin] Cargar el archivo con la superposicin de la respuesta en frecuencia Save overlay [Guardar la superposicin] Guardar la respuesta en frecuencia en un archivo de superposicin Edit [Edicin] Copy [Copiar] copia el mapa de bits del grfico y la etiqueta definida por el usuario al portapapeles Colors and grid styles [Estilos de colores y cuadrcula] activa el cuadro de dilogo de configuracin del color View [Ver] Toolbar [Barra de herramientas] muestra u oculta la barra de herramientas Status Bar [Barra de estado] muestra u oculta la barra de estado B/W background color [Color de fondo B/N] determina el color de fondo en blanco o en negro Fit graph top [Ajustar la parte superior del grfico] Ajusta los mrgenes del grfico al valor mximo de la magnitud Frequency Axis [Eje de frecuencias] Lin - lineal Log - Logartmico Octave Smoothing [Suavizado de octavas] - Respuesta logartmica suavizada en bandas de 1/n octavas) Octave Bands [Bandas de octava] Respuestas en bandas de 1/n octavas Resolution [Resolution] 1/1 octave [1/1 octava] - ancho de banda del filtro de 1/1 octava 1/2 octave [1/2 octava] - 1/2 octava 1/3 octave [1/3 octava] - 1/3 octava 1/6 octave [1/6 octava] - 1/6 octava 1/9 octave [1/9 octava] - 1/9 octava 1/12 octave [1/12 octava] - 1/12 octava 1/24 octave [1/24 octava] - 1/24 octava Scaling [Graduacin] dB FS - relativo a la escala completa (dB) dB V (SPL) - relativo al valor RMS (dB) PSD - Densidad espectral de potencia dB/sqrt(Hz) Voltage units [Unidades de voltaje] dBV - dB re 1V Bu - dB re 0,775V (1mW/600)
Sound pressure units [Unidades de presin de sonido] dB re 20uPa - unidad para el nivel de presin en acstica dB re 1Pa - unidad para el nivel de presin en las recomendaciones ITU_T Weighting [Ponderacin] None [Ninguna] no se usa filtro de ponderacin A filtro A estndar B filtro B estndar C filtro C estndar
115
Mode [Modo] Impulse response [Respuesta de impulsos] - activa la ventana de la respuesta de impulso Spectrum Analyzer [Analizador de espectros] - activa la ventana del anlisis de espectros Dual channel frequency response [Respuesta de frecuencia en canal doble] - activa la ventana de respuesta de frecuencia en canal doble Single channel frequency response [Respuesta de frecuencia en canal sencillo] - activa la ventana de respuesta de frecuencia en canal sencillo Help [Ayuda] About [Acerca de]proporciona informacin sobre el programa Registration [Registro] muestra la informacin de registro / usuario User manual [Manual de usuario] - muestra el archivo de ayuda
116
Gen establece el tipo de generador Fs (Hz) - establece la frecuencia de muestreo FFT establece el tamao de un bloque de adquisicin FFT Wnd - establece la ventana de seal para el anlisis FFT Avg - establece el tipo de compensacin Reset [Borrar] - borra las compensaciones
Acceso rpido con el ratn: Si se pulsa y arrastra el botn izquierdo del ratn se cambia la posicin del cursor Si se hace clic con el botn derecho del ratn en la zona del grfico se activa el cuadro de dilogo de la configuracin del margen del grfico. Si se hace clic con el botn derecho del ratn en la zona del ttulo se activa el cuadro de dilogo de la graduacin del grfico. La rueda de desplazamiento del ratn mueve el margen superior del grfico hacia arriba y abajo. Acceso rpido con el teclado:
117
Men de la ventana de respuesta de frecuencia (2 canales o 1) File [Archivo] Export ASCII [Exportar a ASCII] exporta los valores del grfico a un archivo de formato ASCII. Save as PIR [Guardar como PIR] - transforma la respuesta en frecuencia actual en una respuesta en impulso que la archiva como un archivo .PIR. Exit [Salir] Sale del programa Overlay [Superposicin] Set overlay [Establecer superposicin] establecer la respuesta en frecuencia actual como superposicin Delete overlay [Borrar superposicin] Eliminar la superposicin Load overlay [Cargar superposicin] Cargar el archivo con la superposicin de la respuesta en frecuencia Save overlay [Guardar la superposicin] Guardar la respuesta en frecuencia en un archivo de superposicin Generate target response [Generar respuesta objetivo] Generar una superposicin con respuesta de filtros de cruce estndar Delete target response [Borrar respuesta objetivo] Generar una superposicin con respuesta de filtros de cruce estndar Edit [Edicin] Copy [Copiar] copia el mapa de bits del grfico y el texto definido por el usuario al portapapeles Colors and grid styles [Estilos de colores y cuadrcula] activa el cuadro de dilogo de configuracin del color B/W background color [Color de fondo B/N] determina el color de fondo en blanco o en negro
View [Ver] Toolbar [Barra de herramientas] muestra u oculta la barra de herramientas Status Bar [Barra de estado] muestra u oculta la barra de estado Fit graph top [Ajustar la parte superior del grfico] Ajusta los mrgenes del grfico al valor mximo de la magnitud Frequency Axis [Eje de frecuencias] Linear [Lineal] establece el eje de frecuencias lineal. Log establece el eje de frecuencias logartmico Octave Smoothing [Suavizado de octavas] eje logartmico para magnitudes suavizadas en 1/n octavas Octave Bands [Bandas de octava] barras para bandas de 1/n octavas Resolution [Resolution] 1/1 octave [1/1 octava] - ancho de banda del filtro de 1/1 octava 1/2 octave [1/2 octava] - 1/2 octava 1/3 octave [1/3 octava] - 1/3 octava 1/6 octave [1/6 octava] - 1/6 octava 1/9 octave [1/9 octava] - 1/9 octava 1/12 octave [1/12 octava] - 1/12 octava 1/24 octave [1/24 octava] - 1/24 octava Sound pressure units [Unidades de presin de sonido] dB re 20uPa/1V - unidades para la respuesta de presin en los estndares acsticos dB re 20uPa/2.83V - unidades para la sensibilidad del altavoz dB re 1Pa/1V - unidades para la respuesta de presin en las recomendaciones ITU_T Recorder [Grabadora]
118
119
Acceso rpido con el teclado: Teclas arriba y abajo - cambian el margen del grfico superior Teclas izquierda y derecha - mueven el cursor hacia la izquierda y la derecha Tecla Ctrl + C - copia un grfico de mapa de bits al portapapeles (con las opciones definidas por el usuario) Tecla Ctrl + P copia toda la ventana como mapa de bits al portapapeles Tecla Ctrl + B cambia el color del fondo Tecla Ctrl+A - establece la curva trazada actual como superposicin
Men de la ventana de la respuesta de frecuencia suavizada File [Archivo] Export ASCII [Exportar a ASCII] exporta los valores del grfico a un archivo de formato ASCII. Repeat PIR measurement [Repetir medicin PIR] abre el cuadro de dilogo de la medicin PIR. Si es correcto, calcula la respuesta de frecuencia usando la posicin del cursor y del marcador a partir de la ventana PIR. Save PIR as [Guardar PIR como] guarda la ltima medicin o archivo .PIR Overlay [Superposicin] Set as overlay [Establecer superposicin] se fija la curva de impedancia actual como superposicin Set as overlay Below cursor [Establecer superposicin por debajo del cursor] guarda la parte de la curva actual por debajo del cursor como superposicin Set as overlay Above cursor [Establecer superposicin por encima del cursor] guarda la parte de la curva actual por encima del cursor como superposicin Manage Overlays [Gestor de superposiciones] activa el cuadro de dilogo 'FR Overlay Manager' [Gestor de superposiciones de respuestas de frecuencia] para editar la lista de superposiciones. Delete all [Borrar todo] Borra todas las superposiciones Delete last [Borrar ltima] Borra la ltima superposicin Generate target response [Generar respuesta objetivo] Generar una superposicin con respuesta de filtros de cruce estndar Delete target response [Borrar respuesta objetivo] Generar una superposicin con respuesta de filtros de cruce estndar Edit [Edicin] Copy [Copiar] copia el mapa de bits del grfico y el texto definido por el usuario al portapapeles Colors and grid styles [Estilos de colores y cuadrcula] abre el cuadro de dilogo de configuracin del color del grfico B/W background color [Color de fondo B/N] determina el color de fondo en blanco o en negro Use thick pen [Usar trazo grueso] -usa un trazo grueso para trazar las curvas Cut below cursor [Cortar debajo del cursor] corta los valores del grfico por debajo del cursor Cut above cursor [Cortar encima del cursor] corta los valores del grfico por encima del cursor Scale level [Nivel de escala] escala el nivel con un factor arbitrario (diferencia en dB) LF box diffraction [Difraccin de caja en bajas frecuencias] Escala los niveles con la respuesta de la difraccin en bajas frecuencias de la caja del altavoz Subtract overly [Restar superposicin] resta los valores del nivel de superposicin de la curva actual Subtract from overlay [Restar de la superposicin] - resta los valores del nivel de la curva actual de la superposicin
120
Acceso rpido con el teclado: Teclas arriba y abajo - cambian el margen del grfico superior Teclas izquierda y derecha - mueven el cursor hacia la izquierda y la derecha Tecla Ctrl + C - copia un grfico de mapa de bits al portapapeles (con las opciones definidas por el usuario) Tecla Ctrl + P copia toda la ventana como mapa de bits al portapapeles Tecla Ctrl + B cambia el color del fondo Tecla Ctrl+A - establece la curva trazada actual como superposicin Tecla Ctrl+M - abre el cuadro de dilogo de gestin de superposiciones
Men de la ventana de la respuesta de frecuencia DFT File [Archivo] Export ASCII [Exportar a ASCII] exporta los valores del grfico a un archivo de formato ASCII. Repeat PIR measurement [Repetir medicin PIR] abre el cuadro de dilogo de la medicin PIR. Si es correcto, calcula la respuesta de frecuencia usando la posicin del cursor y del marcador a partir de la ventana PIR. Save PIR as [Guardar PIR como] guarda la ltima medicin o archivo .PIR Edit [Edicin] Copy [Copiar] copia el mapa de bits del grfico y el texto definido por el usuario al portapapeles Colors and grid styles [Estilos de colores y cuadrcula] abre el cuadro de dilogo de configuracin del color del grfico B/W background color [Color de fondo B/N] determina el color de fondo en blanco o en negro
121
Teclas de acceso directo del teclado (tambin para el resto de ventanas): Teclas arriba y abajo - cambian el margen del grfico superior Teclas izquierda y derecha - mueven el cursor hacia la izquierda y la derecha Tecla Ctrl + C - copia un grfico de mapa de bits al portapapeles (con las opciones definidas por el usuario) Tecla Ctrl + P copia toda la ventana como mapa de bits al portapapeles Tecla Ctrl + B cambia el color del fondo
Men de la ventana de la respuesta de etapa File [Archivo] Export ASCII [Exportar a ASCII] exporta los valores del grfico a un archivo de formato ASCII. Edit [Edicin] Copy [Copiar] copia el mapa de bits del grfico y el texto definido por el usuario al portapapeles B/W background color [Color de fondo B/N] determina el color de fondo en blanco o en negro Men de la ventana de deterioro espectral acumulado Edit [Edicin] Copy [Copiar] copia el mapa de bits del grfico y el texto definido por el usuario al portapapeles B/W background color [Color de fondo B/N] determina el color de fondo en blanco o en negro
Ventana de deterioro de energa Edit [Edicin] Copy [Copiar] copia el mapa de bits del grfico y el texto definido por el usuario al portapapeles Set marker [Establecer marcador] - establece un marcador en la posicin del cursor Delete Market [Borrar marcador] Borra el marcador B/W background color [Color de fondo B/N] determina el color de fondo en blanco o en negro
Ventana de la envolvente de respuesta de impulsos (ETC ) Edit [Edicin] Copy [Copiar] copia el mapa de bits del grfico y el texto definido por el usuario al portapapeles B/W background color [Color de fondo B/N] determina el color de fondo en blanco o en negro Set marker [Establecer marcador] - establece un marcador en la posicin del cursor Delete Market [Borrar marcador] Borra el marcador
Men de la ventana MTF Edit [Edicin] Copy [Copiar] copia el mapa de bits del grfico y el texto definido por el usuario al portapapeles
122
123