fft的FPGA实现

首先离散时间傅里叶变换的公式是

令Wn=e^(-j2pi/N)，则X[k]可以表示为

按时间抽取基2的FFT算法的实现不多介绍，自己查资料。下面给出N=8点的示意图

基于时间抽取的FFT算法和按频率抽取的FFT算法的输入数据的顺序是不同的，所以在一开始要在fpga中进行码位倒置。这里简单说一下，第一级的时候每组有2个数，一共有4组，每组要进行1次蝶形运算。第二级每组4个数，一共2组，每组进行2次蝶形运算，所以一共有log2(N)级。本次实验，我用的是1024个点，所以应该是10级。考虑先用matlab生成Wn的实部、虚部和一开始的输入数据。
生成虚部和实部的类似只要把cos改为sin即可。这里需要注意的是，Wn被放大了2^8倍，在后面的代码中需要缩小。