一種基于參數(shù)優(yōu)化的改進廣義s變換的時頻分析方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明屬于數(shù)字信號處理領(lǐng)域，特別涉及到非平穩(wěn)信號的分析方法，基于參數(shù)優(yōu) 化的改進廣義S變換的時頻分析方法，可用于通信、雷達、地震和生物醫(yī)學(xué)信號的分析與處 理。
【背景技術(shù)】
[0002] 在實際工程中廣泛存在著各種類型的非平穩(wěn)信號，非平穩(wěn)信號的分析與處理是數(shù) 字信號處理領(lǐng)域的一個研究熱點。與平穩(wěn)信號不同，非平穩(wěn)信號的頻率隨時間變化，因此對 其分析需要采用時域和頻域聯(lián)合的二維時頻分析方法。
[0003] 常用的時頻分析方法有短時傅里葉變換（STFT)、Wigner-Ville分布（WVD)、小波 變換（WT)和S變換（ST)。短時傅里葉變換的窗函數(shù)固定，在不同頻段的頻率分辨率固定， 時頻聚集性較差；Wigner-Ville分布具有良好的時頻聚集性，但存在較強的交叉項干擾； 小波變換尺度與信號頻率沒有直接的對應(yīng)關(guān)系，在用于信號頻率分析時物理意義不明確； S變換結(jié)合了短時傅里葉變換和小波變換的優(yōu)點，其窗寬度隨頻率呈反向變化，其基本小波 不用滿足容許性條件，與信號頻譜直接對應(yīng)，具有無損可逆性。然而，由于S變換中窗函數(shù) 的形式固定，導(dǎo)致信號高頻部分的時頻聚集性不是很理想，使其在應(yīng)用中受到限制。
[0004] 為了改善S變換的時頻聚集性，許多學(xué)者對其窗函數(shù)進行了改進，提出了不同類 型的廣義S變換。例如Manshiha等人在1997年提出了用f/r替換窗函數(shù)中參數(shù)f，可適當 調(diào)節(jié)窗函數(shù)形式，其計算簡單，但調(diào)節(jié)能力有限；Pinnegar等人在2003年提出了可窗函數(shù) 標準差可調(diào)，且窗口不對稱的雙曲線窗函數(shù)，在其高頻段時窗較窄，選用對稱性較好的窗函 數(shù)，達到提高頻率分辨率的目的；低頻段時窗較寬，選用不對稱的窗函數(shù)，但其中涉及到兩 個參數(shù)的雙曲線方程的計算，較為復(fù)雜；陳學(xué)華在2005年提出將兩個參數(shù)λ和p引入窗函 數(shù)，利用λ|f|ρ替換參數(shù)f，進一步增加了窗函數(shù)調(diào)節(jié)的靈活性，但其中涉及到頻率f的指 數(shù)運算，運算量也較大。
[0005] 上述廣義S變換中的改進窗函數(shù)各有特點，其參數(shù)調(diào)節(jié)越靈活，計算量也越大，且 參數(shù)的選擇受測不準原理的限制，即時間分辨率和頻率分辨率不能同時達到最優(yōu)，需要折 中考慮。因此，如何以較小的計算量對窗函數(shù)進行改進并進行參數(shù)自動優(yōu)化，對于擴展S變 換的應(yīng)用范圍，提高其時頻分析性能具有重要意義。

【發(fā)明內(nèi)容】

[0006] 本發(fā)明要解決的技術(shù)問題是，針對傳統(tǒng)S變換存在的上述不足，提出一種基于參 數(shù)優(yōu)化的改進廣義S變換的時頻分析方法，通過時頻聚集度對一階函數(shù)的參數(shù)進行優(yōu)化， 以較小的計算量改善信號整體的時頻分辨率。
[0007] 本發(fā)明為解決上述技術(shù)問題所采用的技術(shù)方案是：
[0008] -種基于參數(shù)優(yōu)化的改進廣義S變換的時頻分析方法，包括以下步驟：
[0009] S1、輸入非平穩(wěn)信號x(t)，其中t為時間，非平穩(wěn)信號包括單分量和多分量混和信 號；
[0010] S2、對于含有頻率f的一階函數(shù)的改進廣義S變換的窗函數(shù)：
[0011]
CO
[0012] 其中a為乘性調(diào)節(jié)參數(shù)，b為加性調(diào)節(jié)參數(shù)，a> 0,b彡0,確定參數(shù)a和b的取值 范圍RjPRb;
[0013] S3、對RjPRb內(nèi)一組參數(shù)（a，b)，利用下式計算輸入非平穩(wěn)信號x(t)的時頻分布：
[0014]
(2)
[0015] 其中GS^Xr,/)為信號的改進廣義S變換，式中t和τ表示時間，τ為窗函數(shù)的 中心點，控制著窗函數(shù)在時間軸上的位置，f表示頻率，t、τ和f均為實數(shù)；
[0016] S4、對于步驟S3中參數(shù)（a，b)對應(yīng)的GS^r,/)，進行能量歸一化處理：
rπC3)
[0017]
[0018] S5、利用步驟S4歸一化處理后的時頻分布表達式計算整個時頻分布的時頻聚集 度：
「 ^ .(4.).
[0019] -- .χ,
[0020] S6、重復(fù)步驟S3~S5,遍歷RjPRb取值范圍內(nèi)的所有參數(shù)（a，b)，得到一系列的 CM(a，b)，取其中最大值對應(yīng)的一組參數(shù)作為最優(yōu)化參數(shù)：
[0021]
⑴
[0022] S7、將步驟S6得到的最優(yōu)化參數(shù)hb)。#帶入式（2)，計算輸入非平穩(wěn)信號x(t) 的改進廣義S變換：
[0023]
(6)
[0024] 式（6)即為參數(shù)優(yōu)化后的改進廣義S變換。
[0025] 進一步，所述步驟S2中，改進廣義S變換源于基本S變換，信號x(t)的基本S變 換為：
[0026] (7)
[0027] 其窗函數(shù)為：
[0028]
(I)
[0029] 這里將頻率f的一階函數(shù)引入，得到改進廣義S變換的窗函數(shù)如式（1)所述，其中 a和b的取值范圍為：0 <a彡2,0彡b彡20。
[0030] 進一步，所述步驟S3中，令t=kT,τ=jT，
的離散計 算公式如下：
[0031]
[0032] 其中T為信號采樣間隔。
[0033] 進一步，所述步驟S5和S6利用整體時頻聚集度對參數(shù)進行優(yōu)化，即一次利用時頻 面所有時間頻率單元內(nèi)的值進行時頻聚集度的計算，得到的優(yōu)化參數(shù)對應(yīng)整體時頻分布的 最優(yōu)。
[0034] 本發(fā)明的工作原理：在對基本S變換和已有的廣義S變換分析的基礎(chǔ)上，在窗函數(shù) 中用頻率f的一階函數(shù)替換參數(shù)f，通過時頻聚集度對一階函數(shù)的參數(shù)進行優(yōu)化，以較小的 計算量改善信號整體的時頻分辨率。
[0035] 本發(fā)明的有益效果：
[0036] 1、本發(fā)明在窗函數(shù)中用含有頻率f的一階函數(shù)替換f，引入了兩個參數(shù)a和b，用 以調(diào)節(jié)窗函數(shù)的寬度，其中a起主要調(diào)節(jié)作用，b起細微調(diào)節(jié)作用，其計算簡便，調(diào)節(jié)靈活， 通過時頻聚集度對一階函數(shù)的參數(shù)進行優(yōu)化，以較小的計算量改善信號時頻分布的能量聚 集性（整體的時頻分辨率），尤其適用于通信、雷達、地震和生物醫(yī)學(xué)信號的分析與處理；
[0037] 2、本發(fā)明利用時頻聚集度來優(yōu)化參數(shù)a和b的值，每次采用時頻面所有時間頻率 單元的值計算時頻聚集度，以達到信號整體時頻聚集性能的最優(yōu)。
【附圖說明】
[0038] 圖1為本發(fā)明時頻分析方法的實現(xiàn)流程圖；
[0039] 圖2為實施例中信號1的時域波形圖；
[0040] 圖3為實施例中信號1經(jīng)過基本S變換后的時頻分布圖；
[0041] 圖4為實施例中信號1經(jīng)過參數(shù)優(yōu)化的改進廣義S變換后的時頻分布圖；
[0042] 圖5為實施例中信號2的時域波形圖；
[0043] 圖6為實施例中信號2經(jīng)過基本S變換后的時頻分布圖；