matlab 倒位序fft程序,FFT算法设计与实现

FFT算法研究报告

1、程序设计背景(FFT算法理解)

FFT(fast fourier transformation),快速傅里叶变换。是对DFT算法的改进，其利用了WNnk的周期性、共轭对称性和可约性，使得DFT中有些项可以合并，大大减小了计算量。

按输入序列在时间上的次序是属于偶数还是奇数来分解称为“按时间抽取法”(DIT)。另一种是把输出序列X(k)按顺序的奇偶分解为越来越短的序列，称为按频率抽样的FFT算法(DIF)。DIT算法是先作复乘后作加减，而DIF的复乘只出现在减法之后。本次程序采用DIT算法实现FFT。

用c语言实现FFT的难点在于数据倒位序的处理，以及各级蝶形运算的实现。倒位序的实现可以使用“反向进位加法”，即倒位序二进制数的下面一个数是上面一个数在最高位加一并由高位向低位进位而得到的。

对于点数为N = 2^L的FFT运算，可以分解为L阶蝶形图级联，第M阶蝶形图内又分为2^(L-M)个蝶形组，每个蝶形组内包含2^(M-1)个蝶形。而且旋转因子与蝶形阶数和蝶形分组内的蝶形个数存在关联。因此我们就可以构造循环来实现蝶形运算。

2、FFT算法流程图