《数字信号处理》实验指导书课案.docxVIP

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