2数字信号处理56实验指导书20141013.doc

《数字信号处理》实验指导书 实验序号：1 实验名称：利用FFT进行谱分析 1．加深DFT算法原理和基本性质的理解。 2．熟悉FFT算法原理和FFT子程序的应用。 3．学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析 误差及其原因，以便正确应用FFT。 4．加深理解FFT在实现数字滤波（或快速卷积）中的重要作用，更好的利用FFT进行数字信号处理。 5．掌握循环卷积和线性卷积两者之间的关系。 二、实验原理 1．在工程技术的许多分支中，要掌握的基本内容之一就是正确理解时域和频域的关系。对于数字系统来说，就是要精通离散傅立叶变换，因此离散傅立叶变换在数字信号处理中占有十分重要的地位。在实际应用中，有限长序列有相当重要的地位，(DFT)的定义： 快速傅里叶变换（FFT）并不是一种新的变换，而是离散傅里叶变换（DFT）的一种快速算法。 用FFT对信号作频谱分析是学习数字信号处理的重要内容。经常需要进行谱分析的信号是模拟信号和时域离散信号。对信号进行谱分析的重要问题是频谱分辨率D和分析误差。频谱分辨率直接和FFT的变换区间N（即x(n)长度为N）有关，因为FFT能够实现的频率分辨率是，因此要求。可以根据此式选择FFT的变换区间N。误差主要来自于用FFT作频谱分析时，得到的是离散谱，而信号（周期信号除外）是连续谱，只有当N较大时离散谱的包络才能逼近于连续谱，因此N要适当选择大一些。 周期信号的频谱是离散谱，只有用整数倍周期的长度作FFT，得到的离散谱才能代表周期信号的频谱。如果不知道信号周期，可以尽量选择信号的观察时间长一些。 对模拟信号进行谱分析时，首先要按照采样定理将其变成时域离散信号。如果是模拟周期信号，也应该选取整数倍周期的长度，经过采样后形成周期序列，按照周期序列的谱分析进行。 2．数字滤波器根据系统的单位脉冲响应h(n)是有限长还是无限长可分为有限长单位脉冲响应（Finite Impulse Response）系统（简记为FIR系统）和无限长单位脉冲响应（Infinite Impulse Response）系统（简记为IIR系统）。 对于FIR滤波器来说，除了可以通过数字网络来实现外，也可以通过FFT的变换来实现。 首先我们知道，一个信号序列x(n)通过FIR滤波器时，其输出应该是x(n)与h(n)的卷积： 当h(n)是一个有限长序列，即h(n)是FIR滤波器，且时 在数字网络类的FIR滤波器中，普遍使用的横截型结构就是按这个卷积公式构成的。 应用FFT实现数字滤波器实际上就是用FFT来快速计算有限长度序列的线性卷积。 这种方法就是先将输入信号x(n)通过FFT变换为它的频谱采样值X(k)，然后再和FIR滤波器的频响采样值H(k)相乘，H(k)可事先存放在存储器中，最后再将乘积H(k)X(k)通过快速傅里叶变换（简称IFFT）还原为时域序列，即得到输出y(n)。 现以FFT求有限长序列间的卷积及求有限长度序列与较长序列间的卷积为例来讨论FFT的快速卷积方法。 序列x(n)和h(n)的长差不多。设x(n)的长为N1，h(n)的长为N2，要求 用FFT完成这一卷积的具体步骤如下： ①为使两有限长序列的线性卷积可用其循环卷积代替而不发生混叠，必须选择循环卷积长度，若采用基2-FFT完成卷积运算，要求（为整数）。 ②用补零方法使x(n)和h(n)变成列长为N的序列。 ③用FFT计算x(n)和h(n)的N点离散傅里叶变换 ④完成X(k)和H(k)乘积， ⑤用FFT计算Y(k)的离散傅里叶反变换得 三、主要实验仪器及材料 微型计算机、Matlab6.5以上版本的编程环境。 四、实验内容 实验内容 1．编写并运行程序，观察频谱。 （1）对以下序列进行谱分析。 　　　　　　 选择FFT的变换区间N为8和16 两种情况进行频谱分析。分别打印其幅频特性曲线。 并进行对比、分析和讨论。 （2）对周期序列进行谱分析。 选择FFT的变换区间N为8和16 两种情况分别对以上序列进行频谱分析。分别打印其幅频特性曲线。并进行对比、分析和讨论。 2．数字滤波器的脉冲响应为，N2取8。用FFT计算 输入序列x(n)可选下列几种情况 ① ，N1取8。 ② ，N1取8。 当，和时，直接计算线性卷积和利用FFT计算线性卷积的时间。 要求： 1.实验前，预先计算好的值。预先编制一个应用FFT计算频谱和一个实现数字滤波器的通用程序。 2.上机独立调试，并打印或记录实验结果。 3.将实验结果与预先笔算的结果比较，验证其正确性。 五、思考题 1．x(n)长度为N）？（包括非周期信号和周期信号） 3．分析直接计算线性卷积和利用FFT计算线性卷积的时间。 六、实验报告要求 1.简述