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

    聚焦、檢測(cè)、跟蹤、再聚焦:基于經(jīng)典跟蹤前檢測(cè)框架的實(shí)時(shí)小目標(biāo)檢測(cè)技術(shù)

    2019-02-13 08:22:22劉本源宋志勇范紅旗
    航空兵器 2019年6期
    關(guān)鍵詞:弱小加速度距離

    劉本源,宋志勇,范紅旗*

    (1.空軍軍醫(yī)大學(xué) 軍事生物醫(yī)學(xué)工程學(xué)系,西安 710032;2.國(guó)防科技大學(xué) 電子科學(xué)學(xué)院,長(zhǎng)沙 410073)

    0 引 言

    隱身技術(shù)和隱身兵器的發(fā)展和技術(shù)迭代,對(duì)雷達(dá)系統(tǒng)的檢測(cè)和估計(jì)性能提出了挑戰(zhàn)[1-2]。隱身目標(biāo)多為弱小、非合作目標(biāo),典型的截獲過程需要將目標(biāo)在不同時(shí)刻、不同位置的能量累積起來。飛行器的運(yùn)動(dòng)屬于一種物理參數(shù)限制的馬爾可夫過程,因此可利用目標(biāo)的雷達(dá)回波在不同時(shí)刻間的相關(guān)性,通過長(zhǎng)時(shí)間積累使目標(biāo)能量聚焦,進(jìn)而逐步突破檢測(cè)門限生成目標(biāo)的檢測(cè)報(bào)告。

    目前,現(xiàn)有的小目標(biāo)檢測(cè)跟蹤架構(gòu)多為一種聚焦[3]、跟蹤的級(jí)聯(lián)模式。首先,通過聚焦算法累積目標(biāo)能量,隨后利用“恒虛警”(CFAR)檢測(cè)生成過門限的目標(biāo)報(bào)告,最后通過濾波確認(rèn)目標(biāo)航跡。此類傳統(tǒng)方法的核心是聚焦,聚焦方法的好壞直接決定了小目標(biāo)的檢測(cè)跟蹤性能。

    Xu等人[4]提出了Radon Fourier Transform(RFT)方法,其核心思路是

    Li等人[5]提出的Radon Lv變換算法,思路與RFT類似。此類方法需要對(duì)速度進(jìn)行網(wǎng)格搜索,計(jì)算量大,并且速度的離散化精度(步長(zhǎng))有限,影響信號(hào)聚焦效果。同時(shí),該方法未考慮二次項(xiàng)(加速度)對(duì)聚焦性能的影響。在此基礎(chǔ)上,Chen等人[6]提出了RFT結(jié)合數(shù)階FFT(FrFT)變換的RFRFT方法。當(dāng)加速度為0,RFRFT等同于RFT;當(dāng)不存在距離走動(dòng)(Range Cell Migration,RCM)時(shí),等同于利用FrFT估計(jì)加速度,但是RFRFT需要同時(shí)對(duì)速度v、加速度a進(jìn)行網(wǎng)格搜索,計(jì)算量大且實(shí)時(shí)性差。

    Keystone變換(Keystone Transform,KT)是一種無需網(wǎng)格搜索的速度補(bǔ)償方法,通過建立雷達(dá)回波模型,在變換域利用坐標(biāo)伸縮變換對(duì)回波模型中的一次項(xiàng)、二次項(xiàng)進(jìn)行解耦。經(jīng)典的KT可利用插值實(shí)現(xiàn)。Zhu等人[7]提出了一種基于Chirp Z Transform(CZT)無需插值的KT實(shí)現(xiàn)方法。Zheng等人[8]提出了針對(duì)一次、二次誤差的三維變換補(bǔ)償算法,可以看作一種通用的KT(General KT,GKT)算法。Tian等人[9]針對(duì)加速度補(bǔ)償需求,結(jié)合KT和RFT,提出了KT-RFT算法。Jin等人根據(jù)信號(hào)特征,提出了計(jì)算二階相位差[10]的信號(hào)積累算法。Zhang等人[11]提出利用二階KT(Second-order KT,SKT)對(duì)三次相位進(jìn)行補(bǔ)償,實(shí)現(xiàn)弱小機(jī)動(dòng)目標(biāo)的聚焦和估計(jì)。以上方法的主要思路是利用SKT對(duì)距離加速度解耦,隨后估計(jì)相位的chirp rate(CR)及其導(dǎo)數(shù),得到加速度和加加速度,最后進(jìn)行相位補(bǔ)償,提升運(yùn)動(dòng)目標(biāo)能量積累增益。

    值得注意的是,最近提出的檢測(cè)前跟蹤(Track Before Detect,TBD)思路,可以認(rèn)為是聚焦、檢測(cè)、跟蹤這三種模塊的有機(jī)組合。由于此類方法避免了CFAR檢測(cè)對(duì)目標(biāo)信噪比利用率的損失,具備優(yōu)于DBT的性能。但是TBD多采用粒子濾波[12]或隨機(jī)有限集濾波器實(shí)現(xiàn),計(jì)算量[12-13]及實(shí)現(xiàn)難度遠(yuǎn)大于傳統(tǒng)的多目標(biāo)多Kalman濾波器。

    為此,本文針對(duì)實(shí)時(shí)弱小目標(biāo)檢測(cè)和跟蹤的需求,提出一種聚焦、檢測(cè)、跟蹤、再聚焦(Focus、Detect、Track、Re-Focus,F(xiàn)DTR)弱小目標(biāo)檢測(cè)跟蹤方法。該方法是一種基于傳統(tǒng)檢測(cè)后跟蹤架構(gòu)(DBT)的工程實(shí)現(xiàn),所有模塊均采用基本的檢測(cè)跟蹤方法實(shí)現(xiàn),通過精心設(shè)計(jì)的聚焦和恒加加速度(Constant Jerk Model,CJM)模型,利用多目標(biāo)多卡爾曼濾波器(Multiple Targets Multiple Kalman Filter,MTMKF),實(shí)現(xiàn)弱小目標(biāo)的檢測(cè)、跟蹤、航跡維護(hù)等功能。

    本項(xiàng)目的源代碼全部公開,結(jié)果可復(fù)現(xiàn),便于讀者下載評(píng)估。源代碼下載地址為

    https://github.com/liubenyuan/fdtr

    1 數(shù)學(xué)模型

    1.1 變量及定義

    在距離多普勒RD平面內(nèi),快時(shí)間(fast time)維度為319,慢時(shí)間(slow time)以32個(gè)脈沖構(gòu)成一幀(frame)數(shù)據(jù)。雷達(dá)脈沖重復(fù)頻率為32 kHz,幀周期為1 ms。為了檢測(cè)弱小目標(biāo),需要多幀相參積累。本文中相參積累時(shí)間(CPI)為25幀。在一個(gè)CPI內(nèi),完成目標(biāo)的聚焦和檢測(cè);多個(gè)CPI之間,利用多目標(biāo)多卡爾曼濾波器進(jìn)行濾波估計(jì)。

    CPI內(nèi)的原始數(shù)據(jù)定義為X,有X∈C319×32N,其中N為相參積累幀數(shù)。定義Xi為X的第i幀數(shù)據(jù),有Xi∈C319×32。

    令R為CPI內(nèi)的RD圖像序列。RD圖像序列的獲得方法為,依次對(duì)Xi,i∈{1, …,N}的慢時(shí)間進(jìn)行FFT變換,得到第i幀RD圖像Ri∈C319×NFFT,其中NFFT為FFT點(diǎn)數(shù);而R∈C319×NFFT×N即為多幀RD圖像序列,或稱多觀測(cè)分量(Multiple Measurement Vector,MMV)RD模型。R為一個(gè)三維數(shù)據(jù)結(jié)構(gòu),通過對(duì)幀間信號(hào)聚焦,實(shí)現(xiàn)弱小目標(biāo)能量的積累。令F為一個(gè)CPI聚焦后的圖像,F(xiàn)取Focus意,有F∈R319×NFFT。

    1.2 雷達(dá)回波的數(shù)學(xué)模型

    雷達(dá)回波模型可寫為[7]

    (1)

    其中,τ是距離維時(shí)間(快時(shí)間,fast time),t是慢時(shí)間(slow time)。式(1)中第1項(xiàng)是距離回波的沖激響應(yīng),其傅里葉變換為

    (2)

    M(fτ)=FFT(u(τ))

    對(duì)τ進(jìn)行FFT,可得

    (3)

    將r(t)在t=0處展開,可得

    (4)

    (5)

    距離的非線性項(xiàng)(速度、加速度)影響相位和包絡(luò),降低長(zhǎng)相參周期下信號(hào)的積累效果?;夭ㄐ盘?hào)模型可寫作:

    s(fτ,t)=M(fτ)·

    (6)

    本文將在該信號(hào)模型下討論聚焦算法。

    2 FDTR小目標(biāo)檢測(cè)跟蹤方法

    2.1 總體架構(gòu)

    FDTR的總體架構(gòu)如圖1所示。

    針對(duì)距離多普勒(RD)平面內(nèi)的弱小目標(biāo)檢測(cè)問題,通過多幀信號(hào)的累積實(shí)現(xiàn)弱小目標(biāo)的聚焦和參數(shù)提取,利用卡爾曼濾波實(shí)現(xiàn)目標(biāo)狀態(tài)的精確估計(jì)。

    圖1 FDTR算法整體流程圖

    Fig.1 The flowchart of FDTR

    2.2 Focus:信號(hào)聚焦

    2.2.1 Keystone變換

    該模塊主要完成RD信號(hào)Ri的聚焦,利用一階Keystone變換消除距離速度耦合,提升信號(hào)積累的能量。令

    (7)

    則雷達(dá)回波的相位為

    (8)

    通過KT變換,消除了速度v和距離fτ之間的耦合。KT變換無需目標(biāo)速度先驗(yàn)信息或進(jìn)行速度搜索。KT變換可以采用插值法[14]或Chirp-Z Transform(CZT)實(shí)現(xiàn)[15]。CZT實(shí)現(xiàn)[7, 15]基于傅里葉變換,具備優(yōu)異的計(jì)算效率。

    2.2.2 RD平面距離差分目標(biāo)增強(qiáng)方法

    FDTR采用了一種在距離維計(jì)算差分的增強(qiáng)方法。借鑒圖像處理中“微弱線條增強(qiáng)(dim line enhancement)”的思路,對(duì)KT變換后的RD圖像序列R,在距離維上計(jì)算2階差分:

    Renhance=R(r)″

    (9)

    具體實(shí)現(xiàn)中可利用有限差分法逼近,即在距離維計(jì)算與核函數(shù)[1, -2, 1]T的卷積。該方法既可以利用能量變化的特征,還可以保留相位以便進(jìn)行二次補(bǔ)償和糾正;具備相參積累和非相參積累的優(yōu)勢(shì),計(jì)算效率高。

    2.3 Detect:檢測(cè)

    采用二維CFAR檢測(cè),隔離單元gap為1,參考單元ref為2,門限閾值比scale為3.03。

    對(duì)于弱小目標(biāo)的檢測(cè),級(jí)聯(lián)CFAR降低了信噪比利用率。本文采用傳統(tǒng)的級(jí)聯(lián)CFAR方案,目的是探究經(jīng)典方法在弱小目標(biāo)檢測(cè)和估計(jì)中的極限性能。后續(xù)研究可以通過交換模塊順序,或者修改門限閾值比,利用檢測(cè)前跟蹤的思路,提升信噪比的利用率。

    目標(biāo)在距離多普勒維度上具備一定的寬度,較強(qiáng)的目標(biāo)會(huì)跨越幾個(gè)單元。因此,在具體實(shí)現(xiàn)中,并聯(lián)兩個(gè)CFAR,分別針對(duì)差分前和差分后的數(shù)據(jù)進(jìn)行檢測(cè),隨后利用RD聚類方法(代碼:RDScan)生成過門限的目標(biāo)報(bào)告。

    2.4 Track:跟蹤

    采用恒加加速度(Jerk)模型,對(duì)速度和加速度的估計(jì)進(jìn)行平滑濾波xn+1=Fxn,其中:

    (10)

    (11)

    隨后,利用多目標(biāo)多卡爾曼(Multiple Targets Multiple Kalman Filter,MTMKF)濾波器,實(shí)現(xiàn)多目標(biāo)的跟蹤。在該濾波器中,實(shí)現(xiàn)了目標(biāo)生成、消亡、記憶、合并等操作規(guī)則。

    2.5 Re-Focus:再聚焦

    針對(duì)一個(gè)相參周期CPI內(nèi)的數(shù)據(jù)R,F(xiàn)DTR不但估計(jì)距離、速度,同時(shí)利用過門限目標(biāo)所在單元的序列像,估計(jì)加速度。其主要思路是二分法:

    r=R[rcell, vcell]∈CN

    r1,r2=split(r)

    v1,v2=velocity(r1,r2)

    a=2(v2-v1)/T

    (12)

    式中,r為目標(biāo)所在距離多普勒單元的多幀序列像(復(fù)數(shù)),split函數(shù)將r分為前后兩部分,隨后分別估計(jì)這兩部分的速度,利用差分計(jì)算加速度。

    3 實(shí)驗(yàn)結(jié)果及討論

    3.1 聚焦及檢測(cè)

    首先,通過對(duì)Xi的慢時(shí)間做FFT得到單幀RD圖像Ri,見圖2。圖中存在兩個(gè)目標(biāo),弱小目標(biāo)由于信噪比低,在單幀圖像中很難檢測(cè)。

    圖2 單幀RD圖像

    Fig.2 Single frame RD image

    一種有效的策略是利用多幀相參或非相參積累,將弱小目標(biāo)的能量累積起來,使其突破閾值。選定圖2中弱目標(biāo)所在的多普勒單元(Doppler Cell),繪制距離序列像(Range-Frame),見圖3。

    圖3 弱小目標(biāo)的距離序列像

    Fig.3 Range-frame profile of the dim target

    從圖3可見,待檢測(cè)弱目標(biāo)在序列像中是一條微弱的線(Dim Line)。借鑒圖像處理中的“微弱線條增強(qiáng)”方法,即通過高通濾波增強(qiáng)圖像的邊緣。FDTR采用了一種差分聚焦方法:對(duì)距離維計(jì)算2階差分。從處理后的圖像中可看出弱目標(biāo)的序列軌跡(位于距離序列像上半?yún)^(qū)),見圖3(c)。

    圖4為1個(gè)CPI相參積累的結(jié)果。采用Radon傅里葉變換(RFT)和Radon差分(RT-SOD)都可改善目標(biāo)的累積信噪比,但提升效果有限。而本文提出的聚焦方法積累效果較好。

    圖4 強(qiáng)目標(biāo)及弱小目標(biāo)聚焦后的RD圖像

    Fig.4 Focused RD image of the strong and the dim targets

    圖5為弱目標(biāo)速度所在單元的距離像(range profile)。本文提出的聚焦算法可以在1個(gè)CPI中使目標(biāo)能量高于噪聲基底約4.5 dB。

    利用CFAR檢測(cè)并提取FOCUS后圖像中的目標(biāo)。CFAR參數(shù)為:隔離單元為1,參考單元(R軸、D軸)分別為2,門限比3.03。圖6中,利用聚焦后的圖像可有效檢測(cè)強(qiáng)目標(biāo)和弱目標(biāo)。

    采用Keystone變換(KT)或Radon變換(RT)聚焦都可抑制距離走動(dòng)。在較短的CPI內(nèi),若目標(biāo)加速度不致于導(dǎo)致多普勒單元走動(dòng),則兩者積累效果類似,見圖7。

    圖5 弱目標(biāo)聚焦后的距離像

    圖6 聚焦后RD圖像的CFAR檢測(cè)結(jié)果

    圖7 KT及RT弱目標(biāo)積累性能對(duì)比

    Fig.7 Accumulation performance of the dim target using the KT and RT algorithm

    由于KT屬于基于變換的距離多普勒解耦,無需速度先驗(yàn)或者速度搜索,因此相比基于離散網(wǎng)格搜索的RT方法補(bǔ)償精度更高。相比RT算法,KT在目標(biāo)的聚焦效果上有2~3 dB左右的提升。

    3.2 跟蹤及再聚焦

    圖8 二分法估計(jì)加速度

    Fig.8 Estimation of the acceleration using the bisection method

    圖9 目標(biāo)序列像的再聚焦

    Fig.9 Re-focus of the range profile on the target

    圖10為在線公開數(shù)據(jù)集的速度及加速度估計(jì)結(jié)果,可見FDTR具有較高的估計(jì)精度。

    圖10 加速度及速度再聚焦的估計(jì)結(jié)果

    Fig.10 Estimation accuracy of the acceleration and the velocity after re-focus

    實(shí)現(xiàn)了最簡(jiǎn)單的多目標(biāo)多卡爾曼跟蹤器(MTMKF),其優(yōu)點(diǎn)是便于實(shí)現(xiàn)且計(jì)算效率高。MTMKF集成CV/CA/CJ模型,通過設(shè)計(jì)航跡關(guān)聯(lián)算法,并利用距離變化率(Range Rate)排除不符合1重速度模糊假設(shè)的航跡;進(jìn)一步,設(shè)計(jì)航跡生成(激活)、消亡(休眠)、確認(rèn)、關(guān)聯(lián)模塊,實(shí)現(xiàn)多目標(biāo)跟蹤與濾波。

    圖11 弱小目標(biāo)檢測(cè)與估計(jì)得分結(jié)果

    Fig.11 Scores of the detection and estimation of dim targets

    3.3 FDTR算法的特點(diǎn)

    (1)執(zhí)行效率高。完整單一數(shù)據(jù)集(時(shí)長(zhǎng)2 000 ms),運(yùn)行時(shí)間小于56 s。

    (2)序貫檢測(cè)。僅僅使用25 ms數(shù)據(jù)生成距離多普勒RD平面內(nèi)的檢測(cè)結(jié)果,可方便擴(kuò)展至50 ms、100 ms結(jié)構(gòu)。

    (3)集成加速度估計(jì)、濾波與補(bǔ)償。針對(duì)過門限的檢測(cè)數(shù)據(jù),采用二分法估計(jì)加速度;集成恒定加加速度Jerk運(yùn)動(dòng)模型,對(duì)當(dāng)前加速度的估計(jì)值進(jìn)行濾波平滑。

    (4)實(shí)現(xiàn)多目標(biāo)多卡爾曼濾波器。該濾波器集成目標(biāo)的關(guān)聯(lián)、生成、消亡、距離變化率檢測(cè)等功能,實(shí)時(shí)性強(qiáng)。

    (5)采用CZT算法實(shí)現(xiàn)Keystone變換。

    (6)純Python實(shí)現(xiàn),可移植性強(qiáng)。關(guān)鍵模塊可用C++進(jìn)一步優(yōu)化, 采用面向?qū)ο蟮哪K化設(shè)計(jì),可方便評(píng)估不同的聚焦、補(bǔ)償和跟蹤算法。

    3.4 下一步的改進(jìn)方向

    (1)速度估計(jì)方法。FDTR首先在RD平面實(shí)現(xiàn)弱小目標(biāo)的積累和檢測(cè),隨后提取目標(biāo)所在Doppler單元的距離像,進(jìn)行相位補(bǔ)償和速度的精確估計(jì)。但是,若存在Doppler單元走動(dòng)的問題,則會(huì)影響相位補(bǔ)償和速度估計(jì)的精度。后續(xù)可移除聚類模塊,降低檢測(cè)門限,采用檢測(cè)前跟蹤或檢測(cè)前聚焦等方法,準(zhǔn)確提取粗Doppler單元。

    (2)加速度及變加速度補(bǔ)償。若目標(biāo)存在機(jī)動(dòng),則需考慮加速度補(bǔ)償以提升速度估計(jì)精度。可采用二階KT變換方法(Second-order Keystone Transform,SKT),即

    回波相位為

    利用

    (1+fτ/fc)1/2≈1+fτ/(2fc)

    (1+fτ/fc)-1/2≈1-fτ/(2fc)

    可得回波相位為

    (3)級(jí)聯(lián)的檢測(cè)后跟蹤方法損失目標(biāo)信噪比,可考慮使用粒子濾波或PHD等檢測(cè)前跟蹤算法。

    (4)基于隨機(jī)有限集的濾波方法。本文采用最簡(jiǎn)單的多目標(biāo)多Kalman濾波器實(shí)現(xiàn)基于規(guī)則的多目標(biāo)的新生、消亡和數(shù)據(jù)關(guān)聯(lián),后續(xù)可考慮使用PHD或CPHD實(shí)現(xiàn)多目標(biāo)濾波。

    4 總 結(jié)

    本文采用經(jīng)典的DBT框架及多目標(biāo)多Kalman濾波器,通過設(shè)計(jì)聚焦和再聚焦(相位補(bǔ)償)方法,以完整有效的系統(tǒng)集成能力,實(shí)現(xiàn)實(shí)時(shí)小目標(biāo)的檢測(cè)與高精度估計(jì)。本項(xiàng)目屬于面向工程的代碼設(shè)計(jì),采用了數(shù)據(jù)流編程范式,接口規(guī)范、易于修改、擴(kuò)展性強(qiáng)。

    猜你喜歡
    弱小加速度距離
    ??!水手,上岸吧
    “鱉”不住了!從26元/斤飆至38元/斤,2022年甲魚能否再跑出“加速度”?
    天際加速度
    汽車觀察(2018年12期)2018-12-26 01:05:42
    算距離
    創(chuàng)新,動(dòng)能轉(zhuǎn)換的“加速度”
    金橋(2018年4期)2018-09-26 02:24:46
    死亡加速度
    每次失敗都會(huì)距離成功更近一步
    山東青年(2016年3期)2016-02-28 14:25:55
    我有特別的喝水技巧
    柴的嘲笑
    愛的距離
    母子健康(2015年1期)2015-02-28 11:21:33
    一二三四中文在线观看免费高清| 黄色视频在线播放观看不卡| 国内精品宾馆在线| 午夜福利影视在线免费观看| 久久午夜综合久久蜜桃| 亚洲欧洲日产国产| 99久久精品热视频| 一区二区三区精品91| 肉色欧美久久久久久久蜜桃| 亚洲国产日韩一区二区| 亚洲丝袜综合中文字幕| 自线自在国产av| 亚洲成人一二三区av| 国产精品嫩草影院av在线观看| 日本av手机在线免费观看| 欧美精品国产亚洲| 最新中文字幕久久久久| 国产黄色免费在线视频| 天天操日日干夜夜撸| 欧美三级亚洲精品| 亚洲在久久综合| 久久久精品免费免费高清| 国产高清三级在线| 欧美成人精品欧美一级黄| 国产视频内射| 下体分泌物呈黄色| 国产精品嫩草影院av在线观看| 丰满少妇做爰视频| 精品久久久久久电影网| 美女中出高潮动态图| 日日啪夜夜爽| av在线app专区| 亚洲美女搞黄在线观看| 色视频www国产| 少妇 在线观看| 人妻系列 视频| 久久久久精品性色| 日韩免费高清中文字幕av| 最后的刺客免费高清国语| 午夜免费鲁丝| 男女国产视频网站| 国产精品三级大全| 日韩中文字幕视频在线看片| 人人妻人人澡人人看| 亚洲欧洲国产日韩| 国产一区二区在线观看日韩| 午夜日本视频在线| 两个人的视频大全免费| 精品熟女少妇av免费看| 大香蕉久久网| 免费黄频网站在线观看国产| 少妇精品久久久久久久| 国内揄拍国产精品人妻在线| 久久久久久久大尺度免费视频| 亚洲国产精品一区二区三区在线| a级片在线免费高清观看视频| 亚洲人成网站在线观看播放| 熟女av电影| 欧美+日韩+精品| 亚洲性久久影院| 尾随美女入室| 看十八女毛片水多多多| 国产有黄有色有爽视频| 中文字幕人妻丝袜制服| 国产成人精品无人区| 交换朋友夫妻互换小说| 中文天堂在线官网| 蜜臀久久99精品久久宅男| 国产亚洲5aaaaa淫片| 国产精品福利在线免费观看| 国产黄频视频在线观看| 69精品国产乱码久久久| 亚洲欧美成人综合另类久久久| 成人国产麻豆网| 亚洲综合色惰| 色吧在线观看| 黑人高潮一二区| 亚洲内射少妇av| 狂野欧美激情性xxxx在线观看| 欧美日韩视频高清一区二区三区二| 久久综合国产亚洲精品| 人人澡人人妻人| 国产成人免费无遮挡视频| 我要看黄色一级片免费的| 亚洲人与动物交配视频| 天堂俺去俺来也www色官网| 极品教师在线视频| 中文在线观看免费www的网站| a级一级毛片免费在线观看| 午夜福利,免费看| 国产欧美亚洲国产| 成人亚洲欧美一区二区av| 日本-黄色视频高清免费观看| 亚洲一区二区三区欧美精品| 国产伦精品一区二区三区视频9| 亚洲av成人精品一二三区| 97超视频在线观看视频| 色哟哟·www| 噜噜噜噜噜久久久久久91| 91在线精品国自产拍蜜月| 黄色毛片三级朝国网站 | 一级黄片播放器| 成人二区视频| 亚洲人成网站在线播| 精品亚洲成a人片在线观看| 在线观看免费视频网站a站| 99久久精品一区二区三区| 亚洲欧美成人精品一区二区| 久久久久久久久久成人| www.av在线官网国产| 你懂的网址亚洲精品在线观看| 永久网站在线| 欧美精品一区二区免费开放| 国产黄频视频在线观看| 哪个播放器可以免费观看大片| 极品人妻少妇av视频| xxx大片免费视频| 国产高清国产精品国产三级| 男男h啪啪无遮挡| 国产精品久久久久成人av| 国产中年淑女户外野战色| 少妇熟女欧美另类| 欧美成人精品欧美一级黄| 亚洲精品乱码久久久v下载方式| 国产色爽女视频免费观看| 少妇丰满av| 国产精品无大码| 日韩精品有码人妻一区| 久久久久久久久大av| 久久久精品94久久精品| 亚洲国产精品一区三区| 免费观看性生交大片5| 午夜福利影视在线免费观看| 最近2019中文字幕mv第一页| 婷婷色av中文字幕| 日本91视频免费播放| 日韩电影二区| 伦精品一区二区三区| 五月开心婷婷网| 天堂中文最新版在线下载| 亚洲av中文av极速乱| 一区二区av电影网| 日本猛色少妇xxxxx猛交久久| 中国国产av一级| av卡一久久| av女优亚洲男人天堂| 韩国av在线不卡| 日本黄色片子视频| 免费黄网站久久成人精品| 成人美女网站在线观看视频| 亚洲国产精品999| 伦理电影大哥的女人| 三级国产精品欧美在线观看| 日韩一本色道免费dvd| 精华霜和精华液先用哪个| 亚洲精品乱久久久久久| 丰满迷人的少妇在线观看| 亚洲av成人精品一二三区| 亚洲av福利一区| 亚洲欧洲国产日韩| 中文精品一卡2卡3卡4更新| 成人影院久久| 一区二区三区乱码不卡18| 免费黄色在线免费观看| 国产乱人偷精品视频| 另类亚洲欧美激情| 九草在线视频观看| 欧美精品人与动牲交sv欧美| 69精品国产乱码久久久| 国产伦精品一区二区三区视频9| 国产高清国产精品国产三级| 一级毛片我不卡| 99热全是精品| 亚洲国产毛片av蜜桃av| 在线精品无人区一区二区三| 成人亚洲欧美一区二区av| 久久免费观看电影| 26uuu在线亚洲综合色| 午夜激情久久久久久久| 观看美女的网站| 一本一本综合久久| 国产精品蜜桃在线观看| 午夜影院在线不卡| 国产毛片在线视频| 日本免费在线观看一区| 国产精品福利在线免费观看| 亚洲国产精品专区欧美| 欧美少妇被猛烈插入视频| 久久国产亚洲av麻豆专区| 国产精品99久久99久久久不卡 | 波野结衣二区三区在线| 国产精品久久久久久久电影| 深夜a级毛片| 亚州av有码| 国产成人免费无遮挡视频| 国产美女午夜福利| 中文字幕人妻熟人妻熟丝袜美| 日日摸夜夜添夜夜添av毛片| 寂寞人妻少妇视频99o| 99久久中文字幕三级久久日本| 天堂8中文在线网| 亚洲精品国产av成人精品| 国产精品蜜桃在线观看| 亚洲丝袜综合中文字幕| 免费av不卡在线播放| 天天操日日干夜夜撸| 国产伦精品一区二区三区视频9| 久久久久人妻精品一区果冻| 大片免费播放器 马上看| 99热全是精品| 99热这里只有是精品50| 欧美人与善性xxx| 高清黄色对白视频在线免费看 | 国产成人aa在线观看| 欧美成人精品欧美一级黄| 99久久精品一区二区三区| av福利片在线观看| 欧美精品一区二区免费开放| 十分钟在线观看高清视频www | 99热全是精品| 黄色欧美视频在线观看| 亚洲中文av在线| av在线观看视频网站免费| 啦啦啦在线观看免费高清www| 热re99久久国产66热| 天堂俺去俺来也www色官网| 免费少妇av软件| 九草在线视频观看| 中文欧美无线码| 免费av中文字幕在线| 国产在视频线精品| 精品99又大又爽又粗少妇毛片| 成人午夜精彩视频在线观看| 99热这里只有是精品在线观看| 51国产日韩欧美| 少妇精品久久久久久久| 十八禁网站网址无遮挡 | 成人无遮挡网站| 自线自在国产av| 久久av网站| 久久免费观看电影| 亚洲精品色激情综合| 亚洲成人手机| 一级av片app| 久热这里只有精品99| 欧美激情极品国产一区二区三区 | 最近的中文字幕免费完整| 高清在线视频一区二区三区| 尾随美女入室| 日韩av不卡免费在线播放| 国产成人91sexporn| 久久97久久精品| 亚洲精品中文字幕在线视频 | 永久网站在线| 亚洲av.av天堂| av国产精品久久久久影院| 99热全是精品| 久久精品熟女亚洲av麻豆精品| .国产精品久久| 又大又黄又爽视频免费| 日本av手机在线免费观看| 国产精品一二三区在线看| 国产伦精品一区二区三区四那| 九九在线视频观看精品| 五月伊人婷婷丁香| 日本爱情动作片www.在线观看| 国产日韩欧美视频二区| 亚洲第一av免费看| 免费观看在线日韩| 亚洲欧洲国产日韩| 99九九线精品视频在线观看视频| 国产淫片久久久久久久久| 欧美一级a爱片免费观看看| 不卡视频在线观看欧美| av播播在线观看一区| 欧美少妇被猛烈插入视频| 成年人午夜在线观看视频| 国产色爽女视频免费观看| 看十八女毛片水多多多| 国产精品成人在线| 久久影院123| 国产精品99久久99久久久不卡 | 在线观看三级黄色| 亚洲精品国产av成人精品| 伊人久久国产一区二区| 高清午夜精品一区二区三区| 草草在线视频免费看| 久久6这里有精品| 亚洲精品国产av蜜桃| 国产成人精品福利久久| av福利片在线| 校园人妻丝袜中文字幕| 欧美日韩亚洲高清精品| 人人妻人人看人人澡| 日韩av不卡免费在线播放| 男女国产视频网站| 久久久久国产精品人妻一区二区| 国产精品国产av在线观看| 天天躁夜夜躁狠狠久久av| 一边亲一边摸免费视频| 国产片特级美女逼逼视频| 久久狼人影院| av网站免费在线观看视频| 欧美精品国产亚洲| 亚洲国产精品成人久久小说| 成人综合一区亚洲| 久久免费观看电影| 黑人猛操日本美女一级片| 王馨瑶露胸无遮挡在线观看| 一级爰片在线观看| 99视频精品全部免费 在线| 国产成人免费观看mmmm| 国语对白做爰xxxⅹ性视频网站| 久久久国产欧美日韩av| 黑丝袜美女国产一区| 免费久久久久久久精品成人欧美视频 | 中国三级夫妇交换| 99九九在线精品视频 | 国产欧美日韩综合在线一区二区 | 丝瓜视频免费看黄片| 欧美xxⅹ黑人| 色哟哟·www| 日韩 亚洲 欧美在线| 亚洲精品一二三| 2022亚洲国产成人精品| 少妇裸体淫交视频免费看高清| 欧美成人午夜免费资源| 黑丝袜美女国产一区| 在线观看三级黄色| 亚洲四区av| 免费少妇av软件| 这个男人来自地球电影免费观看 | 天堂中文最新版在线下载| 又粗又硬又长又爽又黄的视频| 麻豆成人午夜福利视频| 久久毛片免费看一区二区三区| 亚洲欧美精品专区久久| 一二三四中文在线观看免费高清| 自线自在国产av| 国产成人精品无人区| 黑人猛操日本美女一级片| 国产精品久久久久久久久免| 亚洲成色77777| 日韩av免费高清视频| 少妇熟女欧美另类| 精品人妻熟女毛片av久久网站| 视频区图区小说| 美女cb高潮喷水在线观看| 免费人成在线观看视频色| 午夜av观看不卡| 亚洲av免费高清在线观看| 午夜激情福利司机影院| 伦精品一区二区三区| 欧美bdsm另类| 国产精品熟女久久久久浪| 午夜福利,免费看| 极品人妻少妇av视频| 国产老妇伦熟女老妇高清| 麻豆成人av视频| 国产一区二区在线观看av| 亚洲国产精品成人久久小说| 国产精品国产av在线观看| 日本av免费视频播放| 亚洲成人手机| 国产成人精品一,二区| 91午夜精品亚洲一区二区三区| 少妇裸体淫交视频免费看高清| 晚上一个人看的免费电影| 夫妻性生交免费视频一级片| 国产男人的电影天堂91| 国产乱来视频区| av有码第一页| 国产精品国产三级国产av玫瑰| 国产精品久久久久久精品电影小说| 免费看不卡的av| 日韩中字成人| 精品久久久精品久久久| 久久这里有精品视频免费| 在线观看免费日韩欧美大片 | 精品少妇内射三级| 免费少妇av软件| 国产亚洲5aaaaa淫片| 大片免费播放器 马上看| 80岁老熟妇乱子伦牲交| 国产亚洲5aaaaa淫片| 午夜老司机福利剧场| 午夜视频国产福利| 看非洲黑人一级黄片| 一本一本综合久久| 人妻一区二区av| 日日撸夜夜添| 午夜激情福利司机影院| 五月开心婷婷网| 97在线视频观看| 晚上一个人看的免费电影| 国产精品一区二区性色av| 青青草视频在线视频观看| av线在线观看网站| 久久久亚洲精品成人影院| 亚洲,一卡二卡三卡| 日韩视频在线欧美| av播播在线观看一区| 亚洲欧美清纯卡通| 成人综合一区亚洲| 欧美日韩在线观看h| av国产久精品久网站免费入址| 秋霞伦理黄片| 看非洲黑人一级黄片| 国产精品久久久久久精品电影小说| 在线精品无人区一区二区三| 一本色道久久久久久精品综合| 国产精品一二三区在线看| 一级毛片黄色毛片免费观看视频| 久久综合国产亚洲精品| 亚洲国产毛片av蜜桃av| 国产精品无大码| 女性被躁到高潮视频| 伦理电影大哥的女人| √禁漫天堂资源中文www| 国产欧美日韩一区二区三区在线 | 国产免费一级a男人的天堂| 欧美日韩视频高清一区二区三区二| 国产高清三级在线| 69精品国产乱码久久久| 欧美日韩av久久| .国产精品久久| 国产综合精华液| 丁香六月天网| freevideosex欧美| 久久午夜福利片| 国产精品国产av在线观看| 亚洲国产精品专区欧美| 国产精品伦人一区二区| 国内少妇人妻偷人精品xxx网站| 一区二区三区乱码不卡18| 国产成人freesex在线| 一级a做视频免费观看| 99久久人妻综合| 欧美亚洲 丝袜 人妻 在线| 亚洲精品国产av成人精品| 欧美激情国产日韩精品一区| 有码 亚洲区| 亚洲图色成人| 亚洲精品亚洲一区二区| 最后的刺客免费高清国语| 国产毛片在线视频| 国产精品免费大片| 十八禁高潮呻吟视频 | 99热全是精品| 亚洲成人av在线免费| 成人影院久久| 久久国产精品男人的天堂亚洲 | 精品久久久噜噜| 亚洲怡红院男人天堂| 亚洲欧美一区二区三区国产| 亚洲第一区二区三区不卡| 亚洲精品久久午夜乱码| 国产欧美日韩精品一区二区| 99热网站在线观看| 男人狂女人下面高潮的视频| 日韩欧美 国产精品| 最近手机中文字幕大全| 麻豆成人av视频| 日韩制服骚丝袜av| 亚洲va在线va天堂va国产| 又黄又爽又刺激的免费视频.| 国产欧美亚洲国产| 日产精品乱码卡一卡2卡三| 一区二区av电影网| 日韩成人av中文字幕在线观看| 建设人人有责人人尽责人人享有的| 亚洲av成人精品一区久久| 22中文网久久字幕| 18禁在线无遮挡免费观看视频| 18+在线观看网站| 国产成人精品婷婷| 久久人人爽人人片av| 久久精品国产亚洲网站| 欧美精品一区二区大全| 午夜老司机福利剧场| 亚洲国产日韩一区二区| 多毛熟女@视频| 黄色欧美视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 午夜精品国产一区二区电影| tube8黄色片| 日本wwww免费看| 亚洲av综合色区一区| 这个男人来自地球电影免费观看 | 2022亚洲国产成人精品| 简卡轻食公司| 国产精品一区二区性色av| 一区二区三区乱码不卡18| 看十八女毛片水多多多| 久久人人爽av亚洲精品天堂| 久久这里有精品视频免费| 18禁裸乳无遮挡动漫免费视频| 狠狠精品人妻久久久久久综合| 国产极品天堂在线| 日韩 亚洲 欧美在线| 成人毛片a级毛片在线播放| 一本大道久久a久久精品| 免费av不卡在线播放| 91在线精品国自产拍蜜月| 日日摸夜夜添夜夜添av毛片| 亚洲欧洲日产国产| 国产日韩欧美在线精品| 亚洲av.av天堂| 伊人久久国产一区二区| 国产极品粉嫩免费观看在线 | 日韩欧美精品免费久久| 久久鲁丝午夜福利片| 一级毛片久久久久久久久女| 黄色欧美视频在线观看| 亚洲一级一片aⅴ在线观看| 久久精品国产鲁丝片午夜精品| 看非洲黑人一级黄片| 免费看日本二区| 国产一区二区三区av在线| 99热这里只有是精品在线观看| 男女边摸边吃奶| 五月天丁香电影| 男女免费视频国产| 2021少妇久久久久久久久久久| 久久女婷五月综合色啪小说| 91久久精品国产一区二区成人| 女人久久www免费人成看片| 啦啦啦中文免费视频观看日本| 亚洲在久久综合| 久久精品久久精品一区二区三区| 91精品一卡2卡3卡4卡| 免费av中文字幕在线| 夜夜爽夜夜爽视频| 各种免费的搞黄视频| 久久久a久久爽久久v久久| 亚洲av男天堂| 久久精品国产鲁丝片午夜精品| 亚洲欧美成人综合另类久久久| 久久久久久久大尺度免费视频| 欧美亚洲 丝袜 人妻 在线| 日韩一区二区三区影片| 国产永久视频网站| 亚洲四区av| 女性生殖器流出的白浆| 免费久久久久久久精品成人欧美视频 | 免费大片18禁| 欧美老熟妇乱子伦牲交| 在线观看国产h片| 国产伦精品一区二区三区四那| 亚洲av电影在线观看一区二区三区| 国产在线男女| a级毛片免费高清观看在线播放| 亚洲人成网站在线观看播放| 国产男人的电影天堂91| 成人美女网站在线观看视频| 在线观看三级黄色| 国产免费又黄又爽又色| 黑人高潮一二区| 国产亚洲av片在线观看秒播厂| 日韩制服骚丝袜av| 久久精品国产鲁丝片午夜精品| 久久精品久久久久久噜噜老黄| 精品99又大又爽又粗少妇毛片| 精品久久久久久久久亚洲| 欧美日韩亚洲高清精品| 99热国产这里只有精品6| 国产91av在线免费观看| 亚洲欧美成人综合另类久久久| av线在线观看网站| www.色视频.com| 观看免费一级毛片| 黑人高潮一二区| 美女中出高潮动态图| 国产精品免费大片| av视频免费观看在线观看| 国产高清国产精品国产三级| 日日摸夜夜添夜夜添av毛片| 成人18禁高潮啪啪吃奶动态图 | 啦啦啦啦在线视频资源| 成人黄色视频免费在线看| 久久久久久久国产电影| 成人国产麻豆网| 十八禁高潮呻吟视频 | 观看免费一级毛片| 久久人人爽人人爽人人片va| 在线观看免费视频网站a站| 国产亚洲午夜精品一区二区久久| av专区在线播放| 久久久精品94久久精品| 91精品一卡2卡3卡4卡| 国产黄片视频在线免费观看| 亚洲中文av在线| 高清不卡的av网站| 欧美区成人在线视频| 国产精品嫩草影院av在线观看| 我的老师免费观看完整版| 下体分泌物呈黄色| 汤姆久久久久久久影院中文字幕| 国产视频内射| 一区在线观看完整版| 老熟女久久久| 亚洲精品色激情综合| 国产在线男女| 偷拍熟女少妇极品色| 不卡视频在线观看欧美| 好男人视频免费观看在线| 亚洲精品aⅴ在线观看| 国产精品99久久99久久久不卡 | 精品少妇黑人巨大在线播放| 亚洲,一卡二卡三卡| 黑丝袜美女国产一区| 亚洲四区av| 91午夜精品亚洲一区二区三区| 色视频www国产| av播播在线观看一区|