Ответ: Вот... проверено :)))
#include
#define SQINTWO 1.414213 //Корень квадратный из двойки
#if !defined M_PI //Если PI не определено:
#define M_PI 3.14159265359 //Определяем его :)
#endif
#define FFT -1 //Преобразование Фурье
#define REV_FFT 1 //Обратное преобразование
double y_sig[1000]={0};
double x_sig[1000]={0};
//Быстрое преобразование фурье(БПФ)
void fft(float *x,float *y,int order,int param)
{
unsigned int n,l,e,f,i,j,o,o1,j1,i1,k;
float u,v,z,c,s,p,q,r,t,w,a;
n=1u<