第十章信号时频分析.pptVIP

本章结束 青岛大学机电工程学院 第十章 信号的时域分析 第十章 信号的时频分析 ▲ 问题：随着信号分析应用领域的不断扩大，促使人们在理论上、方法上向更深层次探索，要求研究非高斯、非平稳信号的处理，以及时变、非因果、非最小相位、非线性系统的分析，这已经成为信号分析中的一个热点。 ▲ 什么是时频分析法？ 对非平稳或时变信号的分析方法。 ▲ 时频分析的作用与特点 作用：分析非平稳信号。对于非平稳信号x(t)随时间的任何突变都会传遍X(f)的整个频率轴，如果仍采用FT，在单一域内分析就不够了。FT的不足：多数的时间信息在频域不容易得到；幅值谱只显示任一频率在信号中的强度，通常不能提供有关谱分量的时间局域化的信息。对于平稳信号的分析，用FT亦足够了。 主要特点：由时间坐标和频率坐标组成相平面，即可得到信号在局部时域内的频率组成，也可看出信号各个频带在局部时间上的分布和排列情况。 应用领域：语音处理、地震资料分析、信号检测、数据压缩、状态监测和故障诊断等领域。 ▲ 本章内容 时频分析分线性和非线性两类，线性时频分析中最重要的是短时傅里叶变换和小波变换。本章主要介绍这两种变换的基本概念。 一、短时傅里叶变换（STFT）的基本原理 1、基本原理 第一节 短时傅里叶变换 目的：为了克服FT没有时间分辨率的缺陷。通常采用两种方法来适应分析非平稳信号的需要。方法一：在FT中引入时间相关性而又保持线性不变。其思想是引入一个“局部频率”参数（在某时间内局部）。这样，“局部”FT便是通过一个窗口来观察信号，在这个窗口内信号接近平稳，见下图。通过移动时窗进行分段取样，将整个信号化为若干个局部“平稳”的信号；对这“平稳”信号进行FT，可得到一组原信号的“局部”频谱。方法二：将FT中所用的正弦基函数改为在时间上更集中而在频率上较分散的基函数。 2、定义与性质 信号x(t’)的短时傅里叶变换或短时谱定义为 t x(t’) t’ 以t为中心的“分析窗” ▲ 短时傅里叶变换（或称短时谱）——信号x(t’)在时间t的短时FT就是x(t’)乘上一个以t为中心的“分析窗” （——表示复共轭），所作的FT。x(t’)叫做基信号。 分析窗的选择：乘 等价于取出信号在分析点t’=t附近的一个切片，所以STFT是x(t’)在t附近的“局部谱”。对x(t’)的STFT很大程度上取决于分析窗 。作为分析窗的时窗函数 应该是紧支集（即为有限长度）的或很快趋近于零的函数。满足此条件的“最优”窗函数是任一高斯函数，其他一些函数也可作为窗函数，如矩形窗、三角窗、Hanning窗和Hamming窗等。 为了对 的位置和大小进行定量描述，定义其中心th和宽度半径Δth为 ▲ STFT(t, f)描述了x(t’)在时窗[t’+th-Δth, t’+th+Δth]内的局部信息。 如果x(t’)和 的FT为X(f’)和 ，则可得到STFT的另一种表示式 ej2πt(f ’-f) 可以认为，STFT是“加窗谱” 的IFT，窗谱Г(f)是时间窗γ(t)的FT. STFT的时域表达 式(10-1) STFT的频域表达 式(10-4) 根据用模拟滤波方式求取频率特征的原理，式(10-4)中的STFT可以采用让信号x(t’)通过一频率响应为 的滤波器得出，参见下图。 (a) 用带通滤波器的实现方式。在某个给定的频率f处，x(t’)先通过一中心频率为f的“分析带通滤波器”，然后再将滤波器的输出频移至零频率。由于最后的频移， 作为t的函数，对任意固定的f都是一低通信号。 (b) 采用低通滤波器的实现方法。 它与“带通方式”是等价的。而低通滤波器的冲激响应等于时间反向的分析窗 (a) (b) *共轭 依据FT的时延和频移特性， 对于短时傅里叶变换受基信号变换的影响，根据式(10-1)不难证明： 及 表明STFT保持了信号x(t’)的频移特性 表明STFT在一调制(相位因子)范围内可保存时间移位 二、短时傅里叶变换的时间-频率分辨率 一维信号x(t) STFT 映射成在时间-频率平面(t，f)内的二维函数 根据不确定性原理或Heisenberg不等式，如果用Δt和Δf分别表示时间分辨率和频率分辨率，则它们的乘积满足不等式 面临的问题：高的时间分辨率要求分析窗 尽量窄，而高的频率分辨率要求滤波器的带宽尽可能窄，即所对应的分析窗则尽量宽。二者相矛盾。 表明时间分辨率和频率分辨率不可兼顾 1、理想的时间分辨率 分析窗为无穷窄，可选择δ(t)函数作为窗函数，则 2、理想的频率分辨