周祥鑫，王小敏，楊 揚(yáng)，郭 進(jìn)，王 平

（1．西南交通大學(xué)信息科學(xué)與技術(shù)學(xué)院交通信息工程及控制省重點(diǎn)實(shí)驗(yàn)室，成都 610031；2．西南交通大學(xué)土木工程學(xué)院高速鐵路線路工程教育部重點(diǎn)實(shí)驗(yàn)室，成都 610031）

基于小波閾值的高速道岔振動信號降噪

周祥鑫1，王小敏1，楊 揚(yáng)1，郭 進(jìn)1，王 平2

（1．西南交通大學(xué)信息科學(xué)與技術(shù)學(xué)院交通信息工程及控制省重點(diǎn)實(shí)驗(yàn)室，成都 610031；2．西南交通大學(xué)土木工程學(xué)院高速鐵路線路工程教育部重點(diǎn)實(shí)驗(yàn)室，成都 610031）

在高速道岔傷損監(jiān)測中，道岔振動信號是道岔傷損監(jiān)測的重要信息來源，鑒于該信號在采集和傳輸過程中噪聲干擾嚴(yán)重，影響傷損識別的準(zhǔn)確性，研究了一種基于小波閾值的高速道岔振動信號降噪方法。詳細(xì)討論了道岔振動信號降噪過程中小波基、分解層數(shù)、閾值準(zhǔn)則、閾值函數(shù)等參數(shù)的選擇，并利用頻響函數(shù)、主元分析和平均馬氏距離分析降噪處理對傷損識別的影響。實(shí)驗(yàn)結(jié)果表明，該方法能有效降低噪聲對傷損識別的干擾，為進(jìn)一步對道岔進(jìn)行傷損分析創(chuàng)造了良好的條件。

高速道岔；振動信號；小波閾值；降噪；傷損識別

高速道岔是高速鐵路的基礎(chǔ)設(shè)備，屬于鐵路線路中的活動部件和脆弱環(huán)節(jié)。在輪軌相互作用、外部載荷、環(huán)境和災(zāi)害等因素作用下，道岔易產(chǎn)生水平彎曲、拱腰、剝落掉塊、裂紋甚至斷裂等損傷，結(jié)構(gòu)性能下降，嚴(yán)重影響道岔的使用安全。因此，實(shí)時(shí)掌握道岔的工作狀態(tài)，是保障高速鐵路行車安全和行車效率的關(guān)鍵。

基于結(jié)構(gòu)頻響函數(shù)（FRF）指標(biāo)的傷損監(jiān)測是研究高速道岔結(jié)構(gòu)狀態(tài)的主要手段之一［1］并在故障診斷領(lǐng)域得到廣泛應(yīng)用［2－4］，主要因?yàn)椋孩兕l響函數(shù)反映結(jié)構(gòu)固有特性，與外界載荷無關(guān)，比時(shí)域信號穩(wěn)定；②結(jié)構(gòu)模態(tài)信息豐富，可完備表征結(jié)構(gòu)的動力特性，既可反映結(jié)構(gòu)整體特性，又可反映結(jié)構(gòu)局部信息。環(huán)境隨機(jī)激振下頻響函數(shù)的估計(jì)［5］：

式中，Gyf（ω）表示輸出信號與輸入信號的互功率譜，Gff（ω）表示輸入信號的自功率譜。由式（1）可知，道岔結(jié)構(gòu)頻響函數(shù)估計(jì)的準(zhǔn)確性主要取決于對自譜和互譜的估計(jì)精度，很大程度上取決于實(shí)測信號的信噪比。而實(shí)測信號由不同動量的列車過岔激勵(lì)產(chǎn)生，且受道岔活動部件相互作用影響，其激勵(lì)和外部環(huán)境干擾復(fù)雜，導(dǎo)致道岔振動信號信噪比低，嚴(yán)重影響對頻響函數(shù)的準(zhǔn)確估計(jì)，最終影響傷損識別的準(zhǔn)確率。因此，在傷損識別之前，對道岔振動信號進(jìn)行降噪處理具有重要意義。

目前，針對機(jī)械設(shè)備振動信號，根據(jù)噪聲產(chǎn)生機(jī)理和類型，常用的信號降噪方法有EMD降噪［6－7］、盲源分離［8］、小波閾值降噪［9－12］等。EMD方法對信號分解時(shí)無需進(jìn)行參數(shù)設(shè)置，是一種自適應(yīng)信號分解方法，但存在端點(diǎn)效應(yīng)、模態(tài)混疊和實(shí)時(shí)性較差等不足，雖然Wu等［13］提出的集總平均經(jīng)驗(yàn)?zāi)B(tài)分解（EEMD）算法有效解決了模態(tài)混疊問題，卻導(dǎo)致計(jì)算量增加了1～2個(gè)數(shù)量級，不利于實(shí)時(shí)處理；盲源分離技術(shù)從混合信號中恢復(fù)源信號具有明顯的效果，所需先驗(yàn)知識較少，但考慮道岔監(jiān)測系統(tǒng)的成本因素，傳感器安裝數(shù)量有嚴(yán)格限制，而在欠定條件下準(zhǔn)確提取源信號具有一定的困難。小波分析相對于Fourier分析，具有更好的時(shí)域分辨率和頻域分辨率，被廣泛應(yīng)用于各類結(jié)構(gòu)振動信號處理領(lǐng)域。文獻(xiàn)［9］針對鐵路機(jī)車振動信號，采用小波變換對其進(jìn)行降噪處理，并取得明顯的效果；文獻(xiàn)［10－12］分別提出針對齒輪、發(fā)動機(jī)、軸承等旋轉(zhuǎn)機(jī)械振動信號的小波閾值降噪方法。但對不同振動對象、不同激勵(lì)和不同環(huán)境干擾下的實(shí)測信號進(jìn)行小波降噪時(shí)，在小波參數(shù)的選取上缺乏理論指導(dǎo)，表現(xiàn)在小波基、分解層數(shù)、閾值等參數(shù)選取上存在很大的不確定性，甚至不同小波參數(shù)可能產(chǎn)生相反的去噪效果。

針對高速道岔傷損監(jiān)測對信號降噪的有效性和實(shí)時(shí)性要求，本文提出一種基于小波參數(shù)優(yōu)選的振動信號閾值降噪算法。通過構(gòu)建與實(shí)測振動信號匹配的仿真信號，完成對小波基、分解層數(shù)、閾值等參數(shù)的優(yōu)選，然后基于小波閾值去噪理論，確定優(yōu)選降噪方案。在此基礎(chǔ)上，結(jié)合實(shí)測列車激勵(lì)道岔振動信號和試驗(yàn)平臺道岔振動信號，分析了優(yōu)選降噪方案對道岔傷損識別的影響。實(shí)驗(yàn)表明，本文方法對噪聲的抑制有顯著效果，并且滿足在線監(jiān)測的實(shí)時(shí)性要求；經(jīng)降噪后的道岔振動信號有利于提高道岔傷損識別的準(zhǔn)確性。

1 道岔振動信號的小波閾值降噪方案

1.1 小波閾值降噪的基本原理

由Donoho等［14］提出小波閾值降噪理論，其主要理論依據(jù)是：信號經(jīng)小波分解后，有用成分的小波系數(shù)絕對值較大，而噪聲成分的小波系數(shù)絕對值較小，通過設(shè)定閾值，將小于閾值的小波系數(shù)置零來達(dá)到濾除噪聲的目的。

利用小于閾值降噪理論對一維信號降噪的具體步驟為：

（1）信號的小波分解：確定一個(gè)小波基和適當(dāng)?shù)姆纸鈱訑?shù)，將信號分解得到相應(yīng)的小波系數(shù)；

（2）小波系數(shù)閾值量化：對分解各層設(shè)定一個(gè)合適的閾值，將絕對值小于該閾值的小波系數(shù)置零，將絕對值大于該閾值的小波系數(shù)保留或作收縮處理；

（3）信號重建：將量化后的小波系數(shù)通過小波逆變換重建信號。

1.2 構(gòu)建道岔振動仿真信號

在工程應(yīng)用中，實(shí)測信號所含噪聲未知，不利于對降噪?yún)?shù)的選擇和降噪效果做出客觀評價(jià)。因此，本文首先通過構(gòu)建與實(shí)測振動信號匹配的仿真信號來確定優(yōu)選降噪方案。圖1（a）為實(shí)測列車激勵(lì)的道岔振動信號，虛線域①所示信號為車輪剛好經(jīng)過測點(diǎn)時(shí)的振動信號，將此類信號稱為有效振動信號，將虛線域②所示信號稱為間隙信號。

圖1 實(shí)測列車激勵(lì)的道岔振動信號Fig．1 Turnout vibration signal impelled by train

部分區(qū)域①信號放大如圖1（b）所示，圖示信號具有一定正弦特性，近似圖2所示的拍波信號［15］，其產(chǎn)生函數(shù)為：

式中，A為信號的最大幅值，w1為共振頻率，w2為包絡(luò)的正弦波頻率。w1與w2滿足關(guān)系w2＝w1／2n，n為一個(gè)拍波的共振次數(shù)，將w2＝w1／2n代入（2）式得：

對式（3）乘以一個(gè)包絡(luò)調(diào)節(jié)因子cos（kw1／2n）得到仿真信號：

式中，k為調(diào)節(jié)波形包絡(luò)形狀的系數(shù)，取k＝1．2。從圖1（b）知實(shí)測道岔振動信號最大幅值A(chǔ)約為90，n約為24，經(jīng)計(jì)算w1＝2 048 π，產(chǎn)生的道岔振動仿真信號如圖3所示。由圖可知，仿真信號與實(shí)測振動信號具有較好地相似性。

圖2 拍波信號Fig．2 Beat wave signal

圖3 仿真信號Fig．3 Simulation signal

1.3 選取小波閾值降噪?yún)?shù)

小波基類型的選取是影響降噪效果的重要因素，同一類型的小波基，其消失矩的選取也至關(guān)重要。通常在選取小波基時(shí)，小波形狀與信號形狀的相似度越高降噪效果越好［16］。經(jīng)實(shí)驗(yàn)得到Daubechies小波、Symmlet小波和Coiflet小波通過調(diào)整消失矩能夠與實(shí)測信號達(dá)到較高的相似度。

閾值的選取是影響降噪效果的另一個(gè)重要因素，如果閾值選擇過大，信號中的部分有用信號將會被視為噪聲濾除，導(dǎo)致重構(gòu)信號存在嚴(yán)重失真；相反，如果閾值選擇過小，信號的降噪效果不理想。其中閾值的確定與閾值準(zhǔn)則息息相關(guān)，通常有四種小波閾值準(zhǔn)則［17－19］：固定閾值（sqtwolog）、自適應(yīng)閾值（rigrsure）、啟發(fā)式閾值（heursure）、極大極小閾值（minimaxi）。最后，對小波系數(shù)的量化方法有軟閾值法和硬閾值法［9］。

表1為小波閾值降噪各項(xiàng)參數(shù)，本文對22種小波，4種閾值準(zhǔn)則和2種閾值函數(shù)進(jìn)行對比分析，對信號分別進(jìn)行2～10層分解，因此，上述小波參數(shù)的組合方式共有22×4×2×9＝1 584種。

表1 小波閾值降噪各項(xiàng)參數(shù)Tab.1 Parameters of wavelet threshold denoising.

1.4 確定降噪的優(yōu)選方案

工程應(yīng)用中，由于道岔安裝環(huán)境、傳感器安裝位置、激勵(lì)條件等因素的差異，導(dǎo)致實(shí)測道岔振動信號所含噪聲強(qiáng)度不同，故降噪優(yōu)選方案應(yīng)具備較好的抗噪干擾能力。因此，本文對仿真信號分別添加－10 dB、0 dB、10 dB的高斯白噪聲，分別進(jìn)行降噪，并對降噪結(jié)果進(jìn)行綜合分析來確定優(yōu)選方案。具體實(shí)驗(yàn)步驟為：

（1）確定閾值估計(jì)準(zhǔn)則：對信號分別采用上述三類小波基進(jìn)行10層分解，對分解的小波系數(shù)分別采用上述四類閾值準(zhǔn)則計(jì)算閾值，并用軟閾值函數(shù)對小波系數(shù)進(jìn)行量化；

（2）確定閾值函數(shù)：對信號分別采用上述三類小波基進(jìn)行10層分解，對分解的小波系數(shù)采用步驟（1）確定的閾值準(zhǔn)則計(jì)算閾值，分別采用軟閾值函數(shù)和硬閾值函數(shù)對小波系數(shù)進(jìn)行量化；

（3）確定小波基和分解層數(shù)：根據(jù)步驟（1）和步驟（2），選取降噪效果較好且相當(dāng)?shù)膸追N小波基，并用這幾種小波基分別對信號進(jìn)行2～10層分解，對分解的小波系數(shù)采用步驟（1）、（2）確定的閾值準(zhǔn)則和閾值量化方法進(jìn)行閾值估計(jì)和系數(shù)量化。

對上述各步驟量化后的小波系數(shù)進(jìn)行重構(gòu)得降噪信號，并計(jì)算降噪信號的信噪比。由于加入的白噪聲具有隨機(jī)性，上述三個(gè)步驟的信噪比計(jì)算均采用100次降噪的信噪比統(tǒng)計(jì)平均值，分別對應(yīng)圖4～6。由圖4可知，在不同噪聲強(qiáng)度下，啟發(fā)式閾值準(zhǔn)則降噪效果均優(yōu)于另外三種閾值準(zhǔn)則。由圖5可知，在不同噪聲強(qiáng)度下，經(jīng)軟閾值函數(shù)量化后重構(gòu)信號的信噪比均高于硬閾值函數(shù)的處理，且當(dāng)攜噪信號信噪比較低時(shí)，軟閾值函數(shù)相對于硬閾值的優(yōu)勢表現(xiàn)得更為明顯。

圖4 不同閾值準(zhǔn)則對應(yīng)的信噪比（橫坐標(biāo)：＃1—db2，＃11—sym4，＃18—coif1）Fig．4 SNR of different threshold criteria

圖5 軟硬閾值對應(yīng)的信噪比（橫坐標(biāo)：＃1—db2，＃11—sym4，＃18—coif1）Fig．5 SNR of soft threshold and hard threshold

圖6 不同分解層數(shù)對應(yīng)的信噪比Fig．6 SNR of different decomposition scale

對比分析圖4和圖5可知，上述三類小波隨著消失矩的增加，信噪比均呈上升趨勢，但經(jīng)db12、db14、db16、db18、db20處理后的信號信噪比相對更優(yōu)，且信噪比非常接近。于是，圖6在不同噪聲強(qiáng)度下對這五種較優(yōu)的小波作了進(jìn)一步對比，即采用這五種小波基對信號均進(jìn)行2～10層分解，并計(jì)算軟閾值量化后重構(gòu)信號的信噪比。由圖可知，當(dāng)信號信噪比較低時(shí)db14小波的降噪效果相對較優(yōu)，當(dāng)信號信噪比較高時(shí)，db20小波相對略優(yōu)，但在工程應(yīng)用中，實(shí)測道岔振動信號的信噪比通常較低，因此綜合考慮宜選用db14小波。

據(jù)圖6還可看出，信噪比隨分解層數(shù)的增加而增大，但當(dāng)分解層數(shù)達(dá)到7層時(shí)繼續(xù)增加分解層數(shù)對信噪比的改善較小，因此宜對信號進(jìn)行7層分解。

綜上所述，針對道岔振動信號的降噪，本文確定優(yōu)選降噪方案：采用db14小波，對信號進(jìn)行7層分解，采用啟發(fā)式閾值準(zhǔn)則計(jì)算閾值，并用軟閾值函數(shù)對各層小波系數(shù)進(jìn)行量化。圖7為優(yōu)選方案對含0 dB白噪聲的仿真信號的降噪效果，據(jù)圖可知，在低信噪比的情況下，優(yōu)選降噪方案仍能較準(zhǔn)確重建原始干凈信號，降噪后的信號失真較小。

圖7 仿真信號的降噪Fig．7 De-noising of simulation signal

2 實(shí)測信號降噪效果分析

分析兩類實(shí)測道岔振動信號的降噪效果：

（1）實(shí)測列車激勵(lì)道岔振動信號。該信號采集于成灌線犀浦站18號高速道岔，傳感器在道岔上的安裝位置如圖8所示，每條曲尖軌安裝14個(gè)傳感器，傳感器為MEAS傳感器公司生產(chǎn)的8021－01加速度傳感器。

圖8 傳感器的分布示意圖Fig．8 The distribution diagram of the transducer

（2）高速道岔傷損監(jiān)測試驗(yàn)平臺的道岔振動信號。試驗(yàn)平臺如圖9（a）所示，道岔型號為18號道岔，激勵(lì)源為一臺安裝在鋼架小車上的離心振動機(jī)，其中一條曲尖軌無傷損，另一條曲尖軌尖端存在剝落掉塊和裂紋傷損，如圖9（b）所示。

圖9 高速道岔傷損監(jiān)測試驗(yàn)平臺與傷損位置Fig．9 High-speed rail damage monitoring test platform and the position of damage

圖10為對上述兩類信號降噪的效果，由圖10（a）可以看出，對于實(shí)測列車激勵(lì)道岔振動信號的有效振動信號，降噪前信號的時(shí)域波形含有豐富的噪聲毛刺，經(jīng)降噪后能較好地消除噪聲毛刺，信號的時(shí)域波形變得較為光滑，且降噪處理未對有用信號的能量造成衰減，同時(shí)從頻譜也可以看出，對消除2．5 kHz以上的高頻白噪聲有明顯的效果；對于間隙信號的降噪，如圖10 （b），有用信號被強(qiáng)噪聲淹沒，通過降噪能有效提取出有用信號，高頻噪聲得到抑制；對圖10（c）所示的試驗(yàn)平臺有效振動信號，本文降噪方法對抑制其中的噪聲同樣具有明顯效果。

3 優(yōu)選降噪方案對傷損識別的影響

3.1 道岔的傷損識別策略

道岔振動信號的降噪，最終目的是為傷損識別提供高質(zhì)量的數(shù)據(jù)，降噪效果的好壞直接影響傷損識別的結(jié)果。因此，本文將經(jīng)優(yōu)選方案降噪后的道岔振動信號和未經(jīng)降噪的道岔振動信號分別進(jìn)行傷損識別，并進(jìn)行對比分析。根據(jù)文獻(xiàn)［2－4］提出的基于頻響函數(shù)傷損識別方法，高速道岔的傷損識別可以總結(jié)為如下步驟：

（1）構(gòu)造頻響矩陣：從采集的道岔振動信號中確定其中一個(gè)測點(diǎn)的信號作為虛擬激勵(lì)，分別計(jì)算其他各測點(diǎn)處的頻響函數(shù)，構(gòu)成頻響矩陣；

（2）數(shù)據(jù)降維：對頻響矩陣進(jìn)行主元分析，提取前m階主元形成主元矩陣；

（3）計(jì)算傷損距離：計(jì)算道岔在無傷損狀態(tài)下的參考主元矩陣與待傷損識別主元矩陣間的平均馬氏距離：

式中，Mi為對應(yīng)測點(diǎn)間的馬氏距離，i＝1，2，…，N－1， N為測點(diǎn)個(gè)數(shù)，Ya和Yd分別為參考主元矩陣和待識別的主元矩陣，C－1為矩陣Ya的協(xié)差陣，為平均馬氏距離。

3.2 結(jié)果分析

降噪對頻響函數(shù)的影響：如圖11為成灌線犀浦站第13號測點(diǎn)頻響函數(shù)的幅頻曲線，由圖可知，由未經(jīng)降噪的振動信號計(jì)算所得頻響函數(shù)，在高頻段易形成虛假譜峰，主因?qū)︻l響函數(shù)進(jìn)行估計(jì)時(shí)，自譜和互譜的計(jì)算會因高頻白噪聲的干擾而受到影響。因此，通過降噪處理消除高頻虛假譜峰的干擾，有利于獲得更加準(zhǔn)確的頻響函數(shù)。

圖10 本文優(yōu)選方案的降噪效果Fig．10 De-noising effects of optimal scheme in this paper

表2 降噪前后前10階特征值和累計(jì)貢獻(xiàn)率Tab.2 The fisrt 10 eigenvalue and cumulative contribution rate

圖11 第13號測點(diǎn)頻響函數(shù)Fig．11 The FRF of 13thmeasure point

圖12 實(shí)測列車激勵(lì)道岔振動信號的平均馬氏距離Fig．12 Average Mahalanobis distance of turnout vibration signal impelled by train

圖13 試驗(yàn)平臺道岔振動信號的平均馬氏距離Fig．13 Average Mahalanobis distance of turnout vibration signal from test platform

降噪對數(shù)據(jù)降維的影響：利用主元分析對頻響矩陣降維時(shí)，用前m階主元矩陣取代高維頻響矩陣，通常采用主元分析所得特征值的累計(jì)貢獻(xiàn)率作為控制限。設(shè)控制限為95%，由表2可以看出，降噪前要達(dá)到95%的控制限，需前10階主元，而降噪后前8階主元即可達(dá)到要求，表明降噪處理使主元分析降維效果更好，降維后的主元矩陣所攜帶頻響函數(shù)的主要信息更集中，即提高了數(shù)據(jù)壓縮比。

降噪對傷損判別的影響：對降維后的主元矩陣采用（5）式和（6）式計(jì)算傷損前后的平均馬氏距離。因?yàn)楫?dāng)?shù)啦硎艿揭欢ǔ潭鹊膿p傷時(shí)，道岔的頻響函數(shù)會發(fā)生一定變化，降維后的主元矩陣也會發(fā)生相應(yīng)的變化，這種變化最終通過傷損前后主元矩陣的平均馬氏距離量化，平均馬氏距離越大，表明道岔整體傷損程度越大，反之則傷損程度越小。

文中選取15個(gè)批次的成灌線犀浦站無損傷道岔振動信號，任選其中一個(gè)批次的數(shù)據(jù)作為參考數(shù)據(jù)，對其他14批次數(shù)據(jù)進(jìn)行傷損識別。如圖12為成灌線犀浦站高速道岔傷損識別的平均馬氏距離，由圖可以看出，降噪之前第3和第10批次數(shù)據(jù)的平均馬氏距離波動較大，如果傷損判決閾值設(shè)置不當(dāng)，容易造成誤判；通過降噪之后再進(jìn)行傷損識別，使無傷損狀態(tài)時(shí)平均馬氏距離的值更集中，在一定程度上抑制噪聲對傷損識別的干擾。

對高速道岔監(jiān)測試驗(yàn)平臺的振動信號，采集有傷損和無傷損兩種工況下的數(shù)據(jù)各29個(gè)批次，將無傷損工況的一個(gè)批次數(shù)據(jù)作為參考數(shù)據(jù)。圖13為降噪前后兩種工況下的平均馬氏距離，由圖可知，降噪前兩種工況的平均馬氏距離的平均差值較小，某些批次數(shù)據(jù)的平均馬氏距離存在重疊，不利于判決閾值的設(shè)置；經(jīng)過降噪之后，兩種工況的平均馬氏距離間隔擴(kuò)大，從而提高了兩種工況的辨識度，可在一定程度上減小誤判的概率。

4 結(jié) 論

本文針對道岔振動信號研究了一種基于小波閾值的降噪方法，對三類小波基、閾值準(zhǔn)則、閾值函數(shù)、分解層數(shù)進(jìn)行比較分析，確定一種參數(shù)優(yōu)選降噪方案，并結(jié)合基于頻響函數(shù)的傷損識別分析優(yōu)選降噪方案對傷損識別的影響，結(jié)論如下：

（1）對道岔振動信號的小波閾值降噪，采用db14小波對信號進(jìn)行7層分解、啟發(fā)式閾值估計(jì)準(zhǔn)則、軟閾值量化函數(shù)的降噪方案，在濾除噪聲的同時(shí)能較好地保留有用信號。與其他參數(shù)組合的降噪方案相比，優(yōu)選方案具有更好的抗噪適應(yīng)性。

（2）在傷損識別之前采用優(yōu)選降噪方案對實(shí)測信號進(jìn)行降噪，能更準(zhǔn)確地估計(jì)頻響函數(shù)、降低主元分析階次與傷損識別計(jì)算量。道岔傷損識別結(jié)果表明：優(yōu)選降噪方案能有效提高有損和無損工況間的辨識度，減小誤判的概率，為進(jìn)一步傷損位置定位、工況分類和定量分析提供了數(shù)據(jù)保障。

［1］王平．高速鐵路道岔監(jiān)測技術(shù)研究報(bào)告［R］．西南交通大學(xué)，2013．1．

［2］楊彥芳，宋玉普，紀(jì)衛(wèi)紅．基于實(shí)測頻響函數(shù)主成分的在役網(wǎng)架損傷識別方法［J］．振動與沖擊，2007，29（9）：55 －59．

YANGYan-fang，SONGYu-pu，JIWei-hong．An identification method of existing truss structural damage based on principal component analysis measured frequency response functions［J］．Journal of Vibration and Shock，2007，29 （9）：55－59．

［3］楊杰，李愛群，繆長青．大跨斜拉橋動力特性的主元特征提取［J］．東南大學(xué)學(xué)報(bào)：自然科學(xué)版，2006，36（4）：613 －616．

YANGJie，LIAi-qun，MIAOChang-qing．Principal component feature extraction of dynamic characters in long spancable-stayedbridge［J］．JournalofSoutheast University：Natural science，2006，36（4）：613－616．

［4］Lee U，ShinJ．Afrequencyresponsefunction-based structural damage identification method［J］．Computers and Structures，2002，80：117－132．

［5］李國強(qiáng)，李杰．工程結(jié)構(gòu)動力檢測理論與應(yīng)用［M］．北京：科學(xué)出版社，2002，11－15．

［6］Meignen S，Oberlin T．A new algorithm for multicompon-ent signalsanalysisbasedonSsynchroSqueezing：withan application to signal sampling and denoising［J］．IEEE Transactions on Signal Processing，2012，60（11）：5787 －5798．

［7］Omitaomu O A，Protopopescu V A，Ganguly A R．Empirical modedecompositiontechniquewithconditionalmutual information for denoising operational sensor data［J］．IEEE Sensors Journal，2011，11（10）：2565－2575．

［8］Gelle G，Colas M，Serviere C．Blind source separation：A new pre-processing tool for rotating machines monitoring？［J］．IEEE Transactions on Instrumentation and Measurement，2003，52（3）：790－795．

［9］廖凱，張家棟．小波變換在機(jī)車振動信號分析中的應(yīng)用研究［J］．鐵道學(xué)報(bào)，2003，25（4）：115－118．

LIAO Kai，ZHANG Jia-dong．Application of wavelet transform in the vibrating signal analysis of locomotive bogie［J］．Journal of the China Railway Society，2003，25（4）：115 －118．

［10］Mohanty A R，Kar C．Fault detection in a multistage gearbox by demodulation of motor current waveform［J］．IEEE Transactions on Industrial Electronics，2006，53（4）：1258 －1296．

［11］臧玉萍，張德江．小波分層閾值降噪法及其在發(fā)動機(jī)振動信號分析中的應(yīng)用［J］．振動與沖擊，2009，28（8）：57 －60．

ZANG Yu-ping，ZHANG De-jiang．Per-level threshold denoising method using wavelet and its application in engine vi-bration analysis［J］．Journal of Vibration and Shock，2009，28（8）：57－60．

［12］Li Zhen，He Zheng-jia．Customized wavelet denoising using intra-and inter-scale dependency for bearing fault detection ［J］．Journal of Sound and Vibration，2008，313（1－2）：342－359．

［13］Wu Z，Huang N E．Ensemble empirical mode decomposition：a noise assisted data analysis method［J］．Advances in Adaptive Data Analysis，2009，1（1）：1－41．

［14］Donoho D L．De-noising by soft-thresholding［J］．IEEE Transactions on Information Theory，1995，41（3）：613 －627．

［15］王濟(jì)，胡曉．MATLAB在振動信號處理中的應(yīng)用［M］．北京：中國水利水電出版社，2006，147－149．

［16］Ma X，Zhou C，Kemp I J．Automated wavelet selection and thresholding for PD detection［J］．IEEE Electrical Insulation Magazine，2002，18（2）：37－45．

［17］Donoho D L，Johnstone I M，Adapting to unknown smoothness via wavelet shrinkage［J］．Journal of the American Statistical Association，1995，90（432）：1200－1224．

［18］Donoho D L，Johnstone I M，Minimax estimation via wavelet shrinkage［J］．Annals of Statistic，1998，26（3）：879－921．

［19］Donoho D L，Johnstone I M．Ideal spatial adaptation by avelet shrinkage［J］．Biometrika，1994，81（3）：425－455．

De-noising of high-speed turnout vibration signals based on wavelet threshold

ZHOU Xiang-xin1，WANG Xiao-min1，YANG Yang1，GUO Jin1，WANG Ping2
（1．Key Lab of Traffic Information Engineering and Control，Southwest Jiaotong University，Chengdu 610031，China；2．School of Civil Engineering，Southwest Jiaotong University，Chengdu 610031，China）

Turnout vibration signals are an important information in high-speed turnout damage monitoring．As the signals are interfered by strong noise during the process of field acquisition and transmission，the accuracy of turnout damage identification based on its noisy vibration signals is declined seriously．To solve this problem，a denoising method is generally employed before turnout damage identification．The complex and noisy vibration samples from a site，however，raise the hurdle of denoising．Here，an effective denoising method based on wavelet threshold for turnout vibration signals was proposed．The selection of parameters，such as，wavelet basis，decomposition scale，threshold criteria and threshold function was empirically discussed for wavelet threshold denoising．Then，turnout damage identification analysis was conducted with the principal component analysis（PCA）of frequency response functions（FRF）and average Mahalanobis distance（MD）．The experimental results showed that the proposed method can be used to reduce the noise interference effectively for turnout damage identification，and create better conditions for further damage analysis．

high-speed turnout；vibration signal；wavelet threshold；de-noising；damage identification

TP274

A

10．13465／j．cnki．jvs．2014．23．036

鐵道部科技研究開發(fā)計(jì)劃重大項(xiàng)目（2011G009，2012X004－A，2013X012－A－1，2013X012－A－2）；四川省杰出青年培育基金（2011JQ0027）；中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金（SWJTU11CX041，SWJTU12CX099）

2013－08－01 修改稿收到日期：2013－12－19

周祥鑫男，碩士生，1988年生

王小敏男，博士，教授，1974年生