连续非周期信号频谱分析及Matlab 实现.docVIP

连续非周期信号频谱分析及MATLAB实现 DFT谱分析的原理及谱分析中的相关问题进行了较为深入的探讨，并结合实例用Matlab仿真软件进行了分析和验证。 关键词：DFT；频谱分析；MATLAB 中图分类号：TN911.72 文献标识码：A Continuous-Time Aperiodic Signal Spectrum Analysis Based on MATLAB XIE Haixia1, SUN Zhixiong1 （Qiongzhou University, sanya Hainan 572022, China） Abstract: In order to calculate the spectrum of complex continuous signal with the computer, the DFT method is normally used. DFT can not only reflect the characteristic of the signal in the frequency domain is more convenient for processing by computer. Here do discretization continuous-time aperiodic signal first, then cut short by finite length sequences，finally, DFT transform.This paper has made a deep discuss in the principle and some related questions of DFT in signal spectrum analysis. Matlab simulation software is used to analyze and test the design. Key words: Discrete Fourier transform（DFT）; spectrum analysis;MATLAB 0 引 ( 言 频谱分析在数字信号处理中用途广泛：如滤波、检测等方面，这些都需要DFT（Discrete Fourier Transform）运算。[1-3]信号的Fourier变换建立了信号的时域与频域之间的一一对应关系，信号的Fourier变换可以借助DFT来分析。有限长序列的DFT可以由数字方法直接计算，有DFT存在快速算法，便于用计算机处理。[4-6]本文介绍具体的连续非周期信号结合MATLAB软件来分析其频谱。 1 连续非周期信号谱DFT分析 在已知连续信号数学解析式的情况下，非周期信号的频谱可以根据Fourier变换的定义进行解析计算。实际应用中的多数信号不存在数学解析式，信号的频谱无法利用傅里叶分析公式方法直接计算，一般需采用数值方法进行近似计算分析频谱，在进行数字计算时，需对计算的连续变量进行离散化。由于连续非周期信号x(t)的频谱函数X(jw)是连续函数，因此，需要对其进行离散化处理x[n]以近似分析相应的频谱。有限长序列的DFT可以由数字方法直接计算，且DFT存在快速算法。通过建立序列x[n]有离散Fourier变换X[m]与连续非周期信号x(t)的Fourier变换X(jw)之间的关系，可以利用DFT对连续非周期信号频谱进行近似分析，此近似分析过程中一般将会出现三种现象：混叠现象、泄漏现象和栅栏现象。[7-8]这些现象与应用中信号和DFT的参数选择有关。下面分别讨论近似过程中可能出现的问题及其解决方法。 2 DFT分析过程中出现的若干问题 2.1混叠现象 由DFT计算出的频谱是信号x(t)的频谱X(jw)周期化的抽样值，如果连续信号不是带限信号，或者抽样频率不满足抽样定理，在连续信号离散化时，就会出现信号频谱的混叠。解决连续信号离散化过程中的频谱混叠主要有两种方法：对于带限连续信号，只要提高抽样频率使之满足时域抽样定理：对于非带限连续信号，可根据实际信号对其进行低通滤波，使之成为带限信号。工程实际中的连续信号一般都不是带限信号，连续信号在抽样前通常都经过一个模拟低通滤波器（称为抗混叠滤波器）进行低通滤波，以减少混叠误差，以提高频谱分析精度。 2.2泄漏现象 对连续非周期信号的采样序列x[n]进行DFT运算时，时间长度总是取有限值，在将信号截短即时域加窗处理的过程中，出现了分散的扩展谱线的现象，称为频谱泄漏。对离散序列的加窗实际上是将离散序列与窗函数相乘，加窗后信号的频谱是加窗前信号的频谱与窗函数频谱的卷积，造成截短后信号的频谱与截短前信号的频谱不同，所得的频谱在原来没有频谱的区间出现了频谱。原来比较尖锐