孫飛，張鵬，連德浩

　?。ㄎ髂峡萍即髮W(xué) 信息工程學(xué)院，四川 綿陽 621010）

      摘要：在復(fù)雜環(huán)境數(shù)據(jù)采集過程中，原始數(shù)據(jù)常常混有若干外界噪聲，為提高信號特征提取的準(zhǔn)確性，基于虛擬儀器技術(shù)設(shè)計并實現(xiàn)了經(jīng)驗?zāi)B(tài)分解（Empirical Mode Decomposition，EMD）算法，根據(jù)相關(guān)分析法篩選有效基本模式分量（Intrinsic Mode Function，IMF）完成測量信號的降噪處理，并設(shè)計仿真信號對軟件降噪效果進行驗證。驗證結(jié)果表明，所設(shè)計軟件可有效去除原始信號噪聲，可靠性較好。

　　關(guān)鍵詞：數(shù)據(jù)采集; 經(jīng)驗?zāi)B(tài)分解; 虛擬儀器技術(shù); 降噪

0引言

　　在實際信號采集過程中，原始信號容易受到復(fù)雜環(huán)境噪聲的干擾，增加后期信號特征提取難度。因此，選擇有效降噪算法對測量信號的降噪預(yù)處理顯得尤為重要。

　　經(jīng)驗?zāi)B(tài)分解（Empirical Mode Decomposition，EMD）作為信號時頻分析的重要手段，在非平穩(wěn)信號處理方面相比其他方法具有較好的優(yōu)越性，本文利用LabVIEW開發(fā)工具，根據(jù)EMD降噪原理，設(shè)計針對復(fù)雜非平穩(wěn)信號的降噪處理軟件，并構(gòu)建仿真信號對其降噪效果進行驗證。

1經(jīng)驗?zāi)B(tài)分解

　　EMD是一種與小波變換相類似的信號分解方法，對于非平穩(wěn)信號，它相比其他類信號處理方法有較強的優(yōu)越性。假設(shè)原始信號為x(t)，它采用“篩選”的模式對所分析信號進行分解，具體過程如下［1］：

　　（1）取其上下包絡(luò)局部均值組成的序列m(t)，則

　　h1(t)=x(t)－m1(t)(1)

　?。?）從原始信號中“篩選”出IMF信號，得到剩余信號r1(t)，如式（2）所示。

　　r1(t)=x(t)－ct(t)（2）

　　（3）對剩余信號r1(t)重復(fù)步驟（1）和步驟（2），如式（3）所示。

　　直到rn(t)變化很小，不能“篩選”出新的模式分量為止，則原始信號被分解成多個基本模式分量（Intrinsic Mode Function,IMF)和余項之和，如式（4）所示。

　　不同IMF分量按照從高頻到低頻順序依此排列，在實際應(yīng)用中，利用相關(guān)分析算法選擇有效模式分量，根據(jù)不同噪聲特性，將相應(yīng)模式分量進行組合還原出原始信號，從而達到去噪效果。

2算法實現(xiàn)

　　LabVIEW是一種圖形化編程語言，在測控領(lǐng)域得到廣泛應(yīng)用［23］。本文通過LabVIEW根據(jù)EMD降噪原理設(shè)計并實現(xiàn)復(fù)雜信號降噪軟件，實現(xiàn)原理分為三大模塊，步驟如下：

　?。?）三次樣條插值構(gòu)造包絡(luò)線

　　主要是根據(jù)所分析信號的極值點構(gòu)造包絡(luò)曲線，借助LabVIEW中提供的三次樣條插值功能函數(shù)來實現(xiàn)，構(gòu)造信號包絡(luò)曲線的程序如圖1所示?！　?/p>

　　（2）包絡(luò)曲線求局部均值

　　從步驟（1）獲取信號的上下包絡(luò)曲線，對曲線中各點的極值取均值即可獲取局部均值，程序結(jié)構(gòu)如圖2所示。

　?。?）剩余分量判斷

　　在信號分解過程中，需要對篩分過程進行限制，按照Huang等人的仿柯西收斂準(zhǔn)則［4］，定義閾值SD作為EMD分解的終止條件，如式（4）所示。

　　式（4）中，T為信號的持續(xù)時間，按照Huang的建議，SD取［0.2，0.3］之間，本文根據(jù)試湊法，取值為0.25，整體程序結(jié)構(gòu)如圖3所示。

　　降噪軟件整體界面如圖4所示。

3降噪效果驗證

　　為驗證本文程序?qū)崿F(xiàn)算法的正確性，構(gòu)造式（5）所示的仿真信號S(t)對軟件降噪效果進行驗證。

　　S(t)=Asin(2*π*10*t)+w(t)（5）

　　其中，w(t)為幅值為0.5的高斯白噪聲，采樣率fs=1 000 Hz，采樣數(shù)N=1 000，含噪信號波形如圖5所示。

　　對w(t)進行經(jīng)驗?zāi)B(tài)分解，共產(chǎn)生5個基本模態(tài)函數(shù)，波形如圖6所示。統(tǒng)計各個IMF分量和原始信號的相關(guān)系數(shù)，統(tǒng)計如表1所示。

　　從表1中可以看出，IMF4與原始信號相關(guān)系數(shù)最大，因此，選擇模式IMF4作為有效分量重構(gòu)信號，同時，為了更好顯示降噪效果，與原始無噪聲信號作對比，波形如圖7所示。

　　從圖中可以看出，原始含噪信號經(jīng)過EMD降噪處理，噪聲得到有效控制，較好地還原出真實信號，從而說明本文程序所實現(xiàn)算法的有效性,計算兩類信號的相關(guān)系數(shù)，其結(jié)果為0.951 292，進而驗證本文所設(shè)計軟件的可靠性。

4結(jié)論

　　針對實際數(shù)據(jù)采集過程中容易遭受噪聲干擾的問題，本文根據(jù)EMD算法降噪原理，通過LabVIEW設(shè)計并實現(xiàn)信號降噪處理軟件，同時構(gòu)建仿真信號對軟件降噪效果進行驗證。驗證結(jié)果表明，本文所設(shè)計軟件對復(fù)雜噪聲處理效果較好，在實際數(shù)據(jù)采集過程中具有較好的應(yīng)用價值。

參考文獻

　?。?］ 李貴子，樓軍緯. EMD包絡(luò)譜在虛擬儀器診斷系統(tǒng)中的實現(xiàn)及應(yīng)用［J］.機械制造與自動化,2014(6):205207.

　　［2］ 王水魚,馮曉靖. 一種基于虛擬儀器技術(shù)的任意波形發(fā)生器［J］. 微型機與應(yīng)用,2013,32(18):1719.

　?。?］ 張力,王盈. LabVIEW與三菱PLC串口通信的實現(xiàn)［J］.微型機與應(yīng)用,2013,32(12):5153.

　?。?］ 曹沖鋒. 基于EMD的機械振動分析與振動方法研究［D］. 杭州：浙江大學(xué)，2009.