数字信号处理实FFT变换.docVIP

数字信号处理 实验报告 姓名: 胡剑烽 学 号： 11座机电话号码 学 院: 计算机工程学院 专 业: 通信工程 题 目: FFT变换及其应用 2014 年 11 月 一 实验目的 1. 在理论课学习的基础上，通过本次实验，加深对 DFT 原理的理解，懂得频域DFT 与时域卷积的关系，进一步加深对DFT 基本性质的理解； 2. 研究 FFT 算法的主要途径和编程思路，掌握FFT 算法及其程序的编写过程，掌握最 基本的时域基-2FFT 算法原理及程序框图； 3. 熟悉应用 FFT 实现两个序列的线性卷积的方法，利用FFT 进行卷积，通过实验比较 出快速卷积优越性，掌握循环卷积和线性卷积两者之间的关系； 4. 熟悉应用 FFT 对典型信号进行频谱分析的方法，初步了解用周期图法作随机信号谱 分析的方法，了解应用FFT 进行信号频谱分析过程中可能出现的问题，以便在实际 中正确应用FFT； 5. 掌握使用 MATLAB 等基本开发工具实现对FFT 编程。 二 实验设备 1. Windows 2000 以上操作系统； 2. Visual C++6.0/ Visual Basic6.0/Delphi6.0/MatLab6.5 等以上版本的开发环境； 3. 每人一台 PC 机。 三 实验原理 （一）离散傅里叶变换（DFT） 有限长序列的离散傅里叶变换的定义为： x n 和 X k 是一个有限长序列的离散傅里叶变换对，分别称为 DFT、IDFT。已知其 中的一个，就能惟一地确定另一个。 X k 也可以看作序列x n 的傅里叶变换 X e jω 在区间［0, 2π］上的N 点等间隔采样， 其采样间隔为2 / N ω π N ，这就是 DFT 的物理意义。 DFT 具有以下一些性质：线性、圆周移位、圆周卷积、有限长序列的线性卷积与圆周 卷积及共轭对称性等。 （二）快速傅里叶变换（FFT） 由于 DFT 运算过程复杂及耗时较大，诞生了FFT。这里仅讨论按时间抽取的基-2 FFT 设序列x n 长度为 N，且满足N 2M ，M 为正整数。按n的奇偶把x n 分解为两个 N/2 点的子序列： 则x n 的 DFT 转化为： 其中， 可见，一个N 点DFT 分解成两个N/2 点的DFT。 但由于1X k 、 2X k 只有 N/2 个点，而 X k 却有 N个点，故计算得到的只是 X k 的前一半的结果，要用1X k 、 2X k 来表达全部的 X k 值，还必须应用系数的周期性 即有： 因此得： （三）用FFT 进行谱分析 由上面的分析可知，若信号本身是有限长的序列，计算序列的频谱就是直接对序列进行 FFT 运算求得 X k ， X k 就代表了序列在[0,2π ]之间的频谱值。 · 幅度谱： · 相位谱： 若信号是模拟信号，用FFT 进行谱分析时，首先必须对信号进行采样，使之变成离散 信号，然后就可按照前面的方法用FFT 来对连续信号进行谱分析。按采样定理，采样频率s f 应大于2 倍信号的最高频率，为了满足采样定理，一般在采样之前要设置一个抗混叠低通滤 波器。 用FFT 对模拟信号进行谱分析的方框图如下所示。 在运用FFT 进行频谱分析的时候可能有混淆现象、泄漏现象和栅栏效应。 （五）FFT 求解函数 MATLAB 提供了求解FFT 的函数。参见下表。 四 实验内容 （一）用 FFT 进行谱分析 1．高斯序列： n 0:15; p 8; q 2; x exp -1* n-p .^2/q ; close all; subplot 3,1,1 ; stem fft x ; %利用fft 函数实现傅里叶变换 subplot 3,1,2 ; stem abs fft x ; %绘制幅度谱 subplot 3,1,3 ; stem angle fft x %绘制相位谱 结果如下图所示： P 8,q 2 P 8,q 4 P 8,q 8 q 8,p 13 q 8,p 14 从上面的图中，根据p和q值相应的变化，可以看出： 固定p 8，改变q的值：随着q的增大，经过傅里叶变化后的图像显示值的变化比较缓慢，幅度谱变化随着q的变大而变大，但相位的变化较不同。 当固定q 8，随着p的增大，时域信号幅值变换的比较缓慢。 2．正弦序列 n 0:15; %定义序列长度 a 0.1; f 0.0625; x exp -a*n .*sin 2*pi*f*n ; close all; subplot 2,1,1 ; stem x ; title 衰减正弦序列 ; subplot 2,1,2 ; stem abs fft x ; %绘制幅度谱 title x 信号的频谱 结果如下图所示： a 0.1、f 0.