DFT栅栏效应分析.docVIP

西山大学 课 程 设 计 说 明 书 题目 :DFT栅栏效应分析 目录 一、摘要……………………………………………………4 二、设计目的及分析………………………………………5 三、理论原理知识…………………………………………6 3.1 DFT频谱分析原理…………………………………6 3.2 栅栏效应…………………………………………8 3.3 分辨率……………………………………………10 四、MATLAB编程…………………………………………12 4.1 MATLAB软件简介…………………………………12 4.2 标点符号与语句…………………………………13 4.3 仿真程序…………………………………………13 4.4 仿真结果与分析…………………………………14 五、心得体会……………………………………………17 六、参考文献……………………………………………18 一、摘要 DFT是在时域和频域上都已离散的傅里叶变换，适用于数值计算且有快速算法，是利用计算机实现信号频谱分析的常用工具。本文介绍了利用DFT分析信号频谱的基本流程，重点阐述了频谱分析过程中误差形成的原因及减小分析误差的主要措施。实例例举了MATLAB环境下频谱分析的实现程序。通过与理论分析的对比，解释了利用DFT分析信号频谱时存在的频谱混叠，频谱泄露及栅栏效应，并提出相应的改进方法。在满足时频采样的的条件下，可以通过对连续信号进行采样并进行DFT，来近似的反映连续信号的频谱特性。由于DFT变换时需要进行时频采样和频域采样，因此这种近似必然带来频谱分析的一定误差。本文主要对栅栏效应进行分析和减小。 二、设计目的及分析 技术参数：信号中包含三种频率成分，分别是20HZ，20.5HZ，40HZ，采样频率为100HZ。为了把三种频率分辨出来，对其进行栅栏效应分析。 首先，求出最小记录点数，易知是fs/（20.5-20） 200，因此当频域采样点数N 200时，不出现栅栏效应，而当N 200时，会有栅栏效应误差出现。为了更好的分析DFT栅栏效应，DFT分三种情况：在128点有效数据不补零情况下的分辨率，在128点有效数据且补零至512点情况下分辨率，在512点有效数据下分辨率。然后比较三次仿真结果的异同，进而对其进行比较分析。 三、理论原理知识 3.1 DFT频谱分析原理 所谓的信号的谱分析，就是计算信号的傅里叶变换。工程实际中，经常遇到的是连续信号Xa（t），截取一段进行FT变换后其谱度函数Xa（jΩ）也是连续函数，因此起计算过程不便用于计算机实现。而DFT作为一种时域和频域均离散化的变换，适合数值运算，成为分析离散信号和系统的有力工具。对于连续信号和系统，可以通过Xa（t）时域采样，得到离散时间序列x（n），再对x（n）进行DFT，因此是一种近似的谱分析。 一个连续信号Xa（t）的频谱可以用它的傅里叶变换表示成： 如果对该信号进行理想采样，可以得到采样序列： 同样，可对该序列进行Z变换，其中T为采样周期： 当时，我们就得到了序列的傅里叶变换： 其中为数字频率，它和模拟频率的关系为： 式中为采样频率。上式说明数字频率是模拟频率对的归一化。同模拟域的情况相似，数字频率代表序列值变化的速率，而序列的傅里叶变换称为序列的频谱。序列的傅里叶变换和对应的采样信号频谱有下列对应关系： 即序列的频谱是采样信号频谱的周期性严拓。从式中可以看出，只要分析采样序列的频谱，就可以得到相应的连续信号的频谱。注意：这里的信号必须是带限信号，采样也必须满足Nyquist定理。 在各种信号序列中，有限长序列在信号处理中占有重要地位。无限长序列也往往可以用有限长序列逼近，对于有限长的序列我们可以用离散的傅里叶变换DFT,这一变换可以很好的反应序列的频域特性。并且容易利用快速算法在计算机上实现当序列的长度是N时，我们定义序列离散傅里叶变换为： 其中：， 它的反变换定义为： 令则有： 可以得到：，是Z平面单位圆上幅角为的点 就是将单位圆进行N等份的第K个点，所以，X K 是Z变换在单位圆上的等距采样。时域采样在满足Nyquist定理时，就不会发生频谱混淆：同样，在频域进行采样的时候，只要采样间隔足够小，就不会发生混叠。DFT是对序列傅里叶变换的等距采样，因此可以用于序列的频谱分析。在运用DFT进行频谱分析的时候可能存在栅栏效应误差。 3.2栅栏效应 （1）栅栏效应的定义 由于为非周期的连续信号，它的频谱是连续的，但将采样，截断然后进行DFT分析时，得到的仅仅是连续信号频谱上的有限个点，而有一部分频谱分量将被挡住，好像是通过栅栏观察频谱，这种现象称为栅栏效应。 DFT是对单位圆上Z变换的均匀采样，是一些离散的点，所以它不可能将频谱视为一个连续函数，就一定意义上看，用DFT来观察频谱就好像通过一个