离散信号与系统的频谱分析实验报告.docx

PAGE PAGE # 实验二 离散信号与系统的频谱分析 一、实验目的 掌握离散傅里叶变换（ DFT ）及快速傅里叶变换（ FFT）的计算机实现方 法。 检验序列 DFT 的性质。 掌握利用 DFT （FFT）计算序列线性卷积的方法。 学习用 DFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现 的分析误差，以便在实际中正确应用 DFT。 了解采样频率对谱分析的影响。 了解利用 FFT进行语音信号分析的方法。 二、实验设备 1.计算机 2.Matlab 软件 7.0 以上版本。 三、实验内容 对不同序列进行离散傅里叶变换并进行分析； DFT 共轭对称性质的应用 （通过1次N点FFT计算2个N点实序列的 DFT）。 线性卷积及循环卷积的关系，以及利用 DFT（ FFT）进行线性卷积的方法。 比较计算序列的 DFT 和FFT的运算时间。 利用FFT实现带噪信号检测。 利用FFT计算信号频谱及功率谱。 扩展部分主要是关于离散系统采样频率、时域持续时间、谱分辨率等参数 之间的关系，频谱的内插恢复，对语音信号进行简单分析。 四、实验原理 序列的离散傅里叶变换及性质 离散傅里叶变换的定义：X(k) DFT [ x( n)] 离散傅里叶变换的定义： X(k) DFT [ x( n)] x(n)e j N nk ,0 n0 N1 离散傅里叶变换的性质： （1）DFT的共轭对称性。若 x（n） xep（n） xop （n） ， X（k） DFT x（n） ，则： DFT [ xep( n)] XR(k) ， DFT[xop(n)] jXI(k)。 实序列DFT的性质。若 x(n) 为实序列，则其离散傅里叶变换 X(k)为共 轭对称，即 X(k) X *(N k),0 k N 1。 实偶序列DFT的性质。若 x(n) 为实偶序列，则其离散傅里叶变换 X(k) 为实偶对称，即 X(k) X(N k),0 k N 1。 实奇序列DFT的性质。若 x(n) 为实奇序列，则其离散傅里叶变换 X(k) 为纯虚奇对称，即 X(k) X(N k),0 k N 1。 利用DFT计算线性卷积 时域循环卷积定理：设有限长序列为 x1(n)和 x2 (n) ，它们的N点DFT分别为 X1(k) 和 X2(k) ，如果 X(k) X1(k) X2(k) ，则其 IDFT为两序列的循环卷积 N1 x(n) IDFT[ X(k)] x1(m)x2((n m))N RN(n)。 m0 循环卷积和线性卷积的关系： 循环卷积是线性卷积周期延拓的主值序列。 当 循环卷积的点数大于等于线性卷积的长度， 这时循环卷积等于线性卷积， 就可以 利用 DFT计算线性卷积，方法如图 4-1所示。 x1(n) 补L-N L点 个零点 DFT y(n) L点IDFT x2n) 补L-M L点 个零点 DFT 图 4-1利用 DFT 计算线性卷积流程图 3、利用 DFT对信号进行谱分析 (1)用DFT进行谱分析的参数选择 用DFT进行连续信号谱分析时参数的选择要遵循一定的原则。 fs 2fcTP fs 2fc TP TP min ，T 1fs 1 2fc (2)用DFT进行谱分析的误差问题 由于 DFT计算频谱仅限于离散点上的频谱， 而采样点之间的频谱函数值是不 知道的。因此会引起栅栏效应，解决方法是：对有限长序列在原数据末尾补 0， 对无限长序列增大截取长度和 DFT 变换区间长度。 由于实际中信号可能是无限长 的，而DFT处理的是有限长序列，因此实际中经常要把观测的时域信号截断，相 当于时域加窗， 会引起截断效应， 表现为频谱泄露和谱间干扰， 解决的方法是增 加窗长、改变窗函数形状。 栅栏效应和谱分辨率是不同的概念，对信号补 0，可 以减小栅栏效应，但是由于截断已经使频谱变模糊了，所以补零后虽采样间隔 变小，但得到的频谱包络仍是变模糊的频谱，因此频率分辨率没有提高，要提 高频率分辨率，必须使时域截取长度增加。 DFT和FFT的运算量 离散傅里叶变换在实际应用中非常重要， 利用它可以计算信号的频谱、 功率 谱和线性卷积等。但其计算量太大 (与N的平方成正比), 很难实时地处理问题 , 因 此引出了快速傅里叶变换 (FFT) 。 设序列的长度为 N，则计算 DFT和计算 FFT的运算量如表 4.1所示。 表4.1 DFT 和FFT运算量比较 计算 DFT 计算 FFT 复数乘法次数 N2 N(N 1) 复数加法次数 N 2 log 2 N N log 2 N 周期图法计算信号功率谱 周期图法是一种估计信号功率谱密度的方法。因为序列的 DFT 隐含有周期 性，所以这种功率谱也有周期性， 称为周期图。 该方法的优点是能应用离散傅里 叶变换的快速算法来进行估值。定义如式( 4-1)所