周慶飛 孫常存 李 鵬

（91439部隊 旅順 116041）

基于WPT的魚雷磁引信目標(biāo)信號消噪研究＊

周慶飛 孫常存 李 鵬

（91439部隊 旅順 116041）

目標(biāo)信號消噪處理是確保引信可靠動作的重要手段。為提高魚雷磁引信動作的可靠性，基于小波包變換（WPT）在信號消噪處理方面所具有的良好特性，采用小波包變換對復(fù)雜背景干擾下的魚雷磁引信目標(biāo)信號進行了消噪處理，并利用Matlab軟件對消噪結(jié)果進行了仿真分析。仿真結(jié)果表明，基于WPT的消噪方法能夠有效地實現(xiàn)對魚雷磁引信目標(biāo)信號的噪聲抑制，具有消噪效果好，精度高，失真小的特點，有利于提高魚雷磁引信動作的可靠性。

魚雷；磁引信；目標(biāo)信號；小波包消噪

Class NumberTJ431.7；TN911.7

1 引言

現(xiàn)代海戰(zhàn)的焦點越來越趨向于近海和淺海作戰(zhàn)，淺海極為復(fù)雜的磁和電磁環(huán)境必將成為魚雷磁引信目標(biāo)信號的主要干擾［1～2］，這樣很有可能使引信目標(biāo)信號失真，引發(fā)誤判，從而造成嚴重的后果。因此，魚雷磁引信目標(biāo)信號的消噪處理成為消除背景噪聲、提高引信接收靈敏度并確保引信可靠動作的重要手段。文獻［3～4］分別研究了魚雷電磁引信、磁引信和超聲引信目標(biāo)信號的小波消噪方法，并取得了一定的降噪效果。但由于小波消噪只是將信號的低頻部分進行分解，沒有對高頻部分進行剖分，這樣一來，高頻區(qū)域中一些有用的信號就可能被濾掉［5］。小波包變換（Wave packet transfor－mation，WPT）是小波變換的推廣，它能夠為信號提供一種更加精細的分解方法，將頻帶進行多層次劃分，對離散小波變換沒有細分的高頻部分進一步分解，并能夠根據(jù)被分析信號的特征，自適應(yīng)選擇相應(yīng)的頻帶，使之與信號匹配，從而提高時頻分辨率。鑒于此，本文采用小波包變換的方法，對復(fù)雜背景干擾下的魚雷磁引信目標(biāo)信號進行消噪處理，以期更加有效地保證引信可靠動作。

2 目標(biāo)信號模型

本文以主動磁感應(yīng)引信的目標(biāo)信號為信號模型展開研究。當(dāng)魚雷以弦角攻擊時，目標(biāo)信號的包絡(luò)可近似為一高斯脈沖［6］，可表示為

式中，Hm為包絡(luò)信號的幅度，β為一常數(shù)，大小與魚雷相對目標(biāo)的速度、相遇角、目標(biāo)的等效寬度和過靶距離等參數(shù)有關(guān)，β的經(jīng)驗公式表示為［7］

式中，θk為魚雷與目標(biāo)的相遇角；VTK為魚雷與目標(biāo)的相對速度，VTK＝VT＋VKcos（θk）；Bk為目標(biāo)的等效寬度；h為過靶距離。

3 小波包分解與重構(gòu)的基本理論

對于給定的正交尺度函數(shù)及其對應(yīng)的小波函數(shù)，存在二尺度方程［8］

式中，h0（n）和h1（n）為共軛濾波器，記φ（t）→u0（t），φ（t）→u1（t），h0（n）→hn，h1（n）→gn。

則式（6）可表示為

由遞推關(guān)系定義un（t）為

稱函數(shù)族｛un｝n∈z＋為由尺度函數(shù)u0（t）＝φ（t）所確定的小波包。

信號的小波包分解公式為［7］

式中，hk和gk是小波濾波器組，hk是低通濾波器，gk是高通濾波器。

圖1 3級小波包分解樹

式（6）表明：小波包分解在對信號的低頻分量連續(xù)進行分解的同時對高頻分量也進行連續(xù)分解，這也是小波包分解與小波分解的主要區(qū)別。小波包分解樹［9］是小波包分解的一種較為直觀的表示方法，如圖1所示為尺度是3的小波包分解樹，小波包分解樹允許信號S表示為：S＝A＋AAD＋ DAD＋DD，其中，A表示信號的低頻分量，D表示信號的高頻分量，其下標(biāo)表示小波包分解的層數(shù)。

相反的回歸算法從下而上重構(gòu)上面小波包分解的信號，則小波包重構(gòu)公式為［10］

4 小波包消噪原理

一維信號包含噪聲的模型可表示為

f（k）＝s（k）＋εe（k） k＝0，1，…，n－1（8）

式中，f（k）是包含噪聲的原信號；s（k）是有用信號，在實際中常為低頻或穩(wěn)態(tài)；e（k）是噪聲信號，一般為高頻，包含在高頻細節(jié)里。

信號降噪實質(zhì)上是抑制無用部分而恢復(fù)有用成分，基于小波包的信號消噪過程可分四步進行：

1）信號的小波包分解，選擇一個小波函數(shù)并確定所需分解的層次，然后對信號進行小波包分解；

2）確定最優(yōu)小波包基，對于一個給定的熵標(biāo)準，計算最優(yōu)樹；

3）小波包分解系數(shù)的閾值量化，對于每一個小波包分解系數(shù)，選擇一個恰當(dāng)?shù)拈撝挡ο禂?shù)進行閾值量化；

4）信號的小波包重構(gòu)，根據(jù)最底層的小波包分解系數(shù)和經(jīng)過閾值量化處理的系數(shù)，進行小波包重構(gòu)。

閾值的選擇是小波包降噪方法的一個關(guān)鍵點，從某種程度上說，它直接關(guān)系到信號消噪的質(zhì)量。如果閾值太小，去噪后的信號仍然有噪聲；相反，閾值太大，重要的信號特征又被過濾掉，引起偏差。本文采用默認閾值消噪和調(diào)節(jié)閾值消噪兩種方法。

1）默認閾值消噪

這種方法是消噪處理中最典型的方法。首先獲取在消噪過程中的默認閾值，然后用默認閾值進行降噪。

2）調(diào)節(jié)閾值降噪

根據(jù)前面消噪效果，調(diào)節(jié)閾值大小進行降噪。根據(jù)需要選取不同的閾值，由于閾值的選取關(guān)系到降噪的好壞，所以閾值調(diào)節(jié)很關(guān)鍵。

5 仿真分析

在式（1）中選取典型的魚雷過靶參數(shù)，在高斯白噪聲背景、信噪比SNR＝8的條件下，磁引信目標(biāo)信號及含有噪聲的目標(biāo)信號的仿真波形如圖2所示。

圖2 魚雷磁引信的目標(biāo)信號及含噪信號

選取db8小波函數(shù)對含噪信號進行8層小波包分解，則消噪結(jié)果如圖3所示。

圖3 白噪聲背景下磁引信目標(biāo)信號的消噪結(jié)果

若背景噪聲為有色噪聲，處理方式同上，消噪結(jié)果如圖4所示。

圖4 色噪聲背景下磁引信目標(biāo)信號的消噪結(jié)果

由圖3和圖4可以看出：對于磁引信目標(biāo)信號，當(dāng)信噪比較大時，不論背景噪聲為白噪聲還是色噪聲，經(jīng)過小波包在默認閾值下消噪后的信號波形幾乎沒有失真，只是背景噪聲為色噪聲時，消噪后信號的幅度略有降低，因此不需要做調(diào)節(jié)閾值的處理。

為了檢驗在強噪聲背景下小波包的消噪效果，首先令背景噪聲為高斯白噪聲，仍然選擇db8小波函數(shù)對含噪信號進行8層小波包分解，則在信噪比SNR＝0.2的條件下對含噪信號進行小波包消噪的結(jié)果如圖5所示。

再令背景噪聲為有色噪聲，同樣在信噪比SNR＝0.2的條件下，此時信號完全淹沒在有色背景噪聲中，選擇db8小波函數(shù)進行8層小波包分解，消噪結(jié)果如圖6所示。

由以上仿真圖形可得出如下結(jié)論：當(dāng)信噪比較小，即信號完全淹沒在噪聲中時，不論背景噪聲為白噪聲還是色噪聲，經(jīng)過小波包在默認閾值下消噪后的信號都有一定的失真，并且色噪聲背景下的波形失真比白噪聲背景下的大，這說明色噪聲比白噪聲要難處理一些，這與理論是相符的。但是經(jīng)過閾值調(diào)節(jié)后，波形失真基本消失。在研究中還發(fā)現(xiàn)，用Daubechies小波族進行消噪處理比用其他小波函數(shù)效果明顯。

圖5 強白噪聲背景下磁引信目標(biāo)信號的消噪結(jié)果

圖6 強色噪聲背景下磁引信目標(biāo)信號的消噪結(jié)果

6 結(jié)語

本文采用小波包變換的方法對魚雷磁引信的目標(biāo)信號進行了消噪處理，獲得了與真實信號十分接近的波形。理論和仿真實驗皆表明，基于WPT的魚雷磁引信目標(biāo)信號消噪方法具有消噪精度高，失真小，不產(chǎn)生畸變，保留信號原有特性等優(yōu)點，且可用于強背景噪聲下的弱信號提取，是一種理想的磁引信目標(biāo)信號消噪手段。在進行消噪處理時，如何選取小波包基和進行閾值的再調(diào)節(jié)是影響信號消噪質(zhì)量好壞的關(guān)鍵。隨著科技的發(fā)展，信號消噪將會在更多的領(lǐng)域中應(yīng)用，本文只是利用小波包理論對信號消噪進行了一般性的研究，對具體的信號要結(jié)合信號的特點選擇具體的方法，例如小波包基的選擇、閾值的調(diào)節(jié)、小波包分解層次的選擇等。

［1］陳春玉，肖素娟.魚雷主動電磁引信防誤動技術(shù)［J］.魚雷技術(shù)，2008，16（1）：29－32.

［2］肖素娟.魚雷感應(yīng)引信誤動機理分析［J］.魚雷技術(shù)，2002，10（3）：32－34.

［3］陳光，任志良，李耀波，等.魚雷引信目標(biāo)信號的小波消噪研究［J］.魚雷技術(shù)，2007，15（5）：15－19.

［4］尹熙鵬，樊養(yǎng)余，段哲民，等.基于自適應(yīng)小波神經(jīng)網(wǎng)絡(luò)實現(xiàn)魚雷超聲引信信號消噪［J］.魚雷技術(shù)，2009，17（5）：19－22.

［5］程正軍，張運陶.基于MATLAB的小波包分析在信號降噪中的應(yīng)用［J］.西華師范大學(xué)學(xué)報（自然科學(xué)版），2004，25（1）：48－50.

［6］郭恩全，楊保生.魚雷引信通用目標(biāo)信號仿真技術(shù)［J］.艦船科學(xué)技術(shù)，1994，16（1）：14－17.

［7］王紹卿.魚雷近炸引信原理與設(shè)計［M］.西安：西北工業(yè)大學(xué)出版社，1992.

［8］馮宏亮.基于最優(yōu)小波包基的信號去噪方法研究［J］.西安文理學(xué)院學(xué)報（自然版），2009，12（2）：60－63.

［9］范海寧，郭英，吳劍鋒，等.基于小波包分解的聲信號特征提取方法［J］.現(xiàn)代電子技術(shù)，2005（4）：20－21.

［10］冉啟文.小波變換與分數(shù)傅立葉變換理論及應(yīng)用［M］.哈爾濱：哈爾濱工業(yè)大學(xué)出版社，2001.

Torpedo Magnetic Fuze Signal DenoisingMethod Based on WPT

ZHOU Qingfei SUN Changcun LI Peng
（No.91439Troops of PLA，Lvshun 116041）

Torpedo fuze signal denoising is an important action to ensure reliable operation of fuze.Based on the good characteristics of wavelet packet transform（WPT）in signal denoising，wavelet packet transform is used to denoise the magnetic fuze signal under a complex background interference，and a simulation of the denoising results with Matlab is performed.Simulation result shows that the WPT denoising method can effectively eliminate background noise existed in torpedo magnetic fuze target signal with higher precision and less distortion，leading to the improvement of the reliability of torpedo magnetic fuze operation.

torpedo，magnetic fuze，signal denoising，wavelet packet transform

TJ431.7；TN911.7DOI：10.3969／j.issn.1672－9730.2015.11.043

2015年5月11日，

2015年6月23日

周慶飛，男，工程師，研究方向：魚雷試驗總體。