周期LFMCW的时频分析 .pdfVIP

雷达对抗技术实验报告

实验题目：周期LFMCW的时频分析

院系：电子与信息工程学院

班级：

姓名：

学号：

指导教师：

实验时间：2012年6月

雷达对抗技术实验报告

周期LFMCW时频分析

一、实验目的

通过周期LFMCW时频分析实验，加深对线性调频信号的理解，

对非平稳信号的时频分析的基本思想与实现方法的认识，并掌握

Matlab的基本语法、基本函数的使用以及周期拓展，时频分析的方法。

二、实验原理

1.信号分析

非平稳信号是指信号的统计特征随时间变化的时变信号，其频率也是时间的

函数。线性调频信号是典型的非平稳信号。传统的傅立叶变换可求得信号的频率，

但该方法是基于信号的全局信息，并不能反映信号的局部特征，也不能反映其中

某个频率分量出现的具体时间及其变化趋势，不具备分析信号的瞬时有效性。而

瞬时频率，能给出信号的调制变化规律，具有它独特的优势和瞬时有效性。

瞬时频率作为描绘非平稳信号特征的一个重要物理量，其估计和提取一直是

非平稳信号处理中的研究热点。目前，人们已提出如瞬时自相关法、相位法、过

零点法、时频分析等多种手段和方法。

实验中采用时频分析方法。

在信号的时频分析中用的最多的就是短时傅立叶变换（STFT），短时傅立叶

变换是典型的线性时频表示。这种变换的基本思想就是用一个窗函数乘时间信

号，该窗函数的时宽足够窄，使取出的信号可以看成是平稳的，然后进行傅立叶

变换，可以反映该时宽中的频谱，如果让窗函数沿时间轴移动，可以得到信号频

谱随时间变化的规律。现对短时傅立叶变换及其性质介绍如下。它在傅里叶分析

中通过加窗来观察信号，因此，短时傅里叶变换也称加窗傅里叶变换。其表达式

为：

其中表示的复共轭，是输入信号，是窗函数。在这个变换中，起着频限的作

用，起着时限的作用。随着的变化，所确定的“时间窗”在轴上移动，使“以某

一时间间隔步进”进行分析。因此，往往被称为窗口函数，大致反映了在时刻频

率的“信号成分”相对含量。在实际应用中，有时需要研究信号能量在时频平面

中的二维分布情况，为此将短时傅立叶变换取模平方，得到二次型时频分布，称

第1页

雷达对抗技术实验报告

为短时功率或谱图。

通过谱图我们可以从整体上观测信号的频率范围以及时频分布情况。

可以看出，短时傅立叶变换用线性时频表示，它不存在交叉项：而谱图用二

次型的时频表示，如果两信号的短时傅立叶变换在时频平面的支撑区域不重叠，

仍可认为其谱图满足叠加性。

在短时傅里叶的分析中，窗函数常常起关键的作用。所加的窗函数能否正确

反映信号的时频特性(即窗函数是否具有较高的时间分辨率和频率分辨率)，与待

分析信号的平稳特性有关。为了了解窗函数的影响，假设窗函数取两种极端情况。

第一种极端情况是取，此时信号的STFT可表示为

其中表示傅立叶算子。这种情况下，STFT退化为信号的傅立叶变换，没有任

何的时间分辨率，却有最好的频率分辨率。

第二种极端情况是取，此时

STFT退化为信号，有理想的时间分辨率，但不提供任何频率分辨率。

短时傅立叶变换由于使用了一个可以移动的时间窗，使其具有一定的时间分

辨率。短时傅立叶变换的时间分辨率取决于窗函数的长度，为了提高信号的时间

分辨率，希望的长度愈短愈好。但是频域分辨率取决于窗函数的频域函数宽度，

为了提高频域分辨率，希望尽量加宽的窗口宽度，这样必然又会降低时域分辨率。

所以，时宽和带宽不可能同时达到任意小，既有任意小时宽，又有任意小带宽的

窗函数是不存在的。归根到底，局部谱的正确表示还在于窗函数的宽度与信号的

局部平稳长度相适应。在实际应用中，我们希望选择的窗函数具有很好的时间和

频率聚集性(即能量在时频平面是高度集中的)，使得能够有效地反映信号在时频

附近的“内容”，也就是的宽度应该与信号的局部平稳长度相适应。

利用STFT可以估计信号在每片短时窗内的频率得到信号的瞬时频率，该曲线

由一组时间