本發(fā)明屬于雷達(dá)雜波抑制技術(shù)領(lǐng)域，特別涉及一種非均勻雜波環(huán)境下的雷達(dá)雜波抑制方法，是一種基于非均勻雜波環(huán)境下的級聯(lián)二維空時(shí)處理方法，適用于斜側(cè)陣機(jī)載雷達(dá)雜波抑制或非正側(cè)陣機(jī)載雷達(dá)雜波抑制。

背景技術(shù)：

機(jī)載預(yù)警雷達(dá)已經(jīng)成為現(xiàn)代戰(zhàn)爭的預(yù)警情報(bào)、戰(zhàn)場監(jiān)視和作戰(zhàn)指揮系統(tǒng)的重要組成部分，且能夠顯著增加低空突防目標(biāo)的探測距離，大大提高武器系統(tǒng)的效能，使得在現(xiàn)在戰(zhàn)爭中扮演著越來越重要的角色，并受到各軍事大國的高度重視。機(jī)載預(yù)警雷達(dá)的波束指向一般是水平或稍下方向，使得機(jī)載預(yù)警雷達(dá)不可避免地會(huì)接收到地雜波。由于技術(shù)水平限制，使得陣列機(jī)載雷達(dá)天線的旁瓣電平較高，且陣列機(jī)載雷達(dá)回波中地雜波較強(qiáng)，可達(dá)90dB，遠(yuǎn)強(qiáng)于目標(biāo)回波；此外，由于飛機(jī)平臺(tái)的運(yùn)動(dòng)，地雜波的多普勒范圍大大展寬，使得可供檢測的多普勒范圍大大降低，有很大一部分目標(biāo)都“湮沒”在地雜波中。因此，為了可靠檢測遠(yuǎn)距離目標(biāo)，首先要解決機(jī)載雷達(dá)的雜波抑制問題。

空時(shí)自適應(yīng)處理(STAP)技術(shù)的誕生追溯到20世紀(jì)70年代，在雜波協(xié)方差矩陣和目標(biāo)信號(hào)分別確知的條件下，Brennan和Reed等在期刊Aerospace and Electronic Systems(AES)上，提出了全空時(shí)自適應(yīng)處理(STAP)的概念和理論，其思想是將陣列信號(hào)處理的基本原理推廣到由脈沖和陣元采樣的兩維場中；全空時(shí)自適應(yīng)處理(STAP)能夠取得比較理想的機(jī)載雷達(dá)雜波抑制效果，但機(jī)載雷達(dá)雜波協(xié)方差矩陣精確已知的條件在工程實(shí)現(xiàn)中難以滿足，通常情況下是由參考單元估計(jì)得到；為了使空時(shí)自適應(yīng)處理(STAP)器的性能損失不超過3dB，空時(shí)自適應(yīng)處理(STAP)器訓(xùn)練樣本個(gè)數(shù)至少應(yīng)為機(jī)載雷達(dá)系統(tǒng)自由度的兩倍，而且，機(jī)載雷達(dá)雜波協(xié)方差矩陣須從獨(dú)立同分布的訓(xùn)練樣本估計(jì)而來，所以全空時(shí)自適應(yīng)處理(STAP)的運(yùn)算量和設(shè)備復(fù)雜度令人難以接受。

為了降低運(yùn)算量，德國的R.Klemm博士于1987年在期刊Signal Processing上提出了輔助通道法，將空時(shí)自適應(yīng)處理(STAP)器維數(shù)從NM降至N+M-1，此處N表示機(jī)載雷達(dá)天線陣元個(gè)數(shù)，M表示一個(gè)相干處理間隔內(nèi)機(jī)載雷達(dá)接收的脈沖個(gè)數(shù)；在沒有誤差的理想情況下，該方法能夠取得比較理想的雜波抑制效果，性能接近最優(yōu)空時(shí)自適應(yīng)處理(STAP)處理，并且降低了估計(jì)雜波協(xié)方差矩陣所需要的獨(dú)立同分布參考單元數(shù)目，但在有誤差情況下，該方法的雜波抑制性能并不理想。

國內(nèi)的西安電子科技大學(xué)在降維空時(shí)自適應(yīng)處理(STAP)研究方面也展開了大量工作，并提出了局域化處理法，即因子化方法(1DT)和擴(kuò)展因子化方法(mDT)；該局域化處理法首先利用低旁瓣的多普勒濾波器對機(jī)載雷達(dá)回波進(jìn)行局域化處理，然后在空域自適應(yīng)處理或空時(shí)域聯(lián)合自適應(yīng)處理將雷達(dá)雜波抑制掉；但這種局域化的處理方法依然需要大量滿足獨(dú)立同分布的訓(xùn)練樣本來構(gòu)造協(xié)方差矩陣，在非均勻雜波環(huán)境下，幾乎很難得到大量滿足獨(dú)立同分布要求的樣本，所以在非均勻雜波環(huán)境下的雜波抑制性能并沒有得到大的改善。

技術(shù)實(shí)現(xiàn)要素：

針對上述現(xiàn)有技術(shù)存在的不足，本發(fā)明的目的在于提出一種非均勻雜波環(huán)境下的雷達(dá)雜波抑制方法，該種非均勻雜波環(huán)境下的雷達(dá)雜波抑制方法能夠充分利用待檢樣本附近的少量樣本作為局部訓(xùn)練樣本，采用相干積累的形式保證信號(hào)增益，同時(shí)又保證了足夠多的自由度進(jìn)行雷達(dá)雜波的充分抑制。

為達(dá)到上述技術(shù)目的，本發(fā)明采用如下技術(shù)方案予以實(shí)現(xiàn)。

一種非均勻雜波環(huán)境下的雷達(dá)雜波抑制方法，包括以下步驟：

步驟1，確定機(jī)載雷達(dá)，并獲取機(jī)載雷達(dá)的原始雷達(dá)回波數(shù)據(jù)X，然后確定機(jī)載雷達(dá)的時(shí)域滑窗孔徑K_t；

步驟2，初始化：令k∈{1,2,…,K}，k表示第k個(gè)多普勒通道，K表示機(jī)載雷達(dá)的原始雷達(dá)回波數(shù)據(jù)X在頻率域包含的多普勒總個(gè)數(shù)，k的初始值為1；令l∈{1,2,…,L}，l表示第l個(gè)距離門，L表示機(jī)載雷達(dá)的原始雷達(dá)回波數(shù)據(jù)X包含的距離門總個(gè)數(shù)，l的初始值為1；

步驟3，依次計(jì)算機(jī)載雷達(dá)第k個(gè)多普勒通道進(jìn)行空時(shí)處理的時(shí)域?qū)蚴噶縎_kt，以及機(jī)載雷達(dá)第k個(gè)多普勒通道的最終空時(shí)導(dǎo)向矢量S_k；

步驟4，將第l個(gè)距離門的原始雷達(dá)回波數(shù)據(jù)記為x_l，并計(jì)算第l個(gè)距離門的擴(kuò)展雷達(dá)回波數(shù)據(jù)Q_l，然后依次計(jì)算第l個(gè)距離門對應(yīng)的訓(xùn)練樣本數(shù)據(jù)Train_l、時(shí)域滑窗后第l個(gè)距離門對應(yīng)的擴(kuò)展訓(xùn)練樣本數(shù)據(jù)Train_l′，以及第l個(gè)距離門對應(yīng)的協(xié)方差矩陣R_l；

步驟5，計(jì)算機(jī)載雷達(dá)第k個(gè)多普勒通道、第l個(gè)距離門對應(yīng)的自適應(yīng)權(quán)值w_kl，進(jìn)而依次計(jì)算空時(shí)處理過程后第k個(gè)多普勒通道、第l個(gè)距離門的對應(yīng)輸出矢量y_kl和機(jī)載雷達(dá)第k個(gè)多普勒通道、第l個(gè)距離門對應(yīng)的濾波輸出z_kl；

步驟6，令k加1，依次重復(fù)步驟3至步驟5，直到得到機(jī)載雷達(dá)第K個(gè)多普勒通道、第l個(gè)距離門對應(yīng)的濾波輸出z_Kl，并將此時(shí)得到的機(jī)載雷達(dá)第1個(gè)多普勒通道、第l個(gè)距離門對應(yīng)的濾波輸出z_1l至機(jī)載雷達(dá)第K個(gè)多普勒通道、第l個(gè)距離門對應(yīng)的濾波輸出z_Kl，作為機(jī)載雷達(dá)第l個(gè)距離門對應(yīng)的濾波輸出z_l，z_l＝[z_1l z_2l … z_kl … z_Kl]^H，上標(biāo)H表示共軛轉(zhuǎn)置；

步驟7，令l加1,依次重復(fù)步驟3至步驟6，直到得到機(jī)載雷達(dá)第L個(gè)距離門對應(yīng)的濾波輸出z_L，并將此時(shí)得到的機(jī)載雷達(dá)第l個(gè)距離門對應(yīng)的濾波輸出z_l至機(jī)載雷達(dá)第L個(gè)距離門對應(yīng)的濾波輸出z_L，作為雜波抑制后機(jī)載雷達(dá)的距離—多普勒雷達(dá)回波數(shù)據(jù)Z，Z＝[z₁ z₂… z_l … z_L]。

本發(fā)明與現(xiàn)有技術(shù)相比具有以下優(yōu)點(diǎn)：

第一，本發(fā)明中提出的方法在訓(xùn)練樣本較少的情況下，能夠通過時(shí)域滑窗的方式提供充足的自由度來抑制雷達(dá)雜波。

第二，使用常規(guī)的擴(kuò)展因子化方法(mDT)估計(jì)機(jī)載雷達(dá)雜波協(xié)方差矩陣時(shí)，選取的訓(xùn)練樣本距離維跨度較大，導(dǎo)致在非均勻雜波環(huán)境下，不能很好地估計(jì)當(dāng)前距離門對應(yīng)的協(xié)方差矩陣，而本發(fā)明方法使得訓(xùn)練樣本選擇局域化，修正了擴(kuò)展因子化方法協(xié)方差矩陣估計(jì)不準(zhǔn)的缺陷，m＝3。

附圖說明

下面結(jié)合附圖和具體實(shí)施方式對本發(fā)明作進(jìn)一步詳細(xì)說明。

圖1是本發(fā)明的一種非均勻雜波環(huán)境下的雷達(dá)雜波抑制方法流程圖；

圖2(a)是使用擴(kuò)展因子化方法(mDT)進(jìn)行雜波抑制后得到的距離—多普勒圖，其中橫坐標(biāo)為多普勒通道，縱坐標(biāo)為距離門，m＝3；

圖2(b)是使用本發(fā)明方法進(jìn)行雜波抑制后得到的距離—多普勒圖，其中橫坐標(biāo)為多普勒通道，縱坐標(biāo)為距離門；

圖3(a)是分別使用擴(kuò)展因子化方法(mDT)和本發(fā)明方法進(jìn)行雜波抑制處理后1-150號(hào)距離門對應(yīng)的雜波剩余圖，其中橫坐標(biāo)表示多普勒通道，縱坐標(biāo)表示雜波剩余，單位為dB，m＝3；

圖3(b)是分別使用擴(kuò)展因子化方法(mDT)和本發(fā)明方法進(jìn)行雜波抑制處理后150-349號(hào)距離門對應(yīng)的雜波剩余圖；其中橫坐標(biāo)表示多普勒通道，縱坐標(biāo)表示雜波剩余，單位為dB，m＝3。

具體實(shí)施方式

參考圖1，為本發(fā)明的一種非均勻雜波環(huán)境下的雷達(dá)雜波抑制方法流程圖；本發(fā)明的一種機(jī)載雷達(dá)雜波抑制方法，包括以下步驟：

步驟1，確定機(jī)載雷達(dá)，并獲取機(jī)載雷達(dá)的原始雷達(dá)回波數(shù)據(jù)X，然后確定機(jī)載雷達(dá)的時(shí)域滑窗孔徑K_t。

具體地，確定機(jī)載雷達(dá)，并獲取機(jī)載雷達(dá)的原始雷達(dá)回波數(shù)據(jù)X，所述原始雷達(dá)回波數(shù)據(jù)X為N×M×L維矩陣，N表示機(jī)載雷達(dá)的天線陣面方位向均勻包含的陣元個(gè)數(shù)，M表示機(jī)載雷達(dá)在每個(gè)相干處理周期內(nèi)發(fā)射的脈沖個(gè)數(shù)；L表示機(jī)載雷達(dá)的原始雷達(dá)回波數(shù)據(jù)X包含的距離門總個(gè)數(shù)；定義機(jī)載雷達(dá)的原始雷達(dá)回波數(shù)據(jù)X在頻率域包含的多普勒總個(gè)數(shù)為K，且K得取值為大于M且為2的冪次方整數(shù)；然后確定機(jī)載雷達(dá)的滑窗時(shí)域孔徑K_t，使其滿足下式：

2×(M-K_t+1)≥2×K_t×N

或者K_t為整數(shù)。

步驟2，初始化：令k∈{1,2,…,K}，k表示第k個(gè)多普勒通道，k的初始值為1；K表示機(jī)載雷達(dá)的原始雷達(dá)回波數(shù)據(jù)X在頻率域包含的多普勒總個(gè)數(shù)；令l∈{1,2,…,L}，l表示第l個(gè)距離門，L表示機(jī)載雷達(dá)的原始雷達(dá)回波數(shù)據(jù)X包含的距離門總個(gè)數(shù)，l的初始值為1。

步驟3，依次計(jì)算機(jī)載雷達(dá)第k個(gè)多普勒通道進(jìn)行空時(shí)處理的時(shí)域?qū)蚴噶縎_kt，以及機(jī)載雷達(dá)第k個(gè)多普勒通道的最終空時(shí)導(dǎo)向矢量S_k。

具體地，將機(jī)載雷達(dá)空時(shí)處理過程中的空域?qū)蚴噶坑洖镾_s，其表達(dá)式為：

其中，n∈{0,1,2,…,N-1}，表示第n個(gè)陣元接收機(jī)載雷達(dá)的脈沖回波數(shù)據(jù)時(shí)相對第1個(gè)陣元的相位偏移，f_s表示機(jī)載雷達(dá)的空域頻率，N表示機(jī)載雷達(dá)天線陣面方位向均勻包含的陣元個(gè)數(shù)，上標(biāo)H表示共軛轉(zhuǎn)置。

將機(jī)載雷達(dá)第k個(gè)多普勒通道進(jìn)行空時(shí)處理的時(shí)域?qū)蚴噶坑洖镾_kt，其表達(dá)式為：

其中，上標(biāo)H表示共軛轉(zhuǎn)置，m∈{0,1,2,…,M-1}，M表示機(jī)載雷達(dá)在每個(gè)相干處理周期內(nèi)發(fā)射的脈沖個(gè)數(shù)，表示機(jī)載雷達(dá)第m個(gè)脈沖相對第1個(gè)脈沖的相位偏移，f_dk表示機(jī)載雷達(dá)第k個(gè)多普勒通道歸一化的時(shí)域多普勒頻率，且f_dk＝k-1/K；機(jī)載雷達(dá)的滑窗時(shí)域孔徑K_t遠(yuǎn)遠(yuǎn)小于M，則將機(jī)載雷達(dá)第k個(gè)多普勒通道進(jìn)行空時(shí)處理的時(shí)域?qū)蚴噶縎_kt中第1到第K_t個(gè)元素，作為機(jī)載雷達(dá)第k個(gè)多普勒通道進(jìn)行空時(shí)處理的優(yōu)化時(shí)域?qū)蚴噶縎′_kt，其表達(dá)式為：

S′_kt＝S_kt[1:K_t]

其中，S_kt[1:K_t]表示取機(jī)載雷達(dá)第k個(gè)多普勒通道進(jìn)行空時(shí)處理的時(shí)域?qū)蚴噶縎_kt中第1到第K_t個(gè)元素構(gòu)成的矢量；進(jìn)而計(jì)算得到機(jī)載雷達(dá)第k個(gè)多普勒通道的最終空時(shí)導(dǎo)向矢量S_k，S_s表示機(jī)載雷達(dá)空時(shí)處理過程中的空域?qū)蚴噶浚硎綤ronecker乘積。

步驟4，將第l個(gè)距離門的原始雷達(dá)回波數(shù)據(jù)記為x_l，并計(jì)算第l個(gè)距離門的擴(kuò)展雷達(dá)回波數(shù)據(jù)Q_l，然后依次計(jì)算第l個(gè)距離門對應(yīng)的訓(xùn)練樣本數(shù)據(jù)Train_l、時(shí)域滑窗后第l個(gè)距離門對應(yīng)的擴(kuò)展訓(xùn)練樣本數(shù)據(jù)Train_l′，以及第l個(gè)距離門對應(yīng)的協(xié)方差矩陣R_l。

具體的，定義一種滑窗處理過程(SlideWinProc)如下：

首先，將第l個(gè)距離門的原始雷達(dá)回波數(shù)據(jù)記為x_l，其維度為N×M，表達(dá)式為：

x_l＝[x₁′_l x′_2l … x′_ml … x′_Ml]

其中，x′_ml表示第m個(gè)脈沖、第l個(gè)距離門的原始雷達(dá)回波數(shù)據(jù)，且是N×1維列矢量；m＝1,2,…,M，l＝1,2,…,L；由于機(jī)載雷達(dá)的時(shí)域滑窗孔徑K_t滿足且K_t為整數(shù)，則任意取一個(gè)滿足條件的整數(shù)r，作為機(jī)載雷達(dá)的時(shí)域滑窗孔徑值，進(jìn)而對第l個(gè)距離門的原始雷達(dá)回波數(shù)據(jù)x_l進(jìn)行孔徑值為r的時(shí)域滑窗，得到第l個(gè)距離門的擴(kuò)展雷達(dá)回波數(shù)據(jù)Q_l，維度為N_r×V，其表達(dá)式為：

Q_l＝[q_1l q_2l … q_vl … q_Vl]

其中，N_r＝N×r，V表示第l個(gè)距離門的擴(kuò)展雷達(dá)回波數(shù)據(jù)Q_l包含的樣本個(gè)數(shù)，V＝M-r+1；q_vl表示第v個(gè)樣本、第l個(gè)距離門的擴(kuò)展雷達(dá)回波數(shù)據(jù)，q_vl＝[x′_vl x′_(v+1)l … x′_(v+r-1)l]^H，x′_vl表示第v個(gè)樣本、第l個(gè)距離門的原始雷達(dá)回波數(shù)據(jù)，上標(biāo)H表示共軛轉(zhuǎn)置，v＝1,2,…,V；所述第v個(gè)樣本與對機(jī)載雷達(dá)的原始雷達(dá)回波數(shù)據(jù)X進(jìn)行滑窗的起始脈沖編號(hào)對應(yīng)。

根據(jù)第l個(gè)距離門的原始雷達(dá)回波數(shù)據(jù)x_l，分別將第l-1個(gè)距離門的原始雷達(dá)回波數(shù)據(jù)x_l-1、將第l+1個(gè)距離門的原始雷達(dá)回波數(shù)據(jù)x_l+1，進(jìn)而得到第l個(gè)距離門對應(yīng)的訓(xùn)練樣本數(shù)據(jù)Train_l，Train_l＝[x_l-1 x_l+1]，其維度為N×M×2。

然后，對第l個(gè)距離門對應(yīng)的訓(xùn)練樣本數(shù)據(jù)Train_l中包含的所有元素分別進(jìn)行孔徑值為r的時(shí)域滑窗，得到時(shí)域滑窗后第l個(gè)距離門對應(yīng)的擴(kuò)展訓(xùn)練樣本數(shù)據(jù)Train_l′，維度為N_r×2V，N_r＝N×r；Train_l′＝[Q_l-1 Q_l+1]，Q_l-1表示第l-1個(gè)距離門的擴(kuò)展雷達(dá)回波數(shù)據(jù)，Q_l+1表示第l+1個(gè)距離門的擴(kuò)展雷達(dá)回波數(shù)據(jù)。

最后，根據(jù)時(shí)域滑窗后第l個(gè)距離門對應(yīng)的擴(kuò)展訓(xùn)練樣本數(shù)據(jù)Train_l′，計(jì)算第l個(gè)距離門對應(yīng)的協(xié)方差矩陣R_l，其表達(dá)式為：

其中，上標(biāo)H表示共軛轉(zhuǎn)置。

步驟5，計(jì)算機(jī)載雷達(dá)第k個(gè)多普勒通道、第l個(gè)距離門對應(yīng)的自適應(yīng)權(quán)值w_kl，進(jìn)而依次計(jì)算空時(shí)處理過程后第k個(gè)多普勒通道、第l個(gè)距離門的對應(yīng)輸出矢量y_kl和機(jī)載雷達(dá)第k個(gè)多普勒通道、第l個(gè)距離門對應(yīng)的濾波輸出z_kl。

具體地，所述機(jī)載雷達(dá)第k個(gè)多普勒通道、第l個(gè)距離門對應(yīng)的自適應(yīng)權(quán)值w_kl，其計(jì)算式為：

其中，R_l表示第l個(gè)距離門對應(yīng)的協(xié)方差矩陣，S_k表示機(jī)載雷達(dá)第k個(gè)多普勒通道的最終空時(shí)導(dǎo)向矢量，上標(biāo)H表示共軛轉(zhuǎn)置，上標(biāo)-1表示求逆操作。

將第l個(gè)距離門的擴(kuò)展雷達(dá)回波數(shù)據(jù)Q_l，作為空時(shí)處理過程中第l個(gè)距離門的輸入數(shù)據(jù)，然后將第k個(gè)多普勒通道、第l個(gè)距離門對應(yīng)的自適應(yīng)權(quán)值記為w_kl，進(jìn)而計(jì)算空時(shí)處理過程后第k個(gè)多普勒通道、第l個(gè)距離門的對應(yīng)輸出矢量y_kl，其表達(dá)式為：

其中，上標(biāo)H表示共軛轉(zhuǎn)置。

根據(jù)機(jī)載雷達(dá)第k個(gè)多普勒通道空時(shí)處理的時(shí)域?qū)蚴噶縎_kt和空時(shí)處理過程后第k個(gè)多普勒通道、第l個(gè)距離門的對應(yīng)輸出矢量y_kl，計(jì)算得到機(jī)載雷達(dá)第k個(gè)多普勒通道、第l個(gè)距離門對應(yīng)的濾波輸出z_kl，其表達(dá)式為：

z_kl＝y(tǒng)_kl×conj(S_kt[1:V])

其中，conj[·]表示取共軛操作，S_kt[1:V]表示取機(jī)載雷達(dá)第k個(gè)多普勒通道空時(shí)處理的時(shí)域?qū)蚴噶縎_kt中第1到第V個(gè)元素構(gòu)成的矢量；V表示第l個(gè)距離門的擴(kuò)展雷達(dá)回波數(shù)據(jù)Q_l包含的樣本個(gè)數(shù)。

所述機(jī)載雷達(dá)第k個(gè)多普勒通道、第l個(gè)距離門單元對應(yīng)的濾波輸出z_kl完成了第k個(gè)多普勒通道、第l個(gè)距離門的雜波抑制。

步驟6，令k加1，依次重復(fù)步驟3至步驟5，直到得到機(jī)載雷達(dá)第K個(gè)多普勒通道、第l個(gè)距離門對應(yīng)的濾波輸出z_Kl，并將此時(shí)得到的機(jī)載雷達(dá)第1個(gè)多普勒通道、第l個(gè)距離門對應(yīng)的濾波輸出z_1l至機(jī)載雷達(dá)第K個(gè)多普勒通道、第l個(gè)距離門對應(yīng)的濾波輸出z_Kl，作為機(jī)載雷達(dá)第l個(gè)距離門對應(yīng)的濾波輸出z_l，z_l＝[z_1l z_2l … z_kl … z_Kl]^H，上標(biāo)H表示共軛轉(zhuǎn)置。

步驟7，令l加1,依次重復(fù)步驟3至步驟6，直到得到機(jī)載雷達(dá)第L個(gè)距離門對應(yīng)的濾波輸出z_L，并將此時(shí)得到的機(jī)載雷達(dá)第l個(gè)距離門對應(yīng)的濾波輸出z_l至機(jī)載雷達(dá)第L個(gè)距離門對應(yīng)的濾波輸出z_L，作為雜波抑制后機(jī)載雷達(dá)的距離—多普勒雷達(dá)回波數(shù)據(jù)Z，Z＝[z₁ z₂… z_l … z_L]。

至此，本發(fā)明的一種非均勻雜波環(huán)境下的雷達(dá)雜波抑制方法結(jié)束。

通過以下仿真實(shí)驗(yàn)對本發(fā)明效果作進(jìn)一步驗(yàn)證說明。

(一)實(shí)驗(yàn)參數(shù)及實(shí)驗(yàn)條件

機(jī)載雷達(dá)的天線陣面結(jié)構(gòu)采用橢圓斜側(cè)陣面，并且各個(gè)陣元均勻排列在橢圓斜側(cè)陣面上，機(jī)載雷達(dá)的天線陣面的方位向和俯仰向分別均勻排列56個(gè)陣元和10個(gè)陣元，陣元間距為d＝λ/2，λ為機(jī)載雷達(dá)的載波波長，機(jī)載雷達(dá)在每個(gè)相干處理周期內(nèi)發(fā)射的脈沖個(gè)數(shù)M為220，機(jī)載雷達(dá)的原始雷達(dá)回波數(shù)據(jù)X包含的距離門個(gè)數(shù)為349。

(二)實(shí)驗(yàn)內(nèi)容及結(jié)果分析

分別采用擴(kuò)展因子化方法(mDT)和本發(fā)明方法對原始雷達(dá)回波數(shù)據(jù)進(jìn)行雜波抑制，得到的兩幅距離—多普勒圖分別如圖2(a)、圖2(b)所示；圖2(a)是使用擴(kuò)展因子化方法(mDT)進(jìn)行雜波抑制后得到的距離—多普勒圖，其中橫坐標(biāo)為多普勒通道，縱坐標(biāo)為距離門，m＝3；圖2(b)是使用本發(fā)明方法進(jìn)行雜波抑制后得到的距離—多普勒圖，其中橫坐標(biāo)為多普勒通道，縱坐標(biāo)為距離門；圖2(a)、圖2(b)中的白色區(qū)域分別表示雜波抑制后的雜波剩余功率分布，黑色區(qū)域分別表示雜波抑制后的噪聲剩余功率分布。

從圖2(a)中可以看出，采用擴(kuò)展因子化方法(mDT)對原始雷達(dá)回波數(shù)據(jù)進(jìn)行雜波抑制后，雜波剩余功率較強(qiáng)，尤其在旁瓣雜波區(qū)，m＝3；從圖2(b)中可以看出，本發(fā)明方法能夠有效降低雜波剩余，對旁瓣雜波的抑制性能較好，有利于后續(xù)對目標(biāo)進(jìn)行檢測。

分別使用擴(kuò)展因子化方法(mDT)和本發(fā)明方法對對原始雷達(dá)回波數(shù)據(jù)進(jìn)行雜波抑制，分別得到兩種方法的雜波剩余比較圖，如圖3(a)和圖3(b)所示；圖3(a)是分別使用擴(kuò)展因子化方法(mDT)和本發(fā)明方法進(jìn)行雜波抑制處理后1-150號(hào)距離門對應(yīng)的雜波剩余圖，其中橫坐標(biāo)表示多普勒通道，縱坐標(biāo)表示雜波剩余，單位為dB，m＝3；圖3(b)是分別使用擴(kuò)展因子化方法(mDT)和本發(fā)明方法進(jìn)行雜波抑制處理后150-349號(hào)距離門對應(yīng)的雜波剩余圖；其中橫坐標(biāo)表示多普勒通道，縱坐標(biāo)表示雜波剩余，單位為dB，m＝3。

從圖3(a)和圖3(b)可以看出，與擴(kuò)展因子化方法(mDT)相比，本發(fā)明方法在保證信號(hào)增益不變的前提下有效抑制旁瓣雜波，m＝3。

綜上所述，仿真實(shí)驗(yàn)驗(yàn)證了本發(fā)明的正確性，有效性和可靠性。

顯然，本領(lǐng)域的技術(shù)人員可以對本發(fā)明進(jìn)行各種改動(dòng)和變型而不脫離本發(fā)明的精神和范圍；這樣，倘若本發(fā)明的這些修改和變型屬于本發(fā)明權(quán)利要求及其等同技術(shù)的范圍之內(nèi)，則本發(fā)明也意圖包含這些改動(dòng)和變型在內(nèi)。