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

    一種針對臨近空間滑躍目標(biāo)的軌跡跟蹤方法

    2020-10-17 02:12:08秦武韜陸小科邢曉勇
    中國慣性技術(shù)學(xué)報 2020年3期
    關(guān)鍵詞:臨空滑翔航向

    秦武韜,陸小科,邢曉勇

    (南京電子技術(shù)研究所,南京 210039)

    隨著科學(xué)技術(shù)的飛速發(fā)展,各主要國家的反導(dǎo)系統(tǒng)日臻完善[1,2],傳統(tǒng)彈道導(dǎo)彈的突防能力持續(xù)下降,在此背景下,以臨近空間內(nèi)高速滑翔為標(biāo)志的臨近空間飛行器逐漸成為遠(yuǎn)程戰(zhàn)略打擊武器研發(fā)的重要方向[3,4]。美國于1987年率先開始了高超聲速滑翔式飛行器(Hypersonic Glide Vehicle,HGV)項目的研究,近十年更是多次進(jìn)行了HTV-2、X-51A 等臨近空間飛行器的飛行試驗,俄羅斯則提出了“冷”和“鷹”的研究計劃,德國、日本、印度也開展了相關(guān)研究[5]。因此,研究針對臨空目標(biāo)的軌跡跟蹤方法具有極大的理論和實際意義。與常規(guī)的軌道目標(biāo)不同,臨近空間飛行器不僅擁有10 馬赫以上的飛行速度,還可在大氣層內(nèi)利用氣動力進(jìn)行大范圍的機動飛行,其典型的跳躍式軌跡更是給雷達(dá)的高精度穩(wěn)定跟蹤帶來了嚴(yán)峻挑戰(zhàn)[6]。

    在機動目標(biāo)跟蹤中,準(zhǔn)確的目標(biāo)運動模型是精確穩(wěn)定跟蹤的基礎(chǔ)[7,8]。例如針對勻速運動的常速度模型(Constant Velocity,CV)和針對勻加速度運動的常加速度模型(Constant Acceleration,CA);對于轉(zhuǎn)彎目標(biāo)的跟蹤問題,文獻(xiàn)[9]則給出了協(xié)同轉(zhuǎn)彎模型(Constant Turn,CT)。但CV、CA 和CT 模型均屬于理想化的建模,為了應(yīng)對跟蹤過程中目標(biāo)隨時可能出現(xiàn)的轉(zhuǎn)彎、加減速等突發(fā)機動情況,Singer 把機動控制項作為相關(guān)噪聲進(jìn)行建模,提出了Singer 模型,該模型在處理機動目標(biāo)時具有良好的跟蹤效果。文獻(xiàn)[10]則專門針對臨近空間飛行器周期性跳躍的特點提出了正弦波模型(Sine Wave,SW),較好地描述了臨空目標(biāo)跳躍運動的周期特性,其結(jié)果表明該方法具有比Singer 模型、Jerk模型更高的跟蹤精度,文獻(xiàn)[11]則對SW 模型做了進(jìn)一步的分析和研究,并基于X-51A 軌跡驗證了SW 模型的跟蹤精度。

    然而,以HTV-2 為代表的無動力臨近空間飛行器在進(jìn)行跳躍滑翔時,其能量不斷減少,跳躍幅度逐漸降低,各周期的相關(guān)性持續(xù)減弱,又呈現(xiàn)出衰減性的特點。SW 模型在進(jìn)行建模時僅考慮到了臨近空間飛行器跳躍滑翔的周期性,而忽視了其衰減性,導(dǎo)致其在跟蹤無動力臨空飛行器時跟蹤精度較低。因此,本文基于一種二階時間自相關(guān)的零均值隨機過程,借鑒“當(dāng)前”統(tǒng)計模型的思想,提出了一種自適應(yīng)衰減震蕩模型(Adaptive Damped Oscillation,ADO)。該模型通過衰減參數(shù)和震蕩參數(shù)的合理調(diào)配,實現(xiàn)了周期性和衰減性的有機統(tǒng)一,較好地描述了臨空滑躍目標(biāo)的運動特點。考慮到臨空目標(biāo)機動范圍大的特點,論文在地理坐標(biāo)系下推導(dǎo)建立了目標(biāo)運動模型,利用無跡卡爾曼濾波算法設(shè)計了跟蹤濾波器,并通過數(shù)學(xué)仿真對所提算法進(jìn)行了驗證。

    1 臨近空間跳躍滑翔飛行器運動描述及特性分析

    1.1 目標(biāo)運動描述

    由于臨近空間飛行器具有機動范圍廣的顯著特點,為直觀地對其運動進(jìn)行描述,本文擬在地理坐標(biāo)系下對其進(jìn)行描述,即以經(jīng)度、緯度和高度表示目標(biāo)位置,以速度大小、飛行路徑角和航向角表示目標(biāo)速度。其中,飛行路徑角γ表示速度矢量與當(dāng)?shù)厮矫娴膴A角,航向角ψ表示速度矢量在當(dāng)?shù)厮矫嫔系耐队芭c正北方向的夾角。

    由飛行路徑角及航向角的定義,有

    式中,V為合速度大小,vE、vN和vU分別為東向、北向和天向速度。

    圖1 地理系下位置與速度關(guān)系Fig.1 The relation of position and velocity in GEO

    圖1給出了經(jīng)度、緯度、高度變化與東、北、天速度的關(guān)系,有:

    式中,λ、L和h分別表示經(jīng)度、緯度和高度,Re表示地球半徑。

    將式(1)帶入式(2),可得

    1.2 運動特性分析

    與傳統(tǒng)彈道導(dǎo)彈不同,臨近空間高速滑翔飛行器可以依靠氣動力在大氣層內(nèi)進(jìn)行高速滑翔,根據(jù)滑翔時高度方向的運動特點大致可以劃分為平衡滑翔和跳躍滑翔兩種模式,其中平衡滑翔模式下目標(biāo)高度方向機動較小,跟蹤難度較低,而跳躍滑翔模式的跟蹤難度更大,因此,本文將主要對這一類型的臨空目標(biāo)跟蹤問題進(jìn)行研究。

    圖2 跳躍式滑翔目標(biāo)運動軌跡Fig.2 The trajectory of NSHV

    呈“周期性衰減”式的跳躍滑翔是臨近空間飛行器所特有的飛行樣式,彈道導(dǎo)彈、飛機等傳統(tǒng)飛行器均無法做出類似機動。在這種機動模式下,臨近空間飛行器不斷上下跳躍,呈現(xiàn)出顯著的周期性(如圖2所示),此時,目標(biāo)的加速度和飛行路徑角也呈現(xiàn)出同樣的周期性,各個周期間的目標(biāo)運動狀態(tài)強相關(guān)。但與此同時,目標(biāo)又在氣動阻力的作用下不斷消耗能量,其速度不斷下降,跳躍幅度逐漸降低,各個周期間的相關(guān)性持續(xù)減弱,表現(xiàn)出衰減特性。因此,臨近空間跳躍滑翔目標(biāo)的運動呈現(xiàn)出周期性與衰減性相統(tǒng)一的特點。

    2 基于自適應(yīng)衰減震蕩模型的縱向機動建模

    2.1 零均值衰減震蕩模型

    由于臨空滑躍目標(biāo)飛行路徑角和加速度的變化規(guī)律相似,因此本文以加速度為例推導(dǎo)描述其衰減震蕩過程的運動方程。根據(jù)臨空滑躍目標(biāo)周期性運動的特點,其t時刻的加速度a(t)與一個周期ΔTc后的加速度a(t+ΔTc)呈強相關(guān)。同時,隨著間隔周期數(shù)的增多,目標(biāo)加速度逐漸降低,該相關(guān)性逐漸減弱,表現(xiàn)出了相關(guān)衰減性的特點。為描述這一運動特征,可利用如下二階時間自相關(guān)隨機過程進(jìn)行加速度建模,即

    當(dāng)α和β取極限值時,有:

    可以看出,在α逼近0 時,式(4)給出的衰減震蕩模型不再具備衰減性的特點,而僅表現(xiàn)出周期性,其退化為SW 模型[10];當(dāng)β逼近0 時,衰減震蕩模型則喪失了周期性的特點,而僅表現(xiàn)出衰減性,其退化為Singer 模型。因此,通過參數(shù)α及β的調(diào)配,衰減震蕩模型可以實現(xiàn)周期性和衰減性的統(tǒng)一。

    由式(4)所給出的相關(guān)函數(shù)Ra(τ),可得其譜密度為:

    其中,ω表示角速率。

    其中,W(ωj)表示白噪聲的傅里葉變換。

    將G(ω)代入H(ωj)中,可得

    令s=ωj,式(8)可化簡為

    由式(9)可得衰減震蕩模型下的加速度和加加速度狀態(tài)方程如下:

    其中,wa是均值為零,方差為的高斯白噪聲,d(t)為加加速度。

    2.2 自適應(yīng)衰減震蕩模型

    對于式(10)所建立的加速度變化模型,可看作SW模型和Singer 模型的組合,其中SW 模型認(rèn)為目標(biāo)加速度以零為均值進(jìn)行周期性的跳動,Singer 模型同樣認(rèn)為目標(biāo)機動加速度的均值為零。但事實上,臨空滑躍目標(biāo)在飛行過程中受氣動力的作用不斷減速,其加速度均值并不為零,而是隨其氣動外形、飛行速度和跳躍階段的變化而變化。為更準(zhǔn)確地對目標(biāo)的加速度變化規(guī)律進(jìn)行描述,本節(jié)將對式(10)所示模型做進(jìn)一步改進(jìn),推導(dǎo)出自適應(yīng)均值衰減震蕩模型。

    定義變量b(t)如下:

    則有

    由式(10)和(12)可得

    其中,

    將式(13)展開,有

    同理,可得飛行路徑角的變化模型如下:

    其中,ωγ(t)為飛行路徑角速度,(t)和(t)為上一周期中濾波器對γ及ωγ估計值的平均,wγ是均值為零、方差為的高斯白噪聲。

    3 側(cè)向運動建模及狀態(tài)方程的構(gòu)建

    3.1 側(cè)向運動模型

    當(dāng)目標(biāo)進(jìn)行側(cè)向機動時,其航向角速度會隨之發(fā)生變化,但由于目標(biāo)飛行速度較快,受結(jié)構(gòu)及氣動性能的限制,下一時刻的航向角速度的變化往往有限,僅能在當(dāng)前角速度的領(lǐng)域內(nèi),這與“當(dāng)前”統(tǒng)計模型中關(guān)于機動加速度的假設(shè)相一致。因此,本文借鑒“當(dāng)前”統(tǒng)計模型的思想,建立起關(guān)于航向角速度的自適應(yīng)均值一階馬爾科夫模型如下:

    其中,(t)服從零均值一階馬爾科夫過程,(t)為航向角速度的“當(dāng)前”均值,αψ為機動時間常數(shù),wωψ為零均值高斯噪聲。

    根據(jù)式(16),有

    結(jié)合式(16),(17)和(18),可得航向角及航向角速度的狀態(tài)方程如下:

    其中,航向角速度的方差為

    其中,為ωψ的估計值,ωψmax為最大正航向角速度,ωψ-max為最大負(fù)航向角速度。由式(20)可以看出,方差的大小與ωψmax和ωψ-max密切相關(guān),而又直接影響到過程噪聲方差陣Q 的大小。因此,為了獲得更好的跟蹤精度,在應(yīng)用過程中需要根據(jù)目標(biāo)側(cè)向機動能力的大致情況選取合適的ωψmax和ωψ-max。

    3.2 狀態(tài)方程

    綜合上述分析,本文選取目標(biāo)狀態(tài)量x如下:

    根據(jù)各狀態(tài)量的定義,結(jié)合式(3)(14)(15)和(19),可建立狀態(tài)方程如下:

    式中,0m×n表示m行n列的零矩陣。

    4 跟蹤濾波器設(shè)計

    現(xiàn)階段,可在大范圍內(nèi)對非合作目標(biāo)進(jìn)行連續(xù)探測的設(shè)備主要包括天基紅外預(yù)警衛(wèi)星和雷達(dá),但對于臨空滑躍目標(biāo),其飛行過程中的紅外特征不明顯,無法利用紅外衛(wèi)星對其進(jìn)行穩(wěn)定跟蹤。而雷達(dá)可以通過主動發(fā)射電磁波進(jìn)行探測,且具有作用距離遠(yuǎn)的顯著優(yōu)勢,可以滿足對臨空滑躍目標(biāo)的跟蹤要求。因此,本文以雷達(dá)為探測傳感器建立測量模型。圖3給出了雷達(dá)對目標(biāo)進(jìn)行探測的示意圖。

    圖3 雷達(dá)探測示意圖Fig.3 The schematic diagram of the detection of radar

    雷達(dá)通過獲得目標(biāo)的斜距和回波到達(dá)角對目標(biāo)進(jìn)行定位,因此測量量可定義為:

    式中,γ為高低角,η為方位角,ρ為斜距。

    則測量方程為

    式中,vγ,vη和vρ表示雷達(dá)的測量誤差;xe,ye和ze表示目標(biāo)在雷達(dá)直角坐標(biāo)系中的位置分量,有

    其中,地心系至雷達(dá)直角坐標(biāo)系的轉(zhuǎn)換矩陣,有

    由目標(biāo)在地理系下的位置信息,可得其在地心系下的位置[x,y,z]T為

    同理,雷達(dá)在地心系下的位置[xr,yr,zr]T為

    在式(30)中,[λr,Lr,hr]為雷達(dá)的經(jīng)緯高位置信息,可通過標(biāo)定獲得。

    圖4 無跡卡爾曼濾波算法流程圖Fig.4 The flow chart of UKF

    考慮到該系統(tǒng)的狀態(tài)方程和測量方程均具有較強的非線性,且雅克比矩陣求取困難,因此論文選用無跡卡爾曼濾波算法進(jìn)行濾波器的設(shè)計,圖4給出了無跡卡爾曼濾波的流程圖,其詳細(xì)內(nèi)容可參考文獻(xiàn)[12]。

    5 數(shù)學(xué)仿真分析

    本文以美國X-51A 及HTV-2 的試驗軌跡為參考,設(shè)置跳躍滑翔跟蹤場景如圖5所示,其中,目標(biāo)的初始位置為[100°E,30°N,48000 m],初始速度為4800 m/s,最大飛行高度為51108 m,最小飛行高度為46449 m,跳躍周期約為150 s,飛行時間為800 s,圖6給出了目標(biāo)的高度變化曲線,圖7給出了目標(biāo)的速度和飛行路徑角變化曲線。

    圖5 跳躍滑翔目標(biāo)運動軌跡Fig.5 The trajectory of the NSHV

    圖6 高度變化曲線Fig.6 Altitude-time curve

    圖7 速度和飛行路徑角變化曲線Fig.7 The velocity-time and flight path angle-time curves

    在仿真中,雷達(dá)位于[133 ° E,30 ° N,200 m],其距離測量誤差σr=50m,高低角和方位角測量誤差為σγ=ση=0.1°,測量頻率為10 Hz。

    各模型中的主要參數(shù)如下表所示:

    表1 各模型中的主要參數(shù)Tab.1 Main parameters of three models

    基于上述仿真場景,獲得單次仿真結(jié)果如下:

    圖8 東北天位置估計誤差Fig.8 Position estimated error of east-north-up direction

    在上述仿真結(jié)果中,圖8給出了目標(biāo)在東北天方向的位置估計誤差,圖9則給出了速度、航向角和飛行路徑角估計誤差,可以看出目標(biāo)的位置跟蹤誤差在1000 m 以內(nèi),速度跟蹤誤差在10 m/s 以內(nèi),具有估計精度高,收斂速度快的特點。

    圖9 速度、航向角和飛行路徑角估計誤差Fig.9 The estimation error of velocity,course angle and flight path angle

    為進(jìn)一步對比分析該算法的精度,論文分別利用常速度模型、Singer 模型、SW 模型和本文提出的自適應(yīng)衰減震蕩模型對臨空目標(biāo)進(jìn)行跟蹤,基于100 次獨立的蒙特卡洛試驗,獲得仿真結(jié)果如圖10所示:

    圖10 位置跟蹤誤差對比曲線Fig.10 The comparison of position estimation error between four algorithms

    圖11 速度跟蹤誤差對比曲線Fig.11 The comparison of velocity estimation error between four algorithms

    表2 各算法跟蹤誤差對比Tab.2 Comparison of tracking error of four algorithms

    圖10和圖11分別給出了四種算法的位置和速度跟蹤誤差對比曲線,表2對跟蹤誤差進(jìn)行了匯總??梢钥闯霰疚乃岢龅腁DO 算法具有最高的跟蹤精度,SW 模型次之,而CA 模型的跟蹤精度最差。這是因為CA 模型簡單地認(rèn)為目標(biāo)進(jìn)行勻加速直線運動,這與目標(biāo)的實際運動規(guī)律并不符合,導(dǎo)致其跟蹤精度較低。同理,Singer 模型也難以對跳躍滑翔目標(biāo)的運動進(jìn)行描述。SW 模型認(rèn)為目標(biāo)進(jìn)行周期性的跳躍運動,這一點較為準(zhǔn)確地描述了跳躍滑翔目標(biāo)的運動特性,但其在推導(dǎo)過程中假設(shè)各個周期目標(biāo)的跳躍高度相等,這與臨空目標(biāo)的實際飛行過程不符,使得其跟蹤精度略為下降。而本文所提出的自適應(yīng)衰減震蕩模型,不僅對臨空滑翔目標(biāo)的周期運動特性進(jìn)行了準(zhǔn)確描述,同時也考慮到了目標(biāo)飛行過程中的衰減特性,其對目標(biāo)運動的描述最為準(zhǔn)確,因而獲得了最高的軌跡跟蹤精度。ADO 算法其位置跟蹤精度和速度跟蹤精度較SW 模型分別提高了17.91%和26.33%。

    6 結(jié) 論

    針對臨近空間跳躍滑翔飛行器運動特性復(fù)雜,現(xiàn)有目標(biāo)運動模型描述失準(zhǔn),跟蹤精度較低的問題,本文基于二階時間自相關(guān)隨機過程提出了一種自適應(yīng)衰減震蕩模型,建立了地理系下的目標(biāo)運動模型,準(zhǔn)確描述了臨近空間滑翔跳躍目標(biāo)運動的周期性和衰減性。在此基礎(chǔ)上,論文給出了雷達(dá)的測量模型并利用無跡卡爾曼濾波算法完成了跟蹤濾波器的設(shè)計,數(shù)學(xué)仿真證明基于ADO 模型設(shè)計的跟蹤濾波器具有比CA模型、Singer 模型和SW 模型更高的跟蹤精度。

    猜你喜歡
    臨空滑翔航向
    蘭花螳螂會滑翔
    北京新機場臨空經(jīng)濟區(qū)生態(tài)空間識別與生態(tài)安全格局構(gòu)建
    攻天掠地的先鋒武器——滑翔導(dǎo)彈
    知坐標(biāo),明航向
    袁奇峰 空鐵上的廣州都市圈,花都臨空經(jīng)濟紅利空間巨大
    考慮幾何限制的航向道模式設(shè)計
    臨空高速目標(biāo)多傳感器自主協(xié)同資源調(diào)度算法
    探析臨空經(jīng)濟對區(qū)域物流競爭力的提升
    讓你心跳加速的“滑翔飛板”
    基于干擾觀測器的船舶系統(tǒng)航向Backstepping 控制
    電子制作(2017年24期)2017-02-02 07:14:16
    黄网站色视频无遮挡免费观看| 女人爽到高潮嗷嗷叫在线视频| 高清在线国产一区| 99久久综合免费| 在线十欧美十亚洲十日本专区| 一本大道久久a久久精品| 国产成人免费无遮挡视频| 精品国产乱子伦一区二区三区 | 午夜免费观看性视频| 国产一卡二卡三卡精品| 国产免费av片在线观看野外av| 热99久久久久精品小说推荐| 国产免费现黄频在线看| 母亲3免费完整高清在线观看| 亚洲精品久久久久久婷婷小说| 永久免费av网站大全| 大型av网站在线播放| 国产av精品麻豆| 国产亚洲精品久久久久5区| 国内毛片毛片毛片毛片毛片| 亚洲国产精品一区二区三区在线| 91精品三级在线观看| av网站在线播放免费| 一区二区av电影网| 亚洲激情五月婷婷啪啪| 三级毛片av免费| 两人在一起打扑克的视频| 18禁裸乳无遮挡动漫免费视频| 两人在一起打扑克的视频| 91麻豆av在线| 熟女少妇亚洲综合色aaa.| 日本猛色少妇xxxxx猛交久久| 大片电影免费在线观看免费| 国产亚洲午夜精品一区二区久久| 丰满人妻熟妇乱又伦精品不卡| 久久热在线av| 丰满迷人的少妇在线观看| 丁香六月欧美| netflix在线观看网站| 欧美性长视频在线观看| 9热在线视频观看99| 国产精品免费视频内射| 少妇精品久久久久久久| 久久亚洲国产成人精品v| 亚洲一码二码三码区别大吗| 99久久国产精品久久久| 久久久久视频综合| 国产99久久九九免费精品| 亚洲精品第二区| 少妇精品久久久久久久| 欧美日本中文国产一区发布| 久久99一区二区三区| 热99re8久久精品国产| 欧美另类一区| 又黄又粗又硬又大视频| 亚洲精品国产av蜜桃| 国产91精品成人一区二区三区 | 亚洲国产日韩一区二区| 国产亚洲精品久久久久5区| 成年人黄色毛片网站| 精品一区二区三区四区五区乱码| 大码成人一级视频| 国产一区二区 视频在线| 热99国产精品久久久久久7| 纯流量卡能插随身wifi吗| 欧美黄色片欧美黄色片| 人成视频在线观看免费观看| 91国产中文字幕| 午夜激情久久久久久久| 中文字幕最新亚洲高清| 正在播放国产对白刺激| 欧美 亚洲 国产 日韩一| 狠狠精品人妻久久久久久综合| 成人国产一区最新在线观看| 搡老岳熟女国产| 1024香蕉在线观看| 久久久久久久国产电影| 日韩中文字幕欧美一区二区| 啦啦啦视频在线资源免费观看| 天堂8中文在线网| 最黄视频免费看| 桃红色精品国产亚洲av| 亚洲七黄色美女视频| 少妇被粗大的猛进出69影院| 婷婷成人精品国产| 在线十欧美十亚洲十日本专区| 99九九在线精品视频| 国产男女内射视频| av又黄又爽大尺度在线免费看| 免费看十八禁软件| 肉色欧美久久久久久久蜜桃| 欧美黄色片欧美黄色片| avwww免费| 久久精品熟女亚洲av麻豆精品| 午夜福利一区二区在线看| 美女高潮喷水抽搐中文字幕| 欧美日本中文国产一区发布| 日韩,欧美,国产一区二区三区| 国产精品自产拍在线观看55亚洲 | 热99国产精品久久久久久7| 久久久精品区二区三区| 亚洲少妇的诱惑av| 国产91精品成人一区二区三区 | 99香蕉大伊视频| 日本av免费视频播放| 欧美日韩亚洲高清精品| 青青草视频在线视频观看| 午夜免费观看性视频| 久久久国产成人免费| 美女主播在线视频| 久久中文看片网| 亚洲欧美日韩高清在线视频 | 中文欧美无线码| 亚洲中文av在线| 男女午夜视频在线观看| 青青草视频在线视频观看| 777米奇影视久久| 中文字幕人妻丝袜制服| 精品亚洲成国产av| 久久久久久免费高清国产稀缺| 菩萨蛮人人尽说江南好唐韦庄| 精品乱码久久久久久99久播| 另类精品久久| 丝袜脚勾引网站| 国产成人av激情在线播放| 操美女的视频在线观看| 99香蕉大伊视频| 精品国产乱码久久久久久男人| av在线app专区| 99九九在线精品视频| 99久久精品国产亚洲精品| 国产av精品麻豆| 国产在线观看jvid| 亚洲欧美一区二区三区久久| bbb黄色大片| 久久国产精品大桥未久av| 国产精品麻豆人妻色哟哟久久| 女人爽到高潮嗷嗷叫在线视频| 极品人妻少妇av视频| 少妇裸体淫交视频免费看高清 | 青青草视频在线视频观看| 中文字幕人妻丝袜一区二区| 久久国产精品人妻蜜桃| 久久精品久久久久久噜噜老黄| 免费在线观看完整版高清| 五月开心婷婷网| 国产国语露脸激情在线看| 啦啦啦在线免费观看视频4| 亚洲专区中文字幕在线| 中文字幕另类日韩欧美亚洲嫩草| 日韩一卡2卡3卡4卡2021年| 成年美女黄网站色视频大全免费| 无限看片的www在线观看| a级毛片黄视频| 啦啦啦视频在线资源免费观看| 啪啪无遮挡十八禁网站| 三上悠亚av全集在线观看| 中国美女看黄片| 少妇精品久久久久久久| 窝窝影院91人妻| 午夜久久久在线观看| 精品久久蜜臀av无| 亚洲一卡2卡3卡4卡5卡精品中文| 精品少妇内射三级| 国产精品香港三级国产av潘金莲| 如日韩欧美国产精品一区二区三区| 午夜免费鲁丝| 精品亚洲成a人片在线观看| 免费高清在线观看日韩| 亚洲国产中文字幕在线视频| 蜜桃国产av成人99| 狠狠婷婷综合久久久久久88av| 国产男人的电影天堂91| 妹子高潮喷水视频| 丝袜美腿诱惑在线| 国产人伦9x9x在线观看| 国产深夜福利视频在线观看| 纵有疾风起免费观看全集完整版| 国产成人av激情在线播放| 久久人人97超碰香蕉20202| 国产成人精品在线电影| 亚洲欧洲日产国产| 国产深夜福利视频在线观看| 欧美精品一区二区大全| 久久青草综合色| 亚洲中文日韩欧美视频| 伊人亚洲综合成人网| 久久久久精品国产欧美久久久 | 欧美xxⅹ黑人| 国内毛片毛片毛片毛片毛片| 亚洲一卡2卡3卡4卡5卡精品中文| 国产高清videossex| 亚洲av电影在线观看一区二区三区| 精品熟女少妇八av免费久了| 十分钟在线观看高清视频www| 亚洲av男天堂| 久久ye,这里只有精品| 亚洲欧美日韩另类电影网站| 精品久久久久久久毛片微露脸 | 狠狠精品人妻久久久久久综合| 欧美日韩亚洲综合一区二区三区_| 国产区一区二久久| 亚洲成人免费电影在线观看| 天天添夜夜摸| 亚洲欧美日韩高清在线视频 | 肉色欧美久久久久久久蜜桃| 人人妻人人添人人爽欧美一区卜| 亚洲国产毛片av蜜桃av| 国产一区二区三区av在线| 午夜日韩欧美国产| 伊人久久大香线蕉亚洲五| 亚洲自偷自拍图片 自拍| 亚洲欧美日韩另类电影网站| 日本黄色日本黄色录像| 午夜精品国产一区二区电影| av欧美777| 欧美日本中文国产一区发布| 久久久久精品人妻al黑| 日韩视频一区二区在线观看| 久久人妻熟女aⅴ| 69精品国产乱码久久久| 欧美性长视频在线观看| 国产三级黄色录像| 精品欧美一区二区三区在线| 日韩欧美国产一区二区入口| 免费高清在线观看视频在线观看| 亚洲自偷自拍图片 自拍| 人妻一区二区av| 中文欧美无线码| 老司机深夜福利视频在线观看 | 国产亚洲欧美精品永久| 国产精品一区二区精品视频观看| 精品视频人人做人人爽| 国产成人精品无人区| 成年女人毛片免费观看观看9 | 建设人人有责人人尽责人人享有的| 少妇粗大呻吟视频| 婷婷成人精品国产| 男人舔女人的私密视频| 汤姆久久久久久久影院中文字幕| 不卡av一区二区三区| 一本一本久久a久久精品综合妖精| 日本撒尿小便嘘嘘汇集6| 狂野欧美激情性bbbbbb| 亚洲九九香蕉| 丁香六月欧美| 亚洲成人国产一区在线观看| 人妻久久中文字幕网| 悠悠久久av| 亚洲精品中文字幕一二三四区 | 中文字幕制服av| 天天操日日干夜夜撸| 久久国产亚洲av麻豆专区| 午夜免费观看性视频| 国产一区二区在线观看av| 久久国产精品影院| 亚洲精品中文字幕一二三四区 | 成人国产一区最新在线观看| 桃红色精品国产亚洲av| 美女国产高潮福利片在线看| 99精品欧美一区二区三区四区| 欧美精品一区二区免费开放| 婷婷成人精品国产| 日日夜夜操网爽| 欧美人与性动交α欧美软件| 黄色 视频免费看| 午夜福利,免费看| 1024视频免费在线观看| 久久久久久亚洲精品国产蜜桃av| 99久久99久久久精品蜜桃| 又黄又粗又硬又大视频| 中国国产av一级| 搡老乐熟女国产| 夜夜骑夜夜射夜夜干| 亚洲国产精品成人久久小说| 啦啦啦视频在线资源免费观看| 午夜精品国产一区二区电影| 免费女性裸体啪啪无遮挡网站| 热99久久久久精品小说推荐| 老司机影院毛片| 两性夫妻黄色片| 精品久久久久久电影网| 久久 成人 亚洲| 国产精品99久久99久久久不卡| 啦啦啦啦在线视频资源| 亚洲中文字幕日韩| 日韩欧美国产一区二区入口| 黄色视频不卡| 日本av免费视频播放| 国产精品久久久久久精品电影小说| 老司机影院毛片| 免费在线观看日本一区| 宅男免费午夜| 亚洲欧洲日产国产| 一区福利在线观看| 大片电影免费在线观看免费| 波多野结衣av一区二区av| 中文字幕人妻熟女乱码| 亚洲熟女毛片儿| 亚洲av日韩在线播放| 99国产精品99久久久久| 国产一区二区激情短视频 | 国产又爽黄色视频| 各种免费的搞黄视频| 老熟妇仑乱视频hdxx| 久久精品成人免费网站| 老司机在亚洲福利影院| 国产免费福利视频在线观看| 在线观看舔阴道视频| 日本精品一区二区三区蜜桃| 亚洲情色 制服丝袜| 一区福利在线观看| 国产激情久久老熟女| 老熟女久久久| 99久久国产精品久久久| 亚洲国产成人一精品久久久| 亚洲熟女精品中文字幕| 亚洲国产欧美日韩在线播放| 亚洲美女黄色视频免费看| 欧美日韩av久久| 久久人妻福利社区极品人妻图片| 99九九在线精品视频| av有码第一页| 精品少妇内射三级| 国产三级黄色录像| 国产成人a∨麻豆精品| 亚洲国产日韩一区二区| 亚洲精品av麻豆狂野| 蜜桃在线观看..| svipshipincom国产片| 搡老岳熟女国产| 亚洲国产av影院在线观看| 在线观看www视频免费| 视频在线观看一区二区三区| 亚洲精品一卡2卡三卡4卡5卡 | 女人爽到高潮嗷嗷叫在线视频| 性高湖久久久久久久久免费观看| 日日爽夜夜爽网站| 性高湖久久久久久久久免费观看| 99久久综合免费| 国产在线视频一区二区| 啦啦啦在线免费观看视频4| 99re6热这里在线精品视频| 男女国产视频网站| 另类亚洲欧美激情| 伦理电影免费视频| 女人高潮潮喷娇喘18禁视频| 伊人亚洲综合成人网| 亚洲国产欧美网| 成年女人毛片免费观看观看9 | 精品少妇一区二区三区视频日本电影| 欧美日韩亚洲高清精品| 狠狠婷婷综合久久久久久88av| 91大片在线观看| 日韩电影二区| 久久久国产一区二区| 欧美av亚洲av综合av国产av| 自拍欧美九色日韩亚洲蝌蚪91| 精品人妻一区二区三区麻豆| 亚洲精品乱久久久久久| 99国产极品粉嫩在线观看| 一区二区av电影网| 久久精品国产a三级三级三级| 欧美日韩福利视频一区二区| 狠狠精品人妻久久久久久综合| 精品熟女少妇八av免费久了| 日韩欧美国产一区二区入口| 十分钟在线观看高清视频www| 91成人精品电影| 国产一区二区在线观看av| 黑人巨大精品欧美一区二区蜜桃| 午夜老司机福利片| 91成人精品电影| 一本色道久久久久久精品综合| 成人国产一区最新在线观看| 国产成人免费观看mmmm| 精品亚洲成国产av| 亚洲综合色网址| 搡老岳熟女国产| 亚洲综合色网址| 日韩电影二区| 搡老乐熟女国产| 国产伦理片在线播放av一区| 久久99热这里只频精品6学生| 欧美黄色片欧美黄色片| 亚洲 欧美一区二区三区| www.av在线官网国产| 啦啦啦视频在线资源免费观看| av有码第一页| av电影中文网址| 黄频高清免费视频| 女性被躁到高潮视频| 亚洲av日韩在线播放| 亚洲午夜精品一区,二区,三区| 欧美精品一区二区大全| 国产视频一区二区在线看| 欧美日韩成人在线一区二区| 女人被躁到高潮嗷嗷叫费观| 老司机午夜十八禁免费视频| 91成人精品电影| 中文字幕色久视频| 青春草视频在线免费观看| 夫妻午夜视频| 亚洲精品一二三| 淫妇啪啪啪对白视频 | 大片电影免费在线观看免费| 国产精品国产av在线观看| 丰满人妻熟妇乱又伦精品不卡| 考比视频在线观看| 久久久久精品国产欧美久久久 | 欧美中文综合在线视频| 国产日韩欧美在线精品| www.精华液| 18禁国产床啪视频网站| 日日夜夜操网爽| 人人妻人人澡人人爽人人夜夜| 99久久综合免费| 丝袜美腿诱惑在线| 欧美日韩精品网址| a在线观看视频网站| 午夜两性在线视频| 99热国产这里只有精品6| 飞空精品影院首页| 亚洲国产av影院在线观看| 在线观看人妻少妇| 国产亚洲欧美精品永久| 久久中文看片网| 人人澡人人妻人| 亚洲精品国产精品久久久不卡| 国产成人免费无遮挡视频| 黑丝袜美女国产一区| 亚洲va日本ⅴa欧美va伊人久久 | 欧美日韩国产mv在线观看视频| 狠狠婷婷综合久久久久久88av| 亚洲午夜精品一区,二区,三区| 嫁个100分男人电影在线观看| 一级片免费观看大全| 91精品国产国语对白视频| 国产av一区二区精品久久| 国产欧美日韩综合在线一区二区| 午夜激情av网站| 国产男女超爽视频在线观看| 成人18禁高潮啪啪吃奶动态图| 国产亚洲精品第一综合不卡| 十八禁人妻一区二区| 亚洲专区字幕在线| 99精品欧美一区二区三区四区| 日本猛色少妇xxxxx猛交久久| 老熟女久久久| 国产精品成人在线| 国产精品一区二区精品视频观看| 国产精品av久久久久免费| 淫妇啪啪啪对白视频 | videosex国产| 国产精品欧美亚洲77777| 亚洲第一av免费看| 亚洲一卡2卡3卡4卡5卡精品中文| 动漫黄色视频在线观看| 久久久久久久久久久久大奶| 无限看片的www在线观看| 美女中出高潮动态图| 少妇人妻久久综合中文| 女人爽到高潮嗷嗷叫在线视频| 黑人巨大精品欧美一区二区蜜桃| 久久国产精品男人的天堂亚洲| 一区二区av电影网| 亚洲国产欧美日韩在线播放| 欧美亚洲 丝袜 人妻 在线| 777米奇影视久久| 国产在线一区二区三区精| 国产精品 国内视频| av在线播放精品| 久久亚洲国产成人精品v| 亚洲综合色网址| 狠狠精品人妻久久久久久综合| 汤姆久久久久久久影院中文字幕| 免费日韩欧美在线观看| 亚洲精品粉嫩美女一区| 国产精品一区二区在线不卡| 国产精品香港三级国产av潘金莲| 亚洲免费av在线视频| 飞空精品影院首页| 国产色视频综合| 99热全是精品| 久久免费观看电影| 欧美性长视频在线观看| 国产精品自产拍在线观看55亚洲 | 成年人免费黄色播放视频| av网站在线播放免费| 国产主播在线观看一区二区| 各种免费的搞黄视频| 日本91视频免费播放| 亚洲一区二区三区欧美精品| 欧美精品一区二区免费开放| 欧美日韩av久久| 国产亚洲精品久久久久5区| 国精品久久久久久国模美| 人妻一区二区av| 狠狠精品人妻久久久久久综合| 中国美女看黄片| 丝袜脚勾引网站| 亚洲国产精品一区二区三区在线| 搡老熟女国产l中国老女人| 久久女婷五月综合色啪小说| 日日摸夜夜添夜夜添小说| 亚洲男人天堂网一区| 高清在线国产一区| 精品人妻在线不人妻| 2018国产大陆天天弄谢| 97在线人人人人妻| 欧美日韩一级在线毛片| 久久久久国产精品人妻一区二区| 国产精品 国内视频| av免费在线观看网站| 麻豆av在线久日| 午夜两性在线视频| 18在线观看网站| 亚洲中文字幕日韩| 国产精品偷伦视频观看了| 亚洲精品国产区一区二| 免费黄频网站在线观看国产| 女人精品久久久久毛片| 男女无遮挡免费网站观看| 精品一区在线观看国产| 高清视频免费观看一区二区| 精品一区在线观看国产| 男女床上黄色一级片免费看| 日日夜夜操网爽| 精品久久久久久久毛片微露脸 | 美女视频免费永久观看网站| 国产精品偷伦视频观看了| 国产欧美日韩一区二区三区在线| 99热网站在线观看| 精品人妻一区二区三区麻豆| 亚洲第一欧美日韩一区二区三区 | 老熟妇乱子伦视频在线观看 | 人人澡人人妻人| 国产在线视频一区二区| 99九九在线精品视频| 国产成人精品在线电影| 老司机深夜福利视频在线观看 | 亚洲成人免费电影在线观看| 欧美人与性动交α欧美精品济南到| 午夜精品国产一区二区电影| 国产成人免费观看mmmm| 日本欧美视频一区| 爱豆传媒免费全集在线观看| 国产亚洲av片在线观看秒播厂| 大片电影免费在线观看免费| 久久亚洲国产成人精品v| xxxhd国产人妻xxx| 91成年电影在线观看| 一本大道久久a久久精品| 国产免费福利视频在线观看| 九色亚洲精品在线播放| 国产精品成人在线| 999久久久国产精品视频| 国产成人免费观看mmmm| 男女无遮挡免费网站观看| 久久午夜综合久久蜜桃| 久久久水蜜桃国产精品网| 青青草视频在线视频观看| 国产精品.久久久| 无限看片的www在线观看| 亚洲精华国产精华精| 国产无遮挡羞羞视频在线观看| 亚洲熟女精品中文字幕| av线在线观看网站| 成年女人毛片免费观看观看9 | 国产男女内射视频| 精品一品国产午夜福利视频| 正在播放国产对白刺激| 国产在线视频一区二区| 日韩有码中文字幕| 久久精品国产综合久久久| 日韩欧美免费精品| 国产老妇伦熟女老妇高清| 国产成人精品久久二区二区91| 又紧又爽又黄一区二区| 久久九九热精品免费| 另类精品久久| 免费在线观看黄色视频的| 人人妻,人人澡人人爽秒播| 香蕉国产在线看| 国产欧美日韩一区二区三区在线| e午夜精品久久久久久久| 一区福利在线观看| 成人18禁高潮啪啪吃奶动态图| 日韩熟女老妇一区二区性免费视频| 无限看片的www在线观看| 精品福利永久在线观看| 国产精品一二三区在线看| 最新的欧美精品一区二区| 国产一级毛片在线| 性少妇av在线| 欧美激情久久久久久爽电影 | 中文字幕人妻丝袜制服| 99久久综合免费| 欧美xxⅹ黑人| 午夜激情av网站| 两性夫妻黄色片| 亚洲伊人色综图| 日本猛色少妇xxxxx猛交久久| 久久av网站| 啦啦啦视频在线资源免费观看| 日韩欧美免费精品| 秋霞在线观看毛片| 操出白浆在线播放| 午夜福利视频在线观看免费| 日本猛色少妇xxxxx猛交久久| www.999成人在线观看| 亚洲中文日韩欧美视频| 男女下面插进去视频免费观看| 久久 成人 亚洲| 国产真人三级小视频在线观看|