la semana pasada estuve leyendo unos artículos sobre métodos de análisis basados en la escala ERB y su aplicación a la acústica de recintos:
http://www.tml.tkk.fi/~ktlokki/Publs/lokkidafx.pdf
http://rvl4.ecn.purdue.edu/~malcolm/app ... onsEar.pdf
El caso es que me pareció que los métodos descritos eran facilmente desarrollables a partir de un script en python y los python-bindings que estoy preparando sobre DRC (http://drc-fir.sourceforge.net/doc/drc.html)
Estos bindings están ya casi terminados para hacerlos públicos (la pu.. documentación). Y estaba desarrollando unos scripts de prueba y de demo de la utilidad del asunto (escasa).
El caso es que tengo preparado un script para Windows (¿vaaaale?) para python 2.3 (http://www.python.org/ftp/python/2.3.5/Python-2.3.5.exe) que toma un impulso en formato raw-32 bits floating point y lo descompone en las (creo) 40 bandas ERB y lo integra a saltos de tiempo siguiendo (creo) lo descrito en el artículo primero.
Cada banda filtrada (por convolución) la guarda en un wav, para análisis de la propia implementación del algoritmo.
Y finalmente guarda el resultado de la integración en un fichero (output.dat) que es manejable por gnuplot (para windows).
http://sourceforge.net/project/showfile ... up_id=2055
Como necesito alguna mente pensante más para depurar o validar o sugerir o lo que sea, os enlazo un rar con todo el pack completo.
El script se llama test_erb.py. Si se instala python 2.3 con hacer doble click en el script arranca, siempre y cuando haya un fichero impulse.pcm
en el mismo directorio.
http://rapidshare.com/files/13001942/erb_test.rar.html
Os incluyo un impulse.pcm que fue medido en Molingordo para unas ATC SCM-12 para que os valga como ejemplo de lo que se necesita. Se puede importar con audacity:
http://audacity.sourceforge.net/
En formato RAW, 32-bits-floating-point-endianless, 44100 Hz. Así que la señal impulsiva que os interese analizar la guardáis en el directorio donde este el script con el nombre impulse.pcm y a darle al click.
Básicamente, es un impulso de sweep-scope (creo).
Incluye un script para gnuplot que dibuja algo así:
![Imagen](http://img213.imageshack.us/img213/2086/erbatc2ng.png)
Una vez abrierto gnuplot y dando a File -> Open y abriendo erb_test.plt, que es el dichoso script. Y a partir de ahí con el ratón se puede rotar el dibujo.
Se agradecerán todo tipo de comentarios. Se trata de un alfa-testing.
El que quiera código fuente de los python-bindings para DRC que me envie un privado y se los paso.
En cuanto tenga terminadas pruebas y documentación lo publicaré por aquí.
Gracias.