肖 瑛,馮長建

(大連民族學(xué)院機(jī)電信息工程學(xué)院,遼寧大連 116600)

0 引言

傅里葉變換方法作為一種全局的線性的處理方法[1],無法滿足頻譜結(jié)構(gòu)隨時(shí)間變化的非平穩(wěn)信號的分析,用時(shí)頻聯(lián)合分析的方法來處理信號,可以揭示信號的時(shí)頻細(xì)節(jié),從而更有效地識別信號特征。目前,已經(jīng)廣泛應(yīng)用于語音、聲納、雷達(dá)、機(jī)械震動、生物醫(yī)學(xué)和地球物理信號等處理領(lǐng)域[2]。常用的時(shí)頻分析方法分為線性時(shí)頻表示和二次型時(shí)頻表示兩大類:線性時(shí)頻表示以STFT和小波變換為代表;二次型時(shí)頻表示主要包括譜圖、Cohen類時(shí)頻表示和Affine類時(shí)頻表示等。眾多時(shí)頻表示中,STFT方法計(jì)算簡單,不受二次型時(shí)頻分布固有的交叉項(xiàng)影響,特別適用于跳頻類信號的分析。但是,STFT方法分析中采用固定時(shí)間窗函數(shù),時(shí)間分辨率和頻率分辨率相互制約,影響STFT的時(shí)頻聚集性,為此,文中提出了一種組合窗函數(shù)的STFT時(shí)頻表示方法。

1 短時(shí)傅里葉變換(STFT)

STFT定義為[3]:

式中,g*(t)表示時(shí)間窗函數(shù),上標(biāo)“*”表示復(fù)共軛,從式(1)中可看出:對于特定的時(shí)刻t,STFT是信號x(t′-t)與中心在t的平滑移動窗函數(shù)g*(t′-t)乘積的傅里葉變換。從這個(gè)意義上,STFT可以理解為信號x(t′)在分析時(shí)刻t的局部頻譜。由于分析窗的作用,在分析窗以外的信號得到抑制,對應(yīng)的STFT在頻域可表示為:

比較式(1)和式(2)可看出:STFT的時(shí)域表示和頻域表示很相似。實(shí)際上,頻域中STFT的表示可以看作信號頻譜與時(shí)間窗譜乘積的反傅里葉變換。因此,STFT在時(shí)域和頻域具有同樣的局域性,能實(shí)現(xiàn)對信號的時(shí)間和頻率的局域化分析。STFT的等效低通和帶通濾波器如圖1和圖2所示。

圖1 STFT的低通實(shí)現(xiàn)Fig.1 STFT implementation by low-pass filter

圖2 STFT的帶通實(shí)現(xiàn)Fig.2 STFT implementation by band-pass filter

與其他時(shí)頻表示相比較,小波分析可根據(jù)信號頻率變化,通過時(shí)間尺度變化在時(shí)頻表示區(qū)域具有不同窗,但是小波基的選擇是一大難題。Ville-Wigner分布(WVD)已經(jīng)被證明,時(shí)間帶寬積達(dá)到了Heisenberg不確定性原理的下界[4],對單分量信號具有最佳的時(shí)頻分辨率,但對多分量信號進(jìn)行分析時(shí),由于其二次型時(shí)頻表示的本質(zhì),具有交叉項(xiàng)的影響[5],這使得WVD的時(shí)頻表示可讀性大打折扣。

2 STFT的時(shí)頻聚集性

通過STFT等效低通和帶通濾波器實(shí)現(xiàn)的分析可知,STFT在時(shí)刻t的時(shí)頻表示是由信號x(t′)通過加窗函數(shù)g*(t′-t)作傅里葉變換得到的,如果想獲得好的時(shí)間分辨率,那么希望用短的時(shí)間窗函數(shù)來刻畫時(shí)刻的信號特征,相反,STFT在頻率 f處的STFT可看作信號x(t′)通過帶通濾波器G*(f'-f)得到的,要想具有好的頻率分辨率,則希望窄帶濾波器對應(yīng)的也就是長的時(shí)間窗函數(shù),因?yàn)闀r(shí)域變短對應(yīng)頻域展寬?？梢?STFT的時(shí)頻表示在時(shí)間分辨率和頻率分辨率之間是矛盾的,存在一個(gè)基本的折衷,即得到好的時(shí)間分辨率就要犧牲頻率分辨率,反之亦然[6]。

舉例來說明,如果時(shí)間窗函數(shù)取δ(t),那么

此時(shí)STFT變成了信號x(t),保持了信號的所有時(shí)間特征,具有完美的時(shí)間分辨率,但是頻率分辨率為0,即沒有任何頻率分辨率。另外一種極端情況,取時(shí)間窗函數(shù)無窮長,即g(t)=1,則其傅里葉變換為δ(f),此時(shí)

可見,STFT變成了傅里葉變換,有最好的頻率分辨率,但是沒有任何時(shí)間分辨率。

通常用時(shí)間帶寬積來表征STFT時(shí)頻表示的時(shí)頻聚集性,定義帶通濾波器的帶寬Δf為:

定義時(shí)間窗函數(shù)的時(shí)寬Δt為

把一個(gè)信號的時(shí)寬Δt和帶寬Δf的乘積大于等于一個(gè)常數(shù)稱為不確定性原理,即Heisenberg不等式[7]

在不確定性原理約束下,時(shí)頻表示的時(shí)間分辨率和頻率分辨率不可能同時(shí)任意好。

3 組合窗函數(shù)的STFT

STFT的缺陷是對于窗函數(shù)一旦選定,則在整個(gè)分析過程中都使用相同的窗,其分辨率在時(shí)間-頻率平面上的所有局域都是相同的。選擇最佳的時(shí)間窗函數(shù)(時(shí)間帶寬積等于下限)也是非常困難的,雖然已經(jīng)有專家提出了一些改進(jìn)方法,如自適應(yīng)窗函數(shù)的STFT方法ASTFT[8],但是自適應(yīng)窗函數(shù)是否為最佳窗函數(shù)在理論上也很難證明,并且計(jì)算復(fù)雜度相當(dāng)高。

通過分析也可以發(fā)現(xiàn)STFT具有如下優(yōu)點(diǎn):若信號在給定的時(shí)間間隔[-T,T]和頻率間隔[-F,F]內(nèi)具有大多數(shù)能量,則其STFT將局域化在區(qū)域[-T,T]×[-F,F]上,而在信號沒有多少能量的時(shí)間和頻率間隔處,STFT接近為0,根據(jù)這一性質(zhì),給出一種組合窗函數(shù)的STFT時(shí)頻表示方法。

以離散STFT來說明組合窗函數(shù)的實(shí)現(xiàn)過程,離散STFT的表達(dá)式如式(8):

對應(yīng)于每一個(gè)時(shí)刻n0,STFT都是一組DFT,考慮N點(diǎn)DFT,對應(yīng)的ω可以表示成式(9):

可見:即使采用不同長度的窗函數(shù),如果N相等,最終得到的STFT時(shí)頻表示結(jié)果將有相同的維數(shù)。

設(shè)長、短窗函數(shù)的頻率分辨率和時(shí)間分辨率分別為 Δf 1、Δt1和 Δf 2、Δt2,顯然,Δf 1 ＜Δf 2,Δt1 ＞Δt2,并設(shè)信號x(n)的能量分布在[-T,T]×[-F,F]內(nèi),那么,對于較長窗函數(shù)得到的STFT時(shí)頻表示將局域化在[-T-Δt1/2,-T+Δt1/2]×[-f-Δf1/2,-T+Δf1/2]內(nèi),其他區(qū)域?yàn)榱?同理,對于短窗函數(shù)的STFT時(shí)頻表示局域化在[-T-Δt2/2,-T+Δt2/2]×[-f-Δf2/2,-T+Δf2/2]內(nèi),其他區(qū)域?yàn)榱?將兩組時(shí)頻表示點(diǎn)乘,可以得到最終的組合窗函數(shù)時(shí)頻表示局域化在[-T-Δt2/2,-T+Δt2/2]×[-f-Δf1/2,-T+Δf1/2]內(nèi),具有短時(shí)間窗函數(shù)的時(shí)間分辨率和長時(shí)間窗函數(shù)的頻率分辨率。

4 仿真驗(yàn)證

以式(10)仿真跳頻數(shù)據(jù)為例來對文中提出的方法進(jìn)行驗(yàn)證,其中,f 1=100 Hz,f 2=200 Hz,采樣頻率 fs=1 000 Hz,信號的時(shí)域波形如圖3所示。

圖3 仿真跳頻信號Fig.3 Simulation frequency-hopping signal

加1/2數(shù)據(jù)長度的Hamming窗函數(shù)得到STFT時(shí)頻表示如圖4所示,加1/8數(shù)據(jù)長度的Hamming窗函數(shù)得到STFT時(shí)頻表示如圖5所示,采用組合窗函數(shù)方法得到的時(shí)頻表示如圖6所示,準(zhǔn)確地刻畫了仿真跳頻信號包含的頻率成分以及頻率跳變時(shí)刻。從圖7中可看出WVD在兩分量信號之間明顯存在交叉項(xiàng)干擾,而圖8所示的小波譜的時(shí)頻分辨率不如文中提出的組合窗函數(shù)STFT方法,同時(shí),小波時(shí)頻表示在小波基的選取上存在著困難。

圖4 STFT(加1/2數(shù)據(jù)長度的 Hamming窗)Fig.4 STFT(with 1/2 length of data Hamming window)

圖5 STFT(加1/8數(shù)據(jù)長度的 Hamming窗)Fig.5 STFT(with 1/8 length of data Hamming window)

圖6 組合窗函數(shù)的STFTFig.6 STFT with combination of time window

圖8 小波時(shí)頻表示Fig.8 Wavelet time-frequency representation

5 試驗(yàn)數(shù)據(jù)驗(yàn)證

以采集的遙測高頻振動信號為例來對組合窗函數(shù)STFT時(shí)頻表示方法進(jìn)行驗(yàn)證。該試驗(yàn)數(shù)據(jù)按試驗(yàn)過程中的特征時(shí)刻應(yīng)有兩處明顯的頻率跳變,采用時(shí)頻分析的目的是找出頻率分布范圍和頻率跳變時(shí)刻,以確定試驗(yàn)過程是否正常。其中傳感器采樣頻率為10 k Hz。圖9為信號的時(shí)域波形(經(jīng)過了最小二乘趨項(xiàng)去除和小波軟閾值降噪),圖10和圖11為分別采用1/4和1/8數(shù)據(jù)長度的Hamming窗函數(shù)得到的STFT時(shí)頻表示,圖12為組合窗函數(shù)的STFT時(shí)頻表示,圖13為WVD時(shí)頻表示結(jié)果,圖14為小波時(shí)頻表示,比較幾種時(shí)頻表示結(jié)果可看出,組合窗函數(shù)的STFT時(shí)頻表示有效的組合了長時(shí)間窗函數(shù)和短時(shí)間窗函數(shù)對應(yīng)的頻率分辨率和時(shí)間分辨率的優(yōu)點(diǎn),具有最佳的時(shí)頻聚集性。

圖9 某傳感器采集高頻振動信號時(shí)域波形Fig.9 The high-frequency signal from one sensor

圖10 1/4數(shù)據(jù)長度Hamming窗的STFTFig.10 STFT(with 1/4 length of data Hamming window)

圖11 1/8數(shù)據(jù)長度Hamming窗的STFTFig.11 STFT(with 1/8 length of data Hamming window)

圖12 組合窗函數(shù)的STFTFig.12 STFT with combination of time window function

圖13 信號的WVDFig.13 The WVD of signal

圖14 信號的小波時(shí)頻表示Fig.14 Wavelet time-frequency representation

6 結(jié)論

固定時(shí)間窗函數(shù)的STFT時(shí)間分辨率和頻率分辨率相互制約,影響時(shí)頻聚集性。通過分別采用長、短窗獲得兩組STFT時(shí)頻表示并以乘積結(jié)果作為最終時(shí)頻表示的組合窗函數(shù)方法,可使STFT時(shí)頻表示方法具有長時(shí)間窗函數(shù)STFT的頻率分辨率,同時(shí)具有短時(shí)間窗函數(shù)STFT的時(shí)間分辨率。仿真和試驗(yàn)數(shù)據(jù)處理結(jié)果表明:組合窗函數(shù)有效提高了STFT時(shí)頻表示的時(shí)頻聚集性,對于多分量疊加信號和跳頻類信號具有很好的時(shí)頻表示結(jié)果,而實(shí)際工程上如故障檢測、醫(yī)學(xué)信號處理等領(lǐng)域中,這類信號具有普遍性,因此,該法具有實(shí)際工程應(yīng)用價(jià)值。

[1]趙立昕,蔡志堅(jiān),周正.超寬帶信號的時(shí)頻分析[J].高技術(shù)通訊,2006,16(2):133-135.ZHAO Lixi,CAI Zhijian,ZHOU Zheng.The time-frequency analysis on Ultra-Wide Band signal[J].Chinese High Technology Letters,2006,16(2):133-135.

[2]樊永生,鄭鐵鋼.時(shí)頻分布的弱信號檢測技術(shù)的研究與應(yīng)用[J].振動工程學(xué)報(bào),2005,18(6):324-328.FAN Yongsheng,ZHENG Tiegang.Research on theweak signal detection technique and application based on timefrequency distribution[J].Journal of Vibration Engineering,2005,18(6):324-328.

[3]趙兆,是湘全.基于STFT和時(shí)變?yōu)V波的跳頻干擾抑制方法[J].探測與控制學(xué)報(bào),2009,31(3):46-50.ZHAO Zhao,SHI Xiangquan.FM interference suppression based on STFT and time-varying filtering[J].Journal of Detection&Control,2009,31(3):46-50.

[4]鄒紅星,戴瓊海,李衍達(dá),等.不含交叉項(xiàng)干擾且具有WVD聚集性的時(shí)頻分布之不存在性[J].中國科學(xué)E輯,2001,31(4):348-353.

[5]張曦,杜興民,茹樂.改進(jìn)的快速短時(shí)傅里葉變換算法在跳頻信號分析中的應(yīng)用[J].探測與控制學(xué)報(bào),2007,29(2):30-34.ZHANG Xi,DU Xingmin,RU Le.Application of a modified fast STFT method in frequency hopping signals analysis[J].Journal of Detection&Control,2007,29(2):30-34.

[6]Pinnegar C R,Mansinha L.Time-local Fourier analysis with a scalable,phase-modulated analyzing function:the S transform with a complex window[J].Signal Processing,2004,84(7):1 167-1 176.

[7]Pinnegar C R.A new subclass of complex-valued S-transform windows[J].Signal Processing,2006,86(8):2 051-2 055.

[8]Kwok K,JonesL.Improved instantaneous frequency estimation using an adaptive short-time fourier transform[J].IEEE trans on signal processing,2000,48(10):2 964-2 972.