实验2matlab中基2ditfft的实现.pdf

电 子 科 技 大 学 实 验 报 告 学生姓名： 学 号：2010013080 指导教师 一、实验室名称：数字信号处理实验室 二、实验项目名称： FFT 的实现 三、实验原理： 一． FFT 算法思想： 1．DFT的定义： 对于有限长离散数字信号 {x[n]} ，0 n N-1 ，其离散谱 {x[k]} 可以由离 散付氏变换（ DFT）求得。 DFT的定义为： N 1 2 j nk X [ k] x[ n]e N ，k=0,1, …N-1 n 0 2 j 通常令 e N WN ，称为旋转因子。 2．直接计算 DFT的问题及 FFT的基本思想： 由 DFT的定义可以看出，在 x[n] 为复数序列的情况下，完全直接运算 N 点 2 DFT需要（N-1） 次复数乘法和 N （N-1）次加法。因此，对于一些相当大的 N值 （如1024）来说，直接计算它的 DFT所作的计算量是很大的。 FFT的基本思想在于，将原有的 N 点序列分成两个较短的序列，这些序列的 DFT可以很简单的组合起来得到原序列的 DFT。例如，若 N 为偶数，将原有的 N 2 2 点序列分成两个 （N/2）点序列，那么计算 N 点 DFT将只需要约 [(N/2) ·2]=N /2 2 次复数乘法。即比直接计算少作一半乘法。因子（ N/2 ） 表示直接计算（ N/2 ） 点 DFT所需要的乘法次数，而乘数 2 代表必须完成两个 DFT。上述处理方法可以 反复使用，即（ N/2）点的 DFT计算也可以化成两个（ N/4 ）点的 DFT （假定N/2 为偶数），从而又少作一半的乘法。这样一级一级的划分下去一直到最后就划分 成两点的 FFT运算的情况。 3．基 2 按时间抽取（ DIT）的 FFT算法思想： L 设序列长度为 N 2 ，L 为整数（如果序列长度不满足此条件，通过在后面 补零让其满足）。 1 / 13 将长度为 N 2L 的序列 x[ n]( n 0,1,..., N 1) ，先按 n 的奇偶分成两组： x[2 r] x [ r ] 1