Una FFT muy pedagógica 
He encontrado una implementación muy sencilla y pedagógica (aunque no muy eficiente en términos de memoria) del algoritmo de Cooley-Tukey de FFT para arrays de tamaño potencia de 2. Se pueden obtener del depatramento de informática de la universidad de Princeton, más concretamente en los siguientes enlaces: FFT.java y Complex.java.

He estado haciendo pruebas calculado el error cuadrático medio entre x y IFFT(FFT(x)) y se obtienen valores inferiores a 1E-30 para tamaños de entrada lo suficientemente grandes (65536 o más). El algoritmo no es muy eficiente en términos de memoria ya que hace varios new en cada llamada recursiva, pero eso facilitará la comprensión del mismo.

El siguiente paso será utilizar esta implementación de la FFT para implementar el algoritmo que dejé aparcado a un lado hace algunas semanas: el de separación de fuentes de sonido mediante discriminación por acimut. El objetivo final es extraer la voz de una grabación comercial estéreo. A ver si sale. Ya iré contando mis logros y mis sinsabores :-).

[ 9 comentarios ] ( 2233 visualizaciones )   |  [ 0 trackbacks ]   |  enlace permanente
  |    |    |    |   ( 3 / 2704 )

<< <Anterior | 53 | 54 | 55 | 56 | 57 | 58 | 59 | 60 | 61 | 62 | Siguiente> >>