999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

波形分集陣列雷達(dá)抗干擾進(jìn)展

2021-05-31 11:19:32許京偉朱圣棋廖桂生張玉洪
關(guān)鍵詞:方向信號(hào)方法

蘭 嵐, 許京偉, 朱圣棋, 廖桂生, 張玉洪

(1. 西安電子科技大學(xué)雷達(dá)信號(hào)處理國(guó)家級(jí)重點(diǎn)實(shí)驗(yàn)室, 陜西 西安 710071;2. 西安電子科技大學(xué)電子工程學(xué)院, 陜西 西安 710071)

0 引 言

雷達(dá)系統(tǒng)面臨復(fù)雜電磁環(huán)境的嚴(yán)峻挑戰(zhàn),來(lái)自周圍環(huán)境的有意和無(wú)意干擾、強(qiáng)雜波等嚴(yán)重限制了其信息獲取能力[1-4]。欺騙式干擾是一種極具威脅的電磁干擾形式,其通過(guò)輻射與真實(shí)目標(biāo)回波相似的電磁波,使雷達(dá)系統(tǒng)誤將虛假目標(biāo)當(dāng)成真實(shí)目標(biāo),導(dǎo)致真實(shí)目標(biāo)丟失、資源占用、異??涨榈葒?yán)重后果,雷達(dá)系統(tǒng)性能急劇惡化[5-7]。干擾機(jī)將截獲的雷達(dá)信號(hào)進(jìn)行距離、多普勒調(diào)制后可在一定空域內(nèi)實(shí)施欺騙,掩護(hù)目標(biāo)進(jìn)行入侵、突防等任務(wù)。數(shù)字射頻存儲(chǔ)器(digital radio frequency memory, DRFM)技術(shù)迅速發(fā)展,通過(guò)對(duì)存儲(chǔ)的信號(hào)進(jìn)行讀取、復(fù)制和參數(shù)調(diào)制(信號(hào)延時(shí)及多普勒參數(shù)),干擾機(jī)可轉(zhuǎn)發(fā)生成與目標(biāo)雷達(dá)信號(hào)相干的欺騙干擾樣式,即通過(guò)“截獲-調(diào)制-轉(zhuǎn)發(fā)”的方式產(chǎn)生大量虛假目標(biāo),加劇了雷達(dá)系統(tǒng)對(duì)干擾環(huán)境感知的難度[8-9]。

對(duì)于來(lái)自旁瓣的欺騙式信號(hào),陣列雷達(dá)能夠采用廣義旁瓣相消器(generalized sidelobe canceller,GSC)[10]、超低旁瓣天線[11]、空時(shí)自適應(yīng)處理(space-time adaptive processing,STAP)[12]等技術(shù)有效地進(jìn)行自適應(yīng)干擾對(duì)抗。然而,若假目標(biāo)信號(hào)位于主瓣區(qū)域,則會(huì)影響雷達(dá)系統(tǒng)對(duì)真實(shí)目標(biāo)的探測(cè),嚴(yán)重限制了雷達(dá)系統(tǒng)獲取信息的能力。由于雷達(dá)波束通常具有一定的寬度,因此位于主波束內(nèi)的干擾在空間上具有一定的角度范圍。通常考慮兩類主瓣欺騙式干擾。一般意義上的主瓣干擾是指干擾位于雷達(dá)主波束內(nèi),典型的應(yīng)用如誘餌、吊艙,通過(guò)轉(zhuǎn)發(fā)截獲的雷達(dá)信號(hào),與目標(biāo)直接散射的信號(hào)構(gòu)成相干散射源,起到對(duì)目標(biāo)的保護(hù)作用。目前,國(guó)內(nèi)外的武器系統(tǒng)中普遍配備有源誘餌設(shè)備,如美國(guó)“愛(ài)國(guó)者”反導(dǎo)系統(tǒng)的AN/MPQ-53雷達(dá)系統(tǒng),美國(guó)小型空射誘餌MALD,美國(guó)AN/ALQ-144型紅外干擾吊艙等。另外,主瓣干擾方向可與目標(biāo)角度完全相同,通常由目標(biāo)本身攜帶的干擾機(jī)產(chǎn)生。此時(shí)在空域無(wú)法對(duì)目標(biāo)與干擾進(jìn)行區(qū)分,抗干擾難度加大。因此,對(duì)主瓣欺騙式干擾抑制方法的研究已迫在眉睫。

目前,主要在空域、時(shí)域、頻域、極化域及聯(lián)合域?qū)χ靼旮蓴_進(jìn)行抑制,具體包括:① 頻域抗干擾:如頻率捷變[13-14],其中文獻(xiàn)[15]提出了一種捷變頻聯(lián)合Hough 變換方法,文獻(xiàn)[16]結(jié)合壓縮感知(compressed sensing,CS)理論進(jìn)行成像;文獻(xiàn)[17]通過(guò)構(gòu)造一組具有頻域正交性的發(fā)射波形實(shí)現(xiàn)距離波門拖引下假目標(biāo)的對(duì)消;此外可利用參差重頻進(jìn)行對(duì)抗[18]。② 空域抗干擾:在自適應(yīng)波束形成基礎(chǔ)上,通過(guò)利用阻塞矩陣[19]、正交投影[20]、特征投影矩陣[21]等設(shè)計(jì)數(shù)據(jù)預(yù)處理矩陣克服主波束畸變,此外,可結(jié)合單脈沖技術(shù)[22]、構(gòu)造輔助子波束空間[23]等。③ 時(shí)域抗干擾:可利用盲源分離將干擾和目標(biāo)分開(kāi)[24-26],可結(jié)合脈沖分集技術(shù)進(jìn)行了對(duì)密集假目標(biāo)進(jìn)行抑制[27],此外,可利用凸優(yōu)化時(shí)域稀疏恢復(fù)抑制主瓣欺騙式干擾信號(hào)[28],也可根據(jù)干擾和目標(biāo)波形差異進(jìn)行抑制[29]。④ 極化域抗干擾:通過(guò)研究真、假目標(biāo)的極化信息差異,采用極化分集[30]、極化濾波[31]、聯(lián)合極化-空域[32]、脈沖間極化隨機(jī)捷變[33]等進(jìn)行對(duì)抗。⑤ 分布式抗干擾:通過(guò)組網(wǎng)雷達(dá)協(xié)同抗干擾,基于數(shù)據(jù)級(jí)和信號(hào)級(jí)利用信息融合技術(shù)對(duì)干擾進(jìn)行有效識(shí)別和抑制[34-36]。傳統(tǒng)的抗主瓣欺騙式干擾方法能在一定程度上解決主瓣欺騙式干擾問(wèn)題。然而,面對(duì)主瓣欺騙式干擾密集性強(qiáng)、逼真度高等特點(diǎn),上述方法在實(shí)際中仍存在相當(dāng)?shù)木窒扌?故需更多的雷達(dá)系統(tǒng)自由度進(jìn)行主瓣欺騙式干擾對(duì)抗。

雷達(dá)體制不斷革新,如從早期的機(jī)械掃描雷達(dá)發(fā)展至相控陣?yán)走_(dá),以及近年來(lái)備受關(guān)注的多輸入多輸出(multiple input multiple output, MIMO)雷達(dá),雷達(dá)系統(tǒng)對(duì)目標(biāo)和環(huán)境的信息獲取能力不斷擴(kuò)展,系統(tǒng)可控自由度不斷增加。波形分集陣新體制雷達(dá)是在 MIMO 雷達(dá)基礎(chǔ)上發(fā)展起來(lái)的,通過(guò)在發(fā)射陣元之間采用頻率/時(shí)延/相位調(diào)制,并在接收端采用相應(yīng)的信號(hào)處理技術(shù),達(dá)到提升雷達(dá)系統(tǒng)性能、擴(kuò)展雷達(dá)系統(tǒng)能力的目的[37-39]。其中,以美國(guó)空軍實(shí)驗(yàn)室的高級(jí)研究員 Paul Antonik 等人提出的頻率分集陣(frequency diverse array,FDA)雷達(dá)為代表[40],其各天線發(fā)射信號(hào)的頻率不同,相比于傳統(tǒng)相控陣發(fā)射方向圖僅依賴于角度,FDA 雷達(dá)形成的發(fā)射方向圖不僅是角度的函數(shù),而且是時(shí)間/距離的函數(shù),這一點(diǎn)是陣列發(fā)展歷程中革命性的變化[41-43]。然而,FDA 在固定時(shí)刻的發(fā)射方向圖呈現(xiàn)距離-角度耦合特性,且發(fā)射方向圖非穩(wěn)態(tài),因此,需要聯(lián)合 MIMO技術(shù),通過(guò)分離發(fā)射波形,從而進(jìn)一步提取出目標(biāo)的距離維信息。FDA-MIMO雷達(dá)相比于傳統(tǒng)的 MIMO 雷達(dá)能夠獲取額外的距離維自由度,針對(duì)這一特性,諸多文獻(xiàn)探索了 FDA-MIMO雷達(dá)在目標(biāo)距離-角度聯(lián)合參數(shù)估計(jì)[44-45]、雜波抑制[46-48]、合成孔徑雷達(dá)解模糊[49-51]等中的應(yīng)用。近年來(lái),信號(hào)處理領(lǐng)域重要期刊如《IEEE Journal of Selected Topics in Signal Processing》《雷達(dá)學(xué)報(bào)》等組織了FDA相關(guān)??W(xué)者們就FDA技術(shù)進(jìn)行了綜述性總結(jié)[52-54],并提出了空距頻聚焦方法[55]、雙基地距離模糊雜波抑制方法[56]等。實(shí)際上,通過(guò)發(fā)射陣元間頻率調(diào)制,等價(jià)于引入了一個(gè)脈間時(shí)變、空間通道數(shù)依賴的相位項(xiàng),而獲取 FDA-MIMO的距離維自由度是在空間上對(duì)該相位進(jìn)行合理處理?;诖?文獻(xiàn)[57]提出了陣元-脈沖編碼-MIMO(element phase code MIMO, EPC-MIMO)雷達(dá)的概念,能夠區(qū)分來(lái)自不同距離模糊區(qū)間的回波信號(hào)。

波形分集陣?yán)走_(dá)通過(guò)陣元間調(diào)制主動(dòng)改變發(fā)射波形,制造目標(biāo)和干擾間的多維參數(shù)差異,從而使系統(tǒng)具有更高的可用自由度,為從發(fā)射端解決主瓣欺騙式干擾問(wèn)題提供了一條有效思路,具有了對(duì)抗主瓣欺騙式干擾更高的系統(tǒng)自由度。目前,已有大量文獻(xiàn)開(kāi)展FDA-MIMO體制下的抗主瓣欺騙式干擾研究[58-80]。其中,文獻(xiàn)[58]在FDA-MIMO雷達(dá)體制下利用自適應(yīng)二維波束形成技術(shù)抑制主瓣欺騙式干擾。然而,該模型忽略了干擾信號(hào)在干擾機(jī)中的存儲(chǔ)調(diào)制時(shí)間,認(rèn)為由同一干擾機(jī)產(chǎn)生的所有假目標(biāo)具有相同的發(fā)射導(dǎo)向矢量,僅是實(shí)際中一種特例。對(duì)此,文獻(xiàn)[59-60]對(duì)假目標(biāo)時(shí)延量進(jìn)行分析并做了模型修正。文獻(xiàn)[61-62]利用子空間投影法挑選樣本,實(shí)現(xiàn)有效的主瓣欺騙式干擾對(duì)抗。文獻(xiàn)[63]基于對(duì)數(shù)形式的頻偏,提出了一種自適應(yīng)波束形成方法,克服了FDA-MIMO雷達(dá)方向圖中存在的距離周期性。文獻(xiàn)[64]采用了一種“低秩+低秩+稀疏”分解提取低秩階期望信號(hào)并抑制干擾。文獻(xiàn)[65]對(duì)均勻線陣和均勻面陣的頻控陣拖曳式干擾抑制方法進(jìn)行了研究。文獻(xiàn)[66]分析了上下變頻、匹配濾波、信號(hào)混合等因素的影響,指出FDA能夠通過(guò)提取相鄰陣元相位差信息,實(shí)現(xiàn)對(duì)假目標(biāo)的鑒別與抑制。文獻(xiàn)[67]將FDA發(fā)射陣分成兩個(gè)子陣列,采用兩種不同的頻率步進(jìn)量,聯(lián)合距離、角度和多普勒域在子陣列級(jí)實(shí)現(xiàn)假目標(biāo)的鑒別與對(duì)抗。文獻(xiàn)[68]在先驗(yàn)信息未知情況下,通過(guò)基于最大信噪比(signal to noise ratio, SNR)的盲源分離算法構(gòu)建分離矩陣,將目標(biāo)和干擾分離在不同通道,有效抑制雷達(dá)主瓣欺騙式干擾。文獻(xiàn)[69]提出了一種基于特征向量剔除法(feature vector elimination, FVE)的干擾抑制方法。文獻(xiàn)[70]提出一種FDA極化MIMO(簡(jiǎn)稱為FDA-PMIMO)雷達(dá)新技術(shù),可通過(guò)對(duì)發(fā)射極化和發(fā)射頻率步進(jìn)間隔進(jìn)行聯(lián)合優(yōu)化設(shè)計(jì)來(lái)進(jìn)一步提高其抗干擾性能。文獻(xiàn)[71]提出一種基于直接數(shù)據(jù)域方法挑選訓(xùn)練樣本。文獻(xiàn)[72]提出一種基于協(xié)方差矩陣重構(gòu)的穩(wěn)健抗干擾方法。文獻(xiàn)[73]提出一種欺騙式干擾背景下的自適應(yīng)目標(biāo)檢測(cè)方法。文獻(xiàn)[74]針對(duì)非均勻間隔的FDA雷達(dá),提出了一種基于相位中心的認(rèn)知自適應(yīng)抗干擾方法。文獻(xiàn)[75]基于模擬退火算法優(yōu)化各陣元間頻率步進(jìn)量進(jìn)行自適應(yīng)抗干擾。盡管上述方法能對(duì)抗主瓣欺騙式干擾,但是假設(shè)的理想正交波形在實(shí)際工程中具有一定的困難。對(duì)此,筆者所在團(tuán)隊(duì)對(duì)FDA-MIMO非理想正交波形下的抗干擾方法開(kāi)展了研究[76-79],并進(jìn)一步考慮了實(shí)際誤差影響下的穩(wěn)健抗干擾方法[78]。此外,在文獻(xiàn)[80]中提出了基于EPC-MIMO體制的穩(wěn)健抗干擾方法。

本文對(duì)波形分集陣列雷達(dá)抗干擾原理以及性能進(jìn)行了分析。首先介紹了波形分集陣列雷達(dá)模型,重點(diǎn)對(duì)FDA-MIMO以及EPC-MIMO雷達(dá)收發(fā)信號(hào)模型進(jìn)行介紹;然后對(duì)波形分集陣列雷達(dá)抗干擾原理進(jìn)行詳細(xì)闡述;并從自適應(yīng)波束形成器設(shè)計(jì)、空間投影、方向圖設(shè)計(jì)3個(gè)方面介紹和分析了波形分集陣?yán)走_(dá)空域抗主瓣欺騙式干擾的方法;其次簡(jiǎn)要介紹了波形分集陣原理樣機(jī)的研制以及抗干擾實(shí)驗(yàn)驗(yàn)證。最后對(duì)波形分集陣列雷達(dá)抗干擾技術(shù)的難點(diǎn)與未來(lái)研究趨勢(shì)進(jìn)行了展望。

1 波形分集陣列雷達(dá)模型

1.1 波形分集陣?yán)走_(dá)收發(fā)信號(hào)模型

假設(shè)一個(gè)共置MIMO 雷達(dá)具有M個(gè)發(fā)射陣元、N個(gè)接收陣元,其收發(fā)陣列均為以d為間隔的半波長(zhǎng)等距線陣。FDA-MIMO在發(fā)射信號(hào)各陣元間引入一個(gè)很小的頻率步進(jìn)量(遠(yuǎn)小于載頻)后Δf,其第m(m=1,2,…,M)個(gè)陣元的發(fā)射信號(hào)為

(1)

式中,fm=f0+(m-1)Δf為第m個(gè)陣元載頻,f0為參考載頻(通常選為第一個(gè)發(fā)射陣元);φm(t)為第m個(gè)陣元發(fā)射波形。

圖1 FDA-MIMO雷達(dá)接收匹配濾波處理流程Fig.1 Receiving matched filtering processing procedure of FDA-MIMO radar

最終,來(lái)自N個(gè)接收陣元的經(jīng)過(guò)匹配濾波之后的和信號(hào)可以組成一個(gè)MN×1維的列矢量[78]:

ys=β1r⊙[b(θ0)?a(R0,θ0)]

(2)

其中,β1表示目標(biāo)的復(fù)回波幅度;r∈CMN×1為匹配濾波輸出矢量;R0和θ0分別表示目標(biāo)的距離和角度;a(R0,θ0)∈CM×1和b(θ0)∈CM×1分別為發(fā)射、接收導(dǎo)向矢量,其具體表達(dá)式為

(3)

(4)

其中,λ0表示波長(zhǎng)??梢?jiàn),FDA-MIMO雷達(dá)接收導(dǎo)向矢量?jī)H與角度有關(guān),而發(fā)射導(dǎo)向矢量卻呈現(xiàn)距離和角度的兩維相關(guān)性。不同于傳統(tǒng)意義上的回波時(shí)延測(cè)量的目標(biāo)距離,FDA-MIMO 雷達(dá)經(jīng)過(guò)接收信號(hào)處理后得到獨(dú)立的距離維可控自由度。若Δf=0,則FDA-MIMO雷達(dá)退化為傳統(tǒng)的MIMO雷達(dá),其發(fā)射、接收導(dǎo)向矢量均只與角度有關(guān)。

Ys=α0{r⊙[b(θ0)?a(γs,θ0)]}fT(fds)

(5)

其中,r∈CMN為匹配濾波輸出矢量;f(vs)∈CK表示Doppler矢量;a(γs,θ0)∈CM和b(θ0)∈CN分別表示發(fā)射、接收導(dǎo)向矢量,具體分別寫作:

a(γs,θ0)=[1,e-j2πγps,…,e-j2πγps(M-1)]T⊙

(6)

(7)

圖2 EPC-MIMO雷達(dá)接收匹配濾波處理流程Fig.2 Receiving matched filtering processing procedure of EPC-MIMO radar

1.2 假目標(biāo)產(chǎn)生模型

目標(biāo)告警系統(tǒng)發(fā)現(xiàn)自己被雷達(dá)跟蹤以后,啟動(dòng)干擾機(jī)截獲雷達(dá)信號(hào),并進(jìn)行調(diào)制轉(zhuǎn)發(fā)。通過(guò)控制調(diào)制延遲時(shí)間,干擾機(jī)產(chǎn)生假目標(biāo)包括停拖期、拖引期和休止期3個(gè)過(guò)程,最終雷達(dá)鎖定跟蹤假目標(biāo)而丟失真實(shí)目標(biāo)。對(duì)于距離拖引式干擾,跟蹤制導(dǎo)雷達(dá)采取脈沖前沿跟蹤的方式提取目標(biāo)的距離維信息。盡管真、假目標(biāo)的距離在一個(gè)無(wú)模糊區(qū)間內(nèi)可能存在超前和滯后的關(guān)系,假目標(biāo)在的實(shí)際等效距離一定大于目標(biāo)(干擾機(jī))的真實(shí)距離[61-62]。若假目標(biāo)產(chǎn)生器(false target generator, FTG)位于角度θj和距離Rj的位置,根據(jù)文獻(xiàn)[59],FTG收到的信號(hào)經(jīng)過(guò)下調(diào)頻后可以表示為

(8)

其中,τj=2Rj/c,c為光速;τj,m表示由第m個(gè)陣元的波程差造成的時(shí)延。值得注意的是,對(duì)于主瓣欺騙式干擾θj=θ0,若考慮自衛(wèi)式干擾,有Rj=R0。FTG經(jīng)過(guò)延時(shí)τjam后產(chǎn)生的第q個(gè)(q=1,2,…,Q)假目標(biāo)位于Rq=Rj+Δr且Δr=cτjam/2,則調(diào)制后的信號(hào)可以表示為

(9)

其中,Ajam表示FTG中的幅度調(diào)制;φτ表示由時(shí)延帶來(lái)的相位變化。以FDA-MIMO雷達(dá)為例,由FTG產(chǎn)生的假目標(biāo)信號(hào)被雷達(dá)系統(tǒng)接收并經(jīng)過(guò)處理(包括混頻、匹配濾波等)后,第n個(gè)接收陣元的接收的由第m個(gè)陣元發(fā)射信號(hào)可以表示為

(10)

其中,ξq表示第q個(gè)假目標(biāo)的復(fù)回波幅度;Cn,m表示第n個(gè)接收通道對(duì)第m個(gè)波形的匹配濾波輸出。

2 波形分集陣列雷達(dá)抗干擾原理

根據(jù)FDA-MIMO雷達(dá)中真實(shí)目標(biāo)的發(fā)射導(dǎo)向矢量和接收導(dǎo)向矢量的表達(dá)式,則真實(shí)目標(biāo)和第q個(gè)假目標(biāo)發(fā)射空間頻率表示如下:

(11)

(12)

據(jù)此,真、假目標(biāo)發(fā)射空間頻率中包含有不同的距離信息,因此可在發(fā)射-接收二維空間頻率域內(nèi)區(qū)分。實(shí)際中,需對(duì)接收數(shù)據(jù)進(jìn)行逐個(gè)距離門的距離補(bǔ)償從而消除距離依賴性,則補(bǔ)償后:

(13)

(14)

其中,Ru=c/2fr表示最大無(wú)模糊距離,fr為PRF;p為第q個(gè)假目標(biāo)的相對(duì)于真實(shí)目標(biāo)的延遲脈沖數(shù)??梢?jiàn),補(bǔ)償后,對(duì)于同一個(gè)發(fā)射脈沖內(nèi)任意距離門的目標(biāo),其發(fā)射空間頻率相同。如圖3所示,補(bǔ)償后的真實(shí)目標(biāo)的發(fā)射與接收空間頻率相等,故真實(shí)目標(biāo)位于發(fā)射-接收二維空域平面的對(duì)角線[62],而假目標(biāo)在發(fā)射-接收二維空間頻率域可任意分布,據(jù)此可對(duì)真、假目標(biāo)進(jìn)行鑒別。同理,在EPC-MIMO雷達(dá)中,補(bǔ)償相鄰陣元間相位差e-j2πγs后,真實(shí)目標(biāo)發(fā)射空間頻率與接收空間頻率相等,而第q個(gè)假目標(biāo)發(fā)射空間頻率可表示為

(15)

圖3 真實(shí)目標(biāo)和假目標(biāo)的分布示意圖Fig.3 Distribution of the true and false targets

可見(jiàn),當(dāng)給定系統(tǒng)參數(shù)(頻率步進(jìn)量或編碼系數(shù))時(shí),發(fā)射空間頻率僅與脈沖延遲數(shù)有關(guān)。因此,在發(fā)射空間頻率上可區(qū)分對(duì)應(yīng)于不同的距離模糊區(qū)間目標(biāo)和欺騙式的干擾。在FDA-MIMO雷達(dá)中,假目標(biāo)信號(hào)的距離模糊重?cái)?shù)盡管不同于真實(shí)目標(biāo),但是模糊重?cái)?shù)是個(gè)常量,因此不同脈沖對(duì)應(yīng)的目標(biāo)空間譜保持不變。而EPC-MIMO雷達(dá)發(fā)射陣元編碼在每個(gè)脈沖都變化,等效的發(fā)射角度在不同發(fā)射脈沖間變化。處于不同距離模糊的假目標(biāo)信號(hào)由于對(duì)應(yīng)的EPC編碼不同,可以被辨別出來(lái)。

圖4給出了真、假目標(biāo)在發(fā)射空域-接收空域-脈沖三維(3-D)空間分布。由于假目標(biāo)位于主瓣,則真、假目標(biāo)具有相同的接收空間頻率,因此其分布在發(fā)射空間的頻率-脈沖平面內(nèi)(如中間子圖所示)。最左側(cè)子圖給出了發(fā)射-接收空間域的真、假目標(biāo)分布情況,由于真、假目標(biāo)具有相同的接收空間頻率,則其分布在一條給定接收空間頻率的水平線上。最右側(cè)子圖表示了固定接收空間頻率的發(fā)射空間頻率-脈沖維剖面,其中真、假目標(biāo)以γ為斜率呈斜線分布。值得注意的是,如圖5所示,將假目標(biāo)在雷達(dá)回波中的位置分為I 區(qū)和II區(qū)。其中I 區(qū)的假目標(biāo)經(jīng)干擾機(jī)快速轉(zhuǎn)發(fā),與真實(shí)目標(biāo)位于雷達(dá)同一個(gè)接收脈沖重復(fù)時(shí)間(pulse repetitive time,PRT)內(nèi);II 區(qū)假目標(biāo)經(jīng)過(guò)較大存儲(chǔ)延遲,滯后于真實(shí)目標(biāo)至少一個(gè)PRT。波形分集陣列雷達(dá)空域抗干擾僅考慮II區(qū)假目標(biāo)。

圖4 真、假目標(biāo)在發(fā)射空域-接收空域-脈沖三維空間分布Fig.4 Distributions of the true and false targets in the transmit-receive-pulse space three-dimensional space

圖5 假目標(biāo)示意圖Fig.5 Illustration of false target

3 波形分集陣列雷達(dá)空域抗干擾方法

現(xiàn)有波形分集陣?yán)走_(dá)多從空域波束形成角度設(shè)計(jì)抗干擾方法,對(duì)此,從自適應(yīng)波束形成器設(shè)計(jì)、空間投影、方向圖設(shè)計(jì)3個(gè)方面對(duì)現(xiàn)有波形分集陣?yán)走_(dá)抗干擾方法進(jìn)行總結(jié)與分析。

3.1 自適應(yīng)波束形成器設(shè)計(jì)類抗干擾方法

在FDA-MIMO雷達(dá)中,大量文獻(xiàn)構(gòu)造基于最小方差無(wú)失真響應(yīng)(minimum variance distortionless response, MVDR)的自適應(yīng)波束形成器,以在距離-角度二維空域自適應(yīng)地形成零點(diǎn)[58-60,63,71-73,77]:

s.t.wHu(R0,θ0)=1

(16)

(17)

式中,X表示來(lái)自第l(l=1,2,…,L)個(gè)距離單元、第k個(gè)脈沖的雷達(dá)信號(hào)樣本;S,J,N分別代表真實(shí)目標(biāo)、假目標(biāo)以及噪聲分量。進(jìn)一步,根據(jù)文獻(xiàn)[62],結(jié)合切諾夫界,可得多假設(shè)檢驗(yàn)的誤警率為

(18)

自適應(yīng)波束形成抗干擾的關(guān)鍵在于干擾加噪聲協(xié)方差矩陣的構(gòu)建,而考慮到現(xiàn)實(shí)中難以獲取足夠多的獨(dú)立同分布樣本來(lái)估計(jì)協(xié)方差矩陣。對(duì)此,大量文獻(xiàn)對(duì)協(xié)方差矩陣的估計(jì)方法開(kāi)展了研究[63,69,71-73,75,77]。圖6給出了幾種具有代表性的自適應(yīng)二維匹配濾波器響應(yīng)圖。其中,文獻(xiàn)[77]提出了一種魯棒的非一致樣本檢測(cè)(no-consistent sample detection, NSD)方法,具體分為兩步:① 挑選包括信號(hào)、包括/不包括干擾的非一致樣本(即滿足H1、H2/H3假設(shè));② 根據(jù)目標(biāo)信號(hào)的信息,利用空間平滑法剔除含有目標(biāo)信號(hào)的樣本。圖7給出了距離-角度二維自適應(yīng)匹配濾波器輸出結(jié)果。通過(guò)魯棒的NSD 方法,來(lái)自主瓣方向的假目標(biāo)由于距離維的失配可以被有效地抑制。若未剔除真實(shí)目標(biāo),會(huì)對(duì)匹配濾波結(jié)果造成影響,導(dǎo)致了假目標(biāo)輸出功率的上升。由圖8給出的信干噪比(signal to interference plus noise ratio, SINR)曲線可見(jiàn),經(jīng)過(guò)樣本挑選并且剔除目標(biāo)信號(hào)后,避免了目標(biāo)相消的問(wèn)題,獲得了較好的輸出SINR,且接近理想上界。

圖6 FDA-MIMO自適應(yīng)匹配濾波器響應(yīng)Fig.6 FDA-MIMO adaptive matched filter reponse

圖7 FDA-MIMO雷達(dá)距離-角度自適應(yīng)匹配濾波器輸出結(jié)果Fig.7 Output results of range angle adaptive matched filter for FDA-MIMO radar

圖8 不同樣本挑選策略下的輸出SINR與輸入SNR關(guān)系Fig.8 Relationship between output SINR and input SNR under different sample selection strategies

由于接收數(shù)據(jù)中需要濾除目標(biāo)信號(hào)后才可用來(lái)構(gòu)建協(xié)方差矩陣,文獻(xiàn)[69]提出了一種利用信號(hào)子空間校正失配導(dǎo)向矢量的思路。由于期望的目標(biāo)導(dǎo)向矢量與實(shí)際存在偏差,需要先對(duì)其進(jìn)行約束;隨后,對(duì)包含目標(biāo)信號(hào)的協(xié)方差矩陣進(jìn)行特征分解[69]。在目標(biāo)導(dǎo)向矢量的范圍內(nèi)[69],由于目標(biāo)導(dǎo)向矢量與信號(hào)子空間相關(guān)性最強(qiáng),基于這一特性,在導(dǎo)向矢量中對(duì)目標(biāo)進(jìn)行尋優(yōu)校正,選擇與信號(hào)子空間中特征向量相關(guān)性最強(qiáng)的為校正后的目標(biāo)導(dǎo)向矢量[69]:

(19)

最后,在接收協(xié)方差矩陣中剔除該特征向量后構(gòu)建自適應(yīng)權(quán)矢量。此外,文獻(xiàn)[71]提出一種基于直接數(shù)據(jù)域的方法來(lái)避免訓(xùn)練樣本中的目標(biāo)相消現(xiàn)象。通過(guò)在發(fā)射-接收二維空域構(gòu)建兩個(gè)復(fù)系數(shù):

(20)

(21)

則可通過(guò)空域?yàn)V波的方式獲得三組訓(xùn)練數(shù)據(jù)矩陣:

(22)

(23)

(24)

并由此構(gòu)建聯(lián)合發(fā)射接收空域協(xié)方差矩陣:

vec(X2)H+vec(X3)vec(X3)H}

(25)

如圖9所示,輸出SINR結(jié)果表明,該方法對(duì)訓(xùn)練數(shù)據(jù)中存在的目標(biāo)信號(hào)具有較強(qiáng)的魯棒性。由于目標(biāo)信號(hào)通過(guò)預(yù)濾波方式剔除,避免了目標(biāo)相消。該方法以微小的孔徑損失為代價(jià),使樣本數(shù)增加了3倍。最終得到的自適應(yīng)波束形成器在能夠有效抑制假目標(biāo)。

圖9 不同方法的輸出SINR與輸入SNR關(guān)系Fig.9 Relationship between output SINR and input SNR under different methods

此外,文獻(xiàn)[72]利用Capon譜重構(gòu)信號(hào)協(xié)方差矩陣,并對(duì)目標(biāo)可能出現(xiàn)的區(qū)域進(jìn)行積分。通過(guò)對(duì)信號(hào)協(xié)方差進(jìn)行特征分解,選擇與最大特征值相關(guān)的特征向量來(lái)估計(jì)信號(hào)的導(dǎo)向矢量。最后,利用估計(jì)的目標(biāo)導(dǎo)向矢量計(jì)算自適應(yīng)抗干擾權(quán)值為

(26)

上述方法均基于頻偏線性遞增的FDA-MIMO雷達(dá),文獻(xiàn)[63]和文獻(xiàn)[75]分別基于對(duì)數(shù)遞增和模擬退火算法優(yōu)化頻偏,此類方法消除了自適應(yīng)波束形成聯(lián)合方向圖的距離周期性,避免了柵瓣附近抗干擾能力下降的問(wèn)題,可以有效抑制主瓣欺騙干擾。

3.2 空間投影類抗干擾方法

圖10 發(fā)射-接收二維波束形成器響應(yīng)Fig.10 Responses of transmit-receive two-dimensional beamformers

此外,文獻(xiàn)[64]借助于GoDec方法[81],在FDA-MIMO雷達(dá)中提出了一種交替最小化的兩步GoDec方法來(lái)抑制欺騙式干擾。由于Xs與Xi具有低秩特性,以及Xe的稀疏特性,對(duì)于包含目標(biāo)信號(hào)、干擾信號(hào)、壓制式干擾與高斯噪聲的總快拍數(shù)據(jù)Y=Xs+Xi+Xe+Xn,Y可以被分解為3個(gè)獨(dú)立的矩陣:Ls、Li與Se,其中Ls是對(duì)信號(hào)矩陣Xs的估計(jì)。為了從Y中找到矩陣Ls,可使用Frobenius范數(shù)最小化分解誤差來(lái)構(gòu)建近似分解問(wèn)題:

(27)

其中,rs與ri分別表示Ls與Li的秩;δ表示Se的基數(shù)。對(duì)此,提出兩步GoDec算法來(lái)更新Li、Se與Ls,最終得到Ls的估計(jì)值,并實(shí)現(xiàn)了假目標(biāo)的抑制。文獻(xiàn)[64]進(jìn)一步對(duì)目標(biāo)角度和距離進(jìn)行估計(jì)。

3.3 方向圖設(shè)計(jì)類抗干擾方法

相比與依賴訓(xùn)練數(shù)據(jù)的抗干擾方法,方向圖設(shè)計(jì)法能夠綜合形成期望的空域方向圖,該方法無(wú)需訓(xùn)練樣本,易于實(shí)現(xiàn)且計(jì)算復(fù)雜度更低。實(shí)際上,此類方法需要滿足的前提是假目標(biāo)滯后于真實(shí)目標(biāo)至少一個(gè)脈沖。根據(jù)FDA-MIMO雷達(dá)中補(bǔ)償后的真、假目標(biāo)發(fā)射空間頻率表達(dá)式,真、假目標(biāo)的發(fā)射空間頻率之差為

(28)

圖11給出了聯(lián)合發(fā)射-接收空間頻率域真實(shí)目標(biāo)所對(duì)應(yīng)的等效發(fā)射波束方向圖。由于假目標(biāo)與真實(shí)目標(biāo)存在距離模糊數(shù)的差異,其發(fā)射空間頻率與真實(shí)目標(biāo)存在p2Δf/cRu的頻率差。利用指向目標(biāo)方向的權(quán)矢量疊加接收信號(hào)后,假目標(biāo)由于位于發(fā)射方向圖的零點(diǎn)而被抑制。

圖11 FDA-MIMO雷達(dá)方向圖設(shè)計(jì)抗干擾原理Fig.11 Anti-jamming principle of FDA-MIMO radar pattern design

圖12給出了EPC-MIMO雷達(dá)不同目標(biāo)對(duì)應(yīng)的等效發(fā)射方向圖和接收雙程方向圖。由圖12(a)可見(jiàn),由于模糊相位的存在,來(lái)自不同距離模糊區(qū)間的真、假目標(biāo)回波等效具有不同的發(fā)射波束指向,進(jìn)而真、假目標(biāo)可以在發(fā)射空域上進(jìn)行區(qū)別。由于等效發(fā)射方向圖中,零點(diǎn)之間的頻率間隔為1/M,當(dāng)γ=3/M,延遲了一個(gè)脈沖的假目標(biāo)的等效發(fā)射方向圖主瓣恰好對(duì)準(zhǔn)第3個(gè)零點(diǎn),從圖12(b)可見(jiàn),在接收雙程方向圖中,假目標(biāo)對(duì)應(yīng)的主瓣方向能量下降,且越遠(yuǎn)離真實(shí)目標(biāo)的主瓣衰減越大。

圖12 EPC-MIMO雷達(dá)抗干擾示意圖(γ=3/M)Fig.12 Anti-jamming schematic diagram of EPC-MIMO radar(γ=3/M)

波形分集陣列雷達(dá)基于方向圖設(shè)計(jì)法通過(guò)將干擾對(duì)準(zhǔn)方向圖零點(diǎn)來(lái)進(jìn)行抑制,而這一點(diǎn)需要通過(guò)對(duì)系統(tǒng)參數(shù)進(jìn)行合理設(shè)計(jì)來(lái)實(shí)現(xiàn)。對(duì)于具有M個(gè)陣元的陣列,其發(fā)射方向圖有M-1個(gè)零點(diǎn),延遲脈沖數(shù)對(duì)應(yīng)于1到M-1的假目標(biāo)位于方向圖的第1至第M-1個(gè)零點(diǎn)。因此,由真、假目標(biāo)的發(fā)射空間頻率差,對(duì)于延遲了p個(gè)脈沖的假目標(biāo),Δf和編碼系數(shù)可設(shè)計(jì)為

(29)

需要說(shuō)明的是,當(dāng)欺騙式假目標(biāo)的轉(zhuǎn)發(fā)延遲為整數(shù)倍的M乘以脈沖重復(fù)周期(pulse repetition time, PRT),即M×PRT,2M×PRT,…,nM×PRT,在該情形下,由于真實(shí)目標(biāo)信號(hào)和欺騙式干擾信號(hào)的發(fā)射空域頻率沒(méi)有差異,抗干擾失效。

然而實(shí)際中,目標(biāo)參數(shù)與理論模型存在偏差,則假目標(biāo)可能偏離其在發(fā)射波束形成方向圖上理論的零點(diǎn),真實(shí)目標(biāo)也可能偏離其理論主瓣,導(dǎo)致抗干擾性能變差。其主要影響因素包括距離量化誤差、頻率步進(jìn)量誤差、角度偏差。此外,由于目標(biāo)是運(yùn)動(dòng)的,因此不同脈沖對(duì)應(yīng)的目標(biāo)距離有變化,補(bǔ)償后其等效發(fā)射頻率與理論值存在偏差;陣列幅度相位誤差也會(huì)影響波束形成方向圖的零點(diǎn)位置和深度。針對(duì)上述問(wèn)題,需要對(duì)非自適應(yīng)波束形成方向圖進(jìn)行設(shè)計(jì),展寬方向圖零點(diǎn)和主瓣,采用穩(wěn)健波束形成算法來(lái)應(yīng)對(duì)目標(biāo)參數(shù)與理論模型偏差的影響。

文獻(xiàn)[78]通過(guò)在零點(diǎn)周圍施加具有特定功率值的干擾來(lái)展寬方向圖零點(diǎn),且該功率值的閉式解可以根據(jù)期望的零點(diǎn)方向圖響應(yīng)求得。利用虛擬干擾來(lái)構(gòu)建協(xié)方差矩陣,并通過(guò)迭代算法優(yōu)化權(quán)矢量直到設(shè)計(jì)的方向圖具有理想的零陷(見(jiàn)圖13),并提出了一種預(yù)設(shè)寬零點(diǎn)波束形成(preset broadened nulling beamformer, PBN-BF)算法。圖14給出了采用PBN-BF方法的輸出SINR變化曲線,采用所提的PBN-BF 方法,即使假目標(biāo)偏離了理論的零點(diǎn),仍可處于方向圖的寬零陷而被充分抑制,因此輸出SINR 性能增強(qiáng),且可以在誤差變化時(shí)維持較為穩(wěn)定的輸出SINR。此外,文獻(xiàn)[80]針對(duì)EPC-MIMO雷達(dá)抗干擾實(shí)際情況中真、假目標(biāo)相對(duì)于理論模型同時(shí)存在的角度偏差的問(wèn)題,提出了一種預(yù)設(shè)方向圖綜合(preset beampattern synthesis, PBPS)方法,通過(guò)對(duì)收發(fā)二維方向圖的權(quán)矢量進(jìn)行正交分解,實(shí)現(xiàn)了具有平頂主瓣、寬零點(diǎn)、低副瓣的期望形狀方向圖,在充分抑制假目標(biāo)的同時(shí)擴(kuò)大對(duì)目標(biāo)的接收范圍,其合成后的方向圖如圖15所示。

圖13 展寬零點(diǎn)后的收發(fā)二維方向圖Fig.13 Transmit-receive two dimensional beampattern after null broadening

圖14 輸出SINR變化曲線Fig.14 Output SINR curves

圖15 PBPS方法合成收發(fā)二維方向圖Fig.15 Synthesized transmit-receive two-dimensional beampattern using the PBPS method

圖16給出了角度偏差存在下的抗干擾結(jié)果。由圖16(a)匹配濾波輸出結(jié)果可見(jiàn),采用PBPS方法后假目標(biāo)仍能位于寬零區(qū)而被抑制,且真實(shí)目標(biāo)仍處于寬接收范圍保證其最大輸出。由圖16(b)可知,相比于原始的方法,PBPS的輸出SINR 提高了15 dB。

圖16 PBPS抗干擾結(jié)果Fig.16 Anti-jamming results of PBPS

4 波形分集陣?yán)走_(dá)試驗(yàn)系統(tǒng)研制與原理性試驗(yàn)驗(yàn)證

除了對(duì)波形分集陣?yán)走_(dá)的理論研究,國(guó)內(nèi)外十分重視試驗(yàn)系統(tǒng)的研制。早在2009年,美國(guó)空軍研究實(shí)驗(yàn)室的Antonik等人研制出了FDA雷達(dá)發(fā)射和接收模塊如圖17所示,并進(jìn)行了外場(chǎng)測(cè)試,驗(yàn)證了頻控陣?yán)走_(dá)的發(fā)射波束傳播特性[82]。同年,英國(guó)的Huang等人利用電磁仿真軟件對(duì)頻率分集線性陣列天線的輻射特性進(jìn)行驗(yàn)證,并設(shè)計(jì)了一個(gè)具有4個(gè)陣元的FDA雷達(dá)天線[83,53](見(jiàn)圖18)。來(lái)自土耳其中東理工大學(xué)的Eker、Demir等人于2014年提出了一種有效的FDA天線系統(tǒng)設(shè)計(jì)方法,并基于線性調(diào)頻波信號(hào)研制了FDA雷達(dá)原理樣機(jī)(見(jiàn)圖19),驗(yàn)證其波束形成方向圖的距離和時(shí)間相關(guān)性[84,53]。國(guó)內(nèi)學(xué)者也于近幾年開(kāi)展了FDA系統(tǒng)研制工作。如圖20所示,2017年,來(lái)自電子科技大學(xué)的王文欽教授團(tuán)隊(duì)研制了FDA多通道信號(hào)源和天線[53],并開(kāi)發(fā)了一套“頻控陣?yán)走_(dá)仿真與處理系統(tǒng)”軟件平臺(tái),用來(lái)對(duì)FDA雷達(dá)進(jìn)行系統(tǒng)仿真與數(shù)據(jù)處理[53]。

圖17 Antonik等人研制的FDA雷達(dá)發(fā)射和接收模塊Fig.17 FDA radar transmiting and receiveing elements developed by Antonik etc.

圖18 Huang等人研制的FDA雷達(dá)Fig.18 FDA radar developed by Huang etc.

圖19 Eker等人研制的FDA雷達(dá)Fig.19 FDA radar developed by Eker etc.

圖20 王文欽研究團(tuán)隊(duì)研制的FDA雷達(dá)Fig.20 FDA radar developed by Wang Wenqin’s team

此外,筆者所在課題組較早開(kāi)展了波形分集陣試驗(yàn)系統(tǒng)的研制,并率先進(jìn)行雷達(dá)主瓣干擾抑制的初步實(shí)驗(yàn)驗(yàn)證。FDA雷達(dá)試驗(yàn)系統(tǒng)包括天線與射頻組件、信號(hào)處理機(jī)、數(shù)據(jù)存儲(chǔ)器、雷達(dá)系統(tǒng)顯控等分系統(tǒng)。雷達(dá)錄取真實(shí)運(yùn)動(dòng)目標(biāo)與地物雜波的回波加干擾信號(hào)的實(shí)測(cè)數(shù)據(jù),如圖21所示,實(shí)測(cè)數(shù)據(jù)經(jīng)處理后,來(lái)自主瓣的欺騙式假目標(biāo)干擾能夠抑制,即實(shí)驗(yàn)結(jié)果與理論研究一致。

圖21 筆者所在團(tuán)隊(duì)研制的FDA雷達(dá)系統(tǒng)Fig.21 FDA radar system designed by the author’s research team

5 研究趨勢(shì)

波形分集陣列雷達(dá)增加了雷達(dá)系統(tǒng)發(fā)射端的空間、頻率、時(shí)間等信息和資源,優(yōu)化了雷達(dá)在檢測(cè)、跟蹤、抗干擾等方面的性能,是未來(lái)雷達(dá)發(fā)展的關(guān)鍵技術(shù),具有重要的研究意義。隨著電磁環(huán)境的日趨復(fù)雜,抗干擾的難度不斷增大,仍存在諸多問(wèn)題,需要進(jìn)行更深入的研究。

(1) 波形分集陣發(fā)射波形設(shè)計(jì)與優(yōu)化

由于波形分集陣列雷達(dá)結(jié)合了發(fā)射波形分集和MIMO雷達(dá)雙重優(yōu)勢(shì),需要在接收端對(duì)發(fā)射波形進(jìn)行分離從而提取出距離維的可控自由度。然而,實(shí)際中多通道波形互相關(guān)會(huì)對(duì)雷達(dá)系統(tǒng)性能造成影響??紤]到實(shí)際中的正交波形難以滿足在任意延遲時(shí)間、任意多普勒頻移條件下互相關(guān)為零的條件,波形設(shè)計(jì)與優(yōu)化是一個(gè)理論難題。因此,研究多維(時(shí)間維、頻率維、角度維)空間中的非凸優(yōu)化問(wèn)題具有重要價(jià)值,同時(shí)結(jié)合波形恒模特征等系統(tǒng)約束,需要開(kāi)展波形分集陣?yán)走_(dá)正交波形設(shè)計(jì)。

(2) 復(fù)合干擾樣式多域聯(lián)合對(duì)抗

波形分集陣列雷達(dá)采用波束形成方法抗干擾具有局限性,即只能夠?qū)罐D(zhuǎn)發(fā)大于一個(gè)脈沖的干擾,而如今電子干擾大都能夠完成快速調(diào)制與轉(zhuǎn)發(fā),且假目標(biāo)信號(hào)相干源、偽隨機(jī)分布、目標(biāo)約束不確定性、小樣本、收發(fā)通道誤差等因素均對(duì)抗干擾性能提出挑戰(zhàn)。此外,雷達(dá)系統(tǒng)面臨無(wú)源干擾(箔條、角反組、強(qiáng)雜波等)、壓制式干擾等復(fù)合干擾形式,因此需要聯(lián)合空、時(shí)、頻、極化等多域進(jìn)行對(duì)抗。

(3) 基于環(huán)境感知的認(rèn)知抗干擾

欺騙式干擾與對(duì)抗實(shí)際上是一個(gè)動(dòng)態(tài)博弈的過(guò)程,雙方需要實(shí)時(shí)調(diào)整己方戰(zhàn)略。具備環(huán)境感知能力已成為雷達(dá)系統(tǒng)的發(fā)展趨勢(shì),因此認(rèn)知能力成為未來(lái)電子對(duì)抗中的重要因素。實(shí)際的戰(zhàn)場(chǎng)環(huán)境下,雷達(dá)系統(tǒng)應(yīng)增加干擾環(huán)境感知功能,對(duì)不同類型干擾主動(dòng)認(rèn)知與快速精準(zhǔn)理解,并采取“有的放矢”的抗干擾策略,化被動(dòng)為主動(dòng),實(shí)現(xiàn)環(huán)境感知、精準(zhǔn)鑒別、有效對(duì)抗、反饋調(diào)整的多位一體的多位一體化抗欺騙式干擾能力。

6 結(jié)束語(yǔ)

有源欺騙干擾對(duì)雷達(dá)系統(tǒng)的威力和生存能力都提出了嚴(yán)峻的挑戰(zhàn),干擾和反干擾已經(jīng)成為對(duì)抗雙方必爭(zhēng)的技術(shù)制高點(diǎn)。本文對(duì)波形分集陣列雷達(dá)抗干擾的理論與方法進(jìn)行了總結(jié)和歸納,分析了空域波束形成方法的對(duì)抗準(zhǔn)則和優(yōu)勢(shì),并指出了現(xiàn)有理論方法的局限性,最后指明了波形分集陣列雷達(dá)抗干擾難點(diǎn)與研究趨勢(shì),提供該領(lǐng)域潛在的研究課題與方向。

猜你喜歡
方向信號(hào)方法
2022年組稿方向
2021年組稿方向
信號(hào)
鴨綠江(2021年35期)2021-04-19 12:24:18
2021年組稿方向
完形填空二則
基于FPGA的多功能信號(hào)發(fā)生器的設(shè)計(jì)
電子制作(2018年11期)2018-08-04 03:25:42
用對(duì)方法才能瘦
Coco薇(2016年2期)2016-03-22 02:42:52
基于LabVIEW的力加載信號(hào)采集與PID控制
四大方法 教你不再“坐以待病”!
Coco薇(2015年1期)2015-08-13 02:47:34
捕魚
主站蜘蛛池模板: julia中文字幕久久亚洲| 国产在线98福利播放视频免费 | 天天摸夜夜操| 制服丝袜在线视频香蕉| 欧美日韩中文字幕二区三区| 国产成人亚洲毛片| 免费午夜无码18禁无码影院| 国产精品美女自慰喷水| 亚洲精品卡2卡3卡4卡5卡区| 免费人成在线观看成人片| 国产精品亚洲综合久久小说| 国产地址二永久伊甸园| 有专无码视频| 中文字幕在线不卡视频| 一级爆乳无码av| 欧美中出一区二区| 亚洲成综合人影院在院播放| 欧美第二区| 亚洲精品波多野结衣| 亚洲男人的天堂视频| 国产女人水多毛片18| av一区二区三区在线观看| 国产一区二区影院| 欧美日韩国产一级| 毛片免费高清免费| 国产制服丝袜91在线| 国产成年无码AⅤ片在线| 天堂在线亚洲| 综合网久久| 国产成人一区在线播放| 乱系列中文字幕在线视频| 久草中文网| 欧美国产日韩在线| 亚洲男人天堂2018| 88av在线播放| 免费看美女自慰的网站| 中文字幕乱码二三区免费| 国产精品久久久久鬼色| 国产电话自拍伊人| 一区二区三区成人| 国产系列在线| 三级视频中文字幕| 亚洲国产欧美目韩成人综合| 亚洲成在人线av品善网好看| 久久久久久尹人网香蕉| 国产成人91精品免费网址在线| 欧美一级黄色影院| 中文字幕无码中文字幕有码在线| 欧美天堂在线| 国产精品性| 亚洲欧洲日韩综合| 新SSS无码手机在线观看| 国产一二三区视频| 波多野结衣一级毛片| 亚洲成人免费在线| 精品一区二区三区水蜜桃| 国产在线欧美| 中文字幕亚洲无线码一区女同| 午夜欧美在线| 99久久精品免费看国产免费软件| 国产在线一区二区视频| 日本少妇又色又爽又高潮| 高清无码不卡视频| 欧美黄色a| 亚洲精品欧美日本中文字幕| 97超碰精品成人国产| 亚洲色无码专线精品观看| 亚洲婷婷六月| 在线观看免费黄色网址| 暴力调教一区二区三区| 91精品综合| 中文字幕日韩欧美| 色妞永久免费视频| 国产乱子伦精品视频| 无码免费视频| 色综合综合网| 色婷婷丁香| 又猛又黄又爽无遮挡的视频网站| 天天综合网色| 又猛又黄又爽无遮挡的视频网站| 欧美色图久久| 97av视频在线观看|