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

    基于分布式雷達(dá)的寬帶脈沖三維測(cè)距機(jī)制及方法研究

    2015-07-05 16:46:54巍齊丁赤飚張履謙賀柏森
    電子與信息學(xué)報(bào) 2015年3期
    關(guān)鍵詞:測(cè)距分布式天線

    李 巍齊 巍*丁赤飚張履謙賀柏森

    ①(北京跟蹤與通信技術(shù)研究所 北京 100094)

    ②(中國(guó)科學(xué)院電子學(xué)研究所 北京 100190)

    ③(中國(guó)航天科技集團(tuán)公司 北京 100048)

    基于分布式雷達(dá)的寬帶脈沖三維測(cè)距機(jī)制及方法研究

    李 ?、冽R 巍*①丁赤飚②張履謙③賀柏森①

    ①(北京跟蹤與通信技術(shù)研究所 北京 100094)

    ②(中國(guó)科學(xué)院電子學(xué)研究所 北京 100190)

    ③(中國(guó)航天科技集團(tuán)公司 北京 100048)

    針對(duì)脈沖雷達(dá)相位測(cè)距技術(shù)在相位解模糊和測(cè)量機(jī)理方面存在的局限性,該文提出一種基于分布式雷達(dá)的3維測(cè)距機(jī)制,建立分布式雷達(dá)3維測(cè)距信號(hào)模型。利用天線空間布局引起的序列相位差代替常規(guī)由序列脈沖累積引起的相位差,結(jié)合微波3維成像處理的基本原理,通過(guò)相干積累實(shí)現(xiàn)空中觀測(cè)目標(biāo)的3維測(cè)量和軌跡測(cè)量,并利用觀測(cè)區(qū)域的目標(biāo)空間稀疏分布特性,給出一種基于CLEAN的寬帶脈沖雷達(dá)高精度3維測(cè)距方法。仿真驗(yàn)證和分析結(jié)果表明,該文設(shè)計(jì)的分布式雷達(dá)3維測(cè)距機(jī)制和方法不僅較好地避免了高動(dòng)態(tài)遠(yuǎn)距離測(cè)量環(huán)境中的相位解模糊問(wèn)題,而且還能快速獲取觀測(cè)目標(biāo)的序列3維圖像和位置信息,有效避免了傳統(tǒng)單站脈沖相位測(cè)量體制只能實(shí)現(xiàn)1維測(cè)距的能力,為雷達(dá)測(cè)距應(yīng)用提供一種新的方法與技術(shù)支撐。

    分布式雷達(dá);寬帶脈沖3維測(cè)距;3維成像;CLEAN算法

    1 引言

    寬帶雷達(dá)相位測(cè)距可以實(shí)現(xiàn)傳統(tǒng)脈沖雷達(dá)系統(tǒng)的遠(yuǎn)距離相位測(cè)距和多普勒測(cè)速,并能夠在寬帶信號(hào)形式下實(shí)現(xiàn)目標(biāo)微動(dòng)特征的獲取,為脈沖雷達(dá)高精度測(cè)距和高動(dòng)態(tài)目標(biāo)特征分析提供了一種有效的觀測(cè)方式和技術(shù)途徑,極大地拓展了雷達(dá)相位測(cè)距系統(tǒng)的性能和應(yīng)用領(lǐng)域[1]。目前,較為成熟的相位測(cè)距方法有寬窄帶游標(biāo)測(cè)距技術(shù)[2],多基線干涉測(cè)量技術(shù)[3,4]以及多頻連續(xù)波處理技術(shù)等。2005年,Raytheon公司提出在寬帶雷達(dá)上利用相位導(dǎo)出距離 (Phase Derived Range, PDR)技術(shù)實(shí)現(xiàn)高精度雷達(dá)相位測(cè)距[5],并在彈道中段進(jìn)動(dòng)目標(biāo)上進(jìn)行驗(yàn)證。近年來(lái),國(guó)內(nèi)大量研究機(jī)構(gòu)在脈沖雷達(dá)相位測(cè)距技術(shù)方面也開(kāi)展了廣泛的研究,文獻(xiàn)[6]和文獻(xiàn)[7]分別利用窄帶游標(biāo)測(cè)距原理開(kāi)展了相位模糊距離解算和目標(biāo)微動(dòng)特性提取研究;文獻(xiàn)[8]給出了基于頻率步進(jìn)連續(xù)波寬帶雷達(dá)的目標(biāo)微動(dòng)特征估計(jì)方法。

    上述研究工作多集中在利用目標(biāo)多普勒信息進(jìn)行相位解模糊的問(wèn)題上,且主要是基于單站測(cè)量機(jī)制的相位測(cè)距技術(shù),通過(guò)FFT、匹配濾波等經(jīng)典處理方法實(shí)現(xiàn)多普勒頻率提取和相位解模糊。由于其固有的相位解模糊和測(cè)量機(jī)理限制,常規(guī)單站相位測(cè)距方法主要通過(guò)相位解模糊實(shí)現(xiàn)高精度相對(duì)距離測(cè)量,而對(duì)于動(dòng)態(tài)目標(biāo),特別是伴隨著目標(biāo)高速不規(guī)則機(jī)動(dòng)特性和姿態(tài)起伏特性的影響,通過(guò)序列相位進(jìn)行積分所獲得的目標(biāo)對(duì)象位置精度受到影響較大?,F(xiàn)有測(cè)量體制下,雷達(dá)系統(tǒng)獲取的目標(biāo)回波信息所包含的空間信息單一,加之系統(tǒng)誤差、測(cè)量噪聲等因素均導(dǎo)致現(xiàn)有的相位解模糊方法實(shí)現(xiàn)難度增加,給寬帶雷達(dá)的相位測(cè)距性能和應(yīng)用帶來(lái)一定的局限性,也使得高速機(jī)動(dòng)目標(biāo)的3維測(cè)量與跟蹤顯得更為困難。

    針對(duì)上述影響因素,本文提出基于統(tǒng)一時(shí)間和頻率基準(zhǔn)源的分布式雷達(dá)的3維測(cè)距機(jī)制,采用分布式布局的多副天線實(shí)現(xiàn)發(fā)射和接收多組脈沖信號(hào)的測(cè)距體制獲取觀測(cè)目標(biāo)的回波信號(hào),分析了分布式天線布局方式與信號(hào)模型;在此基礎(chǔ)上,通過(guò)收發(fā)天線間可預(yù)先精確標(biāo)定的物理基線分布關(guān)系實(shí)現(xiàn)觀測(cè)目標(biāo)3維位置實(shí)時(shí)/近實(shí)時(shí)估計(jì)與跟蹤,設(shè)計(jì)了基于雙站幾何的3維極坐標(biāo)格式成像方法,并引入CLEAN方法進(jìn)行圖像模糊和副瓣抑制,給出了基于分布式雷達(dá)的寬帶脈沖3維測(cè)距機(jī)制處理流程,實(shí)現(xiàn)對(duì)觀測(cè)目標(biāo)的高精度3維測(cè)距與跟蹤;最后,通過(guò)仿真實(shí)驗(yàn)驗(yàn)證了本文測(cè)距機(jī)制和處理方法的有效性,并討論了基于分布式雷達(dá)的寬帶脈沖測(cè)距機(jī)制所面臨的技術(shù)難題和需要進(jìn)一步開(kāi)展的研究工作。

    2 分布式雷達(dá)3維測(cè)距機(jī)制

    2.1 測(cè)量幾何與天線構(gòu)型

    設(shè)基于分布式雷達(dá)的3維空間坐標(biāo)系為OXYZ,目標(biāo)飛行高度為H,目標(biāo)沿任意方向以速度V飛行(稱(chēng)為質(zhì)心平動(dòng)),其3維測(cè)距系統(tǒng)的觀測(cè)幾何如圖1(a)所示。該機(jī)制采用多發(fā)多收(MIMO)或單發(fā)多收(SIMO)[9]天線構(gòu)型結(jié)構(gòu),利用地面上分布式布局的多個(gè)同源脈沖雷達(dá),同時(shí)多通道發(fā)射多個(gè)正交脈沖信號(hào)或單通道發(fā)射脈沖信號(hào),并同時(shí)接收同一觀測(cè)目標(biāo)的回波數(shù)據(jù)。常規(guī)的相位測(cè)距主要通過(guò)收發(fā)多組脈沖信號(hào)實(shí)現(xiàn)觀測(cè)目標(biāo)距離測(cè)量,基于分布式雷達(dá)的寬帶脈沖測(cè)距機(jī)制主要通過(guò)增加天線空間分布冗余,且同時(shí)收發(fā)多組相干脈沖信號(hào),在幾個(gè)乃至單個(gè)脈沖重復(fù)周期內(nèi)快速實(shí)現(xiàn)目標(biāo)空間3維定位和位置測(cè)量,在此基礎(chǔ)上,進(jìn)行目標(biāo)跟蹤和軌跡預(yù)測(cè)。

    為便于描述,本文采用SIMO收發(fā)體制。如圖1(b)所示,分布式布局的多副天線分布在不同半徑的圓周上,發(fā)射天線位于圓周中心,通過(guò)收發(fā)天線組合可以虛擬多個(gè)等效相位中心,設(shè)接收天線單元數(shù)為N個(gè),接收天線所在位置的最大直徑為RAM,則單個(gè)脈沖發(fā)射,其它天線同時(shí)接收可以實(shí)現(xiàn)N個(gè)等效相位中心,等效相位中心最大直徑近似為RAM/2。由于天線陣元較寬,目標(biāo)觀測(cè)范圍較大,觀測(cè)目標(biāo)的飛行方向和預(yù)期軌跡均未知,為了避免目標(biāo)測(cè)距和3維定位中存在模糊,收發(fā)天線布局中,盡可能保證至少存在一組圓形陣列能夠?qū)崿F(xiàn)對(duì)觀測(cè)目標(biāo)范圍無(wú)模糊粗精度定位,也即要求至少存在一組圓形布局天線滿足奈奎斯特采樣定理,設(shè)靠近發(fā)射天線單元的最內(nèi)側(cè)接收天線分布半徑為Ra,則相應(yīng)的天線角度間隔Δθ≤λ/(2Ra),則內(nèi)側(cè)接收天線的個(gè)數(shù)下限為表示上取整,外側(cè)接收天線若滿足奈奎斯特采樣定理,則需要放置的天線數(shù)量非常多,實(shí)際應(yīng)用中難以實(shí)現(xiàn),需要稀疏布置,文中將對(duì)這種分布引起的模糊問(wèn)題進(jìn)行討論。因此,在OXY平面上沿半徑方向上布局2個(gè)稀疏圓形天線陣列,內(nèi)側(cè)圓周布局天線決定了觀測(cè)目標(biāo)的范圍,而外側(cè)圓周布局天線則決定了觀測(cè)目標(biāo)的3維位置精度。

    圖1 分布式同源雷達(dá)3維測(cè)距系統(tǒng)幾何示意圖

    2.2 回波信號(hào)建模

    由圖1可知,在確保系統(tǒng)高精度時(shí)間同步和收發(fā)天線相對(duì)位置的毫米級(jí)精度標(biāo)定的前提下,基于分布式天線的3維測(cè)距機(jī)制能夠比利用單站進(jìn)行測(cè)距具有更好的觀測(cè)冗余。設(shè)發(fā)射天線的坐標(biāo)為(xT=0,yT=0,zT=0),接收天線的坐標(biāo)為(xR= RRcosθR,yR=RRsinθR,zR=0),目標(biāo)3維坐標(biāo)為(xn=rncosφncosθn,yn=rncosφnsinθn,zn=rnsinφn+H),rn定義為目標(biāo)到z=H平面原點(diǎn)(0,0,H)的距離,φn為與z=H平面的夾角,θn為與X軸之間夾角。系統(tǒng)采用調(diào)頻脈沖信號(hào)作為發(fā)射信號(hào):

    式中,fc為載頻,Kr為脈沖信號(hào)調(diào)頻率,T為脈沖寬度,t表示以發(fā)射時(shí)刻為起點(diǎn)的快時(shí)間,=mt′tt+表示第m個(gè)脈沖對(duì)應(yīng)的觀測(cè)時(shí)間,tm=m×PRT(脈沖重復(fù)時(shí)間)為慢時(shí)間。假定目標(biāo)散射保持各向同性,相應(yīng)地,接收天線采集到的回波信號(hào)可表示為

    式中,δn(xn,yn,zn)為目標(biāo)散射系數(shù),V為目標(biāo)速度,c為光速,R(tm)為天線等效相位中心與目標(biāo)之間的距離:

    式中,R0T=H,=+H2, cosφR=RR/R0R, sinφR=H/R0R,rT表示發(fā)射天線相位中心到目標(biāo)的距離,rR表示目標(biāo)到各個(gè)接收天線相位中心的距離。常規(guī)相位測(cè)距需要處理多個(gè)脈沖回波數(shù)據(jù),通過(guò)相位解模糊進(jìn)行距離測(cè)量,使得目標(biāo)的機(jī)動(dòng)性不可忽略。本文提出的測(cè)距方法的優(yōu)勢(shì)是在單次脈沖重復(fù)周期內(nèi)實(shí)現(xiàn)目標(biāo)高精度測(cè)距,此時(shí)對(duì)單個(gè)脈沖而言,由目標(biāo)運(yùn)動(dòng)引起的相位變化量相對(duì)每個(gè)接收天線的相位變化情況近似一致,由此引起的位置變化可以不予考慮。因此,回波信號(hào)經(jīng)接收機(jī)下變頻處理后,式(2)可以簡(jiǎn)化為

    3 分布式雷達(dá)3維測(cè)距數(shù)據(jù)處理方法

    分布式雷達(dá)3維數(shù)據(jù)處理流程如圖2所示。

    在通過(guò)常規(guī)手段獲得高度平面2維復(fù)圖像后,首先利用內(nèi)側(cè)接收?qǐng)A周天線陣列采用滿陣構(gòu)型獲取觀測(cè)目標(biāo)的無(wú)模糊低分辨3維率圖像。理論上,通過(guò)增加天線布局的空間范圍一定程度上可以實(shí)現(xiàn)更高精度3維位置測(cè)量,但有限的天線布局單元會(huì)引起較為嚴(yán)重的圖像模糊問(wèn)題,不利于實(shí)際目標(biāo)的位置估計(jì)。為此,在現(xiàn)有頻域和時(shí)域微波3維成像算法基礎(chǔ)上,需進(jìn)一步考慮目標(biāo)區(qū)域大小和處理效率實(shí)現(xiàn)高分辨率3維成像,同時(shí),還需要考慮由于外側(cè)接收?qǐng)A周陣列天線稀疏采樣引起的3維圖像模糊和高副瓣問(wèn)題。本文給出一種基于極坐標(biāo)格式的3維成像算法,并結(jié)合預(yù)先精確標(biāo)定的物理基線分布進(jìn)行3維點(diǎn)擴(kuò)展函數(shù)計(jì)算,引入CLEAN方法實(shí)現(xiàn)高分辨率無(wú)模糊3維圖像獲取,進(jìn)而進(jìn)行目標(biāo)位置提取和距離估計(jì)。

    3.1 目標(biāo)3維成像處理

    由于系統(tǒng)采用調(diào)頻脈沖信號(hào),首先需要對(duì)回波信號(hào)進(jìn)行距離向脈沖壓縮處理,通過(guò)距離向匹配濾波在距離頻域中完成,經(jīng)匹配濾波后,回波信號(hào)可表示為

    圖 2 基于分布式雷達(dá)3維測(cè)距數(shù)據(jù)處理流程

    為方便算法實(shí)現(xiàn)中參數(shù)的計(jì)算和估計(jì),令

    則式(3)中的Tr和Rr可表示為

    式(7)中的距離近似將會(huì)帶來(lái)近似誤差,成像時(shí)一般要求近似誤差的范圍控制在1/4波長(zhǎng)即λ/4。以第4節(jié)的仿真參數(shù)為例,式(8)中R0T=30000,近似引起的誤差在10-7m量級(jí),遠(yuǎn)小于λ/4,可以進(jìn)行上述近似處理。為獲取觀測(cè)目標(biāo)的3維圖像,首先進(jìn)行目標(biāo)距離的粗估計(jì),而后在有效觀測(cè)范圍內(nèi)沿高度向逐層成像,因此,R0T和R0R可以表征與高度相關(guān)的柱坐標(biāo)系中,對(duì)高度平面z(k)=H-kΔz, k=0,1,…,K (觀測(cè)區(qū)域高度Z0=KΔz)上的目標(biāo),以z=z(k)=zn平面中心點(diǎn)(0,0,zn)為參考點(diǎn)建立參考函數(shù)H1(Kω,z(k )),則有:

    根據(jù)式(3)和式(6)對(duì)應(yīng)的關(guān)系可得

    對(duì)式(4)中的SR(Kω)乘以參考函數(shù)H1,則有

    令rnp=rncosφn=rnsin αT,表示rn在成像平面的投影。遠(yuǎn)場(chǎng)條件下,很容易滿足條件R0T?和;同時(shí)分母rR+rR0為關(guān)于Rθ的周期函數(shù),為提高相位補(bǔ)償精度,可進(jìn)行以下近似

    根據(jù)式(10),S1(Kω)可近似表示為

    由于在特定高度平面,式(13)中二次項(xiàng)隨目標(biāo)位置發(fā)生變化,需要進(jìn)一步去除隨目標(biāo)半徑變化的影響。為此,先對(duì)式(13)進(jìn)行適當(dāng)?shù)淖兞刻鎿Q,該替換過(guò)程也即極坐標(biāo)處理過(guò)程的數(shù)學(xué)描述為

    式中,

    此時(shí),通過(guò)2維極坐標(biāo)插值將極坐標(biāo)格式信號(hào)S2(Kωxy,θ;Kx,Ky)轉(zhuǎn)換為直角坐標(biāo)系的信號(hào)S3(Kx, Ky),相應(yīng)的插值過(guò)程如下:

    式中,NIθT和NIKω分別為沿θT和Kω方向上的插值核長(zhǎng)度,由于成像場(chǎng)景區(qū)域rn較小,rn2項(xiàng)的影響可以忽略,極端情況rn=0時(shí)沒(méi)有影響,通過(guò)2維逆傅里葉變換即可獲得高度z(k)平面上的圖像。

    通過(guò)上述處理能夠重建特定高度平面上的目標(biāo),在此基礎(chǔ)上令k=k+1,循環(huán)執(zhí)行上述各步驟,即可獲得整個(gè)空域內(nèi)目標(biāo)的3維圖像。成像時(shí)k的步長(zhǎng)設(shè)置為略小于高程向分辨率,k的范圍根據(jù)回波信號(hào)雙程歷程rT+rR提供初始搜索范圍。當(dāng)需要循環(huán)計(jì)算的次數(shù)很多時(shí),可以先增大k的步長(zhǎng)進(jìn)行成像,根據(jù)偏離目標(biāo)真實(shí)高度引起的環(huán)狀半徑,縮小高度向成像范圍,在此范圍內(nèi)再用小步長(zhǎng)的k成像,以減少計(jì)算量。

    3.2 基于CLEAN的目標(biāo)3維位置提取

    系統(tǒng)采用兩個(gè)圓周陣列獲取低分辨率3維圖像I1(x,y,z)和高分辨率模糊3維圖像I2(x,y,z),其中低分辨率圖像位置精度較差,而高分辨率圖像中又存在模糊和圓周陣列固有的高副瓣問(wèn)題,因此,需要給出模糊和副瓣抑制方法。在圖像模糊抑制方面,文獻(xiàn)[10]最早提出了CLEAN算法用于消除空間頻率覆蓋不完整而產(chǎn)生的影響,改善綜合孔徑射電天文觀測(cè)圖像的質(zhì)量。隨后文獻(xiàn)[11]對(duì)CLEAN算法進(jìn)行了深入的數(shù)學(xué)理論分析與探討。文獻(xiàn)[12]利用CLEAN算法進(jìn)行微波成像中圖像的旁瓣抑制。近年來(lái),有學(xué)者也利用其進(jìn)行機(jī)載毫米波綜合孔徑成像和空間目標(biāo)ISAR成像中的旁瓣抑制[13]。文獻(xiàn)[14]基于最大似然估計(jì)和CLEAN算法對(duì)稀疏陣列的3維ISAR成像進(jìn)行了研究,同時(shí)利用單個(gè)天線下目標(biāo)的ISAR圖像以解決高副瓣的問(wèn)題。

    本文通過(guò)CLEAN消除由于內(nèi)外圓周陣列天線成像中的目標(biāo)高副瓣問(wèn)題、外側(cè)圓周陣列天線成像中的圖像模糊和高副瓣問(wèn)題等對(duì)3維圖像質(zhì)量的影響,其主要過(guò)程是提取當(dāng)前圖像中的強(qiáng)點(diǎn)目標(biāo),計(jì)算收發(fā)系統(tǒng)對(duì)應(yīng)到強(qiáng)點(diǎn)目標(biāo)的響應(yīng)函數(shù),并從圖像中減去相應(yīng)強(qiáng)點(diǎn)目標(biāo)的響應(yīng)函數(shù),去除強(qiáng)點(diǎn)所對(duì)應(yīng)的模糊和旁瓣響應(yīng),通過(guò)多次迭代最終獲得質(zhì)量良好的3維圖像??紤]到CLEAN算法中迭代比例因子γ對(duì)圖像質(zhì)量會(huì)產(chǎn)生影響,在對(duì)外側(cè)圓周陣列天線獲取的3維圖像進(jìn)行迭代的次數(shù)控制上以目標(biāo)位置點(diǎn)不超過(guò)內(nèi)側(cè)圓周3維圖像有效區(qū)域?yàn)闇?zhǔn)。為實(shí)現(xiàn)對(duì)觀測(cè)目標(biāo)的高精度3維測(cè)距,需要聯(lián)合內(nèi)外側(cè)圓周陣列天線獲取的3維圖像進(jìn)行CLEAN處理,具體處理步驟為:

    步驟1 對(duì)低分辨率3維圖像I10(x,y,z)= I1(x,y,z)進(jìn)行最強(qiáng)散射點(diǎn)目標(biāo)提取,獲得其目標(biāo)的幅相信息A1(x1,y1,z1)ejα1,以及空間位置坐標(biāo)(x1,y1,z1);

    步驟 2 根據(jù)預(yù)先精確標(biāo)定的內(nèi)側(cè)圓周陣列布局基線分布關(guān)系進(jìn)行3維點(diǎn)擴(kuò)展函數(shù)計(jì)算,計(jì)算過(guò)程與3維成像處理過(guò)程類(lèi)似,對(duì)應(yīng)的3維點(diǎn)擴(kuò)展函數(shù)記為PSF1(x,y,z);

    步驟 3 對(duì)圖像I1(x,y,z)進(jìn)行CLEAN操作生成新的3維圖像,并分別記錄當(dāng)前強(qiáng)目標(biāo)的散射幅相信息A1(x1,y1,z1)ejα1以及空間位置坐標(biāo)(x1,y1,z1), I11(x, y,z)=I10(x,y,z)-A1(x1,y1,z1)ejα1×PSF1(x1,y1,z1);

    步驟 4 更新迭代輸入圖像I1(x,y,z)=I10(x,y, z),重復(fù)步驟1至步驟3,直到圖像剩下部分的最大值比預(yù)先定義的噪聲電平小,或者達(dá)到某個(gè)特定的迭代數(shù)字,記錄當(dāng)前獲得的目標(biāo)位置信息及I1C(x,y,z);

    步驟5 對(duì)輸入的存在模糊的高分辨率3維圖像I20(x,y,z)=I2(x,y,z )進(jìn)行最強(qiáng)散射點(diǎn)目標(biāo)提取,獲得其目標(biāo)的幅相信息A2(x2,y2,z2)ejα2以及空間位置坐標(biāo)(x2,y2,z2);

    步驟6 根據(jù)預(yù)先精確標(biāo)定的外側(cè)圓周陣列布局基線分布關(guān)系進(jìn)行3維點(diǎn)擴(kuò)展函數(shù)PSF2(x,y,z)計(jì)算;

    步驟7 對(duì)圖像I2(x,y,z)進(jìn)行CLEAN操作生成新的3維圖像,并分別記錄當(dāng)前強(qiáng)目標(biāo)的散射幅相信息A2(x2,y2,z2)ejα2以及位置坐標(biāo)(x2,y2,z2),

    步驟8 更新迭代輸入圖像I2(x,y,z)=I21(x,y, z),重復(fù)步驟5至步驟7,直到(x2,y2,z2)超出I1C(x,y,z)所表示的目標(biāo)范圍,則終止迭代過(guò)程,并記錄當(dāng)前獲得的目標(biāo)位置信息及I2C(x,y,z)。

    在獲取目標(biāo)3維微波圖像的基礎(chǔ)上,通過(guò)CLEAN實(shí)現(xiàn)垂直于波傳播方向平面的目標(biāo)增強(qiáng),進(jìn)而獲取觀測(cè)目標(biāo)對(duì)應(yīng)不同脈沖相對(duì)雷達(dá)坐標(biāo)系原點(diǎn)的3維微波圖像I2C(x,y,z),在此基礎(chǔ)上可以通過(guò)提取目標(biāo)的3維位置坐標(biāo)實(shí)現(xiàn)目標(biāo)測(cè)距。由于空中或空間目標(biāo)分布相對(duì)較為稀疏,目標(biāo)提取的過(guò)程可以采用恒虛警CFAR進(jìn)行,主要分為兩步:第1步,根據(jù)目標(biāo)所處周?chē)尘半s波的統(tǒng)計(jì)特性自適應(yīng)地求取檢測(cè)閾值,目標(biāo)周?chē)尘半s波的統(tǒng)計(jì)特性通常由目標(biāo)像素周?chē)?維參考窗口內(nèi)的像素確定,然后將待檢測(cè)像素和自適應(yīng)閾值進(jìn)行比較,判斷其是否為強(qiáng)目標(biāo);第2步,通過(guò)3維參考窗口滑動(dòng),實(shí)現(xiàn)對(duì)所有像素的自適應(yīng)檢測(cè)。通過(guò)CFAR檢測(cè)確定強(qiáng)目標(biāo)及其相應(yīng)的3維位置(xmn(tm),ymn(tm),zmn(tm)),其中,m表示第m個(gè)脈沖,n表示強(qiáng)目標(biāo)編號(hào),且n=1,2,…,N, N為檢測(cè)到的強(qiáng)目標(biāo)總數(shù)。相應(yīng)地,相對(duì)雷達(dá)中心位置的目標(biāo)距離Rsl可以表示為

    由式(18)可知,目標(biāo)相對(duì)雷達(dá)中心的距離可以表示為發(fā)射脈沖時(shí)間、第m個(gè)脈沖時(shí)間內(nèi)所處的位置(xmn(tm),ymn(tm),zmn(tm))以及雷達(dá)坐標(biāo)位置的函數(shù),單個(gè)脈沖內(nèi)即能實(shí)現(xiàn)目標(biāo)高精度距離解算,從而避免了常規(guī)測(cè)距方法需要的相位積累過(guò)程,在此基礎(chǔ)上利用序列3維圖像的相位信息還可進(jìn)一步提高測(cè)距精度。

    4 算法仿真與分析

    為檢驗(yàn)本文3維測(cè)距機(jī)制的正確性和處理性能,針對(duì)典型動(dòng)態(tài)目標(biāo)進(jìn)行仿真測(cè)試。系統(tǒng)工作頻率為5 GHz,信號(hào)帶寬為500 MHz,采樣率為600 MHz,脈沖寬度為50 sμ,共設(shè)置1個(gè)發(fā)射天線位于圓周中心點(diǎn),內(nèi)側(cè)接收?qǐng)A周陣列天線半徑為10 m,天線單元數(shù)為400;外側(cè)接收?qǐng)A周陣列天線半徑為2000 m,天線單元數(shù)為600,且天線等角度間距分布;單次發(fā)射信號(hào)后,600個(gè)天線單元同時(shí)接收回波信號(hào)。共設(shè)置3個(gè)機(jī)動(dòng)目標(biāo),其中1個(gè)機(jī)動(dòng)目標(biāo)長(zhǎng)度為4 m,有2個(gè)等效散射點(diǎn),速度為220 m/s,目標(biāo)初始距離為30 km,等效散射中心坐標(biāo)分別為(0,0,30000) m和(4,0,30000) m,機(jī)動(dòng)目標(biāo)在運(yùn)動(dòng)過(guò)程中的空間位置分布情況如圖3(a)所示。

    利用3維雙站極坐標(biāo)格式成像處理方法,可以獲得對(duì)應(yīng)內(nèi)側(cè)圓周天線陣列的低分辨率無(wú)模糊3維圖像,還可以獲得對(duì)應(yīng)外側(cè)圓周天線陣列的模糊的高分辨率3維圖像,如圖3(b), 3(c)所示,由于受到高副瓣和模糊的影響,若直接利用這兩幅圖像進(jìn)行目標(biāo)的位置估計(jì),其估計(jì)精度較差。本文首先對(duì)內(nèi)側(cè)天線陣列的低分辨率無(wú)模糊3維圖像進(jìn)行CLEAN處理,消除目標(biāo)副瓣影響,獲得目標(biāo)分布;在此基礎(chǔ)上,對(duì)外側(cè)天線陣列的高分辨率模糊3維圖像進(jìn)行CLEAN處理,獲得目標(biāo)的精確3維圖像及位置信息,如圖4所示。

    圖3 圓周天線陣列獲得某一脈沖內(nèi)機(jī)動(dòng)目標(biāo)的3維圖像

    圖4 通過(guò)CLEAN處理獲得目標(biāo)高分辨率無(wú)模糊3維圖像剖面顯示效果圖

    圖5和圖6分別給出了觀測(cè)目標(biāo)在運(yùn)動(dòng)過(guò)程中,通過(guò)序列3維圖像解算得到的觀測(cè)目標(biāo)9.1 s高動(dòng)態(tài)運(yùn)動(dòng)過(guò)程中測(cè)距結(jié)果和測(cè)量誤差統(tǒng)計(jì)結(jié)果。由于觀測(cè)時(shí)間較短,目標(biāo)軌跡看起來(lái)近似直線。鑒于本文采用“快照式”成像方式實(shí)現(xiàn)測(cè)距,隨著觀測(cè)時(shí)間增加,當(dāng)目標(biāo)軌跡為曲線時(shí)該算法同樣適用。從圖5可以看出,觀測(cè)目標(biāo)在每個(gè)脈沖內(nèi)的測(cè)量得到的距離信息以及目標(biāo)運(yùn)動(dòng)軌跡均與仿真理論值一致。從圖6中的誤差統(tǒng)計(jì)結(jié)果能夠看出,目標(biāo)3維分辨率優(yōu)于(1.125 m, 1.125 m, 0.3 m)。由于空間中目標(biāo)稀疏分布,對(duì)目標(biāo)進(jìn)行CLEAN處理前經(jīng)過(guò)了升采樣處理,其像素間隔最小為0.2 m,提取目標(biāo)的3維位置偏差保持在兩個(gè)像素以?xún)?nèi),相應(yīng)的測(cè)距精度可達(dá)0.7 m。上述仿真驗(yàn)證結(jié)果充分表明,本文方法利用地面分布的多個(gè)同源接收天線節(jié)點(diǎn)進(jìn)行3維成像處理,能夠精確分辨出高動(dòng)態(tài)環(huán)境下的目標(biāo)位置和運(yùn)動(dòng)信息。

    本方法獲得0.7 m測(cè)距精度為絕對(duì)測(cè)距精度,與信號(hào)帶寬、天線位置姿態(tài)、目標(biāo)的稀疏程度和目標(biāo)重建像素大小等因素有關(guān),主要是通過(guò)3維成像首先獲取目標(biāo)強(qiáng)散射點(diǎn)的3維位置,在此基礎(chǔ)上獲得對(duì)應(yīng)強(qiáng)散射點(diǎn)的實(shí)際距離,實(shí)際中目標(biāo)強(qiáng)散射點(diǎn)的位置提取與3維分辨率相關(guān),因此,絕對(duì)測(cè)距精度受限于系統(tǒng)分辨率;同時(shí),由于天線位置姿態(tài)誤差也會(huì)影響到目標(biāo)的響應(yīng)函數(shù),也會(huì)對(duì)測(cè)距精度產(chǎn)生影響;此外,空間中目標(biāo)的稀疏程度會(huì)影響到CLEAN方法的使用,目標(biāo)越稀疏,測(cè)距精度越高。基于上述分析,未來(lái)可結(jié)合目標(biāo)3維圖像所攜帶的相位信息進(jìn)一步提高目標(biāo)測(cè)距精度。

    5 結(jié)論

    圖5 觀測(cè)目標(biāo)動(dòng)態(tài)測(cè)距結(jié)果

    圖 6 觀測(cè)目標(biāo)動(dòng)態(tài)測(cè)距誤差統(tǒng)計(jì)

    本文針對(duì)高動(dòng)態(tài)大機(jī)動(dòng)目標(biāo)的高精度測(cè)量問(wèn)題,設(shè)計(jì)了一種基于分布式雷達(dá)的寬帶脈沖3維測(cè)距機(jī)制和方法。利用分布式同源布局的多個(gè)天線等效相位中心建立基于分布式雷達(dá)測(cè)量機(jī)制的3維測(cè)距信號(hào)模型,提出一種基于極坐標(biāo)格式的3維成像算法;并結(jié)合預(yù)先精確標(biāo)定的物理基線分布關(guān)系進(jìn)行3維點(diǎn)擴(kuò)展函數(shù)計(jì)算,引入CLEAN方法實(shí)現(xiàn)高分辨率無(wú)模糊3維圖像的獲取,進(jìn)而實(shí)現(xiàn)對(duì)觀測(cè)目標(biāo)的高精度3維空間位置測(cè)量。通過(guò)計(jì)算機(jī)仿真驗(yàn)證了本文3維測(cè)距機(jī)制和信號(hào)處理方法的正確性,對(duì)基于單次脈沖信號(hào)的3維測(cè)距精度和效率進(jìn)行分析比對(duì)。仿真測(cè)試結(jié)果表明,該方法能夠?qū)崿F(xiàn)動(dòng)態(tài)目標(biāo)的3維位置與距離估計(jì),在單個(gè)回波脈沖內(nèi)的快速測(cè)距精度能夠達(dá)到亞米級(jí)。未來(lái)將針對(duì)分布式地面天線布局設(shè)計(jì)、通道幅相補(bǔ)償、序列脈沖相位信息處理及其對(duì)3維定位和測(cè)距精度的影響等方面開(kāi)展深入研究,進(jìn)一步驗(yàn)證和提升分布式雷達(dá)3維測(cè)距技術(shù)的解算精度;同時(shí),本文主要基于CLEAN算法進(jìn)行了3維測(cè)量精度分析,未來(lái)可利用目標(biāo)空間稀疏分布特性,結(jié)合稀疏信號(hào)處理進(jìn)一步提升目標(biāo)3維位置測(cè)量精度,發(fā)展高精度快速求解優(yōu)化算法。

    [1] Fitzgerald R J. Range-Doppler coupling and other aspects of reentry vehicle tracking[C]. Symposium on Nonlinear Estimation Theory and Applications, 1973.

    [2] Russell R. Range vernier[R]. ADA-077178, 1979.

    [3] 李金嶺, 喬書(shū)波, 劉鸝, 等. 基于我國(guó)VLBI網(wǎng)航天器定位歸算仿真分析[J]. 宇航學(xué)報(bào), 2011, 32(11): 2333-2338.

    Li Jin-ling, Qiao Shu-bo, Liu Li, et al.. Simulation analysis of positioning reduction of spacecraft as per the Chinese VLBI network[J]. Journal of Astronautics, 2011, 32(11): 2333-2338.

    [4] 吳偉仁, 王廣利, 節(jié)德剛, 等. 基于?DOR信號(hào)的高精度VLBI技術(shù)[J]. 中國(guó)科學(xué): 信息科學(xué), 2013, 43(2): 185-196.

    Wu Wei-ren, Wang Guang-li, Jie De-gang et al.. High-accuracy VLBI technique using DOR signals[J]. SCIENCE CHINA: Information Sciences, 2013, 43(2): 185-196.

    [5] Fritz S. Process for phase-derived-range measurements[P]. United States Patent, US7046190b2, 2006.

    [6] 姚輝偉, 魏璽章, 徐少坤, 等. 彈道中段進(jìn)動(dòng)目標(biāo)非理想散射中心微動(dòng)特性研究[J]. 電子學(xué)報(bào), 2012, 40(9): 1844-1851.

    Yao Hui-wei, Wei Xi-zhang, Xu Shao-kun, et al.. Micro-motion characteristics of non-ideal scattering centers of midcourse targets with precession[J]. Acta Electronica Sinica, 2012, 40(9): 1844-1851.

    [7] 李康樂(lè). 雷達(dá)目標(biāo)微動(dòng)特征提取與估計(jì)技術(shù)研究[D]. [博士論文]國(guó)防科技大學(xué),2010.

    Li Kang-le. Research on feature extraction and parameters estimation for radar targets with micro-motions[D]. [Ph.D. dissertation], National University of Defense Technology, 2010.

    [8] 賀峰, 劉承蘭, 朱國(guó)富, 等. 基于距離和差值序列的復(fù)雜運(yùn)動(dòng)目標(biāo)微動(dòng)特征提取[J]. 電子學(xué)報(bào), 2012, 40(3): 422-428.

    He Feng, Liu Cheng-lan, Zhu Guo-fu, et al.. Micro-motion signature extraction of targets with complex motions based on summation and difference sequence of ranges[J]. Acta Electronica Sinica, 2012, 40(3): 422-428.

    [9] Fisher E, Haimovich A, Blum R S, et al.. Spatial diversity in radars-models and detection performance[J]. IEEE Transactions on Signal Processing, 2006, 54(3): 823-838.

    [10] H?gbom J A. Aperture synthesis with a non-regular distribution of interferometer baselines[J]. Astrophysics Supplement, 1974, 39(15): 417-426.

    [11] Schwarz U J. Mathematical-statistical description of the iterative beam removing technique (method CLEAN)[J]. Astronomy and Astrophysics, 1978, 65(2): 345-356.

    [12] Tao J and Steinberg B D. Reduction of side-lobe and speckle artifacts in microwave imaging: the CLEAN technique[J]. IEEE Transactions on Antennas and Propagation, 1988, 36(4): 543-556.

    [13] 何云濤, 江月松, 鐘宇. CLEAN算法在機(jī)載毫米波綜合孔徑成像中的應(yīng)用[J]. 電子與信息學(xué)報(bào), 2007, 29(7): 1756-1760.

    He Yun-tao, Jiang Yue-song ,and Zhong Yu. An application of the CLEAN algorithm in airborne mm-wave synthetic aperture radar[J]. Journal of Electronics & Information Technology, 2007, 29(7): 1756-1760.

    [14] Ma C Z, Yeo T S, Tan C S, et al.. Sparse Array 3-D ISAR imaging based on maximum likelihood estimation and CLEAN Technique[J]. IEEE Transactions on Image Processing, 2010, 19(8): 2127-2142.

    李 ?。?男,1971年生,研究員,主要研究方向?yàn)楹教鞙y(cè)控總體、雷達(dá)信號(hào)處理、目標(biāo)探測(cè)與識(shí)別等.

    齊 ?。?男,1982年生,博士,工程師,主要研究方向?yàn)楹教鞙y(cè)控、雷達(dá)信號(hào)處理、衛(wèi)星導(dǎo)航.

    丁赤飚: 男,1969年生,研究員,博士生導(dǎo)師,主要研究方向?yàn)楹铣煽讖嚼走_(dá)、雷達(dá)信號(hào)處理、遙感信息處理和應(yīng)用系統(tǒng).

    張履謙: 男,1926年生,院士,主要研究方向?yàn)楹教鞙y(cè)控、雷達(dá)信號(hào)處理、電子對(duì)抗技術(shù).

    Research on Mechanisms and Methods in the 3D Measurement with Wide-band Pulse Based on Distributed Radar

    Li Wei①Q(mào)i Wei①Ding Chi-biao②Zhang Lü-qian③He Bo-sen①①(Beijing Institute of Tracking and Telecommunication Technology, Beijing 100094, China)
    ②(Institute of Electronics, Chinese Academy of Sciences, Beijing 100190, China)
    ③(China Aerospace Science and Technology Corporation, Beijing 100048, China)

    Considering the limitations of phase unwrapping and measurement mechanism in the pulsed radar phase ranging technology, a three-dimension (3D) ranging mechanism based on distributed radar is proposed, and the observation matrix model of distributed radar 3D ranging is established. By using the sequential phase difference caused by the distributed antenna layout instead of the phase difference caused by pulse incubation, a distributed wide-band radar (3D) measurement based on observation matrix is implemented, and in combination with the target sparse distribution characteristics in the observed area, a wide-band pulse radar carrier measurement by phase based on CLEAN method is suggested. As shown by the simulation and analysis results, the distributed radar phase ranging mechanism and method in this paper not only solve the phase unwrapping problem in high dynamic and long distance measurement environment, but also acquire the three-dimension sequential images and position information of the observed target quickly. This proposed method provides a new technical support for the application of the radar ranging technology.

    Distributed radar; 3D measurement with wide-band pulse; 3D imaging; CLEAN algorithm

    TN959.6

    A

    1009-5896(2015)03-0643-08

    10.11999/JEIT140575

    2014-05-04 收到,2014-10-20改回

    國(guó)家自然科學(xué)基金(61471012)和中國(guó)博士后科學(xué)基金(2013M542455)資助課題

    *通信作者:齊巍 robyche@163.com

    猜你喜歡
    測(cè)距分布式天線
    類(lèi)星體的精準(zhǔn)測(cè)距
    科學(xué)(2020年3期)2020-01-06 04:02:51
    分布式光伏熱錢(qián)洶涌
    能源(2017年10期)2017-12-20 05:54:07
    分布式光伏:爆發(fā)還是徘徊
    能源(2017年5期)2017-07-06 09:25:54
    ETC相控陣天線與普通天線應(yīng)用對(duì)比分析
    淺談超聲波測(cè)距
    電子制作(2017年7期)2017-06-05 09:36:13
    ALLESS轉(zhuǎn)動(dòng)天線射頻旋轉(zhuǎn)維護(hù)與改造
    電子制作(2016年1期)2016-11-07 08:42:54
    理論宣講要上接天線、下接地氣
    基于DDS的分布式三維協(xié)同仿真研究
    基于PSOC超聲測(cè)距系統(tǒng)設(shè)計(jì)
    相對(duì)差分單項(xiàng)測(cè)距△DOR
    太空探索(2014年1期)2014-07-10 13:41:50
    欧美日韩瑟瑟在线播放| 免费高清视频大片| 国产伦精品一区二区三区视频9 | 88av欧美| 深夜精品福利| 国产精品日韩av在线免费观看| 亚洲中文av在线| 天堂√8在线中文| 嫩草影院精品99| 一个人免费在线观看的高清视频| 99在线视频只有这里精品首页| 久久久色成人| 三级毛片av免费| 99热这里只有是精品50| 亚洲专区国产一区二区| 日韩 欧美 亚洲 中文字幕| 欧美性猛交黑人性爽| 国产精品精品国产色婷婷| 亚洲成av人片免费观看| 欧洲精品卡2卡3卡4卡5卡区| 国产精品乱码一区二三区的特点| 波多野结衣高清无吗| 亚洲精品美女久久av网站| 中国美女看黄片| 成熟少妇高潮喷水视频| 国产精品99久久久久久久久| ponron亚洲| 国产精品,欧美在线| 亚洲黑人精品在线| 欧美一级毛片孕妇| 日韩欧美在线乱码| 999久久久国产精品视频| 色吧在线观看| 成人无遮挡网站| 亚洲男人的天堂狠狠| 亚洲一区二区三区不卡视频| 精品电影一区二区在线| 久久精品影院6| a在线观看视频网站| 欧美另类亚洲清纯唯美| 欧美不卡视频在线免费观看| 国产一区二区在线观看日韩 | 91久久精品国产一区二区成人 | 色尼玛亚洲综合影院| 两性午夜刺激爽爽歪歪视频在线观看| 国产野战对白在线观看| 九九热线精品视视频播放| 欧美色视频一区免费| 久久午夜综合久久蜜桃| 国产精品亚洲一级av第二区| 国产欧美日韩一区二区三| 欧美激情在线99| 亚洲精品456在线播放app | 成人午夜高清在线视频| 特大巨黑吊av在线直播| 免费看光身美女| 午夜影院日韩av| 听说在线观看完整版免费高清| 99久久国产精品久久久| 亚洲中文日韩欧美视频| 97碰自拍视频| 亚洲狠狠婷婷综合久久图片| 97超级碰碰碰精品色视频在线观看| 波多野结衣巨乳人妻| 一区二区三区高清视频在线| 久久久水蜜桃国产精品网| 国产精华一区二区三区| 亚洲av成人精品一区久久| 女同久久另类99精品国产91| 丰满人妻一区二区三区视频av | 午夜精品久久久久久毛片777| 精品久久久久久久久久免费视频| 精品一区二区三区视频在线观看免费| 亚洲激情在线av| 色精品久久人妻99蜜桃| 精品久久久久久久久久久久久| 欧美不卡视频在线免费观看| 在线永久观看黄色视频| 在线观看午夜福利视频| 久久久精品欧美日韩精品| 这个男人来自地球电影免费观看| av黄色大香蕉| 精品国产亚洲在线| 国产成人福利小说| 精品无人区乱码1区二区| 国产熟女xx| 村上凉子中文字幕在线| 成人午夜高清在线视频| 亚洲天堂国产精品一区在线| 黄色 视频免费看| 麻豆国产97在线/欧美| 午夜精品久久久久久毛片777| 日日摸夜夜添夜夜添小说| 又粗又爽又猛毛片免费看| 久久精品人妻少妇| 国产精品自产拍在线观看55亚洲| 欧美色欧美亚洲另类二区| 波多野结衣巨乳人妻| 亚洲精品一区av在线观看| 午夜亚洲福利在线播放| 欧美中文综合在线视频| 国产极品精品免费视频能看的| 88av欧美| 老司机午夜十八禁免费视频| 亚洲欧美日韩东京热| 国产精品香港三级国产av潘金莲| 男女视频在线观看网站免费| 九色国产91popny在线| 我的老师免费观看完整版| 亚洲中文av在线| 99久久精品国产亚洲精品| 欧美日韩国产亚洲二区| 一a级毛片在线观看| 亚洲精华国产精华精| 国产精品一区二区免费欧美| 黄色 视频免费看| a级毛片a级免费在线| 国产亚洲欧美98| 天天一区二区日本电影三级| 一进一出好大好爽视频| 一级作爱视频免费观看| 久久精品亚洲精品国产色婷小说| 国产高清激情床上av| a在线观看视频网站| 欧美最黄视频在线播放免费| 欧美+亚洲+日韩+国产| 国产精品av久久久久免费| 99久国产av精品| 免费无遮挡裸体视频| 国产精品久久久人人做人人爽| 成年人黄色毛片网站| 小蜜桃在线观看免费完整版高清| 中文字幕人成人乱码亚洲影| 久久久久九九精品影院| 真人做人爱边吃奶动态| 巨乳人妻的诱惑在线观看| 成人午夜高清在线视频| 又爽又黄无遮挡网站| 亚洲第一电影网av| 国内揄拍国产精品人妻在线| 露出奶头的视频| av女优亚洲男人天堂 | 久久久国产欧美日韩av| 麻豆国产97在线/欧美| 这个男人来自地球电影免费观看| 日本撒尿小便嘘嘘汇集6| 日韩大尺度精品在线看网址| 一区二区三区激情视频| www日本黄色视频网| 男女之事视频高清在线观看| 国产淫片久久久久久久久 | 黑人巨大精品欧美一区二区mp4| 亚洲精品在线观看二区| 亚洲国产精品sss在线观看| 黄片小视频在线播放| a级毛片a级免费在线| 欧美成人一区二区免费高清观看 | 日本免费一区二区三区高清不卡| 国产三级在线视频| 美女黄网站色视频| 国产亚洲精品久久久com| 天堂av国产一区二区熟女人妻| 亚洲乱码一区二区免费版| 免费人成视频x8x8入口观看| 一级毛片女人18水好多| 久久国产乱子伦精品免费另类| 亚洲精华国产精华精| 99热精品在线国产| 五月伊人婷婷丁香| 后天国语完整版免费观看| 99久久99久久久精品蜜桃| 美女午夜性视频免费| 日韩高清综合在线| 亚洲成av人片在线播放无| 色噜噜av男人的天堂激情| 久久九九热精品免费| netflix在线观看网站| avwww免费| 亚洲成av人片免费观看| 亚洲最大成人中文| 亚洲国产欧洲综合997久久,| 成人鲁丝片一二三区免费| 99视频精品全部免费 在线 | 日本成人三级电影网站| 国产午夜精品久久久久久| 欧美中文日本在线观看视频| 成人特级黄色片久久久久久久| 亚洲无线在线观看| 校园春色视频在线观看| 午夜视频精品福利| 91久久精品国产一区二区成人 | 一本综合久久免费| 丰满人妻一区二区三区视频av | 久久伊人香网站| 人妻夜夜爽99麻豆av| 少妇裸体淫交视频免费看高清| 99久久成人亚洲精品观看| 欧美一区二区精品小视频在线| 亚洲欧美一区二区三区黑人| 午夜精品一区二区三区免费看| 亚洲欧美日韩无卡精品| 国产亚洲av嫩草精品影院| 欧美日韩国产亚洲二区| 一二三四社区在线视频社区8| 国产单亲对白刺激| 身体一侧抽搐| 欧美zozozo另类| 中文字幕久久专区| 精品国产三级普通话版| 精品一区二区三区四区五区乱码| 黄色成人免费大全| 嫩草影院精品99| 18禁黄网站禁片午夜丰满| 人人妻人人看人人澡| 人人妻人人澡欧美一区二区| 偷拍熟女少妇极品色| 国产久久久一区二区三区| 好男人电影高清在线观看| 国产成人精品久久二区二区91| 亚洲欧美日韩高清在线视频| 久久精品aⅴ一区二区三区四区| 国产成人aa在线观看| 婷婷亚洲欧美| 中文字幕熟女人妻在线| 欧美一级a爱片免费观看看| 18禁黄网站禁片免费观看直播| 亚洲成av人片免费观看| 免费无遮挡裸体视频| 97超视频在线观看视频| 欧美色视频一区免费| 男女那种视频在线观看| av欧美777| 色综合站精品国产| 一个人观看的视频www高清免费观看 | 国产主播在线观看一区二区| 久久天堂一区二区三区四区| а√天堂www在线а√下载| 午夜福利视频1000在线观看| 在线播放国产精品三级| 国产亚洲精品综合一区在线观看| 麻豆久久精品国产亚洲av| www日本黄色视频网| 少妇裸体淫交视频免费看高清| 最近最新免费中文字幕在线| 欧美一级毛片孕妇| 精品久久久久久久人妻蜜臀av| 又紧又爽又黄一区二区| 国产伦精品一区二区三区视频9 | 色综合站精品国产| 精品久久久久久久末码| 最新美女视频免费是黄的| 国产乱人伦免费视频| 巨乳人妻的诱惑在线观看| 制服丝袜大香蕉在线| 男女视频在线观看网站免费| aaaaa片日本免费| 国产精品野战在线观看| 999久久久国产精品视频| 婷婷六月久久综合丁香| 欧美黑人巨大hd| 最新在线观看一区二区三区| 精品国内亚洲2022精品成人| 免费av不卡在线播放| 国产高清激情床上av| cao死你这个sao货| 欧美大码av| 国产乱人伦免费视频| 香蕉av资源在线| 最近最新免费中文字幕在线| 久久人妻av系列| 日韩免费av在线播放| 免费搜索国产男女视频| 一个人看视频在线观看www免费 | 可以在线观看的亚洲视频| 一级黄色大片毛片| 亚洲七黄色美女视频| 波多野结衣巨乳人妻| 97超视频在线观看视频| 国产高清videossex| 搡老岳熟女国产| 女警被强在线播放| 国产99白浆流出| cao死你这个sao货| 日本熟妇午夜| 叶爱在线成人免费视频播放| 国产av麻豆久久久久久久| 宅男免费午夜| 国产高清videossex| 一级a爱片免费观看的视频| 三级毛片av免费| 999久久久精品免费观看国产| 91久久精品国产一区二区成人 | 一二三四在线观看免费中文在| 国内毛片毛片毛片毛片毛片| 少妇的逼水好多| 日韩欧美三级三区| 欧美大码av| 男女下面进入的视频免费午夜| 国产亚洲欧美98| 国内少妇人妻偷人精品xxx网站 | 国产三级在线视频| 日韩有码中文字幕| 1000部很黄的大片| 久久精品91无色码中文字幕| 婷婷精品国产亚洲av| 叶爱在线成人免费视频播放| 黑人操中国人逼视频| 99在线人妻在线中文字幕| 国产 一区 欧美 日韩| 久久亚洲精品不卡| 熟女少妇亚洲综合色aaa.| 久久人人精品亚洲av| 成人性生交大片免费视频hd| 亚洲一区二区三区色噜噜| 国产成人影院久久av| 色老头精品视频在线观看| 亚洲九九香蕉| 国产成+人综合+亚洲专区| 一二三四社区在线视频社区8| 最新在线观看一区二区三区| 看黄色毛片网站| 久久久久性生活片| 国产激情欧美一区二区| 国产高清激情床上av| 国产欧美日韩一区二区精品| 亚洲中文日韩欧美视频| www.自偷自拍.com| 19禁男女啪啪无遮挡网站| 日本撒尿小便嘘嘘汇集6| 免费一级毛片在线播放高清视频| 日韩人妻高清精品专区| 成在线人永久免费视频| 国产精品亚洲av一区麻豆| 欧美激情久久久久久爽电影| 精华霜和精华液先用哪个| 午夜免费激情av| 丁香六月欧美| 久久久久亚洲av毛片大全| 亚洲精品中文字幕一二三四区| 国产高清videossex| 老汉色∧v一级毛片| 欧美日韩乱码在线| 亚洲午夜理论影院| 国产精品女同一区二区软件 | 长腿黑丝高跟| 一级毛片女人18水好多| 国产黄a三级三级三级人| 淫妇啪啪啪对白视频| 男人舔奶头视频| 国产成人精品无人区| 国产激情偷乱视频一区二区| 国产一区二区在线av高清观看| 淫妇啪啪啪对白视频| 日韩欧美三级三区| 人妻丰满熟妇av一区二区三区| 久久中文看片网| 99久久成人亚洲精品观看| 国产伦精品一区二区三区视频9 | 一级毛片高清免费大全| 90打野战视频偷拍视频| 黑人操中国人逼视频| 亚洲狠狠婷婷综合久久图片| 亚洲专区字幕在线| 亚洲欧美日韩卡通动漫| 国产精品一区二区精品视频观看| 丝袜人妻中文字幕| 久久精品91蜜桃| 给我免费播放毛片高清在线观看| 日本黄大片高清| 真人一进一出gif抽搐免费| 国内精品美女久久久久久| 日本一二三区视频观看| 1000部很黄的大片| 午夜激情欧美在线| 无遮挡黄片免费观看| 9191精品国产免费久久| 九九热线精品视视频播放| 少妇的逼水好多| 淫秽高清视频在线观看| 中文字幕高清在线视频| 我要搜黄色片| www日本在线高清视频| 男女之事视频高清在线观看| 网址你懂的国产日韩在线| 一个人观看的视频www高清免费观看 | 国产不卡一卡二| 精品久久久久久成人av| 久久精品国产亚洲av香蕉五月| 亚洲第一电影网av| 99精品欧美一区二区三区四区| 国产在线精品亚洲第一网站| 国产三级黄色录像| 欧美日韩国产亚洲二区| 亚洲aⅴ乱码一区二区在线播放| 欧美日韩精品网址| 亚洲五月婷婷丁香| 午夜视频精品福利| 国产视频内射| 免费看a级黄色片| 亚洲人成伊人成综合网2020| 久久国产乱子伦精品免费另类| 久久久久国产一级毛片高清牌| 国产真人三级小视频在线观看| 亚洲av电影不卡..在线观看| 国产精品久久久久久精品电影| 国产aⅴ精品一区二区三区波| 亚洲精品色激情综合| 免费在线观看日本一区| 露出奶头的视频| 亚洲专区字幕在线| 天天躁日日操中文字幕| 亚洲成人免费电影在线观看| 99精品久久久久人妻精品| 欧美一区二区精品小视频在线| 非洲黑人性xxxx精品又粗又长| 久久人人精品亚洲av| 亚洲欧美精品综合久久99| 国产精品香港三级国产av潘金莲| 国产又黄又爽又无遮挡在线| 国产午夜精品久久久久久| 亚洲中文日韩欧美视频| 亚洲av第一区精品v没综合| 老司机在亚洲福利影院| 日本一二三区视频观看| 欧美+亚洲+日韩+国产| 两人在一起打扑克的视频| 精华霜和精华液先用哪个| 久久伊人香网站| 亚洲av免费在线观看| 免费观看精品视频网站| 国产精品自产拍在线观看55亚洲| 国产精品国产高清国产av| 欧美成人一区二区免费高清观看 | 夜夜夜夜夜久久久久| 精品久久久久久,| 亚洲av电影不卡..在线观看| 99热只有精品国产| 麻豆一二三区av精品| 欧美性猛交╳xxx乱大交人| 熟女人妻精品中文字幕| 91字幕亚洲| 色视频www国产| 午夜福利欧美成人| 天堂影院成人在线观看| 久久久久久国产a免费观看| 国产伦在线观看视频一区| 91久久精品国产一区二区成人 | 日韩欧美国产一区二区入口| 精品日产1卡2卡| 午夜免费观看网址| 免费一级毛片在线播放高清视频| 少妇裸体淫交视频免费看高清| 老司机在亚洲福利影院| tocl精华| 欧美日本亚洲视频在线播放| 色尼玛亚洲综合影院| 国产高清视频在线观看网站| 韩国av一区二区三区四区| cao死你这个sao货| 欧美日韩黄片免| 国内精品久久久久久久电影| 日韩 欧美 亚洲 中文字幕| www.熟女人妻精品国产| 床上黄色一级片| 成人永久免费在线观看视频| 成人无遮挡网站| 午夜视频精品福利| 亚洲第一欧美日韩一区二区三区| 国产私拍福利视频在线观看| 国产成人精品久久二区二区免费| 日日干狠狠操夜夜爽| 麻豆成人av在线观看| 精品熟女少妇八av免费久了| 精华霜和精华液先用哪个| 中文在线观看免费www的网站| 国产成人福利小说| 中文字幕av在线有码专区| 视频区欧美日本亚洲| 天堂av国产一区二区熟女人妻| 国产精品一区二区三区四区免费观看 | 97超级碰碰碰精品色视频在线观看| 日韩欧美国产一区二区入口| 亚洲人成电影免费在线| 欧美一区二区精品小视频在线| 国产欧美日韩一区二区精品| 国产亚洲精品av在线| 18禁观看日本| 丝袜人妻中文字幕| 亚洲中文字幕一区二区三区有码在线看 | 中出人妻视频一区二区| 国产又色又爽无遮挡免费看| 久久久国产成人精品二区| av天堂在线播放| 国内精品久久久久久久电影| 国产男靠女视频免费网站| 国产69精品久久久久777片 | 99riav亚洲国产免费| 亚洲精品一卡2卡三卡4卡5卡| 久久久久久久精品吃奶| 一个人看视频在线观看www免费 | 91av网站免费观看| 人人妻人人看人人澡| 亚洲中文av在线| 国产探花在线观看一区二区| 深夜精品福利| 久久热在线av| 香蕉av资源在线| 精品久久久久久久久久久久久| 国产高清激情床上av| 成人亚洲精品av一区二区| 18禁黄网站禁片免费观看直播| 日本黄色片子视频| 国产av在哪里看| 黄色片一级片一级黄色片| 男人舔女人的私密视频| 最近在线观看免费完整版| 999久久久精品免费观看国产| 国产精华一区二区三区| 999久久久精品免费观看国产| 久久九九热精品免费| 我的老师免费观看完整版| 青草久久国产| 亚洲欧美日韩无卡精品| 久久久久久久久免费视频了| 叶爱在线成人免费视频播放| 国产精品久久久久久久电影 | 国产成+人综合+亚洲专区| 国产成人影院久久av| 国产av在哪里看| 欧美大码av| 久久中文字幕人妻熟女| 亚洲熟女毛片儿| aaaaa片日本免费| 午夜免费观看网址| 香蕉久久夜色| 嫁个100分男人电影在线观看| 亚洲精品一区av在线观看| 国产aⅴ精品一区二区三区波| 国产精品精品国产色婷婷| 国产亚洲精品久久久久久毛片| 久久这里只有精品19| 色综合亚洲欧美另类图片| 亚洲午夜理论影院| 国产成人影院久久av| 国产三级在线视频| 亚洲av成人精品一区久久| 黄频高清免费视频| 黄色 视频免费看| 欧美日韩亚洲国产一区二区在线观看| 女生性感内裤真人,穿戴方法视频| 曰老女人黄片| 三级毛片av免费| 欧美午夜高清在线| 国产伦在线观看视频一区| 两性夫妻黄色片| 亚洲成人中文字幕在线播放| 怎么达到女性高潮| 亚洲黑人精品在线| 国产精品久久久久久亚洲av鲁大| 丰满人妻熟妇乱又伦精品不卡| 精品久久久久久久毛片微露脸| av福利片在线观看| 亚洲欧美激情综合另类| 在线免费观看的www视频| 欧美绝顶高潮抽搐喷水| 久久久色成人| 最新中文字幕久久久久 | 97超视频在线观看视频| 日韩国内少妇激情av| 亚洲美女视频黄频| 啦啦啦韩国在线观看视频| 99国产精品一区二区蜜桃av| 又粗又爽又猛毛片免费看| 日韩欧美国产在线观看| 日本免费a在线| 身体一侧抽搐| 精品99又大又爽又粗少妇毛片 | 亚洲国产精品合色在线| 国产成人aa在线观看| 国产欧美日韩一区二区三| 久久精品亚洲精品国产色婷小说| 最近在线观看免费完整版| 国产主播在线观看一区二区| 欧美日本视频| 亚洲色图 男人天堂 中文字幕| 午夜a级毛片| 中文字幕高清在线视频| 桃红色精品国产亚洲av| av国产免费在线观看| 看片在线看免费视频| 丝袜人妻中文字幕| 国产成人系列免费观看| 国产1区2区3区精品| 欧美xxxx黑人xx丫x性爽| 日韩欧美精品v在线| 国产精品电影一区二区三区| 午夜福利在线观看吧| 日韩欧美精品v在线| 国产伦精品一区二区三区四那| 国产成人精品无人区| 国产精品久久久久久亚洲av鲁大| 国产精品一区二区免费欧美| 欧美在线一区亚洲| 午夜免费成人在线视频| 国产精品一区二区免费欧美| 亚洲黑人精品在线| 女人高潮潮喷娇喘18禁视频| 免费观看精品视频网站| 一夜夜www| 一个人看的www免费观看视频| 99国产综合亚洲精品| 精品熟女少妇八av免费久了| 精品一区二区三区视频在线观看免费| 热99re8久久精品国产| 亚洲美女黄片视频|