基于雜波分布特性的恒虛警率檢測器設計

洪 偉,馬艷艷,程 義

(中國船舶集團有限公司第七二三研究所,江蘇 揚州 225101)

0 引 言

恒虛警率(CFAR)檢測器的主要功能是在包含噪聲、雜波甚至干擾的雷達回波數據中,檢測目標是否存在,其在雷達信號檢測技術中具有重要的地位。Nitzberg將CFAR閾值分為以下4種:固定閾值;基于干擾平均幅度的閾值;基于干擾統計先驗信息的閾值;無干擾統計先驗信息時,基于分布自由的統計假設檢驗的閾值[1]。對固定閾值CFAR檢測器的研究已經非常成熟,它主要應用在高斯白噪聲確知恒定的條件下。然而,實際工程中雷達接收回波中包含了多種類型的環境雜波、海雜波、地雜波、氣象雜波等,這些環境因素導致了雷達目標檢測的背景一般處于非均勻環境,難以達到平穩恒定的要求,因此,自適應閾值CFAR檢測器的研究應運而生。

在待檢測單元所處的背景信息未知時,有關自適應閾值恒虛警率檢測器的研究越來越成熟,文獻[2]研究內容為在背景環境未知時,將均方差和有序差分思想應用到背景功率估計中,提出了一種基于最大類方差的恒虛警率檢測器,基于仿真數據的實驗表明,所提恒虛警率檢測器在多目標環境下有更強的抗干擾能力。文獻[3]針對海雜波情況下的小目標檢測,提出了多幀聯合的檢測思想,分兩步完成:第1步通過設置較低的門限,保證較高的目標檢測概率;第2步利用海雜波的尖峰和目標在時間相關性方面的差異,抑制海雜波虛警。文獻[4]分析了常用的平均類恒虛警檢測器和排序類恒虛警檢測器,并分別針對2類恒虛警檢測器提出了改進算法,仿真結果驗證了改進恒虛警檢測器的有效性。文獻[5]針對檢測前跟蹤(TBD)算法計算量大的問題,改進了二維CFAR的參考滑窗方法,基于仿真實驗驗證了背景噪聲估計的有效性,實驗中在不同背景功率下自適應改變門限的取值,提高了目標檢測的實時性。

有不少文獻研究了在待檢測單元所處的背景信息已知的情況下恒虛警檢測器的設計。文獻[6]針對K分布背景下的海雜波,結合非相參積累檢測方法對脈沖積累數、海雜波模型參數等的依賴特性,提出了基于海雜波白化方法的恒虛警檢測器。文獻[7]通過仿真產生了幅度為瑞利分布、對數正態分布和韋布爾分布3種類型的雜波,針對不同雜波類型,提出了不同的恒虛警檢測器設計。文獻[8]為了解決傳統恒虛警檢測器在非均勻韋布爾雜波環境中性能劇烈下降的問題,提出了一種基于雜波邊緣二元積累算法的恒虛警檢測器,并使用對比實驗驗證了新提出恒虛警檢測器的有效性。文獻[9]通過計算待檢測單元所處雜波的背景功率,判斷待檢測單元位于強雜波區域或弱雜波區域,自適應選擇不同的的恒虛警檢測器。文獻[10]通過計算待檢測單元所處雜波的背景功率,判斷待檢測單元位于強雜波區域或弱雜波區域,自適應選擇不同的恒虛警檢測器。文獻[11]在已知背景分布為對數-正態分布時,設計了一種有序統計恒虛警檢測器,并用數值方法討論了參考單元個數導致的信雜比損失以及最佳序值的選取。

本文基于當前環境確知的前提下,根據檢測單元所處的雜波類型為陸地雜波、海雜波或陸雜波海雜波交界這3種類型,選擇不同的恒虛警檢測器,并通過人工加運動目標的方式,使用實測數據對所提方法進行仿真驗證。

1 雜波特性分析

1.1 雜波時域回波起伏特性

由圖1可知,海雜波背景幅度起伏最大,海陸交叉區域雜波背景幅度起伏較大,陸地雜波背景幅度較為平穩,而且相對于海雜波和海陸交叉區域雜波來說,陸雜波幅度較小;海雜波的幅度分布高振幅回波出現的概率大,時域波形較尖銳,海陸交叉區域雜波次之,陸雜波最平穩。

圖1 3種雜波類別的時域特性

1.2 雜波頻域特性

圖2為對相參處理間隔的數據做16點快速傅里葉變換(FFT)之后的距離-頻譜圖,x軸表示一個脈沖的有效距離單元,y軸為頻譜維,z軸為幅度,(a)所示為海雜波的距離-頻譜圖,(b)表示海陸交叉區域雜波的距離-頻譜圖,(c)表示陸雜波的距離-頻譜圖。由圖2可知,海浪可以在雷達徑向產生相對速度,故海雜波的頻譜中心頻率一般在零頻以外,容易構成虛假目標,而大多數地雜波不能在雷達徑向產生相對速度,其中心頻率一般在零頻,因此海雜波的頻譜幅度與陸雜波相比有明顯的偏移和較長的拖尾。由此可得,海雜波分布范圍較大,在多個距離單元和多個頻率上均有較強的幅度;海陸交叉區域的雜波在較少的距離單元和頻率上有分布;而陸雜波一般在某些距離單元有分布,在更多的距離上沒有分布,且陸雜波幅度相對于其他2種類別雜波的幅度小。

圖2 3種類別雜波的距離-頻譜圖

2 恒虛警檢測器

由第1節分析可知,海雜波背景下,雜波幅度起伏較大,頻譜帶寬較大;海陸交叉區域次之;陸雜波背景幅度起伏平穩,頻譜帶寬較小。故在海雜波區域選擇基于威布爾的恒虛警率檢測器;在海陸交叉區域選擇基于對數-正態分布的恒虛警率檢測器,在陸雜波區域選擇有序恒虛警率檢測器。

2.1 有序恒虛警率檢測器[12]

圖3 有序恒虛警檢測器流程圖

(1)

m=3/4×2L

(2)

與傳統的單元平均恒虛警率檢測器相比,當參考單元中混有干擾時,有序恒虛警檢測器只會引起參考單元排序結果變化,對檢測門限影響相對較少,故其更適應于多目標環境檢測中。

2.2 基于對數-正態分布的恒虛警率檢測器[12]

圖4 對數-正態分布的恒虛警率檢測器流程圖

yi=lgxi

(3)

(4)

(5)

(6)

然后選擇單元平均恒虛警率檢測器對參考單元內雜波電平進行估計,獲得虛警門限。

2.3 基于威布爾分布的恒虛警率檢測器[12]

文獻[12]通過公式推導證明了將威布爾分布的雜波轉換為正態分布后,虛警率只與檢測門限有關,與威布爾雜波的形狀參數和尺度參數均無關,故若海雜波滿足威布爾分布,即可使用該檢測器。首先按照公式(3)～(6)對威布爾分布的雜波進行轉換,將分布轉化為標準正態分布,然后選擇有序恒虛警率檢測器選出參考單元內雜波電平的估計,流程圖如圖5所示。

圖5 威布爾分布的恒虛警率檢測器流程圖

3 實驗結果

本實驗基于實測數據,錄取環境有海和陸地,雜波背景分為陸雜波環境、海陸交叉環境2類。根據不同環境的雜波類別,選擇不同的恒虛警率檢測器,對目標進行檢測。對比的恒虛警率檢測器選擇工程上常用的單元平均恒虛警率檢測器。

通過人工加目標的方式,在雜波較強的多個方位加入1個信噪比為-10 dB的目標,非相參積累點數選擇16,待檢測單元的左右保護單元選擇1,左右兩側的參考單元數目選擇16,比較本文所提恒虛警率檢測器和傳統單元平均恒虛警率檢測器的檢測性能。圖6(a)～(c)所示為上述方位的目標檢測情況,“+”線表示非相參積累后的待輸入恒虛警率檢測器的數據,“×”線表示傳統的單恒虛警率門限,虛線表示本文所提的恒虛警率檢測器的門限,實線為目標。由圖6對比可知,對于雜波較強區域附近的目標,選擇單元平均恒虛警率檢測器時,門限受附近幅度較大的雜波影響較明顯,導致門限較高,往往發生漏檢的情況,目標都未檢測到;而選用本文所提恒虛警率檢測器時,門限受附近幅度較大的雜波影響很小,可以將2個目標檢測出來。

圖6 在雜波較強區域檢測性能對比圖

在雜波相對較平穩的海陸交叉區域和陸雜波區域的多個方位,在距離單元間隔3的位置人工加入2個信噪比均為-10 dB的目標,非相參積累點數選擇16,待檢測單元的左右保護單元選擇1,左右兩側的參考單元數目選擇16,比較本文所提恒虛警率檢測器和傳統單元平均恒虛警率檢測器的檢測性能。圖7(a)～(f)所示為上述方位的目標檢測情況,“+”線表示非相參積累后的待輸入恒虛警率檢測器的數據,“×”線表示傳統的單恒虛警率門限,虛線表示本文所提的恒虛警率檢測器的門限,實線表示目標。圖7(a)、(b)和(d)對比可知,在雜波相對較強區域,選擇單元平均恒虛警率檢測器時,門限受附近幅度較大的雜波影響較大,可能導致漏檢,只能檢測到1個目標或2個目標都未檢測到;而選用本文所提恒虛警率檢測器時,可以同時將2個目標檢測出來。在雜波較平穩區域,本文所提恒虛警率檢測器和單元平均恒虛警率檢測器都可以將2個目標檢測出來,如圖7(c)、(e)、(f)所示。由以上分析可知,本文所提恒虛警率檢測器比傳統的單元平均恒虛警率檢測器檢測性能好。

圖7 在雜波較為平穩區域檢測性能對比圖

4 結束語

CA-CFAR在工程上得到成熟的應用,但當雜波環境不滿足平穩特性的時候,雜波模型已經不服從高斯分布,若依然選擇CA-CFAR,受附近單元雜波幅度的影響,恒虛警率門限會較高,將出現漏警現象,檢測性能大幅下降。因此,本文提出一種基于雜波分布特性的自適應恒虛警率檢測器設計方法。主要目的是當雜波環境不服從高斯分布時,選擇其對應分布模型的恒虛警率檢測器,使恒虛警率得到保障的基礎上,提高目標的檢測性能?；趯崪y數據,通過在不同雜波類別的區域人工加運動目標的方式進行仿真驗證,經過與CA-CFAR檢測結果的對比,驗證了本文所提恒虛警率檢測器的有效性。