{"id":180,"date":"2012-03-02T13:36:47","date_gmt":"2012-03-02T11:36:47","guid":{"rendered":"http:\/\/www.joanillo.org\/?p=180"},"modified":"2020-09-16T12:46:11","modified_gmt":"2020-09-16T10:46:11","slug":"jugando-con-la-fft","status":"publish","type":"post","link":"https:\/\/www.joanillo.org\/?p=180&lang=es","title":{"rendered":"Jugando con la FFT"},"content":{"rendered":"

Hace ya mucho tiempo (realment mucho tiempo) que estudiaba la Transformada de Fourier<\/strong> (http:\/\/en.wikipedia.org\/wiki\/Fourier_transform<\/a>) en la facultad, e incluso utilic\u00e9 el algoritmo de la FFT<\/strong> (Fast Fourier Transform, http:\/\/en.wikipedia.org\/wiki\/Fast_Fourier_transform<\/a>) en un aplicativo escrito en C para procesamiento digital de im\u00e1genes. Desde entonces no me hab\u00eda preocupado de la FFT<\/strong> y realmente lo ten\u00eda muy olvidado.<\/p>\n

Estos d\u00edas he vuelto a la carga con la FFT<\/strong> pues, con la idea de aprender las t\u00e9cnicas de programaci\u00f3n de aplicaciones de audio, quer\u00eda programar un afinador o un analizador de la potencia espectral de una se\u00f1al, que fuese compatible con JACK.<\/p>\n

Antes que nada hay que mirar qu\u00e9 algoritmo de FFT se podr\u00eda utilizar, con licencia GPL, y r\u00e1pidamente he visto que la elecci\u00f3n ser\u00eda FFTW<\/strong> (http:\/\/www.fftw.org\/<\/a>):<\/p>\n

FFTW is a C subroutine library for computing the discrete Fourier transform (DFT) in one or more dimensions, of arbitrary input size, and of both real and complex data (as well as of even\/odd data, i.e. the discrete cosine\/sine transforms or DCT\/DST). We believe that FFTW, which is free software, should become the FFT library of choice for most applications.<\/p><\/blockquote>\n

Me he puesto a compilar unos cuantos ejemplos y a familiarizarme con los c\u00e1lculos. En definitiva yo lo que quiero es calcular la potencia espectral de una se\u00f1al de audio, por tanto el se\u00f1al de entrada ser\u00e1n unas muestras reales (una sola dimensi\u00f3n). La salida de la FFT ser\u00e1n unos valores complejos, y he de tener en cuenta tanto la parte real como la parte imaginaria para hacer los c\u00e1lculos de la potencia espectral.<\/p>\n

fftPlan = fftwf_plan_dft_r2c_1d(fftSize, fftIn, fftOut, FFTW_MEASURE); \/\/r2c_1d: real to complex, one dimension<\/p><\/blockquote>\n

Mirando si hay alg\u00fan afinador de c\u00f3digo libre, que trabaje en modo consola, y compatible con JACK, me he encontrado estos dos proyectos, que precisamente utilizan la librer\u00eda FFTW:<\/p>\n