首先离散时间傅里叶变换的公式是
令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倍,在后面的代码中需要缩小。