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

连续非周期信号频谱分析及 MATLAB 实现 谢海霞 1，孙志雄 1 （琼州学院海南 三亚 572022） 摘要：为了便于计算机辅助计算复杂的连续信号频谱，常常采用DFT方法。DFT不仅能反映信号的频域特征更便于用计算机处理。这里先对连续非周期信号做离散化处理，然后截短得到有限长序列，最后做DFT变换。针对常用信号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 DFTmethod is normally used. DFT can not only reflect the characteristic of the signal in the frequency domain is moreconvenient for processing by computer. Here do discretization continuous-time aperiodic signal first then cut shortby finite length sequences，finally DFT transform.This paper has made a deep discuss in the principle and somerelated questions of DFT in signal spectrum analysis. Matlab simulation software is used to analyze and test thedesign. Key words: Discrete Fourier transform（DFT） spectrum analysisMATLAB 算法，便于用计算机处理。4-6本文介绍具体的连0 引 言 续非周期信号结合 MATLAB 软件来分析其频谱。 频谱分析在数字信号处理中用途广泛：如滤波、检测等方面，这些都需要 DFT（Discrete Fourier 1 连续非周期信号谱 DFT 分析Transform）运算。1-3信号的 Fourier 变换建立了 在已知连续信号数学解析式的情况下，非周期信号的时域与频域之间的一一对应关系，信号的 信号的频谱可以根据 Fourier 变换的定义进行解析Fourier 变换可以借助 DFT 来分析。有限长序列的 计算。实际应用中的多数信号不存在数学解析式，DFT 可以由数字方法直接计算，有 DFT 存在快速 信号的频谱无法利用傅里叶分析公式方法直接计 算，一般需采用数值方法进行近似计算分析频谱，基金项目：海南省自然科学基金项目611133。作者简介：谢海霞（1978-），女，讲师，工学硕士，主要研 在进行数字计算时，需对计算的连续变量进行离散究 EDA 技术在信号处理中的应用。化。由于连续非周期信号 xt的频谱函数 Xjw是连 续信号一般都不是带限信号，连续信号在抽样前通续函数，因此，需要对其进行离散化处理 xn以近 常都经过一个模拟低通滤波器（称为抗混叠滤波似分析相应的频谱。有限长序列的 DFT 可以由数 器）进行低通滤波，以减少混叠误差，以提高频谱字方法直接计算，且 DFT 存在快速算法。通过建 分析精度。 2.2 泄漏现象立序列 xn有离散 Fourier 变换 Xm与连续非周期 对连续非周期信号的采样序列 xn进行 DFT信号 xt的 Fourier 变换 Xjw之间的关系，可以利 运算时，时间长度总是取有限值，在将信号截短即用 DFT 对连续非周期信号频谱进行近似分析，此 时域加窗处理的过程中，出现了分散的扩展谱线的近似分析过程中一般将会出现三种现象：混叠现 现象，称为频谱泄漏。对离散序列的加窗实际上是象、泄漏现象和栅栏现象。7-8这些现象与应用中 将离散序列与窗函数相乘，加窗后信号的频谱是加信号和 DFT 的参数选择有关。下面分别讨论近似 窗前信号的频谱与窗函数频谱的卷积，造成截短后过程中可能出现的问题及其解决方法。 信号的频谱与截短前信号的频谱不同，所得的频谱2 DFT 分析过程中出现的若干问 在原来没有频