[發明專利]一種頻率—頻率平面能量分布圖的實現方法

有效

• 申請號: 201710978474.7
• 申請日: 2017-10-18
• 公開（公告）號: CN107967239B
• 公開（公告）日: 2021-08-31
• 發明（設計）人: 段晨東;徐先峰;張彥寧;高強;梁棟;李光輝
• 申請（專利權）人: 長安大學
• 主分類號: G06F17/14
• 分類號: G06F17/14
• 代理公司: 西安恒泰知識產權代理事務所 61216
• 代理人: 李鄭建

1.一種頻率—頻率平面能量分布圖的實現方法，其特征在于，首先用頻率切片小波變換對信號進行時頻分解，得到信號的時頻分解矩陣，然后，通過對所有變換矩陣元素求模運算，以矩陣元素的模作為信號在對應時間頻率點的幅值，由此得到信號的時頻幅值矩陣；然后，針對頻率空間，在時頻空間上求取每個頻率分量所對應時頻幅值序列的傅里葉變換，以獲取該頻率分量所對應時頻幅值序列的頻譜序列，最后，由時頻空間上所有頻率分量的時頻幅值的頻譜序列構成一個信號的頻率—頻率平面能量分布矩陣，把頻率—頻率平面能量分布矩陣按頻率——頻率——數值的形式投射到2維、3維色度圖上，即得到頻率—頻率平面能量分布圖；

具體按下列步驟實施：

(1)去除信號的直流分量；

(2)選擇頻率切片函數；

(3)對去除直流分量的信號做頻率切片小波變換得到其時頻分解矩陣；

式中κ0，為的共軛函數,是s(t)傅里葉變換；

(4)對時頻分解矩陣的復數元素做模運算，其運算結果作為信號在對應時間頻率點(t,ω)的幅值，得到信號的時頻幅值矩陣；設κ為任意給定值時，由于ω＝2πf，把信號s(t)的頻率切片小波變換時頻分解W(t,ω,κ)寫為W_κ(t,f)；

(5)以時頻幅值矩陣為基礎，按頻率方向，針對每個頻率分量提取它所對應的時間方向上的幅值序列，再對該幅值序列做傅里葉變換，

然后由此計算出該幅值序列的頻譜序列；最后，由時頻空間上所有的頻率分量的頻譜序列構成頻率—頻率平面能量分布矩陣：

FW(θ,f)＝{FW(θ,f),θ＝0～f_s/2,f＝0～f_s/2}

(6)把頻率—頻率平面能量分布矩陣按頻率——頻率——數值的形式投射到2維、3維色度圖上，即得到頻率—頻率平面能量分布圖。