Znalazłem proste i dobrze działające funkcje fft i ifft. Za ich pomocą robie analize 512 próbek. Na wyjściu mam 2 tablice liczb typu double. Spektrum uzyskuje w ten sposób (wyciete ze znalezionego kodu):
double *real, ///czesc rzeczywista
*imag, ///ta druga
*spec; ///widmo
///...
fft(512,real,imag);
double a = 0;
int j = 0;
for(j = 0; j < 512; i++)
{
spec[j] = sqrt(real[j] * real[j] + imag[j] * imag[j]);
if(spec[j] > a) a = spec[j];
}
for(i = 0; i !!!