数字信号处理实验报告1.pdfVIP

数字信号处理实验 快速傅里叶变换与信号频谱分析 一． 实验目的 1. 在理论学习的基础上，通过本实验加深对离散傅里叶变换的理解。 2. 熟悉并掌握按时间抽取编写快速傅里叶变换（FFT）算法的程序。 3. 了解应用FFT 进行信号频谱分析过程中可能出现的问题，例如频谱混淆 二．实验原理，内容与分析 1. 仔细分析教材第六章“时间抽取法 FFT 的 FORTRAN 程序”，编写出相应的使用 FFT 进行信 号频谱分析的Matlab 程序。 (1)时间抽取FFT 算法 有限长序列 x(n)的 N 点DFT 定义为： N 1 X (k ) x (n )W kn N n 0 2 j 式中W e N 。 N 2 若对x (n) 做一次DFT 运算，时间复杂度为 级，其运算效率低下。而利用基于时间（或频率） N 抽取FFT 方法，时间复杂度可降至N log N 级，大大提高了效率。 k 基于时间抽取的FFT 方法将序列x (n) 分解为两组：偶数和技术，利用W 的周期性和对称性， N 可将X (k ) 表达为前后两部分：  X (k ) X (k ) W k X (k ), k 0,1,..., N/ 2 1  1 N 2  X (N / 2 k ) X (k ) W k X (k ), k 0,1,..., N/ 2 1  1 N 2  用蝶形图表示为： 只要 N 是2 的整数次幂，这种分解就可一直进行下去，将DFT 运算转换为.. 级运算，每一级运 算都由N / 2 个碟形运算组成。利用长度为N 的数组存放输入数据和每一级计算的结果。 (2).原位计算 对 M 点的FFT 共进行 级运算，每级由 个蝶形运算组成。在同一级中，每个蝶形的 N 2 M N / 2 输入数据只对本蝶有用，且输出节点与输入节点在同一水平线上，这就意味着每算完一个蝶形块后， 所得数据可立即存入原输入数据所占用的数组元素(存储单元)，这种原位(址)计算的方法可节省大 量内存。 (3). 蝶形运算 蝶形运算是分级进行的；每级的蝶形运算可以按旋转因子的指数大小排序进行；如果指数大小一样则 可从上往下依次计算。对N 2M 点的 FFT 共有 级运算，用 表示从左到右的运算级数 M L L 1 L 1, 2,...M 。第 级共有 个不同指数的旋转因子，用 表示这些不同指数旋转因子从上到