同亞龍 王 彤文 才 吳建新

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

一種穩(wěn)健的機載非正側(cè)視陣?yán)走_(dá)雜波抑制方法

同亞龍 王 彤*文 才 吳建新

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

非正側(cè)視陣列構(gòu)型會導(dǎo)致機載雷達(dá)回波數(shù)據(jù)非均勻，從而嚴(yán)重削弱傳統(tǒng)空時自適應(yīng)處理算法的雜波抑制性能。針對此問題，該文提出一種對陣元誤差穩(wěn)健的非均勻雜波抑制方法。該方法首先根據(jù)雷達(dá)系統(tǒng)參數(shù)等先驗信息構(gòu)造雜波表示基；然后在考慮陣元誤差下對待檢測單元數(shù)據(jù)進(jìn)行迭代的最小二乘擬合，文中推導(dǎo)出了該問題的閉式解；最后對擬合后的剩余數(shù)據(jù)進(jìn)行脈沖多普勒處理以及恒虛警檢測。該方法不需要訓(xùn)練樣本，且能夠在沒有俯仰自由度的情況下對機載非正側(cè)視陣?yán)走_(dá)的非均勻雜波進(jìn)行有效抑制。仿真結(jié)果驗證了該方法的有效性。

機載雷達(dá)；非正側(cè)視陣；非均勻雜波；雜波抑制；空時自適應(yīng)處理

1 引言

機載預(yù)警雷達(dá)在接收目標(biāo)回波信號的同時，會不可避免地受到復(fù)雜強地雜波的影響。如何從極強的雜波背景中對運動目標(biāo)進(jìn)行有效檢測始終是機載預(yù)警雷達(dá)信號處理的核心任務(wù)，因此雜波成為其必須解決的首要問題?？諘r自適應(yīng)處理（Space-Time Adap tive Processing， STAP）技術(shù)能夠從空間和時間兩維上同時區(qū)分運動目標(biāo)和雜波，可以在很大程度上改善機載預(yù)警雷達(dá)的雜波抑制及動目標(biāo)檢測性能，特別是對低速目標(biāo)以及被旁瓣雜波所遮蔽的弱小目標(biāo)的檢測更為有效，因而得到了廣泛的關(guān)注［13］-。然而雜波信號統(tǒng)計特性往往是無法預(yù)知的，STAP需要選取待檢測單元（Cell Under Test， CUT）附近的回波數(shù)據(jù)作為訓(xùn)練樣本對其進(jìn)行估計。為獲得較優(yōu)的STAP性能（相對于最優(yōu)處理的性能損失不超過3 dB），所需要的訓(xùn)練樣本數(shù)不應(yīng)少于系統(tǒng)自由度的兩倍［4］。這里假設(shè)所選用的訓(xùn)練樣本與待檢測單元數(shù)據(jù)滿足獨立同分布（Independent and Identically D istributed， IID）條件，即這些訓(xùn)練樣本是均勻的且能夠反映CUT數(shù)據(jù)的統(tǒng)計特性。然而，由于實際場景中存在多變的地表類型、植被覆蓋以及孤立干擾等情況，機載雷達(dá)通常都工作在非均勻雜波環(huán)境中［5］；另外，當(dāng)雷達(dá)天線非正側(cè)視安置時，近程回波會具有非常明顯的距離依賴性［68］-。這些因素都會導(dǎo)致機載雷達(dá)回波數(shù)據(jù)非均勻，使得上述假設(shè)所要求的IID樣本條件很難滿足，STAP的雜波抑制性能受到很大程度的削弱。

對于非均勻環(huán)境造成的STAP性能下降問題，目前通常采用的解決辦法是在自適應(yīng)信號處理前級聯(lián)有效的非均勻檢測器（Non-Homogeneity Detector， NHD），對數(shù)據(jù)均勻性做出判斷，從而剔除有損自適應(yīng)算法性能的非均勻樣本數(shù)據(jù)［912］-。這類算法能夠在一定程度上提高非均勻環(huán)境下的檢測性能，但當(dāng)環(huán)境非均勻性嚴(yán)重時，NHD需要剔除大量奇異樣本，使得可用來進(jìn)行自適應(yīng)權(quán)值訓(xùn)練的樣本數(shù)變少甚至不足，STAP算法性能得不到有效提高，另外，這類算法大都關(guān)注的是如何更好地去除非均勻的訓(xùn)練樣本，而無法對CUT內(nèi)的孤立雜波進(jìn)行有效抑制。針對雷達(dá)非正側(cè)視安置所帶來的非均勻雜波抑制問題的研究也較多，比較典型的解決方法有多普勒補償［13，14］、角度多普勒補償［15，16］、配準(zhǔn)補償［17，18］等算法，但當(dāng)雷達(dá)發(fā)射中高脈沖重復(fù)頻率信號時，會發(fā)生距離模糊現(xiàn)象，這些補償類算法不能對多次模糊的回波數(shù)據(jù)進(jìn)行一致補償，即在對近程雜波的距離依賴性進(jìn)行補償?shù)耐瑫r會造成模糊的遠(yuǎn)程平穩(wěn)雜波出現(xiàn)新的非均勻性。除此之外，3維STAP方法［19］利用接收系統(tǒng)具有的俯仰維自由度，理論上也可以對非正側(cè)視陣?yán)走_(dá)的近程非均勻雜波進(jìn)行有效抑制，但由于系統(tǒng)自由度的增加，其運算量顯著增大，且對IID樣本數(shù)據(jù)的需求也更大，這在實際應(yīng)用中很難滿足。另外，帶有俯仰預(yù)濾波處理的子陣合成算法［20］在不存在陣元誤差影響下可以有效抑制近程非均勻雜波，但陣元誤差在實際情況中是不可避免的，算法性能因此受到很大程度地削弱。

針對以上問題，本文提出了一種對陣元誤差穩(wěn)健的機載非正側(cè)視陣?yán)走_(dá)非均勻雜波抑制方法。對于均勻線性陣列天線或者列合成后的線性子陣系統(tǒng)，由于缺乏俯仰自由度，不能直接對近程非均勻雜波進(jìn)行有效抑制。在這種情況下，本文方法首先根據(jù)機載雷達(dá)系統(tǒng)參數(shù)以及平臺運動狀態(tài)等先驗信息構(gòu)造雜波表示基，然后在考慮陣元誤差影響下對CUT數(shù)據(jù)進(jìn)行迭代的最小二乘擬合，最后對擬合后的剩余數(shù)據(jù)進(jìn)行常規(guī)的脈沖多普勒處理以及單元平均恒虛警率（Cell-Averaging Constant False A larm Rate， CA-CFAR）檢測。文中推導(dǎo)出了所涉及到的最小二乘問題的閉式解，大大減少了所提算法的運算復(fù)雜度。仿真實驗驗證了本文方法的有效性。

2 信號模型與雜波特性

機載雷達(dá)的幾何構(gòu)型如圖1所示。假定雷達(dá)天線為非正側(cè)視均勻線陣，載機以速度v沿平行于Y軸正方向的水平方向飛行，高度為h；雷達(dá)陣面由間距為d的N個陣元組成，陣面偏置角（陣面軸向與速度矢量的夾角）為α；雷達(dá)工作波長為λ，脈沖重復(fù)頻率為 fr，一個相干處理時間內(nèi)的脈沖數(shù)為M。

圖1 機載雷達(dá)陣面與雜波幾何關(guān)系

若雜波散射體相對于天線陣面的方位角和俯仰角分別為θ和φ，ψ是其對應(yīng)的空間錐角，滿足cos（ψ）= cos（θ） cos （φ），則該散射體回波的歸一化多普勒頻率為

由Ward模型［1］可知，第l個采樣距離單元所接收到的雜波回波信號是由各次模糊距離門內(nèi)不同方位角處散射體回波的線性疊加，因而可以表示為

由于同一距離環(huán)中不同雜波散射體之間是互不相關(guān)的，因此第l個采樣距離單元的回波數(shù)據(jù)所對應(yīng)的協(xié)方差矩陣可表示為

為了防止光束在球面鏡的外表面與內(nèi)表面之間干涉，球面反射鏡的外表面鍍99%的增透膜.若球面反射鏡和微懸臂的反射率分別為R1，R2，球面鏡與微懸臂的距離為L，即半球面F-P腔腔長，則相鄰光束之間的光程差為4L，位相差為

當(dāng)陣列天線不存在誤差時，天線上的每個陣元具有一致的幅相特性，此時由式（2）-式（5）所給出的信號模型可以很好地表征機載雷達(dá)回波數(shù)據(jù)中的雜波信號。但在實際陣列天線中，受機械加工等因素的影響，不可避免地存在著各種誤差，這也就意味著式（4）所采用的理想陣列流型已不再合適，此時需要考慮誤差對陣列天線的影響。本文主要考慮陣元誤差這一典型情況［21，22］，其誤差模型可以表示為

其中diag（·）表示對角化操作， IM表示M維的單位陣。這樣，在存在陣元誤差影響時，雜波的真實空時導(dǎo)向矢量就需要用式（8）的錐削矩陣對式（3）進(jìn)行修正，即

3 非均勻雜波抑制方法

由上述雜波特性分析可以看出，非正側(cè)視陣?yán)走_(dá)系統(tǒng)的雜波2維譜隨距離的變化而變化。傳統(tǒng)STAP方法需要選取CUT附近一段距離范圍內(nèi)的樣本數(shù)據(jù)進(jìn)行自適應(yīng)權(quán)值的訓(xùn)練，而非正側(cè)視陣雜波的這種距離依賴性會導(dǎo)致STAP雜波抑制凹口變寬，降低算法對低速運動目標(biāo)的檢測性能。另外，由于無法估計CUT內(nèi)孤立雜波的統(tǒng)計特性，STAP不能對其進(jìn)行充分抑制從而造成大量虛警。本節(jié)給出一種對陣元誤差穩(wěn)健的機載非正側(cè)視陣?yán)走_(dá)非均勻雜波抑制方法，主要包括雜波表示基的構(gòu)造、雜波與誤差的迭代估計兩部分。

3.1 構(gòu)造雜波表示基

式（6）所給出的雜波協(xié)方差矩陣可以寫成

其中svd（·）表示奇異值分解操作，U和V分別表示Sl的左、右奇異向量矩陣，Σ為奇異值矩陣。

其中η∈［0，1］，越接近于1逼近質(zhì)量越高。這樣可以通過設(shè)置有效秩g為滿足η≥η0的最小整數(shù)來保證雜波子空間的逼近效果，其中η0為接近于1的閾值，例如η0= 0.999等。

在考慮存在陣元誤差影響時，將式（9）中修正的雜波空時導(dǎo)向矢量代入到上述計算可得新的雜波表示基矩陣為

對于雜波表示基的構(gòu)造，這里有兩點需要進(jìn)行說明。一是距離模糊數(shù) Nr的設(shè)定，由于非正側(cè)視陣的遠(yuǎn)程雜波逐漸趨于平穩(wěn)，高次模糊距離門雜波分量的角度多普勒譜幾乎重合，其相應(yīng)的空時信號流型近似一致；另一方面，遠(yuǎn)程雜波信號的回波功率隨距離的衰減非常嚴(yán)重。因此lS表達(dá)式中rN的設(shè)定可根據(jù)實際需要進(jìn)行調(diào)整，這樣就大大降低了高維矩陣奇異值分解所帶來的運算復(fù)雜度。二是本部分雜波表示基的構(gòu)造可以離線完成，即在已知陣列響應(yīng)及平臺運動參數(shù)的情況下，可在實時數(shù)據(jù)處理之前計算并存儲算法所要用到的雜波表示基，有助于算法具有更好的實時性。

3.2 迭代估計雜波與誤差

用式（15）所獲得的雜波表示基矩陣?lΨ對第l個采樣單元數(shù)據(jù)lx進(jìn)行最小二乘擬合，可通過求解式（16）的優(yōu)化問題完成

當(dāng)陣元誤差已知時，式（16）可轉(zhuǎn)化為

當(dāng)擬合系數(shù)已知時，式（16）的目標(biāo)函數(shù)可變形為

也就是說，原N維的優(yōu)化問題可以分解為N個1維的優(yōu)化問題，其中每一個子優(yōu)化問題都是一個超定的最小二乘問題，存在唯一的最小均方誤差解，即

本文算法的詳細(xì)處理流程在圖2中給出。另外需要說明的是，本文采用的是全維處理方式，當(dāng)雷達(dá)系統(tǒng)空時自由度過大時，本文算法也可很容易推廣到空域、時域或空時域降維處理方式。

4 仿真實驗

本文按照表1所給出的雷達(dá)系統(tǒng)參數(shù)進(jìn)行仿真實驗。為了驗證本文算法對非正側(cè)視陣雜波的抑制效果，實驗以前視陣情況（即α=-90o）為例進(jìn)行說明。雜波數(shù)據(jù)按照W ard模型仿真產(chǎn)生，其能量受方向圖及距離聯(lián)合調(diào)制。陣列天線單陣元單脈沖回波數(shù)據(jù)的雜噪比為40 dB，且存在10%的幅度誤差以及6o的相位誤差。另外，為了更好地反映本文算法相對于已有算法在非均勻環(huán)境中的性能優(yōu)勢，仿真數(shù)據(jù)中注入了5個運動目標(biāo)信號以及4個強的旁瓣孤立雜波點信號，這些信號的空域流型也存在著相同的陣元誤差，其具體參數(shù)如表2所示。所注入的強雜波點信號可認(rèn)為來自于具有比背景雜波強得多的雷達(dá)后向散射系數(shù)的地面散射體，如橋梁、鐵塔、角反射器等，有助于增強仿真數(shù)據(jù)的非均勻性。

圖3給出的是仿真數(shù)據(jù)的常規(guī)脈沖多普勒處理結(jié)果，從圖中可以看出近程雜波具有很強的距離非均勻性，這種非均勻性往往使得傳統(tǒng)STAP方法不能獲得足夠的IID訓(xùn)練樣本，從而無法對CUT內(nèi)的雜波統(tǒng)計特性進(jìn)行有效估計，造成非正側(cè)視陣情況下的雜波抑制及運動目標(biāo)檢測性能的嚴(yán)重下降。圖中標(biāo)注出了所注入的運動目標(biāo)及孤立雜波點信號的位置，其中運動目標(biāo)1和目標(biāo)2被近程彎曲雜波所淹沒，運動目標(biāo)4為靠近主瓣雜波多普勒的弱小目標(biāo)。

圖2 算法流程圖

表1 雷達(dá)系統(tǒng)參數(shù)

表2 注入信號參數(shù)

圖3 脈沖多普勒處理結(jié)果

下面通過兩組實驗結(jié)果對本文算法的有效性進(jìn)行說明。實驗中選作性能對比的是兩種比較典型的傳統(tǒng)STAP方法［1］：擴(kuò)展因子化法（Ex tended Factored App roach， EFA）和聯(lián)合域局域化（Joint Domain Localized， JDL）方法。圖4給出了EFA，JDL以及本文算法對仿真數(shù)據(jù)進(jìn)行濾波處理后的距離多普勒結(jié)果，從圖4中可以看出EFA和JDL由于無法充分估計近程非均勻雜波的統(tǒng)計特性因而不能對其進(jìn)行有效的自適應(yīng)抑制，而經(jīng)過本文算法處理后上述非均勻雜波大部分都被有效濾除。圖5給出的是對3種方法的濾波輸出進(jìn)行CA-CFAR檢測的目標(biāo)檢測結(jié)果。在圖5（a）和圖5（b）中可以看到，STAP處理由于無法獲得孤立雜波點信號的訓(xùn)練樣本，因此無法對其進(jìn)行估計和抑制，4個強雜波點信號在CA-CFAR檢測后表現(xiàn)為虛警；強雜波點的存在也抬高了恒虛警檢測門限，造成對其附近運動目標(biāo)（如目標(biāo)3和目標(biāo)5）的漏警；運動目標(biāo)1和目標(biāo)2仍然被湮沒在較強的近程剩余雜波中而無法得到有效檢測；靠近主雜波的弱小目標(biāo)4在兩種STAP方法處理中也同樣沒有得到檢測。從圖5（c）可以看到，經(jīng)本文算法處理后，4個強雜波點信號均被有效濾除，而所有運動目標(biāo)也都得到了檢測。

圖6給出了EFA， JDL以及本文算法的檢測概率曲線比較。實驗選取第21號到30號多普勒通道、第201號到400號距離門范圍內(nèi)的回波數(shù)據(jù)做總共2000次的蒙特卡洛實驗，虛警概率設(shè)置為0.01。從圖6（a）中可以看出，EFA和JDL的檢測性能比較接近，而本文算法的檢測性能則遠(yuǎn)好于這兩種STAP方法，在檢測概率為0.8時，有超過4 dB的性能優(yōu)勢。另外，本文算法對先驗信息的依賴性較強，實驗中也對系統(tǒng)參數(shù)存在誤差的情況進(jìn)行了檢查。圖6（b）和6（c）中分別給出的是實際雷達(dá)系統(tǒng)中較為常見的偏航角誤差及載機速度誤差情況下算法的檢測性能比較，可以看出在偏航角存在0.1°， 0.5°及1°誤差時對算法性能影響不大；在存在較小的載機速度誤差時本文算法性能會受到一定損失，隨著速度誤差增大算法性能損失會加大。但從公開的國內(nèi)外器件性能指標(biāo)及發(fā)展趨勢來看，機載雷達(dá)系統(tǒng)參數(shù)的測量誤差會控制在一定范圍之內(nèi)并且精度會越來越高，因此本文算法仍然具有很高的實用價值。

圖4 雜波抑制后的距離多普勒圖

圖5 目標(biāo)檢測結(jié)果

圖6 檢測曲線比較

5 結(jié)論

本文提出了一種機載雷達(dá)非正側(cè)視陣情況下的非均勻雜波抑制方法。該方法充分利用了雷達(dá)系統(tǒng)參數(shù)及平臺運動狀態(tài)等先驗信息，通過離線構(gòu)造雜波表示基，可以大大提高算法的實時性；無需訓(xùn)練樣本，不用估計雜波統(tǒng)計特性，適用于極度非均勻回波數(shù)據(jù)的雜波抑制問題；對陣元誤差穩(wěn)健，提高對運動目標(biāo)的檢測性能；可以在沒有俯仰自由度的情況下對非正側(cè)視陣雜波進(jìn)行有效抑制。仿真結(jié)果驗證了本文算法的有效性。

［1］ Ward J. Space-time adaptive p rocessing for airborne radar［R］. Lincoln Laboratory， M IT， 1994.

［2］ K lemm R. Princip les of Space-T ime Adaptive Processing［M］. London: the Institution of Electrical Engineers， 2002: 117-149.

［3］ Brennan L E and Reed I S. Theory of adaptive radar［J］. IEEE Transactions on Aerospace and Electronic Systems，1973， 9（1）: 237-252.

［4］ Reed I S， M allett J D， and B rennan L E. Rap id convergence rate in adaptive arrays［J］. IEEE Transactions on Aerospace and Electronic System s， 1974， 10（4）: 853-863.

［5］ Melvin W L. Space-time adaptive radar performance in heterogeneous clutter［J］. IEEE Transactions on Aerospace and Electronic System s， 2000， 36（2）: 621-633.

［6］ Wang Yong-liang， Duan Ke-qing， and Xie Wen-chong. Cross beam STAP for nonstationary clutter supp ression in airborne radar［J］. International Journal of An tennas and Propagation，2013， ID: 276310.

［7］ 文珺， 史林， 廖桂生. 一種機載非正側(cè)視陣?yán)走_(dá)近程雜波對消方法［J］. 系統(tǒng)工程與電子技術(shù)， 2013， 35（10）: 2069-2073. Wen Jun， Shi Lin， and Liao Gui-sheng. Short-range clutter elim ination method for non-side looking airborne radar［J］. System s Engineering and Electronics， 2013， 35（10）: 2069-2073.

［8］ 王杰， 沈明威， 吳迪， 等. 基于主瓣雜波高效配準(zhǔn)的機載非正側(cè)視陣?yán)走_(dá)STAP算法研究［J］. 雷達(dá)學(xué)報， 2014， 3（2）: 235-240. Wang Jie， Shen M ing-wei， Wu Di， et al.. An efficient STAP algorithm for nonsidelooking airborne radar based on m ainlobe clutter com pensation［J］. Journal of Radars， 2014，3（2）: 235-240.

［9］ Yang Xiao-peng， Liu Yong-xu， and Long Teng. Robust non-homogeneity detection algorithm based on prolate spheroidal wave functions for space-tim e adaptive processing［J］. IET Radar， Sonar & Navigation， 2013， 7（1）: 47-54.

［10］ Chen Si-jia， Kong Ling-jiang， and Yang Jian-yu. Target detection for heterogeneous cyclostationary sea clutter［J］. Circuits， Systems， and Signal Processing， 2014， 33（3）: 959-971.

［11］ Aubry A， De M aio A， Pallotta L， et al.. Median matrices and their application to radar training data selection［J］. IET Radar， Sonar & Navigation， 2014， 8（4）: 265-274.

［12］ Zaimbashi A. An adaptive cell averaging-based CFAR detector for interfering targets and clutter-edge situations［J］. D igital Signal Processing， 2014， 31: 59-68.

［13］ Borsari G. M itigating effects on STAP p rocessing caused by an inclined array［C］. Proceedings of the IEEE National Radar Conference， Dallas， USA， 1998: 135-140.

［14］ K reyenkam p O and K lemm R. Doppler com pensation in forward-looking STAP radar［J］. IEE Proceedings of Radar，Sonar & Navigation， 2001， 148（5）: 253-258.

［15］ Himed B， Zhang Yu-hong， and Hajjari A. STAP w ith angle-Doppler com pensation for bistatic airborne radars［C］. Proceedings of the IEEE Radar Conference， California， USA，2002: 311-317.

［16］ Colone F. Spectral slope-based approach for m itigating bistatic space-time adaptive p rocessing clutter dispersion［J］. IET Radar， Sonar & Navigation， 2011， 5（5）: 593-603.

［17］ Ries P， Lapierre F D， Lesturgie M， et al.. Hand ling partially calibrated con formal antenna arrays in range-dependence com pensation for airborne bistatic space-time adaptive processing radar［J］. IET Radar， Sonar & Navigation， 2009，3（4）: 373-383.

［18］ Lapierre F D， Ries P， and Verly J G. Foundation for m itigating range dependence in radar space-time adaptive processing［J］. IET Radar， Sonar & Navigation， 2009， 3（1）: 18-29.

［19］ Corbell P M and Hale T B. 3-dimensional STAP performance analysis using the cross-spectral metric［C］. Proceedings of the IEEE Radar Conference， Philadelphia， USA， 2004: 610-615.

［20］ Meng Xiang-dong， Wang Tong， Wu Jian-xin， et al.. Shortrange clu tter suppression for airborne radar by utilizing prefiltering in elevation［J］. IEEE Geoscience and Rem ote Sensing Letters， 2009， 6（2）: 268-272.

［21］ Li Jian-feng， Zhang Xiao-fei， and Gao Xin. A joint scheme for angle and array gain-phase error estimation in bistatic M IMO radar［J］. IEEE Geoscience and Remote Sensing Letters， 2013，10（6）: 1478-1482.

［22］ Cao Sheng-hong， Ye Zhong-fu， Hu Nan， et al.. DOA estimation based on fourth-order cumulants in the presence of sensor gain-phase errors［J］. Signal Processing， 2013， 93（9）: 2581-2585.

同亞龍： 男，1987年生，博士生，研究方向為陣列信號處理、空時自適應(yīng)信號處理以及地面運動目標(biāo)檢測.

王 彤： 男，1974年生，教授，研究方向為機載雷達(dá)運動目標(biāo)檢測、合成孔徑雷達(dá)成像等.

文 才： 男，1986年生，博士生，研究方向為陣列信號處理、空時自適應(yīng)信號處理以及地面運動目標(biāo)檢測.

吳建新： 男，1982年生，副教授，研究方向為空時自適應(yīng)信號處理、雷達(dá)成像以及動目標(biāo)檢測等.

A Robust Clutter Suppression M ethod for A irborne Non-sidelooking Radar

Tong Ya-long Wang Tong Wen Cai Wu Jian-xin
（National Laboratory of Radar Signal Processing， Xidian University， Xi’an 710071， China）

Non-sidelooking array configuration leads to heterogeneous airborne radar echo data， it degrades seriously the clutter suppression performance of the traditional Space-Time Adap tive Processing （STAP）algorithm s. To solve this p rob lem， a heterogeneous clu tter suppression m ethod， which is robust to array error， is proposed. Firstly， a clutter rep resentation basis is constructed based on a priori information such as system parameters. Then， the test data is fitted in an iterative least square manner with consideration of array error， and the closed-form solution is derived and utilized here. Finally， pulse Doppler p rocessing and constant false alarm detection are conducted on the residual data. The proposed method does not need any training samp le and can effectively suppress the heterogeneous clutter for airborne non-sidelooking radar w ithout elevation degree of freedom. Simu lation results verify the validity of the p roposed m ethod.

Airborne radar； Non-sidelooking array； Heterogeneous clutter； Clutter suppression； Space-Time Adaptive Processing （STAP）

TN959.73

： A

：1009-5896（2015）05-1044-07

10.11999/JEIT141222

2014-09-19收到，2014-12-11改回

國家自然科學(xué)基金（61372133， 61471285）資助課題*通信作者：王彤 twang@mail.xidian.edu.cn