matlab非周期信号的频谱实现,连续非周期信号频谱分析及Matlab 实现

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

连续非周期信号频谱分析及 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 DFT 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 trans.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 trans(DFT) spectrum analysisMATLAB 算法，便于用计算机处理。4-6本文介绍具体的连0 引 言 续非周期信号结合 MATLAB 软件来分析其频谱。 频谱分析在数字信号处理中用途广泛：如滤波、检测等方面，这些都需要 DFT(Discrete Fourier 1 连续非周期信号谱 DFT 分析 Trans)运算。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 分析过程中出现的若干问 在原来没有频谱的区间出现了频谱。原来比较尖锐题 的谱峰变得比较平缓，当两个不同频率的谱峰靠得2.1 混叠现象 比较近时，可能显现不出两个明显的峰值。特 别是 由 DFT 计算出的频谱是信号 xt的频谱 Xjw 强信号谱的旁瓣可能淹没弱信号的主谱或误认为周期化的抽样值，如果连续信号不是带限信号，或 是另一假信号的主谱线。矩形窗的旁瓣幅度大，谱者抽样频率不满足抽样定理，在 连续信号离散化 间干扰严重。频谱泄漏使频谱变模糊，分辨率 变(事时，就会出 现信号频谱的混叠。解决连续信号离散 实上通常规定 DFT 的频率分辨率 为 fs/N， 为采样 fs化过程中的频谱混叠主要有两种方法：对于带限连 频率， N 是指信号 xn的有效长度)差，泄漏程度续信号，只要提高抽 样频率使之满足时域抽样定 与窗函数幅度谱主瓣宽度有关。窗型一定，窗口越理：对于非带限连续信号，可根据实际信号对其进 长，主瓣越窄，频谱泄漏越小，窗口 长度一定，矩行低通滤波，使之成为带限信号。工程实际中的连 形窗主瓣最窄，频谱泄漏最小，但其旁瓣的幅度最大。因此为了尽量减少泄漏现象，应选用旁瓣幅度 在利用 DFT 分析连续时间信号的频谱时，涉小、主瓣窄，即：“ 泄漏” 小的窗函数。相对而言，布 及频谱混叠、频率泄漏及栅栏现象。 频率混叠与连莱克曼窗的旁瓣幅度比矩形窗小，谱间干扰小，但 续信号的时域抽样间隔有关，频率泄漏与信号的时其主瓣过渡带宽，分辨率差。采 样频率或采样周期 域加窗截短的长度及窗型有关，栅栏现象与 DFT是在满足混叠误差前提下选取的，当采样频率或采 的点数有关。在大多数情况下，一般已知待分析连样周期确定后，适当增加窗口 长度有利于减小泄漏 续信号的最高频率，以及希望的 DFT 分析的频率误差。 分辨率。下面根据信号的 Fourier 变换的理论， 讨2.3 栅栏现象 论利用 DFT 进行谱分析的参数(抽样频率、持续 DFT 得到的频谱 Xm只能是连续非周期信号 时间、样点数等)选择的原则。 频谱 Xjw上的有限离散频点采样，由于 Xm是离 (a)首先确定信号抽 样频率 f s ， f s 应满足散序列，因而无法反映抽样点之间的细节，就如同 f m 为待分析 时域抽样定理，即 f s ≥ 2 f m ，其中隔着百叶窗观察窗外的景色，这种现象称为栅栏现 的连续信号的