• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    合成寬帶脈沖多普勒雷達(dá)

    2016-12-24 07:06:21毛二可范花玉
    關(guān)鍵詞:鋼球雜波測(cè)距

    毛二可, 范花玉

    (1. 北京理工大學(xué)嵌入式實(shí)時(shí)信息處理技術(shù)北京市重點(diǎn)實(shí)驗(yàn)室, 北京 100081;2. 北京理工大學(xué)信息與電子學(xué)院雷達(dá)技術(shù)研究所, 北京 100081)

    ?

    合成寬帶脈沖多普勒雷達(dá)

    毛二可1,2, 范花玉1,2

    (1. 北京理工大學(xué)嵌入式實(shí)時(shí)信息處理技術(shù)北京市重點(diǎn)實(shí)驗(yàn)室, 北京 100081;2. 北京理工大學(xué)信息與電子學(xué)院雷達(dá)技術(shù)研究所, 北京 100081)

    合成寬帶脈沖多普勒(pulse Doppler,PD)雷達(dá)可以同時(shí)實(shí)現(xiàn)距離和速度二維高分辨;具備良好的相參特性,利用相推原理能實(shí)現(xiàn)極高的測(cè)距和測(cè)速精度,可用于復(fù)雜目標(biāo)多部件微動(dòng)測(cè)量;具有較強(qiáng)的抗雜波和抗干擾性能。闡述了合成寬帶PD雷達(dá)體制及其主要信號(hào)形式,介紹了合成寬帶PD雷達(dá)的距離速度二維高分辨實(shí)現(xiàn)原理,總結(jié)了合成寬帶PD雷達(dá)的主要特性,最后通過(guò)一系列實(shí)驗(yàn)驗(yàn)證了合成寬帶PD雷達(dá)的微動(dòng)測(cè)量能力、高精度測(cè)距測(cè)速能力和相參ISAR成像能力。

    合成寬帶信號(hào); 脈沖多普勒雷達(dá); 距離速度二維高分辨; 抗雜波; 抗干擾; 微動(dòng)測(cè)量

    0 引 言

    20世紀(jì)60年代,為了解決強(qiáng)雜波下發(fā)現(xiàn)運(yùn)動(dòng)目標(biāo)的難題,人們提出了脈沖多普勒(pulse Doppler,PD)雷達(dá)體制,該體制最早建立在窄帶信號(hào)基礎(chǔ)上,并在國(guó)內(nèi)外得到了廣泛應(yīng)用。PD雷達(dá)不僅能夠在頻域?qū)崿F(xiàn)運(yùn)動(dòng)目標(biāo)與雜波的分離,而且具有單根譜線濾波的能力,能夠提供精確的目標(biāo)速度信息。在預(yù)警探測(cè)和防空反導(dǎo)領(lǐng)域,要求雷達(dá)同時(shí)具備發(fā)現(xiàn)和識(shí)別目標(biāo)的能力,即要求雷達(dá)同時(shí)實(shí)現(xiàn)距離和速度二維高分辨,這就是寬帶脈沖多普勒雷達(dá)。根據(jù)相關(guān)報(bào)道美國(guó)最新的幾種相控陣?yán)走_(dá)均采用了寬帶脈沖多普勒雷達(dá)體制,但相關(guān)寬帶雷達(dá)所采用的波形、參數(shù),以及能夠兼容距離高分辨和速度高分辨的雷達(dá)解決方案等未見(jiàn)公開(kāi)報(bào)導(dǎo)。

    根據(jù)雷達(dá)系統(tǒng)理論,雷達(dá)速度測(cè)量分辨率與相參積累時(shí)長(zhǎng)有關(guān),時(shí)間越長(zhǎng)則分辨率越高,因此PD雷達(dá)采用脈沖串以加大信號(hào)處理時(shí)間;而距離測(cè)量分辨率與信號(hào)帶寬相關(guān),帶寬越大則分辨率越高,然而速度在脈沖串內(nèi)引起的回波包絡(luò)走動(dòng)也越大,造成回波相參積累損失。若采用瞬時(shí)窄帶合成寬帶信號(hào),其中每個(gè)脈沖是窄帶的,對(duì)高速目標(biāo)允許有較長(zhǎng)的相參處理時(shí)間,同時(shí)多個(gè)載頻不同的脈沖可以合成大帶寬信號(hào),提高距離分辨力[1]。綜上可知,合成寬帶信號(hào)可以兼容PD處理,合成寬帶PD雷達(dá)就是基于合成寬帶信號(hào)的寬帶PD雷達(dá)。

    本文首先闡述了合成寬帶PD雷達(dá)兩種主要的信號(hào)形式及其應(yīng)用特點(diǎn),接著介紹了合成寬帶PD雷達(dá)的距離速度二維高分辨實(shí)現(xiàn)原理,然后總結(jié)了合成寬帶PD雷達(dá)的主要特點(diǎn),最后基于自研的S波段實(shí)驗(yàn)雷達(dá),開(kāi)展了一系列微波暗室和外場(chǎng)實(shí)驗(yàn),驗(yàn)證了合成寬帶PD雷達(dá)的微動(dòng)測(cè)量能力、高精度測(cè)距測(cè)速能力和相參ISAR成像能力。

    1 合成寬帶PD雷達(dá)信號(hào)形式

    頻率步進(jìn)信號(hào)是典型的瞬時(shí)窄帶合成寬帶信號(hào)[2],以頻率步進(jìn)信號(hào)為基礎(chǔ)[3],在各脈沖或子脈沖內(nèi)進(jìn)行頻率編碼或相位編碼調(diào)制,可以得到合成寬帶PD雷達(dá)兩種主要的信號(hào)形式:即多周期脈沖間頻率步進(jìn)信號(hào)和多周期脈沖內(nèi)頻率步進(jìn)信號(hào)。

    多周期脈沖間頻率步進(jìn)信號(hào)的時(shí)頻圖如圖1所示,該信號(hào)也稱(chēng)為多周期頻率步進(jìn)脈沖序列信號(hào)。每個(gè)脈沖重復(fù)周期只發(fā)射一個(gè)脈沖,不同脈沖重復(fù)周期的脈沖具有頻率步進(jìn)的載頻,多個(gè)脈沖重復(fù)周期構(gòu)成一個(gè)頻率步進(jìn)周期,多個(gè)頻率步進(jìn)周期構(gòu)成一個(gè)相參處理周期。

    圖1 多周期脈沖間頻率步進(jìn)信號(hào)時(shí)頻圖Fig.1 Time-frequency image of multi-cycle inter-pulse stepped frequency signal

    多周期脈沖內(nèi)頻率步進(jìn)信號(hào)的時(shí)頻圖如圖2所示,該信號(hào)也稱(chēng)為多周期頻率步進(jìn)子脈沖序列信號(hào)。每個(gè)脈沖重復(fù)周期發(fā)射一串由多個(gè)載頻步進(jìn)的子脈沖組成的子脈沖序列,序列內(nèi)的各個(gè)子脈沖相互正交,多個(gè)脈沖重復(fù)周期構(gòu)成一個(gè)相參處理周期。

    圖2 多周期脈沖內(nèi)頻率步進(jìn)信號(hào)時(shí)頻圖Fig.2 Time-frequency image of multi-cycle intra-pulse stepped frequency signal

    上述兩種信號(hào)中,脈沖或子脈沖間載頻步進(jìn)方式可以是線性步進(jìn)或隨機(jī)跳頻,脈沖或子脈沖內(nèi)調(diào)制方式可以采用單頻、線性調(diào)頻或相位編碼等。

    多周期脈沖間頻率步進(jìn)信號(hào)與傳統(tǒng)的頻率步進(jìn)信號(hào)類(lèi)似,其優(yōu)點(diǎn)是信號(hào)產(chǎn)生、發(fā)射和接收簡(jiǎn)單,但是由于頻率步進(jìn)在不同發(fā)射脈沖之間進(jìn)行,相參處理周期相對(duì)較長(zhǎng),因此受目標(biāo)運(yùn)動(dòng)影響較大,適用于探測(cè)低速目標(biāo),且對(duì)數(shù)據(jù)率要求不高的場(chǎng)景。典型的信號(hào)包括簡(jiǎn)單頻率步進(jìn)信號(hào)[4]、調(diào)頻步進(jìn)信號(hào)[5]和頻率步進(jìn)相位編碼信號(hào)[6-7]。這類(lèi)信號(hào)一般采用時(shí)分的方法進(jìn)行同頻點(diǎn)分離。

    多周期脈沖內(nèi)頻率步進(jìn)信號(hào)的產(chǎn)生、發(fā)射和接收與大時(shí)寬帶寬信號(hào)一致,但是相比第一種信號(hào)形式,這種信號(hào)的頻率步進(jìn)發(fā)生在同一個(gè)發(fā)射脈沖的不同子脈沖序列之間,在同樣帶寬條件下,其相參處理周期較短,因此受目標(biāo)運(yùn)動(dòng)影響較小,適用于探測(cè)高速目標(biāo)。典型的信號(hào)包括頻率編碼信號(hào),大時(shí)寬帶寬Chirp信號(hào),頻率相位復(fù)合編碼信號(hào)[8]。為進(jìn)行PD處理,這類(lèi)信號(hào)可以采用頻分的方法對(duì)各子脈沖進(jìn)行同頻點(diǎn)分離,如針對(duì)大時(shí)寬帶寬Chirp信號(hào),提出了通過(guò)多子脈沖切分等效為寬帶PD處理的成像方法[9],該方法不僅能夠?qū)崿F(xiàn)目標(biāo)距離-速度二維高分辨,而且能夠降低信號(hào)采樣率,減少信號(hào)處理量,便于信號(hào)的實(shí)時(shí)處理。

    2 合成寬帶PD雷達(dá)原理

    合成寬帶PD雷達(dá)的目標(biāo)距離和速度二維高分辨處理分兩步進(jìn)行:先對(duì)不同頻點(diǎn)的窄帶脈沖或子脈沖進(jìn)行PD處理,從速度維對(duì)目標(biāo)與雜波進(jìn)行分辨,然后對(duì)目標(biāo)所在速度維進(jìn)行合成脈壓處理,得到目標(biāo)距離-速度二維高分辨成像結(jié)果,從而實(shí)現(xiàn)對(duì)目標(biāo)回波的相參積累。

    具體的合成寬帶PD雷達(dá)原理如圖3所示。合成寬帶PD雷達(dá)發(fā)射多周期頻率步進(jìn)信號(hào),并利用時(shí)分或頻分的方法分離出各個(gè)步進(jìn)頻點(diǎn)的回波;然后對(duì)多周期的同頻點(diǎn)回波脈沖信號(hào)進(jìn)行頻譜分析(PD處理),不僅可以從頻率域?qū)Σ煌俣鹊哪繕?biāo)和雜波進(jìn)行分離,達(dá)到抑制雜波和進(jìn)行多目標(biāo)分辨的目的,而且可測(cè)得各個(gè)目標(biāo)回波的徑向速度;進(jìn)而可對(duì)不同步進(jìn)頻點(diǎn)不同速度通道的回波進(jìn)行運(yùn)動(dòng)補(bǔ)償,減弱運(yùn)動(dòng)對(duì)一維高分辨距離成像的影響;然后對(duì)每個(gè)速度通道的不同步進(jìn)頻點(diǎn)回波信號(hào)進(jìn)行相參合成脈壓處理,進(jìn)一步得到每個(gè)速度通道的一維高分辨距離像,從而得到目標(biāo)的距離-速度二維高分辨成像結(jié)果,實(shí)現(xiàn)目標(biāo)回波的相參積累;對(duì)該相參積累結(jié)果進(jìn)行二維恒虛警檢測(cè),最終實(shí)現(xiàn)目標(biāo)距離和速度二維高分辨測(cè)量。

    圖3 合成寬帶PD雷達(dá)原理框圖Fig.3 Principle diagram of the synthetic wideband PD radar

    3 合成寬帶PD雷達(dá)特點(diǎn)

    合成寬帶PD雷達(dá)具有良好的相參性能和瞬時(shí)窄帶合成寬帶特性,使其易于實(shí)現(xiàn)距離和速度二維高分辨,抗雜波和抗干擾性能顯著增強(qiáng);同時(shí)兩種不同的信號(hào)形式使合成寬帶PD雷達(dá)信號(hào)設(shè)計(jì)更靈活,能夠適應(yīng)不同的應(yīng)用場(chǎng)景需求。下面總結(jié)了合成寬帶PD雷達(dá)的八大特點(diǎn)。

    (1) 幅相一致性校正易于實(shí)現(xiàn),能夠?qū)崿F(xiàn)接近理想的雷達(dá)性能

    大時(shí)寬帶寬Chirp信號(hào)很難進(jìn)行幅相一致性校正,若選定幾個(gè)頻點(diǎn)進(jìn)行幅相校正,校正效果有限,且校正點(diǎn)附近的頻譜也會(huì)受影響;對(duì)多周期頻率步進(jìn)合成寬帶信號(hào),可以采用時(shí)分或頻分方法區(qū)分各個(gè)步進(jìn)頻點(diǎn),這有利于實(shí)現(xiàn)基于每個(gè)頻點(diǎn)的幅相一致性校正處理。

    (2) 同時(shí)具備高距離分辨率和高多普勒分辨率

    頻率步進(jìn)合成寬帶信號(hào)為大時(shí)寬帶寬信號(hào),其模糊函數(shù)近似呈 “圖釘形”,可進(jìn)行理想的匹配濾波處理,能夠?qū)崿F(xiàn)基于寬帶信號(hào)的高距離分辨和高速度分辨。

    (3) 可解速度模糊

    對(duì)于導(dǎo)彈等超高速目標(biāo),可以根據(jù)目標(biāo)速度范圍,遍歷目標(biāo)所有可能的模糊數(shù),對(duì)目標(biāo)進(jìn)行相應(yīng)的速度補(bǔ)償,當(dāng)速度補(bǔ)償結(jié)果最接近目標(biāo)真值時(shí),匹配處理結(jié)果的幅值最大,從而可通過(guò)多通道并行處理求解速度模糊數(shù),即通過(guò)加大硬件處理能力解速度模糊[10],因此寬帶PD雷達(dá)原則上可以不用參差重復(fù)頻率來(lái)解速度模糊,在相同數(shù)據(jù)率條件下提高了雷達(dá)威力。

    (4) 可通過(guò)波形設(shè)計(jì)使雷達(dá)無(wú)距離模糊

    ①對(duì)一般雷達(dá),當(dāng)回波延時(shí)超過(guò)雷達(dá)重復(fù)周期時(shí),目標(biāo)回波會(huì)出現(xiàn)距離模糊;而對(duì)于采用多周期頻率步進(jìn)合成寬帶信號(hào)的寬帶PD雷達(dá),通過(guò)改變本振延時(shí),即利用接收機(jī)處理可以分別得到所需匹配的距離段[11-12],因此,在單個(gè)頻率步進(jìn)周期相參積累時(shí)間內(nèi),不存在距離模糊。

    ②若對(duì)頻率步進(jìn)脈沖進(jìn)行調(diào)制,如脈內(nèi)進(jìn)行相位編碼,且不同頻率步進(jìn)周期內(nèi)脈沖采用不同的相位編碼類(lèi)型,則此時(shí)信號(hào)不存在周期性,即雷達(dá)無(wú)距離模糊。

    (5) 在雜波中檢測(cè)目標(biāo)的能力增強(qiáng)

    ①采用多周期頻率步進(jìn)合成寬帶信號(hào),相對(duì)于其他寬帶信號(hào)形式,對(duì)高速目標(biāo)允許有多個(gè)脈沖重復(fù)周期的相參處理時(shí)間。

    ②對(duì)運(yùn)動(dòng)目標(biāo),PD雷達(dá)體制可以有效去除強(qiáng)雜波[13];對(duì)靜止或低速目標(biāo),由于雷達(dá)的寬帶特性,寬帶PD雷達(dá)的雜波分辨單元面積(或體積)小,雜波強(qiáng)度低,另外由于距離分辨率提高,雜波分布趨于稀疏,雜波區(qū)中出現(xiàn)大量無(wú)雜波或低雜波區(qū),使得雜波所占距離單元數(shù)減少[14]。

    ③距離上無(wú)雜波折疊。對(duì)多周期脈沖間頻率步進(jìn)信號(hào),其雷達(dá)接收機(jī)的本振信號(hào)也為頻率步進(jìn)信號(hào),只有相應(yīng)距離段回波混頻后為中頻,其他距離段回波因偏離中頻被中放的帶通濾波器所抑制,我們把合成寬帶PD雷達(dá)的這種特性稱(chēng)為距離選通性,即能夠通過(guò)改變本振信號(hào)延時(shí),來(lái)選擇性接收不同脈沖重復(fù)周期內(nèi)的目標(biāo)回波[11-12]。當(dāng)雷達(dá)工作在強(qiáng)地/海雜波環(huán)境時(shí),若采用窄帶中高重頻PD雷達(dá),較遠(yuǎn)距離的面雜波可以折疊到近距離回波上,使雜波與目標(biāo)回波在時(shí)間上重合;若采用多周期頻率步進(jìn)信號(hào)的寬帶PD雷達(dá),由于距離上無(wú)雜波折疊,故近距離小目標(biāo)只需與同距離的雜波競(jìng)爭(zhēng),提高了檢測(cè)小目標(biāo)的能力。

    (6) 利用相推原理能實(shí)現(xiàn)極高的測(cè)距精度

    頻率步進(jìn)合成寬帶信號(hào)更容易實(shí)現(xiàn)長(zhǎng)時(shí)間相參處理,有利于實(shí)現(xiàn)更高精度的相推測(cè)距,該特性可以應(yīng)用于微動(dòng)測(cè)量[15-16],為目標(biāo)識(shí)別提供更多的參考信息。

    (7) 可用于高速目標(biāo)探測(cè)

    頻率步進(jìn)合成寬帶PD雷達(dá)能同時(shí)測(cè)速和測(cè)距,再結(jié)合適當(dāng)?shù)男盘?hào)波形選擇和參數(shù)設(shè)計(jì)[17],寬帶PD雷達(dá)能夠具備較強(qiáng)的高速目標(biāo)處理能力。

    (8) 優(yōu)良的抗干擾和低截獲性能

    采用多周期頻率步進(jìn)調(diào)制脈沖和多周期頻率步進(jìn)調(diào)制子脈沖序列信號(hào),不同脈沖重復(fù)周期的雷達(dá)信號(hào)不具有周期性,并且每個(gè)發(fā)射脈沖內(nèi)可以進(jìn)行隨機(jī)的頻率和相位編碼,采用現(xiàn)有基于數(shù)字射頻存儲(chǔ) (DRFM) 技術(shù)的干擾裝備,很難對(duì)其進(jìn)行有效干擾。

    4 合成寬帶PD雷達(dá)實(shí)驗(yàn)驗(yàn)證

    本文基于自研的一套S波段合成寬帶PD雷達(dá),開(kāi)展了一系列微波暗室與外場(chǎng)實(shí)驗(yàn)。首先在微波暗室開(kāi)展了音叉振動(dòng)實(shí)驗(yàn)、角反射器位移實(shí)驗(yàn)和單擺實(shí)驗(yàn),驗(yàn)證了合成寬帶PD雷達(dá)極高的相位穩(wěn)定度和良好的相參性能,展示了其高精度相推測(cè)距和微動(dòng)測(cè)量能力;然后開(kāi)展了外場(chǎng)彈射鋼球?qū)嶒?yàn),分析得到了鋼球的高精度相推測(cè)距和測(cè)速結(jié)果,驗(yàn)證了合成寬帶PD雷達(dá)同時(shí)實(shí)現(xiàn)距離和速度二維高分辨的能力;最后開(kāi)展了民航飛機(jī)實(shí)驗(yàn),獲得了飛機(jī)的高質(zhì)量ISAR圖像,驗(yàn)證了合成寬帶PD雷達(dá)的相參ISAR成像能力。

    S波段實(shí)驗(yàn)雷達(dá)原理框圖如圖4所示[17-18],頻率步進(jìn)信號(hào)源利用DDS產(chǎn)生基帶頻率步進(jìn)信號(hào),經(jīng)過(guò)鎖相倍頻電路產(chǎn)生頻率步進(jìn)發(fā)射脈沖信號(hào)和頻率步進(jìn)本振信號(hào);回波信號(hào)經(jīng)混頻后得中頻信號(hào),放大后經(jīng)中頻采樣得回波I、Q信號(hào),信號(hào)采集完成后即可進(jìn)行非實(shí)時(shí)處理。

    圖4 S波段實(shí)驗(yàn)雷達(dá)原理框圖Fig.4 Principle diagram of the S-band trial radar

    S波段實(shí)驗(yàn)雷達(dá)發(fā)射峰值功率為100 W,信號(hào)中心載頻為3.3 GHz,中頻采樣率為100 MHz,在進(jìn)行信號(hào)處理時(shí)將根據(jù)信號(hào)瞬時(shí)帶寬進(jìn)行抽取操作。實(shí)驗(yàn)采用頻率步進(jìn)合成寬帶信號(hào),根據(jù)目標(biāo)作用距離選用簡(jiǎn)單頻率步進(jìn)信號(hào)或者頻率步進(jìn)相位編碼信號(hào),信號(hào)合成帶寬均為320 MHz。當(dāng)探測(cè)200 m以?xún)?nèi)目標(biāo),如開(kāi)展微波暗室實(shí)驗(yàn)和彈射鋼球?qū)嶒?yàn)時(shí),選用簡(jiǎn)單頻率步進(jìn)信號(hào);當(dāng)探測(cè)10 km以?xún)?nèi)目標(biāo)時(shí),如開(kāi)展民航飛機(jī)ISAR成像實(shí)驗(yàn),選用頻率步進(jìn)相位編碼信號(hào),具體參數(shù)如表 1所示。

    表1 合成寬帶信號(hào)參數(shù)表

    4.1 微動(dòng)測(cè)量能力

    目前微動(dòng)測(cè)量的難點(diǎn)在于提取幅度極小的微動(dòng)特征信號(hào),突破該瓶頸最有效的辦法是提高雷達(dá)對(duì)精細(xì)運(yùn)動(dòng)的刻畫(huà)能力,即提高雷達(dá)精確測(cè)距的能力[15]。利用相位信息進(jìn)行測(cè)距,可以顯著提高測(cè)距精度。下面通過(guò)微波暗室實(shí)驗(yàn),根據(jù)目標(biāo)回波的相位信息,提取目標(biāo)微動(dòng)。

    4.1.1 音叉振動(dòng)實(shí)驗(yàn)

    在微波暗室中,分別用頻率步進(jìn)PD雷達(dá)和麥克風(fēng)音頻處理方法測(cè)定兩個(gè)標(biāo)稱(chēng)值分別為64 Hz和512 Hz的音叉振動(dòng)頻率。采用麥克風(fēng)錄入的音頻處理方法,測(cè)量得到64 Hz音叉頻率為63 Hz,512 Hz音叉頻率為536 Hz。

    實(shí)驗(yàn)測(cè)量64 Hz音叉頻率為63.1 Hz,512 Hz音叉頻率為535.6 Hz,與麥克風(fēng)測(cè)量值基本一致。圖5顯示出了512 Hz音叉回波分析結(jié)果,由音叉測(cè)距結(jié)果可知,音叉振動(dòng)幅度小于0.1 mm。音叉振動(dòng)實(shí)驗(yàn)驗(yàn)證了合成寬帶雷達(dá)測(cè)量極微小運(yùn)動(dòng)的能力,展示了實(shí)驗(yàn)雷達(dá)極高的相位穩(wěn)定度。

    圖5 512 Hz音叉回波測(cè)距及頻譜分析結(jié)果Fig.5 Ranging result and spectrum analysis of 512 Hz tuning fork echo

    4.1.2 角反射器位移實(shí)驗(yàn)

    在微波暗室中,將角反射器置于位移精度可達(dá)0.01 mm的精密位移平臺(tái)上,實(shí)驗(yàn)中通過(guò)精密位移平臺(tái)調(diào)整角反射器相對(duì)于測(cè)量雷達(dá)的徑向距離,并通過(guò)雷達(dá)測(cè)量該距離,從而驗(yàn)證雷達(dá)的測(cè)相和測(cè)距精度。

    實(shí)驗(yàn)中設(shè)置平臺(tái)的位移增量分別為1 mm,0.5 mm,0.2 mm和0.1 mm,且每種情況錄取10個(gè)位置。圖6顯示出了位移增量為0.2 mm時(shí)的實(shí)驗(yàn)結(jié)果。

    圖6 位移增量為0.2 mm時(shí)角反射器回波測(cè)距結(jié)果Fig.6 Ranging results of corner reflector echo when displacement increment is 0.2 mm

    由于積累后角反射器回波信噪比較高,各位置距離測(cè)量均方根誤差約為0.001 6 mm,即相位測(cè)量均方根誤差為0.013°,驗(yàn)證了實(shí)驗(yàn)雷達(dá)具有極高的相位穩(wěn)定度;同時(shí),由測(cè)距結(jié)果可知,其測(cè)距精度達(dá)亞mm量級(jí),驗(yàn)證了合成寬帶PD雷達(dá)的高精度相推測(cè)距能力。

    4.1.3 單擺實(shí)驗(yàn)

    在微波暗室中,利用合成寬帶雷達(dá)觀測(cè)單擺運(yùn)動(dòng)軌跡,并測(cè)量其振動(dòng)頻率。

    當(dāng)θ較小時(shí)(θ<5°),sinθ≈θ,即當(dāng)單擺擺動(dòng)幅度很小時(shí),單擺作簡(jiǎn)諧振動(dòng),振動(dòng)頻率為

    (1)

    當(dāng)θ較大時(shí)(θ≥5°),單擺作非線性運(yùn)動(dòng),其振動(dòng)頻率包括1次、3次和5次諧波:

    (2)

    式中,ω≈ω0。實(shí)驗(yàn)中單擺所用鋼球直徑為3 cm,單擺擺長(zhǎng)為0.776 m,由式(1)可知,鋼球單擺理論頻率為0.565 8 Hz。

    當(dāng)鋼球擺幅較大時(shí),單擺運(yùn)動(dòng)軌跡如7 (a) 所示,由圖可知單擺做類(lèi)正弦運(yùn)動(dòng);對(duì)單擺運(yùn)動(dòng)軌跡進(jìn)行頻譜分析,分析結(jié)果如7(b) 所示,從圖中不僅能夠觀察到0.569 Hz的單擺振動(dòng)頻率,還能觀察到頻率為1.705 Hz的3次諧波和2.843 Hz的5次諧波。

    單擺實(shí)驗(yàn)驗(yàn)證了合成寬帶雷達(dá)測(cè)量目標(biāo)微小運(yùn)動(dòng)的能力,同時(shí)由單擺頻譜測(cè)量值與理論值對(duì)比分析可知,合成寬帶相推測(cè)距的精度很高,該技術(shù)有望應(yīng)用于目標(biāo)微動(dòng)特性反演。

    圖7 鋼球單擺運(yùn)動(dòng)分析Fig.7 Pendulum oscillation analysis of the ejection ball

    4.2 高精度測(cè)距測(cè)速能力

    彈射直徑為5 cm的鋼球,通過(guò)合成寬帶PD處理,即從多普勒域抑制靜止或低速?gòu)?qiáng)雜波,并得到鋼球所在速度維上的一維高分辨距離像,從距離和速度兩維同時(shí)實(shí)現(xiàn)鋼球回波的相參積累,提高回波信噪比;同時(shí)由于寬帶PD雷達(dá)良好的相參性能,可以在所得一維高分辨距離像的基礎(chǔ)上對(duì)鋼球進(jìn)行高精度相推測(cè)距和測(cè)速處理。

    相位測(cè)距存在嚴(yán)重的距離模糊,常采用包絡(luò)測(cè)距解相位測(cè)距模糊的思想[19-20],該方法稱(chēng)為相推測(cè)距。由于實(shí)際中并不確知雷達(dá)初相和目標(biāo)相對(duì)雷達(dá)的真實(shí)距離,故一般通過(guò)測(cè)量相鄰幀距離增量實(shí)現(xiàn)相推測(cè)距[16],而相鄰幀距離變化率即為目標(biāo)速度,故求取相鄰幀距離增量后,即可得到相應(yīng)的高精度相推測(cè)距和測(cè)速結(jié)果。

    在傳統(tǒng)包絡(luò)測(cè)距解相位模糊思想的基礎(chǔ)上,我們提出了一種低信噪比條件下的相位解模糊方法,即進(jìn)一步聯(lián)合航跡平滑和最小熵方法解相位模糊,推進(jìn)了相推技術(shù)的工程化應(yīng)用。下面給出高精度相推測(cè)距和測(cè)速結(jié)果。

    鋼球隨時(shí)間變化的一維高分辨距離像如圖8所示,由于圖像顯示的動(dòng)態(tài)范圍有限,故將第1~2 650幀一維距離像分開(kāi)繪制。

    圖8 鋼球一維高分辨距離像Fig.8 HRRPs of the ejection ball

    鋼球包絡(luò)測(cè)距結(jié)果如圖9所示,圖9(a)中濾波前后結(jié)果之差視為包絡(luò)測(cè)距起伏誤差,如圖9(b)所示,統(tǒng)計(jì)可得包絡(luò)測(cè)距起伏均方根誤差為14.19 mm。

    圖9 包絡(luò)測(cè)距結(jié)果與測(cè)距起伏誤差示意圖Fig.9 Schematic diagrams of envelope range measurement result and fluctuation error

    鋼球相推測(cè)速結(jié)果如圖10所示。由于鋼球是從雷達(dá)左側(cè)斜向右側(cè)彈射,在彈射出去的一小段時(shí)間內(nèi),鋼球相對(duì)雷達(dá)視角變化較大,從而徑向速度變化較大,如圖中速度上升段所示,此后鋼球徑向速度主要受空氣阻力的影響,如圖中減速段所示。由于不知道鋼球真實(shí)速度,故將速度擬合前后的偏差視為相推測(cè)速起伏誤差,如圖10(b)所示,隨著鋼球遠(yuǎn)離雷達(dá),信噪比逐漸降低,約從40 dB以上一直降低到20 dB左右,故相推測(cè)速起伏誤差逐漸增大;當(dāng)脈壓后一維距離像信噪比大于30 dB時(shí),相推測(cè)速起伏均方根誤差為0.057 m/s,當(dāng)信噪比大于20 dB時(shí),相推測(cè)速起伏均方根誤差為0.25 m/s。實(shí)際中將平滑濾波后的結(jié)果作為最終的相推測(cè)速結(jié)果,因?yàn)槠交瑸V波相當(dāng)于對(duì)回波進(jìn)行非相參積累,即等效為提高相推測(cè)速信噪比,故相比于直接測(cè)量結(jié)果,平滑濾波后的測(cè)速精度大約能提高一個(gè)數(shù)量級(jí)。

    圖10 相推測(cè)速結(jié)果與測(cè)速起伏誤差示意圖Fig.10 Schematic diagrams of phase derived velocity measurement result and fluctuation error

    同上,根據(jù)相推原理,可推得相推測(cè)距均方根誤差,即當(dāng)信噪比大于30 dB時(shí),相推測(cè)距均方根誤差為0.048 mm,當(dāng)信噪比大于20 dB時(shí),相推測(cè)距均方根誤差為0.21 mm。

    4.3 相參ISAR成像能力

    基于合成寬帶PD雷達(dá)良好的相參性能,提出了一種基于相推測(cè)速的ISAR成像運(yùn)動(dòng)補(bǔ)償方法,該方法結(jié)合距離像互相關(guān)測(cè)速和相推技術(shù),可以獲得高精度的平動(dòng)速度測(cè)量結(jié)果,其測(cè)量精度同時(shí)滿足ISAR成像包絡(luò)對(duì)齊和相位校正的精度要求,能夠?qū)崿F(xiàn)精確運(yùn)動(dòng)補(bǔ)償,在小轉(zhuǎn)角情況下,再結(jié)合距離-多普勒算法就可以獲得高質(zhì)量的ISAR圖像。

    下面分析波音777(Boeing 777-200)飛機(jī)的回波數(shù)據(jù),該飛機(jī)的機(jī)身長(zhǎng)為63.73 m,翼展為60.93 m,錄取場(chǎng)景為飛機(jī)即將降落,此時(shí)飛機(jī)距機(jī)場(chǎng)約為20 km。飛機(jī)隨時(shí)間變化的一維距離像如圖11所示,圖中共繪出了1 024幀一維距離像,飛機(jī)由遠(yuǎn)及近飛向雷達(dá)。

    圖11 飛機(jī)一維高分辨距離像Fig.11 HRRPs of the airplane

    基于距離像互相關(guān)的相推測(cè)速結(jié)果如圖12所示,濾波前后結(jié)果之差視為相推測(cè)速起伏誤差,如圖12(b)所示,統(tǒng)計(jì)可得相推測(cè)速的起伏均方根誤差為0.061 m/s。同第4.2節(jié)分析,可知平滑濾波相當(dāng)于對(duì)回波進(jìn)行非相參積累,可以進(jìn)一步提高速度測(cè)量精度,故將平滑濾波結(jié)果用于飛機(jī)回波的包絡(luò)對(duì)齊和相位校正處理。

    圖12 相推測(cè)速結(jié)果與測(cè)速起伏誤差示意圖Fig.12 Schematic diagrams of phase derived velocity measurement result and fluctuation error

    飛機(jī)的ISAR成像結(jié)果如圖13所示,其中圖13(a)~圖13(c)分別為飛機(jī)飛向雷達(dá)、與雷達(dá)視線近似垂直和遠(yuǎn)離雷達(dá)時(shí)的ISAR圖像,從圖中可看到清晰的飛機(jī)輪廓,包括飛機(jī)機(jī)頭、機(jī)翼、機(jī)尾以及發(fā)動(dòng)機(jī);同時(shí)由這3張圖可以判斷,該飛機(jī)有兩個(gè)發(fā)動(dòng)機(jī),且當(dāng)飛機(jī)飛向雷達(dá)時(shí),每個(gè)發(fā)動(dòng)機(jī)存在不止一個(gè)散射點(diǎn);當(dāng)飛機(jī)與雷達(dá)視線近乎垂直時(shí),發(fā)動(dòng)機(jī)僅存在一個(gè)強(qiáng)散射點(diǎn);當(dāng)飛機(jī)遠(yuǎn)離雷達(dá)時(shí),發(fā)動(dòng)機(jī)回波被遮擋。

    圖13(b)中飛機(jī)近似與雷達(dá)垂直,故以其為例分析ISAR圖像中飛機(jī)尺寸。飛機(jī)ISAR成像積累時(shí)長(zhǎng)為1.048 6 s,錄取時(shí)刻飛機(jī)速度約為400 km/h,飛機(jī)距雷達(dá)2 280 m,則飛機(jī)轉(zhuǎn)角近似為2.92°,因此橫向分辨率約為0.891 m,由圖13(b)可知其方位向占據(jù)第97~162個(gè)分辨單元,故飛機(jī)橫向尺寸約為58.8 m,又由圖(b)可知飛機(jī)在縱向尺寸約為56 m。飛機(jī)橫向與縱向尺寸分別對(duì)應(yīng)其機(jī)長(zhǎng)和翼展,可知測(cè)量結(jié)果與飛機(jī)真實(shí)尺寸近似相符。

    綜上,ISAR圖像可以用于估計(jì)目標(biāo)尺寸,為目標(biāo)識(shí)別提供參考信息;同時(shí)合成寬帶PD雷達(dá)回波相參,基于距離像互相關(guān)的相推測(cè)速方法可以充分利用ISAR成像時(shí)間內(nèi)全部回波的信息,相當(dāng)于進(jìn)行了整體回波最優(yōu)處理,測(cè)量結(jié)果不存在積累誤差,因此該方法適用于較低信噪比環(huán)境,根據(jù)飛機(jī)回波數(shù)據(jù),可知飛機(jī)各散射點(diǎn)回波幅度相差較大,若定義信噪比為脈壓后一維距離像中最強(qiáng)散射點(diǎn)處的信噪比,則根據(jù)實(shí)測(cè)數(shù)據(jù)統(tǒng)計(jì),可知該方法適用的最低信噪比約為15 dB。

    圖13 飛機(jī)ISAR成像結(jié)果Fig.13 ISAR imaging results of the airplane

    5 結(jié) 論

    合成寬帶PD雷達(dá)可以同時(shí)實(shí)現(xiàn)距離和速度二維高分辨,具有較強(qiáng)的抗雜波和抗干擾性能;此外,微動(dòng)提取和ISAR成像均可為目標(biāo)識(shí)別提供重要的參考信息。本文通過(guò)音叉振動(dòng)實(shí)驗(yàn)、角反射器位移實(shí)驗(yàn)、單擺實(shí)驗(yàn),驗(yàn)證了合成寬帶PD雷達(dá)的微動(dòng)測(cè)量能力;通過(guò)彈射鋼球?qū)嶒?yàn),驗(yàn)證了合成寬帶PD雷達(dá)在強(qiáng)雜波中檢測(cè)運(yùn)動(dòng)小目標(biāo)的能力,以及高精度相推測(cè)距測(cè)速性能,其中相對(duì)距離的測(cè)量精度達(dá)到了亞mm或mm量級(jí),測(cè)速精度達(dá)到了cm/s或亞m/s量級(jí);通過(guò)民航飛機(jī)實(shí)驗(yàn),驗(yàn)證了合成寬帶PD雷達(dá)對(duì)復(fù)雜目標(biāo)的高精度相推測(cè)速能力,以及相參ISAR成像能力;上述實(shí)驗(yàn)均驗(yàn)證了合成寬帶PD雷達(dá)良好的相參性能,為其工程化應(yīng)用奠定了堅(jiān)實(shí)的基礎(chǔ)。

    合成寬帶PD雷達(dá)是下一代新體制雷達(dá)發(fā)展的重要方向,該體制雷達(dá)可應(yīng)用于地面對(duì)空情報(bào)雷達(dá)、預(yù)警機(jī)雷達(dá)、戰(zhàn)場(chǎng)偵察及地面火控雷達(dá)等。

    [1] Liu H B, Lu J D. Target motion compensation algorithm based on keystone transform for wideband pulse Doppler radar[J].TransactionsofBeijingInstituteofTechnology, 2012, 32(6): 625-630.(劉海波,盧俊道.基于Keystone變換的寬帶脈沖多普勒雷達(dá)目標(biāo)運(yùn)動(dòng)補(bǔ)償算法[J].北京理工大學(xué)學(xué)報(bào),2012,32(6):625-630.)

    [2] Wehner D R.High-resolutionradar[M]. 2nd ed. Boston: Artech House, 1995: 197-237.

    [3] Huang M, Long T. Design and realization of the kind of step frequency and pulse Doppler complex processing[J].JournalofProjectiles,Rockets,MissilesandGuidance, 2005, 25(1): 72-75.(黃默,龍騰.一種步進(jìn)頻率和PD復(fù)合體制的信號(hào)處理技術(shù)[J].彈箭與制導(dǎo)學(xué)報(bào),2005,25(1):72-75.)

    [4] Long T, Li D, Wu Q Z. Design methods for step frequency waveform and the target pick-up algorithm[J].SystemsEngineeringandElectronics,2001,23(6):26-31.(龍騰,李眈,吳瓊之.頻率步進(jìn)雷達(dá)參數(shù)設(shè)計(jì)與目標(biāo)抽取算法[J].系統(tǒng)工程與電子技術(shù),2001,23(6):26-31.)

    [5] Lord R T. Aspects of stepped-frequency processing for low-frequency SAR systems[D]. Cape Town: University of Cape Town, 2000.

    [6] Berizzi F, Martorella M, Bernabo M. A range profiling technique for synthetic wideband radar[J].IEEERadarSonarNavigation, 2008, V2(5): 334-350.

    [7] Sitler K L, Temple M A, Novack R C. High range resolution profiling using phase-coded, stepped-frequency waveforms[J].ElectronicsLetters, 2002, 38(1):46-48.

    [8] He Y. Study on the coded multi-modulation signal analysis and processing[D]. Changsha: National University of Defense Technology, 2009. (何元. 編碼復(fù)合調(diào)制信號(hào)分析與處理技術(shù)研究[D]. 長(zhǎng)沙: 國(guó)防科技大學(xué), 2009.)

    [9] Zhang H G, Fan H Y, He S H, et al. Multi-subpulse process of large time-bandwidth product chirp signal[J].JournalofBeijingInstituteofTechnology, 2015, 24(4):501-507.

    [10] Li Y, Zeng T, Long T, et al. Range migration compensation and Doppler ambiguity resolution by keystone transform[C]∥Proc.oftheInternationalConferenceonRadar, 2006:1-4.

    [11] Myers H, Moore R, Butler W, et al. Common aperture multiband radar (CAMBR)[C]∥Proc.oftheRadarProcessing,Technology,andApplications, 1996, V2845:2-13.

    [12] Myers H, Moore R, Butler W, et al. Common aperture multiband radar (CAMBR)[C]∥Proc.oftheIEEENationalRadarConference, 1997: 143-148.

    [13] Bao Y X, Mao E K, He P K, et al. Clutter cancellation and high resolution processing in HPRF stepped-frequency radar[J].ModernRadar, 2009, 31(1): 58-61. (包云霞, 毛二可,何佩琨, 等. 高重頻頻率步進(jìn)雷達(dá)的雜波抑制和高分辨處理[J]. 現(xiàn)代雷達(dá), 2009, 31(1): 58-61.)

    [14] Merrill I S.Introductiontoradarsystems[M]. 3rd ed. Beijing: Publishing House of Electronics Industry, 2007: 427-430.

    [15] Fan H Y, Ren L X, Mao E K. A micro-motion measurement method on wideband radar phase derived ranging[C]∥Proc.oftheIETRadarConference, 2013: 1-4.

    [16] Liu Y X, Zhu D K, Li X, et al. Micromotion characteristic acquisition based on wideband radar phase[J].IEEETrans.onGeoscienceandRemoteSensing, 2014, 52(6): 3650-3657.

    [17] Huang C, Ren L X, Mao E K. Design of a time delayed receiver local oscillator in stepped-frequency radar[C]∥Proc.oftheIETRadarConference, 2009: 1-4.

    [18] Ren L X, Long T, Yuan H P. Signal processing and parameter design of HPRF pulsed-Doppler stepped-frequency radar[J].ActaElectronicaSinica, 2007, 35(9): 1630-1636. (任麗香, 龍騰, 遠(yuǎn)海鵬. HPRF脈沖多普勒頻率步進(jìn)雷達(dá)信號(hào)處理與參數(shù)設(shè)計(jì)[J]. 電子學(xué)報(bào), 2007, 35(9): 1630-1636.)

    [19] Fritz S. An improved process for phase-derived-range measurements[P]. World Intellectual Property Organization Patent, 2005.

    [20] Fritz S.Process for phase-derived-range measurements[P].U.S.Patent,US 7046190 B2, 2005-2-10.

    Synthetic wideband pulse Doppler radar

    MAO Er-ke1,2, FAN Hua-yu1,2

    (1.BeijingKeyLaboratoryofEmbeddedReal-timeInformationProcessingTechnology,BeijingInstituteofTechnology,Beijing100081,China; 2.RadarResearchLab,SchoolofInformationandElectronics,BeijingInstituteofTechnology,Beijing100081,China)

    The synthetic wideband pulse Doppler (PD) radar can simultaneously achieve high range and velocity resolution. And it has good performance of coherence, so high range and velocity measurement precision can be achieved by using the phase derived principle and it can be applied to micro-motion measurement of complex targets with multi-components. Besides, it has good performance of anti-clutter and anti-jamming. In this paper a synthetic wideband PD radar system and its main signal forms are presented; the principle of achieving high range and velocity resolution is introduced; moreover, the main properties of the synthetic wideband PD radar are summarized; finally, the abilities of the synthetic wideband PD radar, including the abilities of micro-motion measurement, precise range and velocity measurement and coherent ISAR imaging, are verified by a series of experiments.

    synthetic wideband signal; pulse Doppler (PD) radar; high range and velocity resolution; anti-clutter; anti-jamming; micro-motion measurement

    2016-04-27;

    2016-07-14;網(wǎng)絡(luò)優(yōu)先出版日期:2016-08-22。

    國(guó)家自然科學(xué)基金(61301189);高等學(xué)校學(xué)科創(chuàng)新引智計(jì)劃(B14010)資助課題

    TN 958

    A

    10.3969/j.issn.1001-506X.2016.12.06

    毛二可(1934-),男,中國(guó)工程院院士,教授,博士研究生導(dǎo)師,主要研究方向?yàn)槔走_(dá)系統(tǒng)、數(shù)字信號(hào)處理。

    E-mail:maoerke@bit.edu.cn

    范花玉(1989-),女,博士研究生,主要研究方向?yàn)楹铣蓪拵Ю走_(dá)信號(hào)處理。

    E-mail:fan_huayu@sina.com

    網(wǎng)絡(luò)優(yōu)先出版地址:http:∥www.cnki.net/kcms/detail/11.2422.TN.20160822.1008.010.html

    猜你喜歡
    鋼球雜波測(cè)距
    STAR2000型空管一次雷達(dá)雜波抑制淺析
    昆鋼球團(tuán)生產(chǎn)降事故率降消耗的生產(chǎn)實(shí)踐
    昆鋼科技(2021年2期)2021-07-22 07:46:54
    8Cr4Mo4V鋼球壓碎載荷不合格的原因及返修
    哈爾濱軸承(2020年4期)2020-03-17 08:13:52
    類(lèi)星體的精準(zhǔn)測(cè)距
    科學(xué)(2020年3期)2020-01-06 04:02:51
    自動(dòng)控制的鋼球分離器設(shè)計(jì)與實(shí)現(xiàn)
    淺談超聲波測(cè)距
    電子制作(2017年7期)2017-06-05 09:36:13
    密集雜波環(huán)境下確定性退火DA-HPMHT跟蹤算法
    相關(guān)廣義復(fù)合分布雷達(dá)海雜波仿真
    基于PSOC超聲測(cè)距系統(tǒng)設(shè)計(jì)
    相對(duì)差分單項(xiàng)測(cè)距△DOR
    太空探索(2014年1期)2014-07-10 13:41:50
    h视频一区二区三区| 黄色怎么调成土黄色| 欧美最新免费一区二区三区| 天天躁夜夜躁狠狠久久av| 美女主播在线视频| 天堂8中文在线网| 曰老女人黄片| 日韩精品免费视频一区二区三区 | 精品人妻偷拍中文字幕| 十分钟在线观看高清视频www| 人人妻人人添人人爽欧美一区卜| 麻豆精品久久久久久蜜桃| 国产xxxxx性猛交| 成人亚洲欧美一区二区av| av在线观看视频网站免费| 中国三级夫妇交换| 亚洲婷婷狠狠爱综合网| 内地一区二区视频在线| 在线天堂中文资源库| 亚洲av中文av极速乱| 黑人高潮一二区| 91精品国产国语对白视频| 久久精品久久久久久久性| 午夜免费观看性视频| 老司机影院毛片| 观看av在线不卡| 成人18禁高潮啪啪吃奶动态图| 五月天丁香电影| 欧美国产精品va在线观看不卡| 日产精品乱码卡一卡2卡三| 中文字幕av电影在线播放| 国产熟女午夜一区二区三区| 欧美最新免费一区二区三区| 少妇高潮的动态图| 日韩制服丝袜自拍偷拍| 国产老妇伦熟女老妇高清| 精品久久久久久电影网| 国产成人一区二区在线| 只有这里有精品99| 久久久久久久精品精品| 在线亚洲精品国产二区图片欧美| 女人久久www免费人成看片| 中文欧美无线码| 成人影院久久| 水蜜桃什么品种好| 国产精品人妻久久久久久| 午夜福利网站1000一区二区三区| 99久久中文字幕三级久久日本| 久久韩国三级中文字幕| 高清在线视频一区二区三区| av.在线天堂| 色婷婷av一区二区三区视频| 国产精品久久久久成人av| av在线app专区| 90打野战视频偷拍视频| 七月丁香在线播放| www.av在线官网国产| 一级片免费观看大全| 最近中文字幕2019免费版| 亚洲国产精品一区三区| 91aial.com中文字幕在线观看| 国产精品嫩草影院av在线观看| 国产精品偷伦视频观看了| 国产在线视频一区二区| 波野结衣二区三区在线| 国产免费一级a男人的天堂| 国产精品国产三级专区第一集| 亚洲欧美一区二区三区黑人 | 久久久精品区二区三区| 啦啦啦啦在线视频资源| 晚上一个人看的免费电影| 欧美成人午夜精品| 男男h啪啪无遮挡| 丝袜美足系列| 啦啦啦中文免费视频观看日本| 日韩av免费高清视频| 国产一区亚洲一区在线观看| 亚洲国产精品专区欧美| 国产精品国产三级国产av玫瑰| 国产一区二区三区av在线| 91午夜精品亚洲一区二区三区| 少妇被粗大猛烈的视频| xxx大片免费视频| 9色porny在线观看| 国产成人免费无遮挡视频| 久久午夜福利片| 国产免费视频播放在线视频| 高清欧美精品videossex| 国产乱来视频区| 亚洲综合精品二区| 婷婷色综合www| 国产精品一二三区在线看| 伊人久久国产一区二区| 美女国产高潮福利片在线看| 日日摸夜夜添夜夜爱| 亚洲欧洲精品一区二区精品久久久 | 在线亚洲精品国产二区图片欧美| 老司机亚洲免费影院| 久久精品国产综合久久久 | videos熟女内射| 国产69精品久久久久777片| 男女高潮啪啪啪动态图| 又黄又爽又刺激的免费视频.| 丝袜喷水一区| 一二三四在线观看免费中文在 | 欧美亚洲日本最大视频资源| 丰满少妇做爰视频| 久久久久国产精品人妻一区二区| 捣出白浆h1v1| 99久久中文字幕三级久久日本| 我要看黄色一级片免费的| videos熟女内射| 在线天堂最新版资源| 亚洲精品,欧美精品| 人妻一区二区av| 考比视频在线观看| 韩国精品一区二区三区 | 午夜福利视频精品| 国产xxxxx性猛交| 母亲3免费完整高清在线观看 | 久久韩国三级中文字幕| 国产又爽黄色视频| 永久免费av网站大全| 精品一区在线观看国产| 国产精品蜜桃在线观看| 你懂的网址亚洲精品在线观看| 亚洲av福利一区| 我要看黄色一级片免费的| 日韩精品免费视频一区二区三区 | 女性被躁到高潮视频| 国产淫语在线视频| 91aial.com中文字幕在线观看| 91在线精品国自产拍蜜月| 久久人人97超碰香蕉20202| 人人妻人人澡人人爽人人夜夜| 亚洲欧洲国产日韩| www日本在线高清视频| 午夜福利乱码中文字幕| av福利片在线| 久久午夜综合久久蜜桃| 午夜福利网站1000一区二区三区| 久久久久久久精品精品| 久久精品国产自在天天线| 精品一区二区三区四区五区乱码 | 一级黄片播放器| 99久久综合免费| 亚洲 欧美一区二区三区| 亚洲人成77777在线视频| 欧美变态另类bdsm刘玥| 考比视频在线观看| 国产精品人妻久久久久久| 欧美变态另类bdsm刘玥| 日本欧美视频一区| 一个人免费看片子| 成人综合一区亚洲| a级片在线免费高清观看视频| 久久女婷五月综合色啪小说| av播播在线观看一区| 色吧在线观看| 成人漫画全彩无遮挡| videossex国产| 丰满饥渴人妻一区二区三| 亚洲国产精品国产精品| 精品国产一区二区三区四区第35| 久久精品熟女亚洲av麻豆精品| 久久鲁丝午夜福利片| 成年美女黄网站色视频大全免费| 国产黄色免费在线视频| 99久久中文字幕三级久久日本| 色婷婷久久久亚洲欧美| 久久久久久久久久成人| 国产男女内射视频| 熟妇人妻不卡中文字幕| 亚洲激情五月婷婷啪啪| 亚洲综合精品二区| 99国产精品免费福利视频| 精品久久久久久电影网| 国产成人aa在线观看| 国产在线视频一区二区| 黑人猛操日本美女一级片| 丝袜在线中文字幕| 人成视频在线观看免费观看| 嫩草影院入口| 日韩成人伦理影院| 在线观看三级黄色| 亚洲av综合色区一区| 99精国产麻豆久久婷婷| 在线观看美女被高潮喷水网站| 亚洲av在线观看美女高潮| 精品福利永久在线观看| 国产成人免费无遮挡视频| 精品一区二区三卡| 精品久久久久久电影网| 三级国产精品片| 亚洲成人av在线免费| 欧美国产精品一级二级三级| 国产亚洲av片在线观看秒播厂| 视频在线观看一区二区三区| 一本大道久久a久久精品| 中国美白少妇内射xxxbb| 亚洲欧美成人综合另类久久久| 国产老妇伦熟女老妇高清| 亚洲精品美女久久av网站| 国产欧美另类精品又又久久亚洲欧美| 成年av动漫网址| 国产福利在线免费观看视频| 亚洲高清免费不卡视频| 欧美激情 高清一区二区三区| 国产麻豆69| 97人妻天天添夜夜摸| 看免费成人av毛片| 少妇的逼好多水| 在线亚洲精品国产二区图片欧美| 超色免费av| 成人无遮挡网站| 国产精品人妻久久久久久| 亚洲综合色惰| 中文字幕制服av| 交换朋友夫妻互换小说| 一区二区三区乱码不卡18| 尾随美女入室| 午夜久久久在线观看| 晚上一个人看的免费电影| 边亲边吃奶的免费视频| 一本久久精品| 久久精品人人爽人人爽视色| 黄色怎么调成土黄色| 亚洲国产看品久久| 国产成人免费观看mmmm| 超碰97精品在线观看| 国产av精品麻豆| 国产精品女同一区二区软件| 精品一品国产午夜福利视频| 99热网站在线观看| 女人被躁到高潮嗷嗷叫费观| 久久99精品国语久久久| 国产精品一区www在线观看| 色94色欧美一区二区| 成人影院久久| 秋霞伦理黄片| 日韩三级伦理在线观看| 乱人伦中国视频| 日韩一本色道免费dvd| 亚洲成人手机| 亚洲精品美女久久av网站| 国国产精品蜜臀av免费| 国产一区有黄有色的免费视频| 欧美人与性动交α欧美精品济南到 | 欧美国产精品一级二级三级| 国产又爽黄色视频| 女人精品久久久久毛片| 国产成人欧美| 久久99蜜桃精品久久| 高清欧美精品videossex| 日韩,欧美,国产一区二区三区| 国产亚洲午夜精品一区二区久久| 亚洲国产毛片av蜜桃av| 极品人妻少妇av视频| 妹子高潮喷水视频| 男女啪啪激烈高潮av片| 亚洲四区av| 女人久久www免费人成看片| 成人影院久久| 国产av一区二区精品久久| 一级毛片 在线播放| 亚洲成人一二三区av| 国产免费一级a男人的天堂| 亚洲婷婷狠狠爱综合网| 人人妻人人爽人人添夜夜欢视频| 丰满乱子伦码专区| 国产乱人偷精品视频| 日本与韩国留学比较| 女性被躁到高潮视频| 免费av不卡在线播放| 免费大片18禁| 日韩一区二区视频免费看| 精品一区二区三区四区五区乱码 | 伊人亚洲综合成人网| 久久午夜福利片| 人妻一区二区av| 国产亚洲欧美精品永久| 日本91视频免费播放| 男人操女人黄网站| 丰满饥渴人妻一区二区三| 交换朋友夫妻互换小说| 日本欧美国产在线视频| 久久国产精品大桥未久av| 黄片播放在线免费| 18禁国产床啪视频网站| 亚洲精品乱码久久久久久按摩| 国产激情久久老熟女| 欧美bdsm另类| 亚洲国产精品专区欧美| 一级黄片播放器| 18+在线观看网站| 哪个播放器可以免费观看大片| 最近手机中文字幕大全| 婷婷色综合www| 新久久久久国产一级毛片| 亚洲av中文av极速乱| 免费观看无遮挡的男女| av在线老鸭窝| 久久99蜜桃精品久久| 啦啦啦视频在线资源免费观看| 国产日韩欧美亚洲二区| 午夜日本视频在线| 成人国产av品久久久| 只有这里有精品99| 91久久精品国产一区二区三区| 中文字幕免费在线视频6| 国产精品一区二区在线观看99| 国产熟女午夜一区二区三区| 国产av国产精品国产| 国产在线一区二区三区精| 中文字幕人妻丝袜制服| 欧美另类一区| 久久久久网色| 国产淫语在线视频| 亚洲精品日韩在线中文字幕| 日韩成人av中文字幕在线观看| 日韩一区二区三区影片| 人成视频在线观看免费观看| 国产精品一区www在线观看| 制服丝袜香蕉在线| 国产精品熟女久久久久浪| 国产熟女欧美一区二区| 婷婷色综合www| 亚洲成人av在线免费| 国产精品久久久久久av不卡| 色94色欧美一区二区| 日日爽夜夜爽网站| 欧美激情国产日韩精品一区| 在线观看美女被高潮喷水网站| 亚洲精品一二三| √禁漫天堂资源中文www| 精品国产一区二区三区四区第35| 久久人人爽人人爽人人片va| 国产日韩欧美在线精品| 国产亚洲欧美精品永久| 日韩一区二区视频免费看| 午夜激情久久久久久久| 夜夜爽夜夜爽视频| 一级片'在线观看视频| 免费少妇av软件| 中国三级夫妇交换| 国产免费又黄又爽又色| 国产一区亚洲一区在线观看| 日韩视频在线欧美| 亚洲国产欧美在线一区| 亚洲国产毛片av蜜桃av| 色婷婷av一区二区三区视频| 久久久国产欧美日韩av| 国产一区亚洲一区在线观看| 日韩电影二区| 大片电影免费在线观看免费| 久久久精品免费免费高清| 九色成人免费人妻av| 天天影视国产精品| 精品久久国产蜜桃| 18+在线观看网站| 狂野欧美激情性xxxx在线观看| 久久毛片免费看一区二区三区| 日韩av免费高清视频| 国产日韩欧美在线精品| 有码 亚洲区| 午夜免费观看性视频| 免费不卡的大黄色大毛片视频在线观看| 国产无遮挡羞羞视频在线观看| 国产淫语在线视频| 国产成人91sexporn| 国产精品久久久久久精品电影小说| 91精品三级在线观看| 日韩伦理黄色片| 国产精品嫩草影院av在线观看| 精品人妻熟女毛片av久久网站| 大码成人一级视频| 中文字幕亚洲精品专区| 美女大奶头黄色视频| 亚洲成人手机| 国产精品不卡视频一区二区| 99久久中文字幕三级久久日本| 啦啦啦啦在线视频资源| 亚洲欧美精品自产自拍| 最近中文字幕高清免费大全6| 国产淫语在线视频| 美女主播在线视频| av女优亚洲男人天堂| 黄色一级大片看看| 欧美另类一区| 亚洲精品一二三| 久久午夜综合久久蜜桃| 国产在线视频一区二区| 丝袜美足系列| 一级毛片黄色毛片免费观看视频| 亚洲国产最新在线播放| 黄色一级大片看看| 亚洲,一卡二卡三卡| 亚洲,欧美精品.| 中文字幕最新亚洲高清| 日本黄色日本黄色录像| 美女内射精品一级片tv| 香蕉丝袜av| 女性被躁到高潮视频| 久久久久久久久久久久大奶| 黄网站色视频无遮挡免费观看| 欧美日韩成人在线一区二区| 亚洲国产成人一精品久久久| 日日撸夜夜添| 成人综合一区亚洲| 国产精品免费大片| 国产精品国产av在线观看| 丰满迷人的少妇在线观看| 国产成人午夜福利电影在线观看| 久久综合国产亚洲精品| 久久精品人人爽人人爽视色| av不卡在线播放| 母亲3免费完整高清在线观看 | 中文字幕精品免费在线观看视频 | 国产精品国产三级国产专区5o| 两性夫妻黄色片 | av播播在线观看一区| 国产精品蜜桃在线观看| 国产成人91sexporn| 亚洲欧美色中文字幕在线| 自拍欧美九色日韩亚洲蝌蚪91| 久久精品熟女亚洲av麻豆精品| 国产片内射在线| 国产精品国产av在线观看| 制服人妻中文乱码| 岛国毛片在线播放| 欧美人与性动交α欧美软件 | a级毛色黄片| 欧美日韩成人在线一区二区| 亚洲精品456在线播放app| 中文欧美无线码| 男人舔女人的私密视频| 亚洲一级一片aⅴ在线观看| 精品视频人人做人人爽| 青春草亚洲视频在线观看| 精品少妇黑人巨大在线播放| 欧美日韩国产mv在线观看视频| 肉色欧美久久久久久久蜜桃| 97人妻天天添夜夜摸| 国产精品偷伦视频观看了| 午夜精品国产一区二区电影| 亚洲一级一片aⅴ在线观看| 午夜福利乱码中文字幕| 天美传媒精品一区二区| 男女免费视频国产| 亚洲综合精品二区| 欧美3d第一页| 搡女人真爽免费视频火全软件| 涩涩av久久男人的天堂| 多毛熟女@视频| 亚洲久久久国产精品| 亚洲成人手机| 交换朋友夫妻互换小说| 成人二区视频| 我的女老师完整版在线观看| 少妇的逼水好多| xxxhd国产人妻xxx| 全区人妻精品视频| 秋霞在线观看毛片| 亚洲欧洲国产日韩| 美女中出高潮动态图| 一级毛片电影观看| 亚洲中文av在线| 99九九在线精品视频| 边亲边吃奶的免费视频| 蜜桃国产av成人99| 美女国产视频在线观看| 午夜视频国产福利| 久久 成人 亚洲| 香蕉丝袜av| 少妇 在线观看| 少妇被粗大猛烈的视频| 巨乳人妻的诱惑在线观看| av天堂久久9| 看十八女毛片水多多多| 精品一区二区三区视频在线| 成人手机av| 久久精品国产a三级三级三级| 又粗又硬又长又爽又黄的视频| 免费在线观看完整版高清| 亚洲 欧美一区二区三区| 大陆偷拍与自拍| 狂野欧美激情性bbbbbb| 欧美bdsm另类| 成人综合一区亚洲| 久久久久久久精品精品| 日本-黄色视频高清免费观看| 最近的中文字幕免费完整| 最黄视频免费看| 国产成人精品一,二区| 亚洲国产精品一区三区| 色网站视频免费| 51国产日韩欧美| 最新中文字幕久久久久| 亚洲精品久久成人aⅴ小说| 午夜免费鲁丝| 国产熟女欧美一区二区| 久久99热6这里只有精品| 少妇人妻久久综合中文| 另类亚洲欧美激情| 美女内射精品一级片tv| 亚洲精品乱码久久久久久按摩| 精品99又大又爽又粗少妇毛片| 国产精品人妻久久久影院| 中国美白少妇内射xxxbb| 好男人视频免费观看在线| 天天躁夜夜躁狠狠久久av| 国产欧美日韩一区二区三区在线| 又粗又硬又长又爽又黄的视频| 99re6热这里在线精品视频| 一本大道久久a久久精品| 日产精品乱码卡一卡2卡三| 亚洲精品一区蜜桃| 搡老乐熟女国产| 色94色欧美一区二区| 黄色视频在线播放观看不卡| 日韩精品有码人妻一区| 精品人妻在线不人妻| 亚洲精品av麻豆狂野| 亚洲av日韩在线播放| 999精品在线视频| 国产精品人妻久久久影院| 欧美激情国产日韩精品一区| 国内精品宾馆在线| 日本午夜av视频| 久久久久久久亚洲中文字幕| 人人妻人人添人人爽欧美一区卜| 中文字幕另类日韩欧美亚洲嫩草| 啦啦啦啦在线视频资源| 国产精品人妻久久久影院| 亚洲在久久综合| 久久国内精品自在自线图片| 国产 一区精品| 国产免费又黄又爽又色| 国产av国产精品国产| 纵有疾风起免费观看全集完整版| 欧美成人精品欧美一级黄| 大片电影免费在线观看免费| 国产一区二区激情短视频 | 中文欧美无线码| 免费观看性生交大片5| 午夜免费观看性视频| 一本色道久久久久久精品综合| 久久久国产欧美日韩av| 国产精品成人在线| 久久久久精品性色| 国产精品国产三级专区第一集| 美女视频免费永久观看网站| 久久精品久久久久久噜噜老黄| 中文字幕亚洲精品专区| 久久精品国产自在天天线| 亚洲伊人色综图| 美女国产视频在线观看| 亚洲精品一区蜜桃| 亚洲成国产人片在线观看| 男人添女人高潮全过程视频| 777米奇影视久久| 成年动漫av网址| 激情五月婷婷亚洲| 午夜福利,免费看| 久久精品久久精品一区二区三区| 伦理电影大哥的女人| 久久国产亚洲av麻豆专区| 亚洲成色77777| av不卡在线播放| 午夜福利在线观看免费完整高清在| 国产片特级美女逼逼视频| 天天操日日干夜夜撸| 国产日韩欧美视频二区| 婷婷色综合大香蕉| 18禁观看日本| 精品久久久久久电影网| 免费黄频网站在线观看国产| 有码 亚洲区| 午夜福利,免费看| 久久 成人 亚洲| 九九在线视频观看精品| av女优亚洲男人天堂| 18在线观看网站| 国产成人欧美| 国产精品国产三级专区第一集| 熟女人妻精品中文字幕| 边亲边吃奶的免费视频| 国产黄色视频一区二区在线观看| 日本与韩国留学比较| 国产成人欧美| 亚洲av.av天堂| 秋霞伦理黄片| 久久青草综合色| 人体艺术视频欧美日本| av女优亚洲男人天堂| 国产极品天堂在线| 欧美日韩亚洲高清精品| av免费在线看不卡| h视频一区二区三区| 两个人看的免费小视频| 伦理电影大哥的女人| av电影中文网址| 1024视频免费在线观看| 99热这里只有是精品在线观看| 亚洲色图 男人天堂 中文字幕 | 亚洲伊人色综图| 国产成人精品婷婷| 免费在线观看完整版高清| 久久精品久久精品一区二区三区| 99视频精品全部免费 在线| av女优亚洲男人天堂| 大香蕉久久成人网| 国产欧美日韩一区二区三区在线| 黑人高潮一二区| 国产白丝娇喘喷水9色精品|