摘  要： 針對UHF波段運動目標檢測中信號弱、雜波展寬嚴重等問題，提出一種基于子孔徑SAR技術的雜波抑制方法。首先通過SAR成像實現運動目標的粗聚焦，再通過距離-方位拓展局部STAP實現雜波抑制，避免了運動目標越距離單元和方位散焦帶來的信號損失。最后，通過仿真結果驗證了算法的有效性。

　　關鍵詞： 地面運動目標指示；空時自適應處理；子孔經；雜波抑制；

0 引言

　　具有同時獲得靜止目標成像和運動目標指示的多通道SAR/GMTI系統受到了各軍事強國的青睞[1-5]。相對于高波段（C/X/Ku等）SAR/GMTI系統，低頻段（UHF）系統具有較強的穿透能力，可穿透植被探測目標，是全透明偵察的重要組成部分。然而，葉簇等植被的穿透衰減使得隱蔽目標的回波強度非常弱，同時受天線尺寸的限制波束角非常寬，從而使得回波中雜波強度大，且雜波的多普勒譜展寬嚴重，直接采用常規的運動目標檢測方法難以獲得足夠的信雜噪比改善。

　　長相干STAP方法[1-4]將傳統的空間—方位慢時間二維STAP在距離快時間維進行擴展，具有處理越距離單元徙動的能力，從而利用更多的脈沖進行相干處理，提高輸入信號質量，實現慢速運動目標的檢測。但直接應用的計算量和協方差矩陣估計所需的樣本數都難以滿足研究要求，需要進行降維處理。

　　本文提出一種基于子孔徑SAR的長相干STAP算法，該算法首先通過子孔徑SAR成像實現運動目標的粗聚焦，將在快時間-慢時間維分散的運動目標信號凝聚，提高運動目標的信號質量，再對粗聚焦的像素進行方位和時間聯合，進而利用多通道空間信息實現UHF波段的運動目標檢測。

1運動目標子孔徑成像

　　根據SAR成像原理，當成像的參數選擇與運動目標參數匹配時可以實現運動目標的聚焦成像[4]。然而，在實際的應用中同時存在多個不同速度的運動目標，不可能同時獲得各個速度目標的完美聚焦。另一方面，在UHF波段的SAR/GMTI中，由于受到植被等遮蔽體的電磁波衰減影響，目標信號衰減嚴重，需要進行較長時間的相干積累來獲得足夠的信號質量。因此，在UHF波段的運動目標檢測中，子孔徑長度一般取的較長。此時，目標運動引入的運動誤差將積累，直接采用SAR成像方法將出現一定的散焦。以表1所列的系統參數進行仿真，對不同速度的運動目標進行子孔徑（2048脈沖）SAR成像，成像結果如圖1所示。

　　由子孔徑成像結果可見，由于子孔徑時間較短，低速的運動目標可以實現較好的聚焦成像；速度較高的運動目標出現了越距離單元的距離校正殘余量，但遠距離單元的徙動量較小，主要集中在鄰近的距離單元。另一方面，運動目標的方位向運動速度分量，使得運動目標在方位向也出現一定的散焦，但同樣地主要集中在鄰近的幾個方位單元內。

　　通過以上分析，可見當運動目標的速度較低時可以實現運動目標的子孔徑聚焦，此時可通過對同一位置的單一分辨單元建立多通道觀測模型，并通過空時處理實現雜波的抑制；當運動目標的速度較高時，運動目標的能量主要集中在相鄰的幾個分辨單元內，對單一分辨單元建立多通道觀測模型，不能包含全部的運動目標信息。因此，本文將單一分辨單元模型在方位向和距離向進行擴展，建立聯合分辨單元的STAP模型。

2 局部STAP建模

　　局部STAP的信號樣本示意圖如圖2所示。

　　將N個通道、M個方位單元、L個距離單元的多通道子孔徑SAR成像數據排列成N×M×L維的列向量：

　　Z（i，j）=[ZT（i1，j1）…ZT（im，j1）…ZT（iM，jl），…ZT（i1，jl）…ZT（im，jl）…ZT（iM，jl），…ZT（i1，jL）…ZT（im，jL）…ZT（iM，jL）]T（1）

　　其中，ZT（im，jl）代表成像位置（im，jl）處的N維列矢量：

　　ZT（im，il）=[z1（im，il），z2（im，il），…，zN（im，il）]T（2）

　　對觀測矢量Z（i，j）建立假設檢驗：

　　Z（i，j）=bS（i，j）+C（i，j）+Q（i，j）  H1假設C（i，j）+Q（i，j）          H0假設（3）

　　其中，b為運動目標的幅度，C為雜波信號，Q為噪聲信號，S為運動目標信號的導向矢量：

　　S（i，j）=[ST（i1，j1），ST（im，jl），…，ST（im，jL）]T（4）

　　那么，滿足輸出信雜噪比最高準則的最優權系數Wopt（i，j）可表示為：

　　其中代表H0假設下觀測樣本的協方差矩陣，即觀測樣本中雜波和噪聲的協方差矩陣，為恒虛警常數。

　　那么，雜波抑制后的輸出可表示為：

　　將輸出檢測統計量的絕對值與檢測門限進行比較，即可判定運動目標是否存在。

3 仿真驗證

　　基于系統的天線參數計算得到天線的收發方向圖，基于系統的PRF和載機的速度得到每個脈沖的天線方位位置，同時根據多通道的間隔得到多個通道多脈沖的等效相位位置。雜波場景的仿真通過高斯分布產生，同時通過距離和方位的間隔距離產生場景的散射點分布。由雜波得到原始回波信號后，添加高斯復噪聲，實現雜噪比因素的仿真。

　　由于在UHF系統中多為遮蔽目標的檢測，此時信號受到嚴重的衰減，信號處理的輸入信雜噪比非常低，仿真中通過衰減目標回波信號實現。仿真中采用50 dB的雜噪比，信雜噪比為-57 dB，脈沖點數為       2 048。運動目標位于天線方向圖的中心，速度分布于-9 m/s~ +9 m/s之間，以3 m/s為間隔。原始回波的距離多普勒如圖3所示。對3通道回波檢測局部STAP模型，并進行雜波抑制后的結果如圖4所示。從圖4中可見到雜波被有效地抑制，可看到清晰的運動目標的分布。

4 結論

　　本文在分析了運動目標子孔徑成像的基礎上，提出了一種適用于UHF波段多通道SAR/GMTI系統慢速運動目標檢測的方法，該方法通過子孔徑SAR實現運動目標的粗聚焦，降低運動目標的分布范圍，再通過方位-距離維的局部STAP實現雜波的抑制，進而實現UHF波段下運動目標的檢測。通過仿真實驗初步驗證了算法的有效性，下一步將進行實測數據的處理和算法的進一步優化。

參考文獻

　　[1] STEVEN E G. Upper bounds on processing loss for wideband long-CPI space-time adaptive processing[C]. IEEE Pacific Rim Conference on Communication， Computers and Signal Processing， 2009：251-256.

　　[2] 常玉林，周紅，黃曉濤，等.多通道SAR頻率多普勒域寬帶長CPI STAP方法[J].電子學報，2011，39（6）：1245-1252.

　　[3] 周紅，黃曉濤，常玉林，等.子帶子孔徑ATI地面運動目標檢測及參數估計方法[J].電子與信息學報，2010，32（1）：62-68.

　　[4] Zhu Shengqi， Liao Guisheng. Ground moving targets imaging algorithm for synthetic aperture radar[J]. IEEE Transaction on Geoscience and Remote Sensing， 2011， 49（1）： 462-477.

　　[5] 殷戴乾，田錚.多區域SAR圖像分割的改進水平集方法[J].微型機與應用，2013，32（20）：38-41.