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 ] ( 2373 visualizaciones ) | [ 0 trackbacks ] | enlace permanente | ( 3 / 2909 )