{"id":178,"date":"2012-03-02T13:35:21","date_gmt":"2012-03-02T11:35:21","guid":{"rendered":"http:\/\/www.joanillo.org\/?p=178"},"modified":"2020-09-14T13:36:33","modified_gmt":"2020-09-14T11:36:33","slug":"jugant-amb-la-fft","status":"publish","type":"post","link":"https:\/\/www.joanillo.org\/?p=178","title":{"rendered":"Jugant amb la FFT"},"content":{"rendered":"

>Fa ja molt de temps (realment molt de temps) que estudiava la Transformada de Fourier<\/strong> (http:\/\/en.wikipedia.org\/wiki\/Fourier_transform<\/a>) a la facultat, i fins i tot vaig utilitzar l’algorisme de la FFT<\/strong> (Fast Fourier Transform, http:\/\/en.wikipedia.org\/wiki\/Fast_Fourier_transform<\/a>) en un aplicatiu escrit amb C per a processament digital d’imatges. Des d’aleshores en\u00e7\u00e0 que no m’havia preocupat de la FFT<\/strong> i realment ho tenia molt oblidat.<\/p>\n

Aquests dies he tornat a la c\u00e0rrega amb la FFT<\/strong> doncs, amb la idea d’aprendre les t\u00e8cniques de programaci\u00f3 d’aplicacions d’audio, volia fer o un afinador o un analitzador de la pot\u00e8ncia espectral d’un senyal, que fos compatible amb JACK.<\/p>\n

Primer de tot he mirat quin algorisme de FFT es podria utilitzar, amb llic\u00e8ncia GPL, i r\u00e0pidament he vist que la tria seria 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

M’he posat a compilar uns quants exemples i familiaritzar-me amb els c\u00e0lculs. En definitiva jo el que vull \u00e9s calcular la pot\u00e8ncia espectral d’un senyal d’audio, per tant el senyal d’entrada seran unes mostres reals (una sola dimensi\u00f3). La sortida de la FFT seran uns valors complexos, i he de tenir en compte tant la part real com la part imagin\u00e0ria per tal de calcular la pot\u00e8ncia espectral.<\/p>\n

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

Mirant si hi ha algun afinador de codi lliure, que treballi per la consola, i compatible amb JACK, m’he trobat amb aquests dos projectes, que precisament utilitzen la llibreria FFTW:<\/p>\n