资源描述：
《信号时频分析作业》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、短时傅里叶变换(STFT)算法研究与仿真实现学号:姓名:短时傅里叶变换(STFT)算法研究与仿真实现摘要：本文首先介绍了时频分析的发展，然后主要介绍了一种时频变换技术——短时傅里叶变换(STFT)的基本原理和特点，并应用短时傅里叶变换的方法对仿真信号进行了时频分析。最后，介绍了时频变换在雷达方面的具体应用。关键词:时频分析，STFT引言傅里叶变换是应用最广泛的信号分析工具Z-,其基本观点是一个任意信号总是可以分解成一组不同频率的正弦信号，即实质上是将信号投影为一组基两数的过程，每一个基函数是频率固定的正弦波
2、，投影的结果形成了原始信号的傅里叶变换，它在一个特定频率的值是信号与该频率正弦基相似性的度量，因此，信号的频率特性可以通过傅里叶变换表现出来。现实世界中许多信号的频率是随时间变换的，在这种情况下，利用简单的正弦波作为基函数并且通过频谱來描述信号不总是最好的办法，时频变换就是为了描述信号的时变频率分量而发展起来的。时间信号的时频表示开始TGabor,称为短时傅里叶变换(STFT)o它是一个移动窗口傅里叶变换，通过移动时间窗口来分析信号频率分量，这样得到一个二维的时频分布，称为谱图，谱图包含了信号在不同时间的频
3、率信息。时频变换主要分为两类：线性时频变换和双线性变换。本文主要讨论的STFT是线性时频变换，而双线性变换的典型算法是Wigner-Ville分布(WVD)。本文主要研究了短时傅里叶变换(STFT)的基木原理和特点，并应用短时傅里叶变换的方法对仿真信号进行了时频分析。最后，介绍了时频分析在雷达方面的具体应用。2短时傅里叶变换(STFT)2.1连续信号的STFT分析时变频率分量信号的一种标准方法是把时间信号分成许多段，然后对每傅里叶变换，即为STFT操作，信号兀(J的短时傅里叶变换定义如下：S777(/,O
4、)=耳心)〃*(一/)严氏(1-1)短时傅里叶变换与傅里叶变换唯一的区别就是给泄了一个窗函数“(r)去截取x(0,对截下来的局部信号做傅里叶变换，即可得到/时刻的该段信号的傅里叶变换。由于窗函数〃(f)的存在使短时傅里叶变换具有了局部特性，它既是时间的因数，又是频率的函数。对给定的时间/,5TFTv(r,O)nJ以看作是该时刻的频谱。为了提高短时傅里叶变换的时间分辨率，需要选择尽可能短的窗两数〃(小另一方面要得到高的频率分辨率，要求选择的时间的窗函数"(/)的时间宽度尽可能的长，因此与时间分辨率的提高相矛盾
5、。对于非平稳信号，利用短时傅里叶变换方法很难找到一个合适的时间窗口来适应不同的时间段，这是它最大的不足Z处。与其它的时频分布(如Wigner分布)的方法相比，基于短时傅里叶变换的微多普勒有速度快、算法简单、易实现等特点。2.2离散信号的STFT实际应用小要实现一个信号的STFT,必须对该信号进行离散化，且为有限长。设采样后的信号为%(n),斤=0,1,…，厶-1,对应式((1-1)有S7TT.Y仙”")二工x(^)g(n-mN卜-伽=、g(n_mN)卫®、(1-2)n式子屮N是在时间窗函数移动的步长；cd=
6、Q>Ts是圆周频率，Ts为由尤(『)得到兀(n)的采样间隔。式(1-2)对应傅里叶变换中的DTFT,即时间是离散的，频率是连续的。为了在工程屮实现，还应将血离散化，令C0=—k,则5m*_・—knSTF7（加®）=》x（町g（n-mN）e"（1-3）n上式将频率域的一个周期2龙分成了M个点，显然，它是一个M个点的DFTo若函数g（〃）的宽度正好也是M点，那么上式可以写成M-1拿=f〃）g（〃—wN）W；W，k=0,l,・・・,M—l（1-4）rt=O若g（町的宽度小于M,那么可以将其补零，使之变为M；若g
7、（町的宽度大于M,则应增人M使之等于窗函数的宽度。总之，上式为一标准的DFT,时域频域的长度都是M。上式屮，N的大小决定了窗函数沿时间轴移动的间距，N越小，加的取值越多，得到的时频曲线越密。若N=,即窗函数在x（n）的方向上每个一个点移动一次，这样按照上式共做厶/N个M点的DFTo上式的反变换是1A/-1x（n）=工工S77T0/a）wh（1・5）式中，血的求和取值范用取决于数据的长度厶及窗函数移动的步长2。3仿真实现以线性调频信号和正弦调频信号为例，仿真结果如下所示。3.1线性调频（LFM）信号3.1.
8、1单个信号分量信号的参数为：调频斜率K=500Hz/s,时长Is,带宽500Hz,屮心频率200Hz,采样频率1600HZSpectrumofSignal4035302520等151050・5-10-800-600-400-2000200400600800frequency(Hz)图2-1单个LFM信号的频谱-800STFToftheSignalo40-o60-oo20■2040其STFT时频分布图如下：图2-2
