数字信号处理,中山大学Chapter 5.ppt

第五章 时 频 分 析 ;5.1 引 言;1、引言;;2、解析信号;Hilbert变换器的传输函数为 ;Z(f)=S(f)+jH(f)S(f)=S(f)［1+jH(f)］ ;3、瞬时频率;4、不确定原理;5.2 短时傅里叶变换;图2.1.3 窗函数无限宽时STFT缺少时域定位功能 注：见胡广书《现代信号处理教程》图2.1.3;图2.1.4窗函数无限窄时STFT缺少频域定位功能 注：见胡广书《现代信号处理教程》图2.1.4; 由于受不定原理的制约，窗函数的有效时宽和带宽不可能同时任意小，窗宽应该与信号的局域平稳长度相适应。 对时间分辨率和频率分辨率只能取一个折中，一个提高了，另一个就必然要降低，反之亦然。;谱图：一般把短时傅里叶变换模的平方称为谱图，它是一种能量分布函数，不服从线性叠加原理，两个信号之和的谱图并不等于它们分别的谱图的和，还存在第三项即交叉项。;5.3 小波变换;1、 引 言 ; 短时傅立叶变换在时频平面各处的分辨率都相同，可以用时频平面的相等网格表示。; 小波基函数的包络随尺度参数的变化而变化，可以实现时频平面的多分辨率分析。;2、连续小波变换（CWT） 连续小波变换的定义 设x(t)是平方可积函数，记作 ，则x(t)的连续小波变换可以定义为：?; 定义式的说明：  （1） 基小波函数可能为复函数，例如Morlet小波的表达式为 ;（4） 在ψab前面所加的因子 的作用是保证在不同的尺度因子下的小波函数的能量保持一致。 设E=∫|ψ(t)|2 dt作为基本小波的能量，则对基本小波进行移位和伸缩后得到的ψab(t)的能量为 ; 连续小波变换的频率域表达式 在定义了连续小波变换后， 对该表达式进行傅里叶变换, ; 小波变换的特点 小波变换的时频关系受不确定原理的制约，在时频平面上的分析窗是可调的，但分析窗的面积保持不变。; 当用较小的a对信号作高频分析时，实际上是用高频小波对信号作细致观察;当用较大的a对信号作低频分析时，实际上是用低频小波对信号作概貌观察。 ;注：见张贤达《现代信号处理》图6.5.3;;5.4 Wigner-Ville分布(WVD);引言;时频分布的定义;二次叠加原理 设 ; 对于有p个分量的信号， 二次叠加原理用下式表示： ;Wigner-Ville分布的定义 ; z(t)在频率域的WVD分布定义如下: ;5.5 Cohen类时频分布;1、 模糊函数 对瞬时相关函数kz(t,τ)=z(t+τ/2)z*(t-τ/2) 关于时间t作傅里叶反变换，则得到模糊函数的时域定义为 ; 模糊函数和WVD之间的关系：; WVD是能量化的时频表示，存在时间边缘特性Pz(t)和频率边缘特性Pz(w)，公式重写如下: ; 模糊函数是相关化的时频表示，将模糊函数的定义重写如下: ;; WVD中交叉项的抑制： 对信号求模糊函数，由于模糊函数的自项始终在 平面的原点处，而交叉项远离原点，故可以设计一个二维低通滤波器，来抑制模糊函数中的交叉项； 对滤波后的模糊函数作二维傅立叶变换，得到信号的维格纳变换，此时的WVD即是抑制了交叉项的新WVD。;2、 Cohen类时频分布; 当核函数φ(τ,v)=1 时， Cohen 类时频分布将转换成WVD。 ;由于Wigner分布的核函数是全通函数，它对AF的互项无抑制作用，因此，其WD也就存在着较大的交叉项。 消除干扰项的方法：应该选择 平面上的二维低通函数来作为核函数。