李宗吉，張西勇

(海軍工程大學(xué) 新兵器技術(shù)與應(yīng)用研究所，湖北 武漢430033)

0 引 言

在傳統(tǒng)的水聲信號(hào)檢測(cè)理論中，干擾背景多被假設(shè)為高斯分布。但在實(shí)際應(yīng)用中，尤其是主動(dòng)聲吶條件下，由于海洋噪聲背景和混響的聯(lián)合作用，干擾背景明顯偏離了高斯分布，具有顯著尖峰脈沖特性，其概率密度函數(shù)的衰減過(guò)程比高斯分布要慢，即呈現(xiàn)代數(shù)拖尾。

廣義中心極限定理指出，無(wú)限個(gè)無(wú)限方差的獨(dú)立分布同隨機(jī)變量的極限分布是一種穩(wěn)定分布，而α 穩(wěn)定分布是唯一的滿(mǎn)足廣義中心極限定理的分布族，比高斯分布有更廣泛的適用性。α 穩(wěn)定分布的統(tǒng)計(jì)特性由其特征函數(shù)的4 個(gè)參數(shù)決定，其中最重要的參數(shù)為特征指數(shù)α ∈(0，2］，它決定著α 穩(wěn)定分布概率密度函數(shù)拖尾的厚度。α 的取值越小，概率密度函數(shù)拖尾越厚，表明隨機(jī)信號(hào)的尖峰脈沖特性越明顯。因此，α 穩(wěn)定分布可以很好地對(duì)水聲干擾背景噪聲信號(hào)進(jìn)行描述。此外，α 能夠在(0，2］之間任意取值，因此α 穩(wěn)定分布可以更加靈活地對(duì)水聲干擾背景噪聲信號(hào)進(jìn)行描述。特別地，當(dāng)α=2 時(shí)，α 穩(wěn)定分布與高斯分布完全相同，即高斯分布是α 穩(wěn)定分布的特例。α 穩(wěn)定分布都可以由SαS 分布變換得到，不失一般性，本文將主要討論基于SαS 分布噪聲建模的白化濾波方法。

SαS 分布沒(méi)有顯性的概率密度函數(shù)表達(dá)式，但其特征函數(shù)可由式(1)明確表達(dá):

1 共變譜與穩(wěn)定白噪聲

對(duì)于SαS 分布過(guò)程，由于不存在二階矩，其功率譜不存在。為了進(jìn)行頻域處理，引入共變譜的概念。

對(duì)于一個(gè)嚴(yán)平穩(wěn)復(fù)SαS 過(guò)程{x(t);－∞＜t＜∞}，總存在如下一種譜表示:

其中ζ(ω)的獨(dú)立增量滿(mǎn)足如下關(guān)系:

式中:C(p，α)為依賴(lài)于p 與α 的常量;Φ(ω)為x(t)的一種稱(chēng)作譜密度的非負(fù)函數(shù)。當(dāng)α=2 時(shí)，譜密度Φ(ω)就是對(duì)應(yīng)高斯過(guò)程的功率譜;0 ＜α ＜2 時(shí)，譜密度Φ(ω)并不表示x(t)的功率譜。但在一些如線性預(yù)測(cè)與濾波的應(yīng)用中，它起著和功率譜一樣類(lèi)似的作用，也就是說(shuō)，譜密度Φ(ω)完全可以描述x(t)的一種頻率分布。當(dāng)且僅當(dāng)dζ(ω)具有全向或者旋轉(zhuǎn)不變?cè)隽繒r(shí)，過(guò)程{x(t)}為嚴(yán)平穩(wěn)過(guò)程，并且有如下穩(wěn)定過(guò)程的共變函數(shù)的表達(dá)式:

因此，穩(wěn)定分布過(guò)程的自共變函數(shù)與二階過(guò)程的自協(xié)方差函數(shù)有著同樣的意義，它與譜密度Φ(ω)構(gòu)成一個(gè)傅里葉變換對(duì)。根據(jù)自共變函數(shù)的定義:

可以建立如下準(zhǔn)則:當(dāng)且僅當(dāng)自共變函數(shù)Γx(τ)在理論上滿(mǎn)足Γx(τ)=γxδ(τ)時(shí)，{x(t)}為穩(wěn)定分布白噪聲過(guò)程，其中δ(τ)是沖激函數(shù)。自共變函數(shù)是對(duì)二階過(guò)程的自相關(guān)函數(shù)的廣義化，對(duì)自共變函數(shù)Γx(τ)進(jìn)行傅里葉變換可以得到一種分?jǐn)?shù)階譜——共變譜(Covariation Spectrum，CS):

若將自共變函數(shù)換成自相關(guān)函數(shù)，則得到功率譜密度函數(shù)。由于自相關(guān)函數(shù)是自共變函數(shù)當(dāng)α=2時(shí)的特例，因此，二階過(guò)程的功率譜是共變譜的特例，共變譜是功率譜的廣義化形式。在離散時(shí)域，一個(gè)嚴(yán)平穩(wěn)復(fù)SαS 序列{x(t)}也總存在如下一種譜表示:

因此，一個(gè)穩(wěn)定分布白噪聲序列{x(n)}的自共變序列Γx(n)=［x(m)，x(m－n)］α應(yīng)為序列γxδ(n)，其共變譜理論上應(yīng)是平坦均勻直線。

2 基于廣義Yule-Walker 方程的白化濾波

對(duì)SαSAR(P)序列x(n)作如下改寫(xiě):

其中u(n)是一個(gè)特征指數(shù)為α，分散系數(shù)為γ 的SαS 分布白噪聲序列。

x(n)可以看成是無(wú)窮多個(gè)獨(dú)立同分布穩(wěn)定白噪聲變量的線性組合，x(n)與u(n)具有相同的特征指數(shù)α。

對(duì)于嚴(yán)平穩(wěn)穩(wěn)定分布序列x(n)，其離差應(yīng)該滿(mǎn)足:

可以得到x(n)的離差如下式所示:

因此，嚴(yán)平穩(wěn)穩(wěn)定分布序列x(n)是與u(n)具有相同的特征指數(shù)α、分散系數(shù)為γx(n)的穩(wěn)定有色噪聲。

當(dāng)滿(mǎn)足條件1 ＜α ≤2 時(shí)，在差分方程式(8)的兩端分別取關(guān)于x(n－m)的共變函數(shù)，可以得到:

可以得到:

即

式中:λu(n)x(n－m)為u(n)和x(n－m)的共變系數(shù)，記為λux(m);λx(n)x(n－m)為x(n)的自共變系數(shù)，記為λ(m)。

共變系數(shù)λux(m)和濾波器的脈沖響應(yīng)相關(guān)，可得

上式最后一步利用了共變的偽線性。u(n)是穩(wěn)定白噪聲，滿(mǎn)足如下條件:

結(jié)合式(12)、式(14)～式(16)可得:

因此，對(duì)于一個(gè)SαSAR(P)過(guò)程，自共變系數(shù)存在如下關(guān)系:

SαSAR(P)參數(shù)矢量a 可以通過(guò)求解Yule-Walker 方程組得到。

將使用自共變系數(shù)描述的Yule-Walker 方程稱(chēng)為廣義(或分?jǐn)?shù)低階)Yule-Walker 方程，記為Ca=p。穩(wěn)定白噪聲離差γu(n)的估計(jì)值^γu(n)可以用下面的等式計(jì)算得到:

對(duì)于高斯過(guò)程而言，自共變系數(shù)矩陣C 是對(duì)稱(chēng)和正定的。但對(duì)于非高斯過(guò)程而言，C 通常是非對(duì)稱(chēng)，甚至是奇異的。

3 自共變系數(shù)的估計(jì)

基于廣義Yule-Walker 方程頻譜估計(jì)算法的性能除了取決于算法本身，自共變系數(shù)λ(i)的估計(jì)也至關(guān)重要。在實(shí)際應(yīng)用中，需要根據(jù)觀測(cè)樣本值估計(jì)自共變系數(shù)矩陣C和加權(quán)矢量p，即需要估計(jì)自共變系數(shù)λ(i)［2］。

對(duì)于任何1 ≤p ＜α，聯(lián)合α 穩(wěn)定分布隨機(jī)變量X和Y 的共變系數(shù)可定義如下:

對(duì)于獨(dú)立觀測(cè)值(X1，Y1)，…，(XN，YN)，共變系數(shù)的估計(jì)值如下式所示:

為了簡(jiǎn)化計(jì)算量，可以令p=1，于是得到:

如果觀測(cè)值獨(dú)立，則上式一致收斂。然而，對(duì)于AR 過(guò)程，其輸出序列并不滿(mǎn)足獨(dú)立性，但FLOM估計(jì)依然滿(mǎn)足一致收斂［1］。在實(shí)際應(yīng)用中，對(duì)于FLOM 估計(jì)，只需將式(22)和式(23)分別作如下改寫(xiě):

4 預(yù)白化濾波及實(shí)驗(yàn)數(shù)據(jù)檢驗(yàn)

在實(shí)際的主動(dòng)信號(hào)檢測(cè)中，往往采用滑動(dòng)窗技術(shù)，對(duì)每次的窗內(nèi)數(shù)據(jù)進(jìn)行“譜估計(jì)→預(yù)白化→檢測(cè)處理”。預(yù)白化濾波器的構(gòu)建基于譜估計(jì)，得到AR 參數(shù)估計(jì)后，即可基于該估計(jì)建立起預(yù)白化濾波器。

顯然，這是具有MA 結(jié)構(gòu)的濾波器。它可以從觀測(cè)樣本序列{xn}估計(jì)出激勵(lì)序列{un}，式(27)亦被稱(chēng)為激勵(lì)估計(jì)器。一般來(lái)說(shuō)，AR 濾波有時(shí)亦被看作基于前P 個(gè)的觀測(cè)樣點(diǎn)預(yù)測(cè)當(dāng)前樣點(diǎn)的線性預(yù)測(cè)器，這樣{un}亦是預(yù)測(cè)誤差，該式又被稱(chēng)為線性預(yù)測(cè)誤差估計(jì)器。

使用［0601 三亞海試］中一段基元級(jí)混響數(shù)據(jù)，發(fā)射信號(hào)為800 ～1 000 Hz 的HFM 脈沖，脈寬為4 s，采樣頻率為6 kHz。帶通濾波已向發(fā)射通帶兩端各擴(kuò)展50 Hz，即為750 ～1 050 Hz。數(shù)據(jù)長(zhǎng)度與脈寬同長(zhǎng)。使用廣義Yule-Walker 方程算法對(duì)其進(jìn)行AR(256)譜估計(jì)和預(yù)白化處理，設(shè)定p=1，處理結(jié)果如圖1所示。由圖1(b)可以看到，在全帶內(nèi)，觀測(cè)樣本擬合CS 曲線與觀測(cè)樣本統(tǒng)計(jì)CS 曲線重合得很好，經(jīng)過(guò)預(yù)白濾波后，通帶內(nèi)成功實(shí)現(xiàn)了白化。

圖1 一段海試數(shù)據(jù)預(yù)白化前后的CSFig.1 Covariation spectrum of a segment of sea trial data before and after prewhitened

需要指出的是，對(duì)于AR 模型階數(shù)P 的選擇目前尚無(wú)有效的定階方法，因此在實(shí)際應(yīng)用中，一般是憑借經(jīng)驗(yàn)來(lái)指定AR 模型的階數(shù)。在數(shù)據(jù)有效長(zhǎng)度損失和運(yùn)算時(shí)長(zhǎng)允許的范圍內(nèi)，為了達(dá)到理想的預(yù)白效果，通常將階數(shù)指定的偏高一些。

5 結(jié) 語(yǔ)

本文主要研究基于分?jǐn)?shù)低階統(tǒng)計(jì)量AR 模型參數(shù)的估計(jì)方法。首先，介紹共變譜、α 譜的概念以及穩(wěn)定白噪聲的判別準(zhǔn)則，并采用實(shí)際海試數(shù)據(jù)檢驗(yàn)分析了共變譜的性能;討論一種基于α 譜的頻域廣義白化濾波方法。然后，從α 譜的角度出發(fā)，推導(dǎo)出基于自共變系數(shù)的廣義Yule-Walker 方程。通過(guò)仿真深入討論了基于廣義Yule-Walker 方程頻譜估計(jì)算法中，各參數(shù)對(duì)算法性能的影響。

本文第2 個(gè)主要研究?jī)?nèi)容是基于譜估計(jì)的預(yù)白化技術(shù)。對(duì)于AR 模型擬合的有色數(shù)據(jù)而言，預(yù)白實(shí)際上就是一個(gè)基于各AR 參數(shù)估計(jì)建立起來(lái)的MA濾波，亦即AR 激勵(lì)估計(jì)，只要模型參數(shù)估計(jì)足夠準(zhǔn)確，就能得到理想的預(yù)白效果。采用大量水池實(shí)驗(yàn)、湖試、海試實(shí)際數(shù)據(jù)，考察了分?jǐn)?shù)低階白化濾波器的性能。

［1］SHAO M，NIKIAS C L.Signal processing with fractional lower order moments:stable processes and their applications［J］.Proceedings of the IEEE，1993，81(7):986-1010.

［2］查代奉.基于穩(wěn)定分布白噪聲的信號(hào)處理新方法研究［D］.大連:大連理工大學(xué)，2006:56-67.

ZHA Dai-feng.Study on new methods of signal processing based on stable white nolse［D］.Dalian:Dalian University of Technology，2006:56-67.

［3］WANG Z，HE Z，CHEN J.Robust time delay estimation of bioelectric signals using least absolute deviation neural Network［J］.IEEE Transacions on Biomedical Engineering，2005，52(3):454-462.

［4］MAKAR A.Estimation of the time delay of hydroacoustic signals for passive location of underwater objects［J］.Archives of Acoustics，2004，29(3):435-445.

［5］SO C.Analysis of an adaptive algorithm for unbiased multipath time delay estimation［J］.IEEE Transactions on Aerospace and Electronic Systems，2003，39(3):776-780.