信號處理中,頻率是信號最重要的表示。傳統(tǒng)的傅里葉變換分析方法并不能分析出信號的某一頻率在甚么時刻出現(xiàn),為此產(chǎn)生了能同時在時間和頻率上表示信號密度和強度的時頻分析,如短時傅里葉變換和小波變換等,但其基本思想都是根據(jù)傅里葉分析理論,對非線性非平穩(wěn)信號的分析能力不足,受限于Heisenberg不確定原理。HHT ( Hilbert Huang Transform)是由N. E.Huang 等人在1998 年提出的一種嶄新的時頻分析方法,能夠?qū)Ψ蔷性非平穩(wěn)的信號進行分析,同時具有良好自適應(yīng)性的特點。其本質(zhì)是對信號進行平穩(wěn)化處理,將具有不同時間尺度的信號逐級分解開來。
HHT 方法在各領(lǐng)域已得到了廣泛應(yīng)用,但依然存在一些不足,例如易產(chǎn)生虛假分量和模態(tài)混疊等。針對傳統(tǒng)經(jīng)驗?zāi)J? Empirical Mode Decomposit iON,EMD)分解方法所導致的模態(tài)混疊現(xiàn)象,法國以Flandrin 為首的EMD 算法研究小組和Huang 本人的研究小組通過對EMD 分解白噪聲結(jié)果統(tǒng)計特性的大量研究,提出通過加噪聲輔助分析( NA DA ) 的EEMD ( EnsembleEmpirical Mode Decomposition) 方法,將白噪聲加入信號來補充一些缺失的尺度,在信號分解中具有良好的表現(xiàn)。
EEMD仿真系統(tǒng)的實現(xiàn)利用了Matlab 平臺,通過GUI 控件實現(xiàn)了系統(tǒng)設(shè)計,能直觀方便地進行比較分析,驗證了EEMD 在抗混疊方面較原有方法的改進。
1 經(jīng)驗?zāi)J椒纸? EMD) 和IMF
HHT 方法包含兩個主要步驟:
( 1) 對原始數(shù)據(jù)進行經(jīng)驗?zāi)J椒纸? EMD) ,把數(shù)據(jù)分解為滿足Hilbert 變換要求的n 階本征模式函數(shù)( IMF) 和殘余函數(shù)之和。
( 2) 對每一階IMF 進行Hilbert 變換,得到瞬時頻率,從而求得時頻圖。
函數(shù)必須關(guān)于時間軸局部對稱,且其過零點與極值點個數(shù)相同。此類函數(shù)被稱為固有模態(tài)函數(shù)( Int rinsicMode Function,IMF) 。
經(jīng)驗?zāi)J椒纸夥椒馨逊瞧椒(wěn)、非線性信號分解成一組穩(wěn)態(tài)和線性的序列集,即本征模式函數(shù)。根據(jù)Huang 的定義,每一階的IMF 應(yīng)滿足兩個條件:
( 1) 數(shù)據(jù)的極值點和過零點交替出現(xiàn),且數(shù)目相等或最多相差一個任何點上;
( 2) 在任何點上,有局部最大值和局部最小值定義的包絡(luò)的均值必須是零。
其篩選算法如下:
( 1) 對于輸入信號x ( t) ,確定x ( t) 所有極值點。
( 2) 用三次樣條函數(shù)對極大點和極小點分別進行擬合得到x ( t) 的上下包絡(luò)線。
( 3) 用原始數(shù)據(jù)序列減去上下包絡(luò)線的均值。
平均曲線:
細節(jié)信號:
( 4) 通常s( t ) 還不滿足IMF 的條件,需重復(fù)進行以上步驟,進行迭代處理,H uang 給出的迭代停止準則為:
SD 是篩選門限值,一般取值為0. 2~ 0. 3,若計算SD 小于這個門限值,篩選迭代將會結(jié)束。
經(jīng)過n 次迭代滿足停止準則后得到的sn ( t) 即為有效IMF,剩余信號則進入下一輪篩選過程。
經(jīng)過多次篩選后,原始數(shù)據(jù)序列被分解為一組IMF 分量和一個殘余量,得到的IMF 都是平穩(wěn)的,通過Hilbert 變換得到的結(jié)果能夠很好地分析非線性非平穩(wěn)的信號。
2 傳統(tǒng)EMD 的不足與缺陷
當信號的時間尺度存在跳躍性變化時,對信號進行EMD 分解,會出現(xiàn)一個IMF 分量包含不同時間尺度特征成分的情況,稱之為模態(tài)混疊。
模態(tài)混疊的出現(xiàn)一方面和EMD 的算法有關(guān),另一方面也受原始信號頻率特征的影響。
Huang 曾經(jīng)提出了中斷檢測的方法來解決模態(tài)混疊現(xiàn)象,即直接對結(jié)果進行觀察,如果出現(xiàn)混疊則重新分解,這種方法需要人為后驗判斷。
重慶大學的譚善文提出了多分辨率的EMD 思想,對每一個IMF 規(guī)定一個尺度范圍來解決模態(tài)混疊,但是這種方法犧牲了EMD 良好的自適應(yīng)性。
3 引入正態(tài)分布白噪聲的EEMD
為了更好地解決模態(tài)混疊問題,Huang 提出了EEMD,這是一種噪聲輔助信號處理方法。
降噪技術(shù)的目的是將噪聲從信號中去除,不過在一些情況下,可以通過加入噪聲的方法來進行輔助分析,這鐘方法就稱為噪聲輔助信號處理( NADA) ,噪聲輔助信號處理方法最常見的就是預(yù)白化。在信號中加入白噪聲來平滑脈沖干擾,被廣泛用于各種信號分析領(lǐng)域。
在EMD 方法中,得到合理IMF 的能力取決于信號極值點的分布情況,如果信號極值點分布不均勻,會出現(xiàn)模態(tài)混疊的情況。為此,Huang 將白噪聲加入待分解信號,利用白噪聲頻譜的均勻分布,當信號加在遍布整個時頻空間分布一致的白噪聲背景上時,不同時間尺度的信號會自動分布到合適的參考尺度上,并且由于零均值噪聲的特性,經(jīng)過多次平均后,噪聲將相互抵消,集成均值的結(jié)果就可作為最終結(jié)果。
EEMD 步驟如下:
( 1) 向信號加入正態(tài)分布白噪聲。
( 2) 將加入白噪聲的信號分解成各IMF 分量。
( 3) 重復(fù)步驟( 1) ,( 2) ,每次加入新的白噪聲序列。
( 4) 將每次得到的IMF 集成均值作為最終結(jié)果。
EMMD 算法流程如圖1 所示。
圖1 EEMD 算法流程圖
4 系統(tǒng)功能介紹和仿真實驗分析
為了驗證EEMD 方法的改進之處,利用Mat lab 的GU I 工具設(shè)計了簡單直觀的仿真系統(tǒng)。
此系統(tǒng)實現(xiàn)的功能是,對輸入信號進行傳統(tǒng)EMD分解和EEMD 分解,可顯示信號分解后的各個模態(tài)函數(shù)IMF 分量及其瞬時頻率,并能對Hilbert 時頻譜進行刻畫。
系統(tǒng)界面如圖2 所示。
圖2 仿真系統(tǒng)界面
參數(shù)設(shè)置功能 可自由設(shè)置加入白噪聲的方差和噪聲組數(shù)目( 范圍1~ 500) ,當方差設(shè)置為0,噪聲組數(shù)目選擇為1 時,該系統(tǒng)實現(xiàn)傳統(tǒng)EMD 分解的功能。
EEMD 分解功能 對信號進行加入上述設(shè)定白噪聲EEMD 分解,并刻畫出輸入信號的Hilbert 時頻譜。
顯示IMFs 功能 可通過彈出FIG 的形式顯示對信號分解后的各IMF 分量及瞬時頻率。
仿真實驗結(jié)果如下:
首先對多分量理想樣本信號進行分解,信號構(gòu)成如下:
其中,歸一化頻率為:
其中,歸一化頻率為:
EMD 分解方法應(yīng)將包含4 個頻率分量的信號分解為4 個包含單一頻率信息的IMF 分量。
分解結(jié)果如圖3 所示。
圖3 傳統(tǒng)EMD 對理想信號H ilber t 譜圖
可以看到,對于無干擾的理想信號,傳統(tǒng)EMD 分解方法具有非常好的效果,清晰地將4 個頻率分量在Hilbert 譜上顯示了出來。
對一組存在中斷干擾的實際信號進行分解,結(jié)果如圖4~ 圖6 所示。
圖4 實際信號時域圖
圖5 傳統(tǒng)EMD 對信號的分解
圖6 傳統(tǒng)EM D 對信號的H ilber t 譜刻畫
通過頻譜圖可以看到,低頻分量混雜在一起,難以分辨。
對EEMD 分解方法進行分析,加入了100 組標準差為0. 2 的高斯白噪聲,結(jié)果如圖7,圖8 所示。
通過Hilbert 譜的比較可以看出,分解結(jié)果有了較大改進。
圖7 EEMD 對信號的分解
圖8 EEMD 對信號的H ilber t 譜刻畫
5 結(jié) 語
EEMD 以噪聲輔助信號處理原理為基礎(chǔ),通過加入小幅度的白噪聲來均衡信號,有效地解決了模態(tài)混疊現(xiàn)象,利用高斯白噪聲零均值的特性,使真實信號得到了保留,是對傳統(tǒng)EMD 分析方法的巨大改進。
來源:維庫