技術(shù)特征：

1.一種非均勻雜波環(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è)距離門對(duì)應(yīng)的訓(xùn)練樣本數(shù)據(jù)Train_l、時(shí)域滑窗后第l個(gè)距離門對(duì)應(yīng)的擴(kuò)展訓(xùn)練樣本數(shù)據(jù)Train_l′，以及第l個(gè)距離門對(duì)應(yīng)的協(xié)方差矩陣R_l；

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

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

2.如權(quán)利要求1所述的一種非均勻雜波環(huán)境下的雷達(dá)雜波抑制方法，其特征在于，在步驟1中，所述原始雷達(dá)回波數(shù)據(jù)X為N×M×L維矩陣，N表示機(jī)載雷達(dá)的天線陣面方位向均勻包含的陣元個(gè)數(shù)，L表示機(jī)載雷達(dá)的原始雷達(dá)回波數(shù)據(jù)X包含的距離門個(gè)數(shù)，M表示機(jī)載雷達(dá)在每個(gè)相干處理周期內(nèi)發(fā)射的脈沖個(gè)數(shù)；

所述機(jī)載雷達(dá)的時(shí)域滑窗孔徑K_t，其滿足條件為：

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

或者K_t為整數(shù)。

3.如權(quán)利要求1所述的一種非均勻雜波環(huán)境下的雷達(dá)雜波抑制方法，其特征在于，在步驟3中，所述機(jī)載雷達(dá)第k個(gè)多普勒通道進(jìn)行空時(shí)處理的時(shí)域?qū)蚴噶縎_kt，以及所述機(jī)載雷達(dá)第k個(gè)多普勒通道的最終空時(shí)導(dǎo)向矢量S_k，其表達(dá)式分別為：

$S_{kt}={\left[\begin{array}{cccccc}1& e^{j2{\mathrm{\πf}}_{dk}}& ...& e^{j2{\mathrm{\πf}}_{dk}m}& ...& e^{j2{\mathrm{\πf}}_{dk}(M-1)}\end{array}\right]}^H$

$S_k={S_{kt}}^{\′}\⊗S_s$

其中，S_s表示機(jī)載雷達(dá)空時(shí)處理過程中的空域?qū)蚴噶浚?img id="icf0004" file="FDA0001161887370000024.GIF" wi="43" he="47" img-content="drawing" img-format="GIF" orientation="portrait" inline="no" />表示Kronecker乘積，n∈{0,1,2,…,N-1}，表示第n個(gè)陣元接收機(jī)載雷達(dá)的脈沖回波數(shù)據(jù)時(shí)相對(duì)第1個(gè)陣元的相位偏移，f_s表示機(jī)載雷達(dá)的空域頻率，N表示機(jī)載雷達(dá)天線陣面方位向均勻包含的陣元個(gè)數(shù)，上標(biāo)H表示共軛轉(zhuǎn)置，

m∈{0,1,2,…,M-1}，M表示機(jī)載雷達(dá)在每個(gè)相干處理周期內(nèi)發(fā)射的脈沖個(gè)數(shù)，表示機(jī)載雷達(dá)第m個(gè)脈沖相對(duì)第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，S′_kt＝S_kt[1:K_t]。

4.如權(quán)利要求1所述的一種非均勻雜波環(huán)境下的雷達(dá)雜波抑制方法，其特征在于，步驟4的具體過程為：

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

x_l＝[x′_1l x′_2l … x′_ml … x′_Ml]

其中，x′_ml表示第m個(gè)脈沖、第l個(gè)距離門的原始雷達(dá)回波數(shù)據(jù)，且是N×1維列矢量；在滿足且K_t為整數(shù)的條件下任意取一個(gè)整數(shù)r，作為機(jī)載雷達(dá)的時(shí)域滑窗孔徑值，進(jìn)而對(duì)第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è)樣本與對(duì)機(jī)載雷達(dá)的原始雷達(dá)回波數(shù)據(jù)X進(jìn)行滑窗的起始脈沖編號(hào)對(duì)應(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è)距離門對(duì)應(yīng)的訓(xùn)練樣本數(shù)據(jù)Train_l，Train_l＝[x_l-1 x_l+1]，其維度為N×M×2；

然后，對(duì)第l個(gè)距離門對(duì)應(yīng)的訓(xùn)練樣本數(shù)據(jù)Train_l中包含的所有元素分別進(jìn)行孔徑值為r的時(shí)域滑窗，得到時(shí)域滑窗后第l個(gè)距離門對(duì)應(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è)距離門對(duì)應(yīng)的擴(kuò)展訓(xùn)練樣本數(shù)據(jù)Train_l′，計(jì)算第l個(gè)距離門對(duì)應(yīng)的協(xié)方差矩陣R_l，其表達(dá)式為：

5.如權(quán)利要求1所述的一種非均勻雜波環(huán)境下的雷達(dá)雜波抑制方法，其特征在于，在步驟5中，所述機(jī)載雷達(dá)第k個(gè)多普勒通道、第l個(gè)距離門對(duì)應(yīng)的自適應(yīng)權(quán)值w_kl，其計(jì)算式為：

$w_{kl}=\frac{{R_l}^{-1}{S}_k}{S_k^H{R_l}^{-1}{S}_k}$

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

根據(jù)機(jī)載雷達(dá)第k個(gè)多普勒通道空時(shí)處理的時(shí)域?qū)蚴噶縎_kt和空時(shí)處理過程后第k個(gè)多普勒通道、第l個(gè)距離門的對(duì)應(yīng)輸出矢量y_kl，計(jì)算得到機(jī)載雷達(dá)第k個(gè)多普勒通道、第l個(gè)距離門對(duì)應(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ù)；y_kl表示空時(shí)處理過程后第k個(gè)多普勒通道、第l個(gè)距離門的對(duì)應(yīng)輸出矢量，Q_l表示第l個(gè)距離門的擴(kuò)展雷達(dá)回波數(shù)據(jù)，同時(shí)將第l個(gè)距離門的擴(kuò)展雷達(dá)回波數(shù)據(jù)作為空時(shí)處理過程中第l個(gè)距離門的輸入數(shù)據(jù)。