21.12.07

Some Gigs

Aqui unas fotitos de las tocatas que he ido. Warehouse project en el cual toco Dj Marky (y caleta de otros locos), simplemente espectacular... la ultima fue un lanzamiento que me pegue con mis phd-mates ingleses, era bien similar a ver trainspotting... el grupo re bueno Jeffrey Lewis



Dj Marky






El otro grupito fue Jeffrey Lewis

No hay muchos comentarios, solo decir que fueron espectaculares lanzamientos y que se vienen mas iguales y/o mejores
Cheers!








26.11.07

This is England

You've thought more than ten times that the car you have just seen was driven by nobody...lol

19.6.07

Bone Thugs and Harmony en Chile

Tuvieron que pasar casi 14 años para ver a mi grupo favorito en vivo. Crecí escuchando a Bone Thugs and Harmony, fueron la inspiración del grupo musical -Punto Rojo, aquí el myspace de uno de mis partner en esa aventura- que tuve en todo el sentido de la palabra y ahi estaban el dia 16 de junio en la Estacion Mapocho de Santiago de Chile, agradablemente ese mismo día mi equipo de futbol favorito, Colo Colo, salía tricampeón.
Como habia que hacerla toda ese dia nos fuimos tempranito con mi polola y mi hermano a almorzar al Mercado Central, su pescado frito con arroz y de ahi caminamos unos pasos para llegar a la Estación Mapocho. Ahi me encontré con la Feña y su hermana, ya eramos 5, un buen numero.

Hicimos la fila, peluseamos un rato y cuando empezaba a chispear abrieron la puerta de la Estación y corriendo pa' quedar lo mas adelante posible, lo cual significaba estar a unos 20 metros -al lado de la consola de audio, como a mi me encantan... tenia que estar cerca- del escenario ya q no tenia 10 lucas mas pa pagar la wea VIP. Asi de paradojico, VIP en un evento hip hop, como han cambiado las cosas...
Otro sintoma de cambio fue que antes me encontraba con varios amigos en este tipo de eventos, ahora no vi a niuno, vi a harto pendejito que con cueva sabia hablar cuando Bone Thugs se llamaba Bone Enterprise, pero mal que mal, esta bien que asi sea. Ya pasaron mis tiempos de artista estrella jajajaj.

Hubieron varios Djs, todos sin criterio. Si viene Bone Thugs coloca algo similar al menos pos cabro. el unico descriteriado pero con buen gusto al menos era un tal DJ supi o suri o una wea asi.
Y ahi esperando pa Bone, empezaron esos Djs que no tenian niun brillo de ahi cantó una banda nacional llamada Yomando Crew, sinceramente lamentables los socios, de ahi vino Tea Time, el ex-de los tetas, fome la wea también, algunas reminiscencias de las bases de Busta Rhymes pero Tea-Time esta lejos de tener la rima y el estilo de Busta pos loco, si el weon tiene cara de niñito de bien, es un especimen común de La Reina o Providencia y ante eso poco y nada se puede hacer.

Yo siempre he odiado a los grupos españoles de Hip Hop, y estos aweonados de los organizadores trajeron 2 pa ma recacha... pa q telonearan a Bone los sacoweas. Cero criterio.
El primero que apareció fue Juani CACA, lo mataron al loco si todos ibamos a ver a Bone.. y como todo español (salvo Nash Scratch que tiene uno q otro tema) era pésimo, si el mismo tipo preguntó si querían q siguiera y la respuesta fue un No rotundo, el tipo dijo.. bueno la sinceridad ante todo y el desgraciado siguió cantando, ya me estaba poniendo de mal humor, ahi se empezaron a pasar desde la entrada general a la parte Vip, en la general estabamos pa la caga de apretados y los niñitos VIP hasta tenian espacio pa sentarse y bailar (si weon suena a resentido pero me da lo mismo), se llenó al toque la wea de adelante lo cual me alivio bastante, ya que me quedé en la parte general no se si fue de weon, por que cada dia acato mas las reglas o porq sabia que cuando todos se fueran al Vip iba a estar todo mas relajado en la zona general. Es que hasta una pendeja rucieta del barrio alto pica a flayte me dijo oiga USTED me podría dejar pasar, con todas las S bien marcaditas (como debe ser), buena actriz la cabra. Sólo atiné a reirme y confié que como todos se iban pa adelante iba a estar relajado, a la altura del sonidista, iba a escuchar bien y ver decentemente, todo ello se cumplió. En resumen todo español vale callampa y el weon canto como 1 hora = 1 hora odiando

De ahi vino Frank T, otro "español" este loco era un poco mejor, pero en realidad el tipo es un inmigrante con acento pepe. El Dj del compadre notable, la intro muy buena, dije al menos no será tan tortuoso, se tiró sus buenas rimas Frank y a seguir aguantando, ya al tercer tema estaba chato con el loco. Ni siquiera bajaría un tema del loco pero lo respeto. Ya al quinto tema estaba odiando en mala, al menos ya se habia todo dispersado y me pude fumar un cigarro sin quemar a nadie. Al fin terminó de cantar y...

Apareció al fin Bone Thugs and Harmony, salieron de a uno: primero Layzie de ahi Krayzie y luego Wish quedó la pata, un sueño hecho realidad, ya con sólo ese primer tema (Flow Motion del Strength and Loyalty 2007) las lucas estaban mas que bien gastadas. No así como cuando vino Lauryn Hill que tocó como 4 temas del unplugged y niun brillo.

Aqui no cabe ese asunto de saltar por saltar y pegarle al de al lado, la musica no da pa eso, aqui bailamos y disfrutamos compadre, cantamos las canciones y, tal como el nombre indica, es todo armonía.. felicidad socio. Asi fueron tocando cada uno de los exitos .

Bone preguntaba si estabamos listos y sonaba. "...We're not against rap. We're not against rappers. But we are against those thugs (thugs thugs)" y Thuggish Ruggish Bone en todo su esplendor, temazo del primer album de ellos ya como Bone Thugs and Harmony.
Asi, de ese mismo album sonó también Foe tha Love of $ y No Surrender, justo antes de ese tema Wish estaba pa la caga porque queria contar la historia de que pasa cuando te persiguen los pacos y te estay fumando un pito, lamentablemente el ingles no es muy comprendido y el chicano que trajeron fue de gran ayuda pa comunicarse con el publico, era una lata igual por que Layzie decia Somebody scream y la bulla aparecia como al minuto despues, por la actitud de los locos al principio creyeron que el público no estaba apañando pero menos mal q se dieron cuenta que era por que no les cachaban a la primera lo que querian decir.
Del orden de los temas no me acuerdo mucho pero si claramente recuerdo que tocaron
1st of tha month (con una version especial con Felecia cantando también), el tema con Notorious BIG, dayz of our livez, y absolutamente casi todos los exitos de los locos (creo que faltó Paper del Resurrection, a ver si completamos la lista en orden algun dia)

Mención especial es que pusieron el comienzo del album Makaveli de 2 Pac y de ahi cantaron el temita que tienen con 2 Pac (Thug Luv), nunca crei que lo fueran a tocar... 2 Pac esta muerto al parecer...y bueno lo hicieron revivir compadre pa que cantara ese dia. Simplemente espectacular.

Sin duda el tema que mas disfruté fue The Crossroads, Wish indicó que esto se lo dedicaban a toda la gente que ya no está con nosotros y empezó a sonar la base , ahi vi pasar millones de imagenes por segundo en mi cerebro, crecí con ese tema, es lejos mi favorito no se si es porq estoy mas viejo o que chucha pero con esas weas me emociono, me caian las lagrimas cantando la wea. Es algo indescriptible el bombardeo de imagenes y recuerdos que sentí escuchando el temazo. Ya me podia ir tranquilo pa la casa, sin embargo estos weones siguieron y creo que ahi vino la parte mas relax del concierto, cuando empezaron tocar los temas dedicados a la marihuana, Weed Song sonó por cierto y se sintió mas fuerte que nunca el olorcito característico de los recitales con la consecuente hambre y sed, i.e. su sorbo de jugo y su pancito. Ahi hubo otra tremenda sorpresa porq colocaron el comienzo del temita ese de Dr Dre que dice smoke weed every day no cache lo que significaba la frase jajaja pero igual la grité pa apañar a Layzie que preguntaba desesperado que es la wea q tenian q hacer a diario, cuando escuchó al publico gritar la frasesita ahi se quedó tranquilo. Entremedio de todo eso Layzie tiro su gorrito y resulto que el cabro culiao del primo de mi polola lo agarró, Grande Biach gato e campo y la conchetumare jajajaja te lo voy a robar de una jajajaj.
De ahi tocaron el temazo I tried, todos coreando pos loco, el disco salio hace poquito, de ese mismo tocaron tambien Bump In The Trunk, y creo que Lil Love y Order My Steps as well.
Estos locos se viraron y weiamos harto pa que aparecieron de nuevo, tocaron como 3 o 4 temas mas, pa suerte mia tocaron el temita Ecstacy del Resurrection, que gran tema ese, otro de mis favoritos

Cuando tocaron Crossroad yo habia empezado a cantar la parte de Bizzy y lamentablemente esa parte ya no está, si bien Bone Thugs and Harmony en estos momentos es sólo Krayzie, Layzie y Wish y es espectacular; Bizzy le daba tremendo toque al grupo, ojala las diferencias entre ellos se solucionen un dia, y poder verlos a los 5, sí cuando salga Flesh de la cana: Free Flesh (escuchen el tema World So Cruel ó NorthCoast de Flesh-n-Bone del disco T.H.U.G.S pura clase)

Este concierto junto con el de Mogwai han sido lejos los mejores a los cuales he asistido. Gracias Bone Thugs and Harmony por darme la oportunidad de disfrutar de un evento tan pulento. Las 7 horas esperando valieron mas q la pena, valieron recuerdos, emociones y una alegría inmensa, mas encima fue el primer concierto que asistí con mi hermano asi que mejor aun.

Eternal E 1999 Bone Thugs and Harmony 4eva

Flaytech aka Lord G

22.3.07

Are you NERD?

Sincerely, This survey is outstanding. Moreover, I did not know that I am so nerd. Do you try?. Click on the image

I am nerdier than 91% of all people. Are you nerdier? Click here to find out!

They suggest me: Apply for a professorship at MIT now!!!. ("..It is a great pleasure to accept this professorship..." Si ohh... algun dia)

Ya cabros, todos los que visitan esto quizas tengan algo de Nerd, sino no me lo explicaría. Notable encuesta, aun me estoy riendo de las preguntas.

Saludos a todos,
Flaytech

8.3.07

Silencer Modeling on Comsol. Part I

Dear Sir or Madam

This is the first post on silencer modeling by using COMSOL MULTIPHYSICS. The goal is to present a tutorial series on how to modelling silencers under different conditions and the influence of the materials involved in the construction process of these acoustics devices.
All information presented here is product of the leisure and ... Bingo!! I don't have job. Is necessary more justification??.

Other important objective is write in English, If there are some mistakes, please let me know what, when and how to solve them. Is clear that the mathematics is more easy??

This and subsequent posts assumes some basic knowledge on acoustics, silencers, finite element method, COMSOL Multiphysics and common sense.

Here we go...

As the title says: I am going to present a brief tutorial on how to modelling a simple expansion chamber (without absortive material) by using Comsol v 3.2 package. In this post, an axial symetric 2D simple expansion chamber model is presented.

The parameters of the simple expansion chamber are:
r01=10 cm; L01=30 cm
r02=30 cm; L02=40 cm
r03=10 cm; L03=30 cm

1) Open Comsol Multiphysics and select Axial Symmetric 2D, Acoustics->Time Harmonic Analysis



2) Create the geometry



3) Set the constants (Alt+O+C)



where
rho_air: Air density
c_air: Sound velocity in the air
p0: Arbitrary pressure

4) Set the scalar expressions (Alt+O+E+c)



Where
Dt: Transmission Loss
Pin: Pressure at the input of the simple expansion chamber
Pout: Pressure at the output of the simple expansion chamber

5) Set Sub-domain settings (Alt+ y+ Enter)



6) Set the boundary conditions (Alt+y+B)









7) Calculate the Intensity at the input of the simple expansion chamber (Alt +O+E+B)



8) Calculate the Intensity at the output of the simple expansion chamber (Alt +O+E+B)



9) Integrate the intensity on the input boundary in order to obtain Pin (Alt+O+r+B)



10) Integrate the intensity on the output boundary in order to obtain Pout (Alt+O+r+B)



11) Create and refine the mesh



12) Set the frequency on Scalar variables (Alt+y+V)



13) Set the solve parameters (Alt +S+ P)



The list of parameters values is arbitrary. In this case the frequency varies from 50 Hz to 2000 Hz in steps of 50 Hz. If you have a powerfull computer with enough RAM, you can calculate the transmission loss with more frequency resolution.

14) Press the “=” button

15) This figure presents the sound propagation at 500 Hz



16) Finally, the following steps are necessary in order to obtain the Transmission Loss
a) Go to the Postprocessing menu and select Domain Plot Parameters (Alt+P+D)
b) Select all frequencies from the Solutions to plot list in the Domain Plot Parameters dialog box.
c) Click the Keep current plot page.
d) Click the Title/Axis button and type Transmission loss in the Title edit field.
e) Click OK.
f) Select Point plot in Domain Plot Parameters dialog box and click on the Point tab.
g) Select point 1 from the Point selection list.
h) Type Dt in the Expression edit field.
i) Click the Line Settings button and select Triangle in the Line marker list.
j) Click OK.
k) Click Apply in the Domain Plot Parameters dialog box.





16) The Transmission Loss is:



Easy and funny, Your comments, suggestions and clicks on the publicity are welcomes (remember that I dont have job, help me pleeeease jajajja)

The next post will cover some Porous Material models

Yours sincerely,
Flaytech

16.1.07

Empirical Mode Decomposition (EMD) and Head Related Transfer Function (HRTF) - Si, es verdad, primer post tematico

Bien, es hora de escribir sobre ciencia

Principlamente son tres los objetivos:
1) Divulgar
2) Que lo que se presente sea criticado
3) Sugerencias del lector

Como no tengo mucho remordimiento, estos post seran mezcla entre español, ingles y los idiomas que sean necesarias, incluyendo obviamente la matematica dentro de esa categoria.

En este post abordaremos rapidamente la tecnica llamada Empirical Mode Decomposition (EMD) porpuesta por Huang et al [1]. Esta tecnica descompe cualquier señal en sus denominados modos intrinsecos, estos son dependientes de la señal y sirven para reconstruir la señal original. El principio es similar a la serie (transformada) de Fourier, la cual utiliza una base de funciones ortogonales (senos y cosenos) para escribir como combinacion lineal cualquier señal que cumpla ciertas caracteristicas, si ese serie de Fourier que sea periodica, si es transformada de fourier se hace periodica a la fuerza y otra importante es que la señal a descomponer tenga un numero finito de discontinuidades para asegurar la convergencia (ver series de fourier y transformada de Fourier).

Cual es el "problema": Transformada de Fourier asume linealidad y estacionariedad de la señal.

El combinar EMD con la transformada de Hilbert da paso a la transformada de Hilbert-Huang la cual no asume absolutamente nada acerca de la naturaleza de la señal y por tanto es ideal para trabajar con señales no lineales y no estacionarias [2], como por ejemplo AUDIO.

Mi objetivo es aplicar EMD al analisis de HRTF. Siendo especifico a cada HRIR, la idea es descomponer las HRIR con dicha técnica en sus modos intrinsecos (IMF) y correlacionar los IMF con data antropomorfica de cada persona, es decir, asociar cada IMF con la influencia del torso, pinna, cabeza, etc...

Vamos al temita en sí:

Brief Introduction

The empirical Mode Decomposition is a relatively new method for analysing nonlinear and non-stationary data. The key part of the method is the `empirical mode decomposition' method with which any complicated data set can be decomposed into a finite and often small number of `intrinsic mode functions' (IMF) that admit well-behaved Hilbert transforms. In this work, the first approach considers the aplication of the Fourier Transform on each Intrinsic Mode Functions in order to analyze the decomposition in the frequency domain and correlated these spectrums with anthropomorphic data.

Given such a signal, the method adaptively decomposes it into a number of modes (IMFs) that are topologically equivalent to amplitude and frequency modulated, sinusoidal signals. In the analytic signal representation, the modes correspond to proper rotations [3]. Thus, the EMD method naturally yields estimates of the significant instantaneous frequencies embedded in the signal, by performing the Hilbert transform on each IMF, in this work the first approach consist in applying the Fourier Transform on each IMFs.

Let x(t) a time dependent signal, the EMD method focuses on the level of local oscillations and decomposes the signal into a finite and often a small number of fundamental oscillatory modes. The bases (IMFs) into which the signal is decomposed are obtained from the signal itself, and they are defined in the time domain. They have the same length as the original signal and preserve the frequency variations with time. The base modes can be made approximately complete and nearly orthogonal with respect to each other. Here, completeness implies that the original signal can be reconstructed without any loss of data by simply adding up the IMFs. Therefore the IMFs can be viewed as linear components of the original or source signal x(t)

In order to achieve this, two conditions need to be satisfied:

(a) The total number of extrema of IMF(t) must be equal to the number of zero crossings (or only differing by at most one)

(b) The mean of the upper envelope IMFu(t) and the lower envelope IMFl(t) be zero.

The IMFs calculation process is called sifting process and is composed by the following steps:

1) Identify all the extrema of x(t)

2) Interpolate between successive maxima and minima (say cubic splines), respectively, to obtain upper and lower envelopes (say, by using cubic splines), IMFl(t) y IMFu(t) respectively.

3) Calculate the local mean m(t) between the envelopes

4) Subtraction of the average from the original to yield d(t)=x(t)-m(t)

5) Repite steps (1-3) until d(t) satisfies the two conditions for being an IMF. Once an IMF is generated, the residual r(t)=x(t)-IMF1(t) is regarded as the original signal, and steps (1-4) are repeated to yield the second IMF, and so on

Huang [2] propose two stopping criteria:

a) The S number: S is defined as the consecutive number of siftings, in which the numbers of zero-crossing and extrema are the same for these S siftings.

b) SD is smaller than a pre-set value, where

Another criteria is proposed by Rilling et al [3] [4]. This criteria is based on two thresholds, teta1 and teta2, and aimed at guaranteeing globally small fluctuations in the mean while taking into account locally large excursions. This amounts to introduce the mode amplitude

And the evaluation function

Thus, the sifting is iterated until

for some preescribed fraction (1-alfa) of the total duration, while

One can tipically set alfa=0.05; teta1=0.05 and teta2=10*teta1 (default values in emd.m at [4])

The procedure is complete when either the residual function becomes monotonic, or when the amplitude of the residue falls below a pre-determined small value so that further sifting would not yield any useful components. These features guarantee the computation of a finite number of IMFs within a finite number of iterations. The outcome of the EMD procedure is the following decomposition of the original signal:

Applying the Fourier Transform (The Fourier Transform is the first stage of this paper. If there are bad results, the second stage considers the Hilbert Transform and consequently to apply the Huang-Hilbert Transform to the HRTF analysis)

Then, the hypotheses are:

1) It is possible to recognize the contributions of the different human organs, such as pinnae, torso, head, etc; to the spectral changes of the HRTF using the empirical mode decomposition method

2) It is possible to label each IMF to morphological data.


Method and Preliminary Results

This work uses the CIPIC HRTF Database [5] and the EMD implementation by Rilling et al [3] [4]

The emd options were setted as:

opt.display=2;
opt.maxiterations=2000;
opt.stop=[0.01,0.1,0.02];

The figure 1 shows the left ear Hrir for 0º azimuth and 0º elevation, the figure 2 for the right ear. (Sampling frequency was 44.1 kHz)

Figure 1 - HRIR (0,0) - Left ear - Subject003 CIPIC HRTF Database

Figure 2 - HRIR (0,0) - Right ear - Subject003 CIPIC HRTF Database

The figures 3 and 4 show the respective HRTFs.

Figure 3 - HRTF (0,0) - Left ear - Subject003 CIPIC HRTF Database

Figure 4 - HRTF (0,0) - Right ear - Subject003 CIPIC HRTF Database

The figure 5 shows the first five IMF of the HRIR (0,0) related to the left ear, the figure 6 show the IMF6 to IMF9 and the residual.

Figure 5 - IMF1 to IMF5 of the HRTF (0,0) - Left ear - Subject003 CIPIC HRTF Database

Figure 6 - IMF6 to IMF9 and the residual of the HRTF (0,0) - Left ear - Subject003 CIPIC HRTF Database

The figure 7 and 8 displays the samen information for HRIR (0,0) related to the rigth ear.

Figure 7 - IMF1 to IMF5 of the HRTF (0,0) -Right ear - Subject003 CIPIC HRTF Database

Figure 8 - IMF6 to IMF9 and the residual of the HRTF (0,0) - Right ear - Subject003 CIPIC HRTF Database

The figure 9 shows the absolute difference between the original signal and the signal constructed by mean of the equation

Figure 9 - Absolute differences between the original signal and the reconstructed signal of the HRIR (0,0) Left ear (above) and HRIR (0,0) Right ear (below)- Subject003 CIPIC HRTF Database

The difference are negligible (order 1e-16). With this result, it is possible to assure that the decomposition process is at least suitable.

The figure 10 shows the HRTF calculated from the original signal (above left = left ear; above right = right ear) and form the reconstructed signal by using the IMF and the residual (below left = left ear; below rigth = rigth ear). Again, the difference are negligible (see figure 1 in the time domain)

Figure 10 - HRTF from the original data (above) and the reconstructed data (below)

The figures 11 and 12, show the fourier transform of the IMF1 to IMF3 for the HRIR (0,0) of the left ear and right ear respectively
Figure 11 - Fourier Transform of the IMF1 (blue), IMF2(green) and IMF3 (red) of the decomposed signal (HRIR (0,0) left ear - Subject 003 - CIPIC HRTF Database)

Figure 12 - Fourier Transform of the IMF1 (blue), IMF2(green) and IMF3 (red) of the decomposed signal (HRIR (0,0) right ear - Subject 003 - CIPIC HRTF Database)

The results indicate that the interpretation of each IMF spetrum to the contribution to the respective HRTF is not direct which can be associated to the nonlinear process involved in the magnitude spectrum calculation process (absolute value and logarithm). It is necessary to carry out more test in order to interpretate each IMF spectrums better.

Finally, the figure 13 and 14 show the Hilbert spectrum of the IMFs (Hilbert-Huang Transform)

Figure 13 - Hilbert-Huang spectrum of the HRIR (0,0) Left ear - Subject 003 - CIPIC HRTF Database)
Figure 14 - Hilbert-Huang spectrum of the HRIR (0,0) Right ear - Subject 003 - CIPIC HRTF Database)

Preliminary Conclussions

1) Es posible descomponer las HRIR en sus modos intinsecos IMF por medio de la tecnica EMD con gran precisión.
2) La intepretacion de los espectros de cada IMF no es inmediata
3) Posiblemente seria adecuado evitar el proceso de calculo de magnitud del espectro (eludir la no linealidad de ese proceso en si) y trabajar con graficas de cada espectro (Fourier o Huang-Hilbert) de cada IMF en las cuales el eje x corresponde a la parte real y el eje y a la parte imaginaria. Salen monos rebonitos, como los que vienen a continuacion:

Los 8 primeros espectros de Hilbert-Huang de los IMF del oido derecho (HRIR(0,0) )


BUENO, ESO ES TO ESO ES TO ESO ES TODO AMIGOS, CUALQUIER SUGERENCIA, CRITICA, COMENTARIO, AYUDA SERA BIENVENID@

Suggestions and comments are always welcome, see you later

By The "serious" FLAYTECH

Se me habia olvidado la bibliografía
[1] Huang, Norden E. et al, The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis, Royal Society of London Proceedings Series A, 1998.
[2] Huang, Norden E., Beyond the Fourier Transform :Coping with Nonlinear, Nonstationary Time Series, Johns Hopkins University,1998.
[3] G. Rilling, P. Flandrin and P. Gonçalvès, On empirical mode decomposition and its algorithms, IEEE-EURASIP Workshop on Nonlinear Signal and Image Processing NSIP-03, Grado (I), 2003
[4] Flandrin, P, Empirical Mode Decomposition, Matlab codes
[5] V. R. Algazi, R. O. Duda, D. M. Thompson and C. Avendano, "The CIPIC HRTF Database," Proc. 2001 IEEE Workshop on Applications of Signal Processing to Audio and Electroacoustics, pp. 99-102, Mohonk Mountain House, New Paltz, NY, Oct. 21-24, 2001.