快速傅利葉（FFT）

快速傅利葉：

更加形象的理解傅利葉變換：

大概了解之後：

從傅利葉級數到傅利葉變換：

（**太大，只能裁剪為兩張）

刨根問底的同學：

雷德演算法：

輸出序列是按自然順序排列的，而輸入序列的順序則是「位元反轉」方式排列的。也就是說，將序號用二進位制表示，然後將二進位制數以相反方向排列，再以這個數作為序號。如011變成110，那麼第3個輸入值和第六個輸入值就要交換位置了。

詳見思路：

#include #include #include #include #define pi 3.1415926   //定義圓周率值

#define n 16 //變換點數
using namespace std;
typedef struct compx
compx;
//複數乘法
compx multi(compx a, compx b)
//複數加法
compx add(compx a, compx b)
//複數減法
compx sub(compx a, compx b)
//快速傅利葉
void fft(compx data)
k = n/2;
//實現進製
while(k <= j) //左邊最高位為1
j += k;
}//計算蝶形級數m
int f = n, m;
for(m = 1; (f=f/2) != 1; m++);
//控制級數
for(int n = 1; n <= m; n++)
u = multi(u, w);}}
}int main()
fft(data);
for(int i = 0; i < n; i++)
printf("%lf \t%lf\n", data[i].real, data[i].img);
return 0;
}