数字信号处理-快速傅里叶变换fft实验报告.pdfVIP

数字信号处理 实验报告 姓名：______ ______ 学号：________ 班级：_____ _____ 实验日期：2014 年 10 月 15 日 提交日期：2014 年 10 月 23 日 实验一 快速傅里叶变换与信号频谱分析 一、 实验目的 1. 在理论学习的基础上，通过本实验加深对离散傅里叶变换的理解； 2. 熟悉并掌握按时间抽取编写快速傅里叶变换（FFT ）算法的程序； 3. 了解应用FFT 进行信号频谱分析过程中可能出现的问题，例如频谱混淆、泄漏、栅 栏效应等，以便在实际中正确使用 FFT 算法进行信号处理。 二、 实验内容 1. 仔细分析教材第六章“时间抽取法 FFT 的FORTRAN 程序”，编写出相应的使用 FFT 进行信号频谱分析的 Matlab 程序。 2. 用 FFT 程序分析正弦信号，分别在以下情况进行分析，并讨论所得的结果： a) 信号频率 F ＝50Hz，采样点数N=32 ，采样间隔T=0.000625s ； b) 信号频率 F ＝50Hz，采样点数N=32 ，采样间隔T=0.005s ； c) 信号频率 F ＝50Hz，采样点数N=32 ，采样间隔T=0.0046875s ； d) 信号频率 F ＝50Hz，采样点数N=32 ，采样间隔T=0.004s ； e) 信号频率 F ＝50Hz，采样点数N=64 ，采样间隔T=0.000625s ； f) 信号频率 F ＝250Hz ，采样点数N=32 ，采样间隔T=0.005s ； g) 将 c) 中信号后补32 个 0 ，做 64 点 FFT ，并与直接采样 64 个点做 FFT 的结果进 行对比。 3. 思考题： 1) 在实验 a)、b) 、c)和 d) 中，正弦信号的初始相位对频谱图中的幅度特性是否有影 响？为什么？ 2) 信号补零后做 FFT 是否可以提高信号频谱的分辨率？为什么？ 三、 实验要求 1. 记录实验内容中各种情况下的 X (k)值，做出频谱图并深入讨论结果，说明参数的变 化对信号频谱产生哪些影响。频谱只做幅度特性，并按照幅度最大值进行归一化。 2. 打印出编写的 FFT 源程序，并且在每一小段处加上详细的注释说明。 3. 用 Matlab 编写 FFT 程序时，要求采用人机界面形式：N ，T ，F 以及初始相位等变 量均由键盘输入，补零或不补零要求设置一开关来选择。 4. 回答思考题，撰写实验报告。 四、 设计思路 4.1 码位倒置 按实验要求，需要对正弦信号进行采样，采样后信号按顺序储存，但此刻 DFT 运算结 果是乱序的，因此需要根据结果所需要的顺序进行码位倒置的变址操作。 变址的实现方法是：将原址用二进制表示后左右翻转，再转回十进制。需要注意的是 1 MATLAB 中数组首元素序号为 1，所以在十进制转二进制之前地址需要减 1，二进制转十 进制之后需要加 1。 4.2 蝶形运算 蝶形运算使用三层循环结构，外层循环控制级数 m 的递推（1 ~ log N ），中层循环体 2 进行不同组（运算系数相同的蝶形运算对视作一组），内层循环对每组内各运算对进行蝶形 运算。同级运算的相邻两组运算系数满足以下表达式： ‒ ‒ 1 1 = 1， = ‒ 1 ∗ 2 五、 实验结果与分析 a) 信号频率 F ＝50Hz，采样点数N=32 ，采样间隔T=0.000625s ；