韓建平,, 錢 炯, 董小軍

(1蘭州理工大學(xué)防震減災(zāi)研究所 蘭州,730050) (2同濟(jì)大學(xué)土木工程防災(zāi)國(guó)家重點(diǎn)實(shí)驗(yàn)室 上海,200092)

引 言

希爾伯特-黃變換 (Hilbert-Huang transform,簡(jiǎn)稱 HHT)在海洋、地震、生物和結(jié)構(gòu)健康監(jiān)測(cè)等領(lǐng)域的應(yīng)用中顯示出了獨(dú)特的優(yōu)勢(shì)[1-4]。 HHT包含兩個(gè)基本過(guò)程:經(jīng)驗(yàn)?zāi)B(tài)分解(empirical mode decomposition,簡(jiǎn)稱 EMD)和Hilbert變換。通過(guò) EMD可將復(fù)雜信號(hào)分解成有限個(gè)固有模態(tài)函數(shù)(intrinsic mode function,簡(jiǎn)稱 IMF)之和,然后對(duì)各個(gè) IMF進(jìn)行 Hilbert變換,研究信號(hào)的時(shí)頻能量分布,得到的Hilbert譜在聯(lián)合的頻率-時(shí)間域中描述原始信號(hào)具有較高的時(shí)頻分辨率[5]。國(guó)內(nèi)外學(xué)者對(duì)EMD邊界處理方法作了許多研究,主要有鏡像延拓法[6]、神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)法[7-8]、AR預(yù)測(cè)法[9]、多項(xiàng)式外延法[10]及波形特征匹配延拓法[11]等。

本文選擇鏡像延拓法和徑向基函數(shù)(radial basis function,簡(jiǎn)稱 RBF)神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)法作為 EMD中邊界處理的方法,對(duì)一組簡(jiǎn)單數(shù)值模擬信號(hào)和一組 12層鋼筋混凝土框架模型振動(dòng)臺(tái)試驗(yàn)實(shí)測(cè)加速度數(shù)據(jù)進(jìn)行了處理和分解,并對(duì)這兩種邊界處理方法的效果進(jìn)行了討論。

1 端點(diǎn)效應(yīng)的產(chǎn)生及影響

經(jīng)驗(yàn)?zāi)B(tài)分解中,求包絡(luò)平均是通過(guò)對(duì)原數(shù)據(jù)中的上極值點(diǎn)和下極值點(diǎn)分別進(jìn)行樣條插值擬合,然后再平均。在樣條插值時(shí),由于端點(diǎn)處極值的不確定性,每一次樣條插值都有擬合誤差,誤差不斷累積,分解出來(lái)的第l個(gè) IMF端點(diǎn)處會(huì)有較大的誤差。第 2個(gè) IM F的分解是建立在原始數(shù)據(jù)減去第 1個(gè)IMF的殘余項(xiàng)的基礎(chǔ)上進(jìn)行的,由于第 1個(gè) IM F的誤差,使殘余項(xiàng)也產(chǎn)生誤差,導(dǎo)致分解的第 2個(gè) IM F產(chǎn)生更大的誤差。依此類推,隨著分解的進(jìn)行,誤差就會(huì)由端點(diǎn)處向內(nèi)逐漸傳播,最后在嚴(yán)重的情況下會(huì)使分解的數(shù)據(jù)完全失去意義[8]。數(shù)值仿真信號(hào)的分解可以明顯地反映端點(diǎn)效應(yīng)的影響,如下

如圖1所示,該數(shù)據(jù)序列只存在3個(gè)極大值點(diǎn)和4個(gè)極小值點(diǎn),所得上、下包絡(luò)線均出現(xiàn)失真,尤其是上包絡(luò)線在信號(hào)兩端出現(xiàn)巨大的失真。同時(shí),由于數(shù)據(jù)序列很短,由端點(diǎn)處造成的包絡(luò)誤差已經(jīng)“污染”到整個(gè)數(shù)據(jù)序列。

圖2是未對(duì)原始信號(hào)x(t)進(jìn)行邊界處理時(shí),通過(guò) EMD得到的各階 IM F分量與各自對(duì)應(yīng)的真實(shí)信號(hào)的對(duì)比。分解得到的各階IMF分量與真實(shí)信號(hào)相比,均出現(xiàn)了偏差,尤其是低頻分量,偏差已經(jīng)非常大,未分解出第3階分量。由于三次樣條插值時(shí)需要用到前后各兩個(gè)臨近點(diǎn),利用三次樣條函數(shù)進(jìn)行曲線擬合出現(xiàn)這樣的問(wèn)題是很自然的,解決的惟一途徑是在數(shù)據(jù)序列的兩端各增加一個(gè)或多個(gè)極大值和極小值。

圖1 原始數(shù)值模擬信號(hào)及其上、下包絡(luò)線

圖2 無(wú)邊界處理的數(shù)值模擬信號(hào) EMD分解

2 鏡像閉合延拓

根據(jù)鏡像的對(duì)稱映射特點(diǎn),為盡量減少鏡像的副作用,將鏡子放置在信號(hào)具有對(duì)稱性極值所在位置,首先繪出信號(hào)左、右兩端的一段曲線及相應(yīng)的極值,根據(jù)曲線的分布特征,決定放置鏡子的位置。如在從信號(hào)左邊起向右的第l個(gè)極值處和從信號(hào)右邊起向左的第r個(gè)極值處分別放置兩面平面鏡,把鏡內(nèi)的信號(hào)向外進(jìn)行映射,得到序列長(zhǎng)度為兩倍于鏡內(nèi)信號(hào)長(zhǎng)度的周期性信號(hào)。經(jīng)鏡像延拓后的信號(hào)不含端點(diǎn)(具有周期性)。由于信號(hào)的閉合性,在后續(xù)的EMD過(guò)程中則不會(huì)遇到端點(diǎn)的延拓問(wèn)題。因此,鏡像閉合延拓法可用于處理 EMD的端點(diǎn)效應(yīng)。

圖3為數(shù)值仿真信號(hào)經(jīng)過(guò)鏡像延拓法處理后的EMD結(jié)果?？梢钥闯?采用鏡像延拓法對(duì)原始信號(hào)進(jìn)行處理后可以分解得出3階 IM F。其中,第 1階為高頻,IMF分量與真實(shí)信號(hào)偏差很小,基本上與原始信號(hào)重合;而后 2階低頻的 IM F分量仍存在一定的偏差,但分解效果明顯好于無(wú)邊界處理時(shí)的分解。

圖3 采用鏡像延拓處理后的數(shù)值模擬信號(hào)的 EM D

3 基于RBF神經(jīng)網(wǎng)絡(luò)的數(shù)據(jù)序列延拓

采用 RBF神經(jīng)網(wǎng)絡(luò),可以根據(jù)整個(gè)時(shí)間序列的表現(xiàn)特性以及時(shí)間序列在端點(diǎn)的具體形式對(duì)其加以延拓,在數(shù)據(jù)兩端延拓出極值點(diǎn)。

首先,把相鄰的 50個(gè)數(shù)據(jù)作為輸入向量,與之相連的1個(gè)數(shù)據(jù)作為輸出向量,形成訓(xùn)練樣本。本文一共選擇了 850個(gè)訓(xùn)練樣本,將其輸入到網(wǎng)絡(luò)中以訓(xùn)練網(wǎng)絡(luò),同時(shí)選取預(yù)測(cè)誤差最小的擴(kuò)展常數(shù)來(lái)調(diào)試網(wǎng)絡(luò)。圖4表明 851個(gè)內(nèi)推值與原始信號(hào)基本重合,然后每次將前一次預(yù)測(cè)的值重新加入網(wǎng)絡(luò)學(xué)習(xí),在此基礎(chǔ)上預(yù)測(cè)后一個(gè)值,不斷地重復(fù)學(xué)習(xí)預(yù)測(cè),最終可以在數(shù)據(jù)序列的兩端延拓出極大值點(diǎn)或極小值點(diǎn),如圖 5所示。

圖4 RBF神經(jīng)網(wǎng)絡(luò)對(duì)數(shù)值模擬信號(hào)的預(yù)測(cè)

圖5 RBF神經(jīng)網(wǎng)絡(luò)延拓后的數(shù)值模擬信號(hào)

利用基于 RBF神經(jīng)網(wǎng)絡(luò)延拓得到的數(shù)據(jù)序列與原始信號(hào)構(gòu)成一個(gè)新的信號(hào),EMD分解過(guò)程中得到的包絡(luò)線如圖 6所示。與原信號(hào)包絡(luò)圖的對(duì)比表明,經(jīng) RBF神經(jīng)網(wǎng)絡(luò)延拓后,三次樣條插值得到的包絡(luò)圖失真現(xiàn)象得到了明顯的改善。利用鏡像法處理后的分解結(jié)果如圖 7所示,與圖 3相比,較好地分解出了全部 3個(gè) IMF分量,同時(shí)各階 IMF分量與真實(shí)信號(hào)的偏差明顯減小,分解精度大為提高。

圖6 經(jīng) RBF神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)和鏡像延拓處理后信號(hào)的上、下包絡(luò)線

圖7 經(jīng) RBF神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)和鏡像延拓處理信號(hào)的 EM D分解

4 結(jié)構(gòu)實(shí)測(cè)振動(dòng)信號(hào)的處理與分解

選用同濟(jì)大學(xué)土木工程防災(zāi)國(guó)家重點(diǎn)實(shí)驗(yàn)室完成的單跨 12層鋼筋混凝土框架模型振動(dòng)臺(tái)試驗(yàn)所得的加速度數(shù)據(jù)作為復(fù)雜信號(hào),來(lái)驗(yàn)證本文提出方法的有效性。某工況下結(jié)構(gòu)頂層X(jué)方向一測(cè)點(diǎn)在 El Centro地震波激勵(lì)下的實(shí)測(cè)加速度信號(hào)如圖 8所示[12]。根據(jù) Fourier分析結(jié)果,選取合適的帶通濾波頻率對(duì)原始數(shù)據(jù)進(jìn)行濾波。濾波后的加速度信號(hào)如圖 9所示。

圖8 某測(cè)點(diǎn)加速度響應(yīng)原始信號(hào)

圖9 濾波后的加速度響應(yīng)信號(hào)

利用 RBF神經(jīng)網(wǎng)絡(luò)對(duì)濾波后的信號(hào)進(jìn)行延拓,把相鄰的 50個(gè)數(shù)據(jù)作為輸入向量,與之相連的1個(gè)數(shù)據(jù)作為輸出向量,形成訓(xùn)練樣本;選取預(yù)測(cè)誤差最小的擴(kuò)展常數(shù)來(lái)調(diào)試網(wǎng)絡(luò),選擇得到2 638個(gè)內(nèi)推值與濾波后的信號(hào)基本重合,如圖 10所示?；谕瑯拥姆椒A(yù)測(cè)出信號(hào)前后未知的100個(gè)數(shù)據(jù),圖 11為延拓后的信號(hào)。未進(jìn)行邊界處理、直接對(duì)濾波后的加速度響應(yīng)信號(hào)的 EMD結(jié)果如圖12所示。采用鏡像法,對(duì)經(jīng) RBF神經(jīng)網(wǎng)絡(luò)延拓后信號(hào)的 EMD結(jié)果如圖13所示。

圖10 RBF神經(jīng)網(wǎng)絡(luò)延拓后的加速度響應(yīng)信號(hào)

圖11 RBF神經(jīng)網(wǎng)絡(luò)對(duì)濾波后加速度響應(yīng)信號(hào)的預(yù)測(cè)

圖12和圖13的對(duì)比表明,信號(hào)未進(jìn)行邊界處理而分解出來(lái)的 IM F分量只有 4個(gè),且各個(gè) IMF分量的端點(diǎn)都有失真的現(xiàn)象;而經(jīng)過(guò) RBF對(duì)信號(hào)延拓、再采用鏡像法處理后的分解比較完全,一共分解出了7個(gè) IMF分量,端點(diǎn)失真現(xiàn)象明顯減弱。

5 結(jié) 論

本文在探討端點(diǎn)效應(yīng)出現(xiàn)原因的基礎(chǔ)上,提出采用鏡像法和 RBF神經(jīng)網(wǎng)絡(luò)延拓法處理 EMD中的端點(diǎn)效應(yīng)問(wèn)題,并利用這些方法對(duì)一組數(shù)值仿真信號(hào)和 12層鋼筋混凝土框架模型振動(dòng)臺(tái)試驗(yàn)實(shí)測(cè)加速度記錄進(jìn)行了處理及分解。結(jié)果表明,這兩種方法基于邊界兩端預(yù)測(cè)數(shù)據(jù),可以抑制邊界效應(yīng)對(duì)分析信號(hào)的影響,在一定程度上都提高了 EMD的分解效果。在處理簡(jiǎn)單周期信號(hào)中,利用鏡像延拓、RBF神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)或先用 RBF神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)、再采用鏡像延拓處理信號(hào)都是有效的。對(duì)于復(fù)雜信號(hào),僅采用RBF神經(jīng)網(wǎng)絡(luò)延拓原始信號(hào)對(duì)減弱端點(diǎn)效應(yīng)的影響沒(méi)有明顯改進(jìn);而對(duì)信號(hào)經(jīng)濾波后先利用 RBF神經(jīng)網(wǎng)絡(luò)延拓,再利用鏡像法進(jìn)行處理,可以明顯減弱端點(diǎn)效應(yīng)的影響。

圖12 無(wú)邊界處理的加速度響應(yīng)信號(hào) EM D分解

圖13 經(jīng) RBF神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)和鏡像延拓處理的濾波后加速度響應(yīng)信號(hào)的 EM D分解

[1] Huang N E,Shen Z,Long S R.The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis[J].Proceedings of the Royal Society A,1998(454):903-995.

[2] Han Jianping,Li Dawen,Li Hui.Modal parameter identification of civil engineering structures based on Hilbert-Huang transform[C]∥Engineering Structural Integrity: Research,Development and Application.Beijing:China Machine Press,2007:366-368.

[3] Yang J N,Lei Y,Lin S,et al.Hilbert-Huang based approach for structural damage detection[J].ASCE Journal of Engineering Mechanics,2004,130(1):85-95.

[4] Rilling G,Flandrin P,Goncalves P.On empirical mode decomposition and its algorithms[C]∥IEEEEurasip Workshop on NonlinearSignal and Image Processing.Grado-Trieste:N SIP,2003.

[5] 胡勁松,楊世錫,吳昭同,等.基于 EMD和 HHT的旋轉(zhuǎn)機(jī)械振動(dòng)信號(hào)時(shí)頻分析 [J].振動(dòng)、測(cè)試與診斷,2004,24(2):106-110.

[6] 黃大吉,趙進(jìn)平,蘇紀(jì)蘭,等.希爾伯特-黃變換的端點(diǎn)延拓[J].海洋學(xué)報(bào),2003,25(1):1-11.

[7] 鄧擁軍,王偉,錢成春,等.EMD方法及 Hilbert變換中邊界問(wèn)題的處理 [J].科學(xué)通報(bào),2001,46(3):257-263.

[8] 瞿偉廉,程磊.應(yīng)用徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)處理EM D方法中的邊界問(wèn)題 [J].華中科技大學(xué)學(xué)報(bào):城市科學(xué)版,2006,23(4):1-4.

[9] 張郁山,梁建文,胡聿賢.應(yīng)用自回歸模型處理 EMD方法中的邊界問(wèn)題 [J].自然科學(xué)進(jìn)展,2003,13(10):1054-1059.

[11]胡愛(ài)軍,安連鎖,唐貴基.Hilbert-Huang變換端點(diǎn)效應(yīng)處理新方法 [J].機(jī)械工程學(xué)報(bào),2008,44(4):154-158.

[12]呂西林,李培振,陳躍慶.12層鋼筋混凝土標(biāo)準(zhǔn)框架振動(dòng)臺(tái)模型試驗(yàn)的完整數(shù)據(jù) [R].上海:同濟(jì)大學(xué)土木工程防災(zāi)國(guó)家重點(diǎn)實(shí)驗(yàn)室振動(dòng)臺(tái)試驗(yàn)室,2004.