用DFT进行频谱分析及其误差问题研究.docx

. 目 录 引言 1 利用 DFT 对有限长序列进行谱分析 1 谱分析原理 1 实验结果及分析 2 利用 DFT 对周期序列进行谱分析 2 谱分析原理 2 实验结果及分析 3 利用 DFT 对连续时间非周期信号进行谱分析 4 谱分析原理 4 实验结果及分析 5 利用 DFS 对连续时间周期信号进行谱分析 5 谱分析原理 5 实验结果及分析 6 利用 DFT 进行谱分析的误差问题及其参数选择 7 谱分析的误差分析 7 谱分析的近似性问题 7 谱分析的参数选择 8 利用 DFT 进行谱分析的误差仿真 9 混叠效应仿真 9 栅栏效应仿真 9 频谱泄露效应仿真 10 结束语 14 参考文献 15 致谢 16 . 引言 随着信息时代和数字世界的到来，数字信号处理己成为当今一门极其重要的学科和技术领域，数字信号处理在通信、语音、图像、自动控制、医疗和家用电器等众多领域得到了广泛的应用。任意一个信号都具有时域与频域特性，信号的频谱完全代表了信号， 因而研究信号的频谱就等于研究信号本身。通常从频域角度对信号进行分析与处理，容易对信号的特性获得深入的了解。因此，信号的频谱分析是数字信号处理技术中的一种较为重要的工具。[1] 众所周知，傅里叶变换和 Z 变换是信号处理中常用的重要数学变换。对于有限长序列，还有一种更加重要的数学变换即离散傅里叶变换（Discrete Fourier Transform，DFT）。 DFT[2]之所以重要，是因为其实质是有限长序列傅里叶变换的有限点离散采样，从而实现了频域离散化，使得数字处理可以在频域采用数值运算的方法进行，这样就大大加大了数字信号处理的灵活性。 信号的频谱分析的实质，就是通过信号的傅立叶变换(FT)来分析信号的频谱结构， 信号的 FT 可以借助于 DFT 用计算机仿真方法实现。一般地，信号按时间是否连续可分为连续时间信号和离散时间信号，按周期性可分为周期信号和非周期信号，在时域内信号可分为 4 大类：离散非周期信号(有限长序列)、离散周期信号(周期序列)、连续非周期信号(一般模拟信号)、连续周期信号。 利用 DFT 对有限长序列进行谱分析 谱分析原理 假设 x(n)为长度为 L 的有限长序列，其 FT 和 N 点 DFT 分别为 ? ?L ?1 ? ?X (e j ) ? ? ?L ?1 ? ? n ? 0 NX (k) ? ?N ?1 x(n)e? j 2?kn N （k=0，1…，-1） (2) n ?0 1 对比式(1)，(2)可得， X (k) = X (e j? ) | ?? 2? k N ，即 X (k ) 是在区间[0，2? ]上对 X (e j?) 的 N 点等间隔采样。因此序列的 FT 可以通过 DFT 近似得到。对于有限长序列，可知其 FT 是周期为2? 的连续谱，其DFT 是在区间[0， 2? ]上对其 FT 进行 N 点等间隔采样得到的离散谱。因此对于不同的变换长度 N，同一个序列的 DFT 也不同。随着 N 的增大， 其 DFT 的包络越来越接近 FT，对其频谱分析也越精确。需要注意的是在选择 DFT 的参数时，应满足 N≥L.[1] 实验结果及分析 R (n) 长度为8，前4个单位为1的有限长序列，对序列 R (n) 进行频谱分析，绘制出 4 4 其幅频特性曲线。其谱分析可以通过直接计算其 N点DFT 来近似。使用matlab仿真[3]的程序和结果如下。图 1(a)和(b)分别为R (n)的8点和64点DFT，是离散谱线。 （a） R (n) 的 8 点 DFT 频谱 （b） R (n) 的 64 点 DFT 频谱 4 4 图 1 R (n) 的 DFT 频谱 4 由仿真结果可知，比较图 1（a）、(b)随着 DFT 点数 N 的增加，其包络越来越接近序 列的 FT，即由离散频谱转换为连续谱。在对有限长序列进行谱分析时，通过适当选取DFT 的长度，达到用 DFT 进行谱分析的目的，尤其需要注意的是第 k(0≤k≤N - 1)条离散 谱线对应的 FT 的频率为 2? k 。 N 利用 DFT 对周期序列进行谱分析 谱分析原理 设 ~x N 为周期为 N 的周期序列，对于周期序列的频谱分析可分 3 步进行： 截取其一个周期对应的主值序列 x N N ?1 (n) ，对主值序列进行 N 点 DFT 得到其 2? 离散谱 x N (k ) ，即 DFT[ x N (n) ] = x N (k ) = ? x n (n)e? j N kn ，k=0，1，…，N-1。 由周期序列的离散傅里叶级数n?(0D FS)及其主值序列的 DFT 之间的对应关系， 2 可得周期序列 DFS 对应的 ~x N (n) 是 x N (k