丁元明，柳力嘉，劉蘇睿，楊 陽

(1. 大連大學 通信與網(wǎng)絡(luò)重點實驗室，遼寧 大連 116622；2. 大連大學 信息工程學院，遼寧 大連 116622)

0 引 言

隨著海洋科學技術(shù)的發(fā)展，對艦船輻射噪聲特征提取與識別技術(shù)展開深入研究，因水下環(huán)境、地勢等多種因素，特征提取與識別分類工作尤為困難。艦船輻射噪聲產(chǎn)生的機理較為復雜，并且受到海洋環(huán)境的影響，具有非線性、非平穩(wěn)和非高斯特性，傳統(tǒng)的處理水聲信號方法比較有局限性，已經(jīng)不能有效的進行處理。

近幾年來，經(jīng)驗?zāi)B(tài)分解(empirical mode decomposition, EMD)的出現(xiàn)，艦船輻射噪聲“三非”的問題得到了解決。EMD[1]可以將非平穩(wěn)信號分解成多個相對平穩(wěn)具有不同尺度的振動信號特性的固有模態(tài)函數(shù)(intrinsic mode function, IMF)，但EMD 具有因極值點分布不均勻?qū)е履B(tài)混疊等不足問題[2]。WU 等[3]提出了集合經(jīng)驗?zāi)B(tài)分解(ensemble empirical mode decomposition, EEMD)方法解決了EMD 模態(tài)混疊的問題，該方法利用白噪聲的特性，多次加入不同幅值白噪聲改變信號極值點特性進行分解。李余興等[4-5]提出利用EEMD 對艦船輻射噪聲分解，采用中心頻率或能量差的方式進行特征提取。劉千里[6]提出EEMD 的自適應(yīng)線譜及連續(xù)譜提取方法，將EEMD 分解的IMF 進行線譜的提取，使用余量和剩余的IMF 進行連續(xù)譜的準確估計。但EEMD 屬于遞歸分解，受到端點效應(yīng)影響，也會存在一定的分解誤差。Dragomiretskiy 等[7]提出了變分模態(tài)分解方法，VMD 通過迭代搜尋變分模型最優(yōu)解來確定所知的IMF 及其對應(yīng)的中心頻率和帶寬。文獻[8] 提出的VMD 和中心頻率相結(jié)合的特征提取方法，獲取信號的特性，有效進行識別分類。在故障診斷中，文獻[9]采用VMD 與譜峭度法結(jié)合對齒輪箱故障信號進行特征提取，該方法成功地檢測出齒輪故障。VMD 與EMD，EEMD 相比，具有良好的魯棒性和堅實的數(shù)學理論基礎(chǔ)。

一般來說，艦船輻射噪聲可以看作是具有非線性的時間序列，可以從時間序列復雜度的角度進行分析。有許多表示時間序列復雜度特征的算法：如樣本熵[10]、近似熵[11]、模糊熵[12]等。其中由Bandt 等提出的排列熵(permutation entropy, PE)[13]是檢測時間序列復雜度的新方法，PE 算法具有計算簡單、抗噪能力強等特點，故在艦船噪聲特征提取有所應(yīng)用[14]。但排列熵沒有考慮幅值問題，因此XIA 等[15]在排列熵的基礎(chǔ)上提出加權(quán)排列熵（weighted permutation entropy, WPE），然而WPE 只考慮單一尺度，忽略其他尺度的有用信息。因此YIN 等[16]提出了MWPE 描述不同尺度下時間序列的復雜性，該算法已經(jīng)廣泛應(yīng)用于故障診斷和管道側(cè)漏領(lǐng)域[17-19]。但MWPE 在粗?；^程中忽略其他粗粒化的特征信息，影響熵值的準確性。

本文提出一種基于VMD 與IMWPE 的艦船輻射噪聲特征提取方法。IMWPE 是對MWPE 的一種改進，通過平移均值法彌補MWPE 的不足。首先將艦船輻射噪聲進行VMD 分解成多個IMF，選取最優(yōu)的IMF 作為研究對象，并通過IMWPE 算法進行特征提取。最后將所提取的特征輸入PSO-SVM 進行識別分類驗證，實驗結(jié)果表明了該方法的穩(wěn)定性和優(yōu)越性。

2 基礎(chǔ)理論

2.1 VMD 理論

VMD 是繼EMD 和EEMD 方法之后的新一種方法，具體步驟如下：

1）構(gòu)造變分問題

變分模態(tài)分解算法將輸入信號分解成不同中心頻率和有限帶寬 {uk(t)}(k=1,2,...,K)的子序列。通過求解變分問題進行分解，約束變分模型如式：

式中： {ωk} 為模態(tài)函數(shù)的中心頻率； {uk}表示各個模式；uk(t) 為第k階調(diào)頻信號； δ (t)為單位脈沖函數(shù)；f(t)為輸入信號。

2）求解變分問題

使用懲罰因子 α和拉格朗日乘算子 θ (t)，再使用二次懲罰因子重建信號，使用貝葉斯推導出懲罰權(quán)重，噪聲越大懲罰權(quán)重越小。 θ(t)用于負責條件的嚴格約束，定義增廣拉格朗日L公式如下：

式(2)中采用一種稱為交替方向乘子法(ADMM)的迭代優(yōu)化方法求取式(2)中的“鞍點”，則式(1)最優(yōu)解求解步驟為：

步驟1初始化

步驟2循環(huán)n=n+1；

步 驟3對所有 ω>0 ，更新，k∈{1,2,...,K}; 更新 ωk；更新λ ；

步驟4重復步驟2 和步驟3 操作，直到滿足式(6)停止迭代，即得到K個IMF 分量

2.2 加權(quán)排列熵

排列熵算法描述了時間序列的復雜性，它只考慮時間序列的順序結(jié)構(gòu)而忽略序列的幅值信息，WPE 解決了對此不足的問題，算法如下：

1）原始時間序列X={x(i),i=1,2,...,N}進行相空間重構(gòu)

其中：m為嵌入維數(shù)；τ 為延遲時間。

2）計算每個子序列的權(quán)重值

3）計算每一種排列出現(xiàn)的概率

4）計算加權(quán)排列熵的值

2.3 多尺度加權(quán)排列熵

多尺度加權(quán)排列熵是多尺度熵與加權(quán)排列熵的結(jié)合，是用來描述多尺度時間序列復雜度的方法。MWPE在不同尺度下求取粗?；瘯r間序列的熵值結(jié)果，研究各尺度下熵值的特征，從而反映出艦船輻射噪聲的變化情況，可以比單一尺度的加權(quán)排列熵獲取更多狀態(tài)量。計算步驟如下:

1）進行數(shù)據(jù)化處理

原始時間序列X={x(i),i=1,2,...,N}進行粗?；幚?，得到序列，即

其中：s為尺度因子；[N/S] 向下取整。

2）計算每個粗?；蛄械呐帕徐?/p>

3 改進多尺度加權(quán)排列熵

MWPE 通過粗?；绞浇鉀Q了WPE 尺度單一的問題，但該算法在艦船輻射噪聲特征提取上仍存在缺陷：當尺度因子較大時，粗?；蛄袝兌?，會使熵值產(chǎn)生誤差；在粗?；^程中，忽略了時間序列上的自相關(guān)性。自相關(guān)度是指信號在同一個時間序列中，某一采樣點的瞬時值與另一個采樣點瞬時值之間的依賴程度。自相關(guān)系數(shù)的絕對值接近1 就代表相關(guān)性很強，絕對值接近0 代表相關(guān)性很弱。自相關(guān)函數(shù)公式為：

式中：n為樣本容量；h為階數(shù)；μ為該時間序列的均值；xi代表各采樣點。

通過對時間序列的分析，任意一個采樣點與其前2 個滯后變量呈強自相關(guān)度，因此MWPE 以s為尺度因子分段時，未考慮分段之間前后的自相關(guān)度，從而導致不能夠有效地進行特征提取。

針對MWPE 存在的不足，提出IMWPE 算法解決MWPE 單一粗?；瘑栴}。為了考慮相鄰時間序列之間的有用信息，本文采用均值移動法計算同一尺度下多個粗粒化時間序列的WPE 值。在較大的尺度因子下，IMWPE 可以最大程度保留原始時間序列的信息特征。由于在采集艦船輻射噪聲信號時會在采集前后給出預(yù)留值，因此艦船噪聲信號的中間部分最能反映噪聲的真實特征，在時間序列上，IMWPE 采用累加求均值法將計算資源集中在噪聲的中間部分，多次計算時間序列中段來獲取重要信息，具體操作步驟如圖1 所示。

圖1 尺度為3 時 IMWPE 粗?；绞紽ig. 1 Coarse graining mode of imwpe when the scale is 3

IMWPE 整體計算步驟如下：

將信號進行改進粗粒化過程處理，得到s組新序列：

計算s個時間序列熵的平均值，如下式：

4 特征提取識別分類方法

基于VMD 與改進多尺度加權(quán)排列熵的特征提取和識別分類流程如圖2 所示。

圖2 艦船輻射噪聲特征提取流程框圖Fig. 2 Flow chart of feature extraction of ship radiated noise

1）對獲取的三類艦船輻射噪聲仿真信號進行EMD分解，確定VMD 的分解階數(shù)；

2）VMD 分解后得到一定數(shù)量的IMF，計算各IMF 的能量，確定能量最大的IMF 為最優(yōu)IMF；

3）提取最優(yōu)IMF 的IMWPE；

4）采用PSO-SVM 方法對多類艦船噪聲VMD 分解后最優(yōu)IMF 的IMWPE 特征參數(shù)進行分類識別并得出結(jié)果。

5 實驗驗證

5.1 艦船噪聲模態(tài)分解

采用艦船輻射噪聲仿真信號模擬真實艦船信號數(shù)據(jù)，按照線譜與連續(xù)譜幅值不同分為A、B、C 三類艦船，分別進行EMD 分解，確定分解階數(shù)為8，因此設(shè)置參數(shù)K=8，采樣頻率為5 120 Hz，數(shù)據(jù)長5 120 點， α=2 000。根據(jù)設(shè)置的參數(shù)，將三類仿真信號再進行VMD分解，圖3 為3 種艦船輻射噪聲VMD 分解圖。

由圖3 看出，艦船輻射噪聲信號分解出8 個IMF和一個殘差信號，在一定數(shù)量的IMF 中，至少有一個IMF 信號幅值明顯高于其他分量，一般將幅值最高的IMF 稱為信號的主能量，可以反映艦船輻射噪聲信號的主要特征。根據(jù)文獻[4]提出的相關(guān)公式，計算每個IMF。設(shè)第k階模態(tài)有N個采樣點，第i個采樣點瞬時頻率為fki，瞬時振幅為bki，瞬時強度為Qki=b2ki，IMF 強度為：

圖3 三種艦船輻射噪聲VMD 分解圖Fig. 3 VMD decomposition of three types of ship radiated noise

則能量強度最大的IMF 為：

根據(jù)式(17) 和式(18) 求出能量最強的IMF，用EIMF 表示，則A 類與B 類的EIMF 階數(shù)都為3，C 類的EIMF 階數(shù)為1。

5.2 IMWPE 特征提取

3 種艦船輻射噪聲信號通過VMD 分解后，根據(jù)式(17)和式(18)提取能量最大的EIMF 作為研究對象，然后提取A 類、B 類、C 類艦船輻射噪聲的加權(quán)排列熵，3 種艦船噪聲的樣本數(shù)各為150 個，樣本總量為450 個，實驗結(jié)果如圖4 和表1 所示。

圖4 三類艦船輻射噪聲加權(quán)排列熵Fig. 4 Weighted permutation entropy of radiated noise of three types of ships

表1 三類艦船輻射噪聲加權(quán)排列熵Tab. 1 Weighted permutation entropy of radiated noise of three types of ships

由圖4 和表1 可知，加權(quán)排列熵特征分布均勻，表明了加權(quán)排列熵算法具有較強的穩(wěn)定性。由于艦船輻射噪聲產(chǎn)生的機理不同，導致3 種艦船噪聲的時間復雜度存在一定的差異。根據(jù)平均加權(quán)排列熵值得出，B 類的噪聲信號復雜度最高，其次是A 類與C 類，從加權(quán)排列熵特征分布來看，雖然A 類與B 類的加權(quán)排列熵特征分布較為重疊，但加權(quán)排列熵特征分布還是可以很好地區(qū)分3 種艦船輻射噪聲。

加權(quán)排列熵是從單尺度方面考慮時間序列復雜度問題，在此引入多尺度加權(quán)排列熵從多個尺度了解艦船輻射噪聲特性情況。由于多尺度加權(quán)排列熵粗?；^程忽略其他序列信息問題，本文提出改進多尺度加權(quán)排列熵算法，提取出的EIMF 分別使用MWPE 和IMWPE 計算后進行對比。一般來說，嵌入維數(shù)m通常取值為～7，m取值太小，重構(gòu)時間序列中的重要信息量會減少，m值太大，計算時間過長，不能揭示信號變化，本文設(shè)置m=5。信號數(shù)據(jù)長度也會影響熵值的精確度，信號越長，熵值精確度越高，太多的數(shù)據(jù)點會消耗大量的計算資源，信號過短，則不能準確反映信號的特征變化，文獻[19]建議取值為N≥5m!，因此本文取N=5 120。延遲時間對熵值計算影響較小，則設(shè)τ=1，根據(jù)研究需要設(shè)s=10。圖5 為提取的EIMF 中三類艦船噪聲的隨機一種熵值對比圖。

圖5 隨機一種三類艦船輻射噪聲熵值對比圖Fig. 5 Comparison diagram of entropy value of radiated noise of three types of ships at random

可以看出，3 種艦船在各尺度因子上分布均勻，但提取的MWPE 熵值曲線比IMWPE 波動較大。這是因為IMWPE 提取熵值采用改進的同一尺度多個粗?；椒ǎ萂WPE 算法更為穩(wěn)定。表2 給出2 種熵值的標準差對比，說明IMWPE 比MWPE 更具有穩(wěn)定性和優(yōu)越性。

表2 三類艦船輻射噪聲不同尺度熵值標準差Tab. 2 Standard deviation of entropy values at different scales of radiated noise of three types of ships

5.3 分類識別檢測

為了進一步比較3 種艦船輻射噪聲特征提取方法，將隨機提取的150 個樣本中70%作為訓練樣本，30%作為測試樣本，并把提取的WPE、MWPE 和IMWPE 采用PSO-SVM 進行分類識別驗證。仿真實驗中，設(shè)置粒子群算法的種群規(guī)模為20，最高迭代50 次，學習因子C1和C2均為1.5[20]，分類結(jié)果如表3 所示。

表3 特征提取算法的PSO-SVM 分類識別結(jié)果Tab. 3 PSO-SVM classification and recognition results of feature extraction algorithm

可以看出，WPE 特征提取算法總識別率較低，由圖4 可知，A 類與B 類的WPE 值略有重疊，并且只從單一尺度方面考慮，A 類與B 類在識別分類中產(chǎn)生誤差導致總識別率降低。從圖5 和表3 看出，MWPE 與IMWPE 都是從多尺度角度分析， C 類艦船基本上可以完全識別。因為IMWPE 采用平移均值法解決了單一尺度的粗?；瘑栴}，所以總識別率高于MWPE 和WPE 特征提取算法。

6 結(jié) 語

本文從時間復雜度的角度去分析艦船噪聲，提出了一種基于VMD 與IMWPE 相結(jié)合的艦船輻射噪聲特征提取方法，先將信號通過VMD 進行分解，選取能量最大的EIMF 作為研究對象，并計算IMWPE 作為特征向量，采用PSO-SVM 進行分類識別。經(jīng)實驗數(shù)據(jù)表明：由于不同種艦船輻射噪聲產(chǎn)生的機理不同，它們的WPE 也大不相同；MWPE 忽略了單一粗粒化相鄰序列之間的關(guān)系，IMWPE 算法解決了MWPE 不足的問題，具有良好的穩(wěn)定性和一致性；在識別分類方面中，提取的IMWPE 識別率效果最好。