劉 明水鵬朗

(西安電子科技大學(xué)雷達信號處理國家重點實驗室 西安 710071)

基于功率中值和歸一化采樣協(xié)方差矩陣的自適應(yīng)匹配濾波檢測器

劉 明*水鵬朗

(西安電子科技大學(xué)雷達信號處理國家重點實驗室 西安 710071)

在非均勻海雜波環(huán)境中，參考單元中的異常單元限制了采樣協(xié)方差矩陣(SCM)的估計性能，從而影響了傳統(tǒng)自適應(yīng)匹配濾波(AMF)檢測器的檢測性能。而考慮刪除異常單元的方法在參考單元數(shù)目有限時可能導(dǎo)致矩陣的奇異性。鑒于此，在保持參考單元數(shù)目不變條件下，該文設(shè)計一種基于功率中值與歸一化采樣協(xié)方差矩陣(MNSCM)的估計方法，并將其與匹配濾波器(MF)相結(jié)合，構(gòu)造一種新型的自適應(yīng)匹配濾波檢測器。與傳統(tǒng)的自適應(yīng)匹配濾波相比，該文設(shè)計的檢測器在實測和仿真海雜波數(shù)據(jù)條件下均具有明顯的性能優(yōu)勢。

目標(biāo)檢測；海雜波；自適應(yīng)匹配濾波器；采樣協(xié)方差矩陣；功率中值和歸一化采樣協(xié)方差矩陣

1 引言

海雜波背景下的動目標(biāo)檢測問題因其具有重要的戰(zhàn)略意義和民用價值，一直以來備受國內(nèi)外學(xué)者關(guān)注[1?4]。由于海雜波具有平均功率大，時空相關(guān)性強和統(tǒng)計特性復(fù)雜等特點，使其成為目標(biāo)檢測問題中的一個研究難點。雖然海雜波在短時間內(nèi)可認(rèn)為是平穩(wěn)的，但從長時間來看，海雜波具有明顯的空時非平穩(wěn)性[5,6]。因此，對于海雜波背景下的動目標(biāo)檢測問題，通常會選用自適應(yīng)檢測算法處理此類檢測問題。

自適應(yīng)匹配濾波[7]作為一種常見自適應(yīng)檢測算法。它使用了各個距離單元的多脈沖匹配結(jié)果作為檢測統(tǒng)計量，從而實現(xiàn)目標(biāo)檢測。但使用該方法之前首先要對待測單元的雜波協(xié)方差矩陣進行估計。在雜波均勻假設(shè)下，即在待測單元與其鄰近的距離單元(參考單元)共用同一個協(xié)方差矩陣條件下，可以使用由參考單元構(gòu)成的采樣協(xié)方差矩陣(Sample Covariance Matrix, SCM)估計待測單元的雜波協(xié)方差矩陣。然而，均勻雜波環(huán)境假設(shè)在實際檢測下往往并不滿足。

在雷達高分辨條件下或低擦地角觀測下，雜波數(shù)據(jù)的呈現(xiàn)出很強的非高斯特性，雜波中含有大量類似目標(biāo)的尖峰信號[8,9]，這種尖峰信號的存在會影響SCM的估計性能[10]。關(guān)于高分辨率條件下雜波的統(tǒng)計特性，國內(nèi)外諸多學(xué)者對其進行了詳細(xì)地研究。其中，文獻[5]分析了X波段不同分辨率條件下雜波的統(tǒng)計特性；文獻[11]分別采用幾種經(jīng)典統(tǒng)計模型給出了L波段雜波在不同環(huán)境參數(shù)條件下的幅度最優(yōu)分布模型。雖然不同條件下獲得的雜波幅度分布模型有所不同，但這些分布幾乎都可以歸類為同一個分布，即復(fù)合高斯分布(在常見的雜波模型中，Lognormal分布不屬于此類)。然而，實際工作環(huán)境中的復(fù)合高斯分布的統(tǒng)計特性是事先未知的，尤其是雜波的協(xié)方差矩陣，這就需要我們從數(shù)據(jù)中進行估計。由于高分辨條件下的雜波不滿足均勻性假設(shè)，并且實際檢測中并不清楚哪個距離單元中含有目標(biāo)信號。因此，選取的參考單元中可能含有目標(biāo)或者雜波特性明顯不同于待測單元的單元，本文稱這樣的單元為異常單元。當(dāng)使用含有異常單元的參考單元估計協(xié)方差矩陣時，勢必使SCM的出現(xiàn)性能損失。

為了避免異常單元對SCM估計性能的影響，文獻[12]使用一種數(shù)據(jù)自適應(yīng)篩選(Data-Adaptive Selection, DAS)方法獲取均勻雜波數(shù)據(jù)。文獻[13]采用兩步自動刪除平均水平的方法獲得均勻雜波數(shù)據(jù)。上述方法的共同點是對參考單元作預(yù)處理，剔除其中含有的異常單元，以獲取相對均勻雜波數(shù)據(jù)。然而，當(dāng)參考單元數(shù)目較少時，刪除其中一些異常單元也會使SCM的估計性能下降甚至出現(xiàn)協(xié)方差矩陣的奇異性。在不改變參考單元數(shù)目條件下，本文設(shè)計一種基于功率中值和歸一化采樣協(xié)方差矩陣(Median Normalized Sample Covariance Matrix, MNSCM)的協(xié)方差矩陣估計方法。它能夠有效抑制非均勻雜波環(huán)境下由于不同距離單元間功率起伏對協(xié)方差矩陣估計性能帶來的影響，與此同時，它保留雜波的功率信息，使其能夠適用于匹配濾波(Matched Filter, MF)檢測器。

實驗部分分別采用實測和仿真的雜波數(shù)據(jù)對基于中值和歸一化采樣協(xié)方差矩陣的自適應(yīng)匹配濾波(Median Normalized Sample Covariance Matrix-Adaptive Matched Filter, MNSCM-AMF)檢測器的有效性作了驗證。對于實測數(shù)據(jù)，選取參考單元時人為將一些功率大的異常單元考慮在內(nèi)；對于仿真數(shù)據(jù)，通過調(diào)高參考單元中的某些距離單元的功率，以便獲得一種含有異常單元的雜波數(shù)據(jù)。對于這兩種數(shù)據(jù)，MNSCM-AMF檢測器均獲得了比基于采樣協(xié)方差矩陣的自適應(yīng)匹配濾波(SCM-AMF)檢測器更好的檢測結(jié)果。同時，通過調(diào)整仿真的雜波環(huán)境的非均勻程度，即異常單元的功率大小和所占的比例，定量地分析MNSCM-AMF檢測器在不同程度的非均勻雜波環(huán)境相比SCM-AMF檢測器具有更優(yōu)檢測結(jié)果的原因。

文章其余部分安排如下，第2節(jié)對檢測問題進行簡要描述；第3節(jié)首先引入傳統(tǒng)的自適應(yīng)匹配濾波檢測器，隨后介紹基于中值和歸一化采樣協(xié)方差矩陣的自適應(yīng)匹配濾波檢測器；第4節(jié)通過實驗方式驗證MNSCM-ANMF檢測器相比SCM-ANMF檢測器的檢測性能優(yōu)勢，并在仿真數(shù)據(jù)條件下，定量分析了兩種檢測器性能差異產(chǎn)生的原因；第5節(jié)總結(jié)全文。

2 檢測問題描述

如果復(fù)信號z(t)表示待測單元在觀測時間內(nèi)的接收到的海雜波的基帶形式，對這一基帶信號進行采樣得到接收向量z=zI+jzQ，其中，zI和zQ分別代表接收向量z的同相部分和正交部分。要判斷采樣所得的向量z中是否含有目標(biāo)信號相當(dāng)于式(1)所示的二元假設(shè)檢驗：

其中，z, s和c分別代表待測單元中的N維接收向量，感興趣的目標(biāo)信號和待測單元中的雜波向量；zk或者ck,k=1,2,…,K ，是來自參考單元中不含目標(biāo)信號的一組接收向量，并假定在H0假設(shè)下zk與z具有相同的協(xié)方差矩陣。在此假定下，當(dāng)待測單元接收向量z中的雜波協(xié)方差矩陣未知時，可以使用參考單元數(shù)據(jù)zk估計待測單元中的雜波協(xié)方差矩陣。這里假定K＞2N以滿足輸出的信雜比的損失小于3 dB。

在相干處理中，目標(biāo)信號s可以看作一個相干脈沖序列，即s=αp，其中，α是與傳輸信道和雷達橫截面積有關(guān)的常量，p是多普勒導(dǎo)向矢量，p=[1 exp(j2πfdTr) … exp(j2π(N?1)fdTr)]T，其中，fd是目標(biāo)的多普勒頻率，Tr是雷達脈沖重復(fù)間隔。對于雜波向量c，采用復(fù)合高斯模型進行建模，即雜波可以看作是復(fù)高斯分布的散斑分量u(n)和非負(fù)的紋理分量τ(n)的乘積，其中，相比散斑分量，紋理分量具有更長的去相關(guān)時間[9,14]。因此，當(dāng)觀測時間在幾個相干處理時間內(nèi)時，復(fù)合高斯分布的紋理分量τ(n)可以近似為某一常數(shù)，此時的復(fù)合高斯模型退化為球不變隨機向量(Spherically Invariant Random Vector SIRV)模型，即c=u。

3 兩種自適應(yīng)匹配濾波檢測器

針對式(1)中的檢測問題，自適應(yīng)匹配濾波(AMF)檢測器給出如下檢測方案。當(dāng)雜波協(xié)方差矩陣已知時，由似然比不難得到AMF檢測器的檢測統(tǒng)計量為

其中，M是真實的雜波協(xié)方差矩陣, η是滿足指定虛警概率條件下的判決門限。然而，在實際工作環(huán)

將式(2)中的M替換為式(3)，從而得到均勻雜波環(huán)境中傳統(tǒng)的AMF檢測器。境中，雜波的協(xié)方差矩陣往往是未知的，需要從參考單元數(shù)據(jù)中估計得到。在雜波均勻假設(shè)前提下，可以使用SCM得到真實的雜波協(xié)方差矩陣的估計，即式(3)。

式(4)中的檢測器具有歸一化匹配濾波器的形式，而分母上的歸一化因子使得AMF檢測器具有恒虛警的特性。文獻[7]指出：AMF的檢測統(tǒng)計量在H0假設(shè)下與真實的雜波協(xié)方差矩陣相互獨立，是一種恒虛警率檢測方法。

然而，在高分辨雷達觀測下，絕對均勻的雜波環(huán)境幾乎不存在，這就意味采用SCM估計真實的雜波協(xié)方差矩陣會對檢測器的性能造成損失。實測數(shù)據(jù)顯示：高分辨率低視角條件下獲得海雜波數(shù)據(jù)具有空時非平穩(wěn)性，即海雜波在短的時間內(nèi)是平穩(wěn)的，長的時間內(nèi)是非平穩(wěn)的；短的距離上是均勻的，長的距離上是非均勻的[15?17]。因此，在非均勻雜波環(huán)境中，選取的用于估計雜波協(xié)方差矩陣的參考單元數(shù)據(jù)中可能含有類似目標(biāo)或者雜波特性明顯異于待測單元的單元(統(tǒng)稱這樣的單元為異常單元)，這就使得SCM出現(xiàn)性能損失，從而影響了AMF檢測器的檢測結(jié)果。

為了避免這種異常單元對估計的雜波協(xié)方差矩陣的影響，文獻[12,13]均考慮采用刪除異常單元的方法，這類方法在一定程度取得相對較好的檢測結(jié)果，但它們均未考慮在有限參數(shù)單元數(shù)目時，刪除異常單元可能導(dǎo)致矩陣的奇異性的問題。在保持參考單元數(shù)目條件下，該文設(shè)計了一種基于功率中值和歸一化采樣協(xié)方差矩陣(MNSCM)的估計方法，即

其中，符號Med{·}表示取中值運算，式(5)中的歸一化處理去除了雜波紋理分量的影響，很好地抑制了因為各個參考單元數(shù)據(jù)功率起伏對協(xié)方差矩陣的估計性能的影響；與此同時，它又保留了雜波的功率信息，使其能夠與自適應(yīng)匹配濾波檢測器相結(jié)合。將式(2)中的真實的雜波協(xié)方差矩陣替換為式(5)估計的雜波協(xié)方差矩陣，得到基于功率中值和歸一化采樣協(xié)方差矩陣的自適應(yīng)匹配濾波檢測器

不難證明，與式(4)中的AMF檢測器相同，式(6)對應(yīng)的AMF檢測器也是一種恒虛警率檢測方法。如果從形式上分析兩種自適應(yīng)檢測器，可以看出：在任何類型的雜波環(huán)境中，SCM估計器認(rèn)為每個參考單元對最終估計的雜波協(xié)方差矩陣的貢獻是均等，即所有的參考單元的對應(yīng)的權(quán)值均為1/K，顯然，這種做法只有在均勻雜波環(huán)境中才是合理的；而MNSCM估計器在含有功率大的異常單元的雜波環(huán)境中，會自適應(yīng)地賦予那些功率大的異常單元一個較小的權(quán)值，這相當(dāng)于對那些異常單元作相應(yīng)的懲罰，從而減弱了異常單元對協(xié)方差矩陣估計的影響。因此，在含有功率大的異常單元的參考單元中，相比SCM-AMF檢測器，MNSCM-AMF檢測器會帶來更好的檢測性能。而對于均勻的雜波環(huán)境，MNSCM-AMF檢測器與SCM-AMF檢測器的形式相同，兩種檢測器應(yīng)該具有同等的檢測性能。

4 實驗結(jié)果與性能分析

本節(jié)分別采用實測數(shù)據(jù)和仿真數(shù)據(jù)對SCM-AMF和MNSCM-AMF的檢測性能進行評價。實測數(shù)據(jù)(1990227_215640_ANTSTEP.mat,VV極化)由(Ice Multi-parameter Imaging X)IPIX雷達[18]于1998在加拿大的安大略湖上采集得到，并且每個距離單元對應(yīng)的接收數(shù)據(jù)為一組零均值、單位方差的復(fù)序列。表1總結(jié)了IPIX雷達參數(shù)和所使用數(shù)據(jù)的詳細(xì)信息。通過對實驗數(shù)據(jù)作分析可知，雜波數(shù)據(jù)的幅度分布可以采用不同形狀參數(shù)的K分布進行描述。圖1(a)給出了除了目標(biāo)單元及其周圍4個距離單元以外的其余距離單元采用K分布擬合時對應(yīng)的形狀參數(shù)。為了說明K分布的有效性，圖1(b)中給出了相應(yīng)距離單元對應(yīng)的Kolmogorov-Smirnov (KS) 距離。由圖1可知，一方面，K分布能很好地描述雜波的幅度統(tǒng)計特性；另一方面，K分布的形狀參數(shù)并不是距離單元的常函數(shù)，各個距離單元的幅度須用不同形狀參數(shù)的K分布描述，即雜波呈現(xiàn)非均勻特性。而仿真數(shù)據(jù)根據(jù)實測數(shù)據(jù)的一階、二階統(tǒng)計特性產(chǎn)生，其中，仿真數(shù)據(jù)幅度服從K分布，其形狀參數(shù)等于圖1(a)中各個距離單元的形狀參數(shù)的均值；協(xié)方差矩陣等于實測數(shù)據(jù)不同時間段上的協(xié)方差矩陣的平均。K分布的形狀參數(shù)的估計采用經(jīng)典的矩估計方法，而雜波協(xié)方差矩陣?yán)盟c自相關(guān)函數(shù)的關(guān)系，它可以根據(jù)式(7)的自相關(guān)函數(shù)得到。

表1 IPIX 雷達參數(shù)列表[18]

其中，zk,l(n)為第l個時間段上第k個距離單元上回波數(shù)據(jù)的第n個樣本，q是同一距離單元數(shù)據(jù)中的接收信號的時間差，L由每個距離單元中的脈沖數(shù)目除以脈沖累積數(shù)N得到，星號“*”表示復(fù)數(shù)的共軛。

實驗中，取虛警概率為1.0e-3，分別添加仿真目標(biāo)信號到實測數(shù)據(jù)和仿真數(shù)據(jù)，觀察兩種檢測器的檢測性能，其中，實驗門限采用Monte Carlo 方法計算獲得。對于實測數(shù)據(jù)，保留原有的目標(biāo)信號及其周圍4個距離單元數(shù)據(jù)，而仿真目標(biāo)添加到其余的距離單元中，這樣，對于仿真目標(biāo)所在距離單元而言，真實目標(biāo)所在距離單元數(shù)據(jù)即為異常單元數(shù)據(jù)。以脈沖累積數(shù)N=8為長度，將每個距離單元劃分為L=7500的不重合的短向量。每個短向量都會作為H0假設(shè)下的待測向量，而同一時間段上的其余距離單元中短向量會被看作是用來估計雜波協(xié)方差矩陣的參考單元。因此，用來估計門限的實驗次數(shù)為7500×28，這等于H0假設(shè)下檢測統(tǒng)計量數(shù)目的總和。將所有的檢測統(tǒng)計量按降序重新排列，相應(yīng)的門限值將是重新排列后的檢測統(tǒng)計量的第210個(7500 ×28 × 1.0e-3)。仿真數(shù)據(jù)采用相同的方法獲得虛警概率為1.0e-3時的門限值。

實驗中所加仿真目標(biāo)信號的形式為其中，fd是目標(biāo)的多普勒頻率，初相?是服從區(qū)間?[π,π]的均勻分布，a是與傳播路徑和雷達橫截面積有關(guān)的常量。定義信號的平均信雜比為

4.1 實測數(shù)據(jù)條件下兩種檢測器性能對比

當(dāng)N=8, K=28，平均信雜比為-5 dB時，兩種檢測器檢測結(jié)果如圖2所示。由圖可知，采用MNSCM的自適應(yīng)匹配濾波檢測器相比傳統(tǒng)的自適應(yīng)匹配濾波檢測器具有明顯的性能優(yōu)勢。這是因為MNSCM估計方法很好地抑制了參考單元中的異常單元對雜波協(xié)方差矩陣估計的影響，從而改善了自適應(yīng)匹配濾波檢測器的檢測性能。圖3給出了實驗所用數(shù)據(jù)的歸一化功率譜與頻率之間的關(guān)系。

圖1 各距離單元K分布擬合形狀參數(shù)和K分布擬合時的KS距離

圖2 平均信雜比為-5 dB時，MNSCMAMF與SCM-AMF性能對比

圖3 實驗所用實測數(shù)據(jù)的 歸一化功率譜

對于指定的多普勒頻率處，兩種自適應(yīng)檢測器的檢測概率隨信雜比的變化關(guān)系如圖4所示。圖中給出了3個代表性目標(biāo)多普勒頻率，即頻率-200 Hz, 0 Hz 和80 Hz時，兩種檢測器的檢測結(jié)果。在這3個頻率處，MNSCM-AMF檢測器的性能均優(yōu)于SCM-AMF檢測器，而3個頻率處的檢測性能曲線的差異是因為雜波功率在頻域的非均勻分布引起的。由圖4可知，當(dāng)目標(biāo)多普勒頻率靠近功率譜的峰值所在頻率時，由于此處雜波的功率密度較大，兩種檢測器的檢測性能相比其他頻率處下降了許多。

4.2 仿真數(shù)據(jù)條件下兩種檢測器性能對比

為了說明當(dāng)參考單元中存在功率大的異常樣本時，即雜波環(huán)境不均勻時，SCM確實會因為估計性能差引起自適應(yīng)匹配濾波檢測器性能損失，而MNSCM能夠改善AMF檢測器的檢測性能。仿真一組大小為40000×28相對均勻的雜波數(shù)據(jù)，即數(shù)據(jù)由28個距離單元40000個脈沖組成。選擇其中任意一個距離單元作為待測單元，并添加仿真目標(biāo)到該距離單元，其余的27個距離單元上的數(shù)據(jù)作為輔助數(shù)據(jù)。當(dāng)以N=8為長度作相干檢測時，27個距離單元中的數(shù)據(jù)被劃分為135000個短向量，即正常單元的數(shù)目。然后，隨機地將這27個距離單元中一些正常單元的功率調(diào)高，但它依然與距離鄰近的其他正常單元共享同一個歸一化協(xié)方差矩陣。由于這些被調(diào)整過的正常單元相比未調(diào)整過正常單元具有功率大的特點，這些單元可被看作是異常單元。通過控制異常單元的在原有正常單元中百分比以及異常單元相對正常單元的功率比值，構(gòu)造一系列不同程度的非均勻雜波環(huán)境，觀察兩種不同協(xié)方差矩陣估計方法的自適應(yīng)檢測器在這些環(huán)境中的檢測性能。

當(dāng)異常單元相對同一時間段上所有正常單元平均功率的比值固定時(簡稱功率比)，通過調(diào)整異常單元在所有單元中的比例，觀察由于異常單元的數(shù)目的增加引起的兩種檢測器之間的性能差異。設(shè)平均信雜比為0 dB，當(dāng)異常單元功率比為5，百分比(即異常單元在所有單元中的比例)依次為5%，15%和20%，圖5給了兩種檢測器的檢測結(jié)果；當(dāng)異常單元功率比為20，圖6給了兩種檢測器的檢測結(jié)果。

圖4 3個頻率處，MNSCM-AMF與SCM-AMF性能對比

圖5 功率比為5, 異常單元所占百分比為5%, 15%和20%時，MNSCM-AMF與SCM-AMF性能對比

圖6 功率比為20，異常單元所占百分比為5%, 15%和20%時, MNSCM-AMF與SCM-AMF性能對比

由圖5和圖6可知，當(dāng)異常單元數(shù)目較少時，MNSCM-AMF檢測器相比SCM-AMF檢測器的檢測性能優(yōu)勢并不明顯，二者具有非常相似的檢測結(jié)果；當(dāng)異常單元數(shù)目增大時，兩種檢測器的性能出現(xiàn)不同程度的下降，并且異常單元數(shù)目越多，兩種檢測器的性能差異越顯著。這是由于異常單元數(shù)目的增多不僅會導(dǎo)致兩種檢測器對應(yīng)的門限值增大，而且異常單元數(shù)目越多，采用SCM方法估計協(xié)方差矩陣產(chǎn)生的損失越大。通過對比圖5和圖6可知，相同異常單元數(shù)目條件下，功率比越大，一方面，兩種檢測器的性能都出現(xiàn)了不同程度的下降趨勢，這是由AMF檢測器的檢測機理決定的。由于AMF檢測器的檢測統(tǒng)計量與雜波的能量相聯(lián)系，異常單元功率的增大勢必提高了相應(yīng)的檢測門限。另一方面，兩種檢測器的性能差異越明顯，這是因為功率比的增大加劇了雜波的非均勻程度，使得SCM估計方法出現(xiàn)更大性能損失，從而影響了檢測器的檢測性能。

當(dāng)多普勒頻率為80 Hz，即目標(biāo)位于強雜波區(qū)，異常單元功率比分別設(shè)為5和20，異常單元所占的百分比依次為5%, 15%和20%，兩種檢測器的檢測結(jié)果如圖7和圖8所示。由圖7和圖8可知，功率比增大或異常單元數(shù)目增多時，兩種檢測器的檢測曲線右移，即檢測性能出現(xiàn)下降。這是因為：在相同數(shù)目的異常單元條件下，異常單元功率的增大直接提高了檢測門限，從而使得兩種檢測器的性能呈現(xiàn)下降趨勢；而相同功率比條件下，異常單元數(shù)目的增大，一方面會增大檢測門限，另一方面會引起SCM估計方法更多的估計損失，使得兩種檢測器的性能差異變得更顯著。對于其他多普頻率處，具有與上述類似的結(jié)果，限于篇幅，這里沒有給出詳細(xì)結(jié)果。

5 結(jié)束語

本文分析了在高分辨海雜波背景下的動目標(biāo)檢測的兩種自適應(yīng)匹配濾波檢測算法，即基于采樣協(xié)方差矩陣(SCM)和基于功率中值和歸一化采樣協(xié)方差矩陣(MNSCM)的自適應(yīng)匹配濾波(AMF)檢測器。實驗結(jié)果表明，在非均勻雜波環(huán)境中，由于異常單元的存在，SCM的估計性能受阻，限制了傳統(tǒng)AMF檢測器的檢測性能；而MNSCM能很好地抑制異常單元對估計的協(xié)方差矩陣的影響，提高了AMF檢測器在非均勻雜波環(huán)境中的檢測能力。

圖7 頻率為80 Hz，功率比為5，異常單元所占百分比為5%, 15%和20%時，MNSCM-AMF與SCM-AMF性能對比

圖8 頻率為80 Hz，功率比為20，異常單元所占百分比為5%, 15%和20%時，MNSCM-AMF與SCM-AMF性能對比

[1] Javier M C, Javier M G, Alvaro B D C, et al.. Small target detection in high-resolution heterogeneous sea-clutter: an empirical analysis[J]. IEEE Transactions on Aerospace and Electronic Systems, 2011, 47(3): 1880-1898.

[2] Chen Xiao-long, Cai Yong, and Cai Fu-qing. Application of the sparse decomposition to micro-motion target detection embedded in sea clutter[C]. Proceedings of the 2013 International Conference on Radar, Adelaide, Canada, 2013:163-166.

[3] Chen Si-jia, Kong Ling-jiang, and Yang Jian-yu. Adaptive detection in compound-Gaussian clutter with partially correlated texture[C]. Proceedings of the 2013 IEEE Radar Conference, Ottawa, Canada, 2013: 1-5.

[4] Chen Xiao-long, Guan Jian, He You, et al.. Detection and extraction of target with micro-motion in spiky sea clutter via short-time fractional fourier transform[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(2), 1002-1018. [5] Maria G, Fulvio G, and Muralidhar R. Statistical analysis of measured polarimetric clutter data at different resolutions [J]. IEE Proceedings-Radar, Sonar and Navigation, 2006, 153(6): 473-481.

[6] Maria G, Fulvio G, Muralidhar R, et al.. Impact of the sea clutter nonstationarity on disturbance covariance matrix estimation and CFAR detector performance[J]. IEEE Transactions on Aerospace and Electronic Systems, 2010, 46(6): 1502-1513.

[7] Frank R C, Daniel F R, Edward K J, et al.. A CFAR adaptive matched filter detector[J]. IEEE Transactions on Aerospace and Electronic Systems, 1992, 28(1): 208-216.

[8] Chan H C, Radar sea-clutter at low grazing angles[J]. IEE Proceedings F Radar and Signal Processing, 1990, 137(2): 102-112.

[9] Ward, K D, Baker, C J, and Watts S. Maritime surveillance radar. part I: radar scattering from the ocean surface[J]. IEE Proceedings, F Radar and Signal Processing, 1990, 137(2): 51-62.

[10] Fulvio G and James H M. Performance analysis of two covariance matrix estimators in compound-Gaussian clutter [J]. IEE Proceedings-Radar, Sonar and Navigation, 1999, 146(3): 133-140.

[11] 張玉石, 許心瑜, 尹雅磊, 等. L波段小擦地角海雜波幅度統(tǒng)計特性研究[J]. 電子與信息學(xué)報, 2014, 36(5): 1044-1048. Zhang Yu-shi, Xu Xin-yu, Yin Ya-lei, et al.. Research on amplitude statistics of L-band low grazing angle sea clutter [J]. Journal of Electronics & Information Technology, 2014, 36(5): 1044-1048.

[12] Wei Na. Sea clutter covariance estimation with data-adaptive selection[C]. Proceedings of the IET International Radar Conference, Xi'an, China, 2013: 1-5.

[13] Aghaabdellahian N and Hashemi M M. Improving performance of adaptive detectors in nonhomogeneous environment[C]. Proceedings of the 21st Iranian Conference on Electrical Engineering, Mashhad, Iran, 2013: 1-5.

[14] Moya J C, Menoyo J G, Lopez A A, et al.. Statistical analysis of a high-resolution sea-clutter database[J]. IEEE Transactions on Geoscience and Remote Sensing, 2010, 48(4): 2024-2037.

[15] Fulvio G and Maria G. Texture modeling, estimation, and validation using measured sea clutter data[J]. IEE Proceedings-Radar, Sonar and Navigation, 2002, 149(3): 115-124.

[16] Maria G, Federica B, and Fulvio G. X-band sea-clutter non-stationarity: influence of long waves[J]. IEEE Journal of Oceanic Engineering, 2004, 28(2): 269-283.

[17] Haykin S. Adaptive Radar Signal Processing[M]. Hoboken, NJ: John Wiley and Sons, Inc., 2007: 159-192.

[18] Cognitive System Laboratory, McMaster University, Canada. The IPIX radar database [OL]. http://soma.mcmaster.ca// ipix.php.2012.10.11.

劉 明： 男，1986年生，博士生，研究方向為海面微弱動目標(biāo)檢測.

水鵬朗： 男，1967年生，博士，教授，研究方向為多速率濾波器理論及應(yīng)用、圖像處理和雷達目標(biāo)檢測.

Adaptive Matched Filter Detector Based on Power Median and Normalized Sample Covariance Matrix

Liu Ming Shui Peng-lang
(National Key Laboratory of Radar Signal Processing, Xidian University, Xi'an 710071, China)

In nonhomogeneous sea clutter, abnormal cells included reference cells constrain the performance of the Sample Covariance Matrix (SCM), and then influence the detection performance of the traditional Adaptive Matched Filter (AMF) detector, while censoring abnormal cells may cause singularity of the covariance matrix in the case of limited reference cells. Without changing number of the reference cells, this paper devises the median and normalized covariance matrix estimator and uses in the detection scheme of the AMF. Compared with the traditional AMF, the newly devised AMF obtains better performance in both measured and simulated clutter.

Target detection; Sea clutter; Adaptive Matched Filter (AMF); Sample Covariance Matrix (SCM); Median and Normalized Sample Covariance Matrix (MNSCM)

TN957.51

： A

：1009-5896(2015)06-1395-07

10.11999/JEIT140900

2014-07-09收到，2014-10-08改回

國家自然科學(xué)基金(61271295)資助課題

*通信作者：劉明 mingliu1608@163.com