張 林，王 捷，黃 勇，關(guān) 鍵

（海軍航空工程學院 a.研究生管理大隊；b.電子信息工程系，山東 煙臺 264001）

0 引言

CFAR處理方法從總體上可粗略分為兩大類：參量和非參量的方法[1]。參量CFAR方法是假設(shè)雜波包絡(luò)的分布類型已知的條件下，通過估計一些未知參數(shù)，從而獲得CFAR的一類檢測算法。均值(ML)類和有序統(tǒng)計量(OS)CFAR方法就屬于這一類[2]。在雜波類型已知的場合，一般來說，參量方法的檢測性能優(yōu)于非參量方法[3]，但是由于實際雷達工作環(huán)境中，雜波類型往往無法準確已知并且是時變的，這時傳統(tǒng)的參量方法就失去了CFAR能力。因此，與雜波分布類型無關(guān)的非參量方法的研究受到越來越廣泛的關(guān)注[4-8]。

對非參量檢測器的研究始于上世紀五六十年代，比較經(jīng)典的如1945年Frank Wilcoxon提出了Wilcoxon檢測器[1]和1964年Carlyle J.W提出了符號檢測器[9]。但這兩類檢測器要求信號相位和雜噪聲的中值已知[1]的條件在實際雷達設(shè)計中往往難以實現(xiàn)，目前廣泛應用于實際雷達系統(tǒng)的是克服上述限制的兩樣本非參量檢測器。文獻[10-12]闡述了在瑞利雜波中兩樣本非參量檢測器的檢測性能，對韋伯爾分布(Weibull)、對數(shù)正態(tài)(log-normal)分布和K分布條件下兩樣本非參量檢測器的檢測性能可在文獻[13-14]中找到，另外在文獻[15]中對兩樣本非參量檢測器在實測海雜波數(shù)據(jù)的基礎(chǔ)上疊加起伏目標的檢測性能進行了研究。

本文在實測海雜波數(shù)據(jù)及真實目標環(huán)境中，利用蒙特卡洛方法對3種兩樣本非參量檢測器、廣義符號GS[10](Generalized sign)、MW[11](Mann-Whitney)、和S[12](Savage) CFAR檢測器的目標檢測性能和海雜波抑制能力進行了分析，并與參量CA-CFAR檢測器的檢測性能進行了比較。

1 非參量CFAR檢測器模型

圖1給出了兩樣本非參量檢測器的原理圖，在N個檢測周期中，檢測單元的采樣(t0時刻的采樣)用yj(j=1,2,…,M)表示，參考單元的采樣用xki(i=1,…,N,k=1,…,M)表示，其中檢測單元兩側(cè)分別有一個保護單元以防止檢測單元中可能的目標信號泄漏到鄰近的參考單元中影響檢測。下面給出前面提到的3種非參量檢測器及其檢驗統(tǒng)計量的數(shù)學表達式：

圖1 兩樣本非參量CFAR檢測器原理框圖

2 參量CFAR檢測器模型

參量檢測器的種類很多，本文主要介紹了CA-CFAR檢測器在實測海雜波數(shù)據(jù)中的檢測性能，其原理框圖如圖2所示。

圖2 CA-CFAR檢測器原理框圖

除了檢測單元采樣的表示方法外，其余按照大多數(shù)CFAR 文獻的表示方法：分別用xi(i=1,…,n)和yi(i=1,…,n)表示兩側(cè)參考單元(也稱作參考滑窗)采樣，參考單元為M=2n；n為前沿和后沿參考滑窗長度，X和Y分別是前沿和后沿滑窗中的局部估計。此時自適應判決準則為

式中：H1表示有目標；H0表示沒有目標；Z是干擾功率水平估計，代表參考滑窗中的平均雜波包絡(luò)估計；T是標稱化因子；D表示檢測單元中的檢測統(tǒng)計量。與檢測單元鄰近的是2個保護單元，主要用在單目標情況，防止目標能量泄漏到參考單元影響ML類檢測器兩側(cè)參考滑窗的局部估計值。對于CA-CFAR檢測器有

式中：

3 實測海雜波數(shù)據(jù)分析

本實測數(shù)據(jù)取自強弱海雜波背景下海上漁船目標實驗，該數(shù)據(jù)可以用來分析CFAR檢測器在強弱海雜波背景下對海上單一目標的分辨能力以及目標周圍海雜波的抑制能力。

圖3給出了強弱海雜波背景下原始數(shù)據(jù)回放結(jié)果，其中橫坐標為距離信息，縱坐標為脈沖信息，幅度值可由右側(cè)顏色條刻度讀出，圖4為CA-CFAR檢測器在虛警概率 Pfa=10?4的條件下對圖3數(shù)據(jù)的處理結(jié)果。

圖3 雷達原始視頻數(shù)據(jù)回放

圖4 CA-CFAR檢測器處理結(jié)果

從圖3原始視頻數(shù)據(jù)回放可以看出，在強海雜波背景下，由于目標周圍的海雜波較強，漁船目標湮沒于周圍的雜波背景之中，無法判定其位置。在弱海雜波背景下，漁船目標可以從原始視頻數(shù)據(jù)圖像中讀出，但是周圍存在一定的雜波干擾。圖4為CA-CFAR檢測器對該原始數(shù)據(jù)處理結(jié)果。從圖中可以看出，在強海雜波背景下，相比于原始視頻數(shù)據(jù)CA-CFAR檢測器對雜波有一定的抑制能力，但是抑制作用并不明顯，從圖像中仍無法檢測出漁船目標所在位置。在弱海雜波背景下，由于漁船目標相比于周圍雜波較強，CA-CFAR檢測器成功地檢測出了目標的位置，并對周圍的雜波有一定抑制能力。

分析原因：

1)實際雷達工作環(huán)境中的雜波包絡(luò)的分布類型往往無法準確已知，這就需要對雜波包絡(luò)的分布類型進行假設(shè)從而推算出門限因子，但是當假定的雜波包絡(luò)類型與實際的雜波包絡(luò)類型不匹配時，CA-CFAR檢測器就會喪失恒虛警能力，從圖4可以看出強海雜波背景下CA-CFAR檢測仍有較多的虛警點。

2)CA-CFAR檢測器的雜波功率水平是由周圍參考單元的采樣值形成，在弱海雜波背景下目標周圍的雜波較弱，雜波功率水平的估計值較低。這樣目標往往容易超過該值，從而被成功檢測，但是在強海雜波背景下目標周圍的雜波較強，雜波功率水平的估計值較高，目標往往無法超過門限，超過門限的個別點，又因為CA-CFAR的失配，周圍存在較多的雜波點而無法準確檢測。

圖5～7分別給出 GS-CFAR、MW-CFAR和Savege-CFAR檢測器在虛警概率 Pfa=10?4，參考單元個數(shù)M=16，非參量檢測器積累脈沖個數(shù)N=10的條件下對圖3所示原始視頻數(shù)據(jù)的處理結(jié)果。這里所取的假設(shè)條件與CA-CFAR檢測器相同，其中橫坐標表示距離信息，縱坐標表示脈沖信息，幅度值可由右側(cè)顏色條刻度讀出。

圖5 GS-CFAR檢測器處理結(jié)果

圖6 MW-CFAR檢測器處理結(jié)果

圖7 Savage-CFAR檢測器處理結(jié)果

從圖5～7可以看出：

1)3種非參量CFAR檢測器均可準確地檢測出2批漁船目標，其中強海雜波背景下的目標位于159號距離單元，弱海雜波背景下的目標位于76號距離單元。相比于雷達原始視頻數(shù)據(jù)圖像經(jīng)過非參量CFAR檢測器處理后的顯示效果改善明顯，顯示畫面較為“干凈”，特別是3種非參量檢測器均檢測出了在強海雜波背景下CA-CFAR檢測器并未檢測出的目標。另外從圖像可以看出3種非參量CFAR檢測器的虛警點較少，分析原因主要是非參量檢測器與雜波分布類型無關(guān)。因此，其在分布不確定或變化的背景中仍能保持恒虛警。

2)利用monte-carlo仿真對上述3種非參量檢測器的檢測概率進行計算發(fā)現(xiàn)：在弱海雜波背景下，GS-CFAR檢測器為68.5%，MW-CFAR檢測器為72.3%，Savege-CFAR檢測器為83.9%；在強海雜波背景下，GS-CFAR檢測器為19.3%，MW-CFAR檢測器為17.4%，Savege-CFAR檢測器為19.0%。從上述數(shù)據(jù)可以發(fā)現(xiàn)，在強海雜波條件下，GS-CFAR檢測器的檢測概率最高，在弱海雜波條件下，Savage-CFAR檢測器的檢測概率最高。分析原因主要是因為在強海雜波條件下，目標周圍的海雜波較強，目標的幅值與周圍雜波的幅值較為接近，在利用多個脈沖回波信息形成統(tǒng)計量的MW-CFAR和Savege-CFAR檢測器由于目標的優(yōu)勢不明顯，往往使產(chǎn)生的檢驗統(tǒng)計量低于門限值。因此，效果不如GS-CFAR檢測器，但是在弱海雜波條件下，目標的幅值相對于周圍的雜波的幅值有明顯優(yōu)勢，在利用多個脈沖信息形成的檢驗統(tǒng)計量往往比利用單個脈沖回波信息形成的檢驗統(tǒng)計量的優(yōu)勢積累更明顯，因此MW-CFAR、Savage-CFAR檢測器的檢測概率更高。

4 總結(jié)

本文基于實測雷達海雜波數(shù)據(jù)分別以 GSCFAR、MW-CFAR和Savege-CFAR檢測器為例研究了兩樣本非參量CFAR檢測算法在實際雷達中的目標檢測性能，并與參量CA-CFAR檢測器進行了對比。研究表明，相比于參量CA-CFAR檢測器，3種非參量檢測器在強弱海雜波背景下均有更好的雜波抑制能力和不錯的檢測性能，特別是在強海雜波背景下，3種非參量檢測器成功地檢測出了 CACFAR檢測器并未檢測出的位于159號距離單元的漁船目標。在強海雜波條件下，GS-CFAR檢測器的檢測性能最優(yōu)，在弱海雜波條件下，Savage-CFAR檢測器的檢測性能最優(yōu)，但是綜合考慮，在實際工作中更關(guān)心強海雜波條件下目標檢測能力。另外GS-CFAR檢測器的復雜度相對于Savage-CFAR檢測器更低，因此建議使用GS-CFAR檢測器。

[1]何友,關(guān)鍵.雷達自動檢測與恒虛警處理[M].北京:清華大學出版社,1999:230.

[2]祝本玉,畢大平,王正.兩類典型的CFAR檢測器性能仿真[J].艦船電子對抗,2008,31(2):64-68.

[3]孟祥偉.韋布爾雜波下非參數(shù)量化秩檢測器的性能[J].電子學報,2009,37(9):2030-2034.

[4]SEYFE B,SHARAFAT A R.Signed-rank nonparametric multiuser detection in non-Gaussian channels[J].IEEE transactions on information theory,2005,51(4):1478-1486.

[5]HAO CHEN.Noise enhanced nonparametric detection[J].IEEE transactions on information theory,2009,55(2):499-506.

[6]BROWN C L.A nonparametric approach to signal detection in impulsive interference[J].IEEE transactions on signal processing,2000,48(9):2665-2669.

[7]ORTON M.Nonparametric rank detectors under k-distributed clutter in radar application[J].IEEE transactions on aerospace and electronic systems,2005,41(2):702-710.

[8]SEYFE B,SHARAFAT A R.Nonparametric multiuser detection in non-Gaussian channels[J].IEEE transactions on signal processing,2006,54(1):23-33.

[9]CARLYLE J W.On nonparametric signal detectors[J].IEEE Trans.On IT,1964,10(2):146-152.

[10]HANSEN V G.Nonparametric radar extraction using a generalized sign test[J].IEEE Trans.On AES,1971,7(5):942-950.

[11]ZEOLI G W.Performance of a two-sample mannwhitney nonparametric detector in a radar application[J].IEEE Trans.on AES,1971,7(5):951- 959.

[12]SAVAGE I R.Contributions to the theory of rank order statistics-the two-sample case[J].Ann Math Statistics,1956,27(2):590-615.

[13]吳靖巍,王小謨,曹晨.兩種非參量檢測器在非瑞利雜波中的檢測性能[J].計算機仿真,2007,24(11):5-16.

[14]SANZ-GONZáLEZ J L,ALVAREZ-VAQUERO F.Nonparametric rank detectors under K distributed clutter in radar applications[J].IEEE Trans on AES,2005,41(2):702-710.

[15]趙志堅.海雜波中非參量恒虛警檢測器性能分析[J].雷達科學與技術(shù),2010,8(1):65-73.