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

    潛航器高度和載體姿態(tài)測(cè)量誤差對(duì)地形匹配定位的影響研究

    2022-08-08 12:39:22宋解元范昊鵬周義博
    測(cè)繪工程 2022年4期
    關(guān)鍵詞:底圖定位精度姿態(tài)

    宋解元,范昊鵬,周義博,賈 真

    (1.61365部隊(duì),天津 300143;2.信息工程大學(xué),鄭州 450001;3.智慧中原地理信息技術(shù)河南省協(xié)同創(chuàng)新中心,鄭州 450001;4.時(shí)空感知與智能處理自然資源部重點(diǎn)實(shí)驗(yàn)室,鄭州 450001;5.海南大學(xué),???570228)

    水下載體定位導(dǎo)航一直是海洋環(huán)境探測(cè)中的一項(xiàng)技術(shù)難題。由于電磁波在水下傳輸距離十分受限,因而目前水下航行器定位過程多采用慣性主動(dòng)導(dǎo)航技術(shù)或聲學(xué)導(dǎo)航技術(shù)。其中,慣性導(dǎo)航技術(shù)因誤差積累問題,長航時(shí)的水下航行器一般采用“慣導(dǎo)+”的方式定位,而慣導(dǎo)+地形匹配因匹配背景是地形數(shù)據(jù),分辨率高,且可采用面匹配策略,故匹配精度更為理想。挪威國防研究組織成功將自主水下潛航器(Autonomous Underwater Vehicle,AUV)同時(shí)應(yīng)用于商業(yè)和軍事領(lǐng)域,其研發(fā)的地形匹配定位系統(tǒng)在實(shí)驗(yàn)中獲得了較高的精度[1]:英國南安普頓大學(xué),以條帶狀數(shù)據(jù)作為海底地形的底圖數(shù)據(jù),通過水下地形匹配實(shí)驗(yàn),研究傳感器系統(tǒng)偏差對(duì)匹配定位的影響;日本東京大學(xué)及日本海洋工程研究所等單位聯(lián)合研制“TUNA-SAND”自主水下潛航器,利用多波束聲吶對(duì)海底地形進(jìn)行實(shí)時(shí)測(cè)量,進(jìn)而匹配定位,取得了較好的效果。

    哈爾濱工程大學(xué)對(duì)經(jīng)典的地形輪廓匹配算法 (TERCOM)進(jìn)行研究,并對(duì)其導(dǎo)航定位的可靠性進(jìn)行驗(yàn)證,提出一種基于Hu 矩的二維地形匹配算法,取得較好的抗噪聲效果[2];同時(shí),相關(guān)研究人員對(duì)常用的ICCP(Iterative Closest Contour Point)匹配算法進(jìn)行改進(jìn),減小迭代次數(shù)和最近點(diǎn)計(jì)算量,提高計(jì)算效率[3]。東南大學(xué)針對(duì)低成本低分辨率底圖情況,提出并行ICCP地形匹配導(dǎo)航方法[4-6]。另外,大連艦艇學(xué)院開展的基于三角格網(wǎng)底圖的地形匹配算法[7],南京工程學(xué)院提出的ICCP-KF融合地形匹配算法[8],武漢大學(xué)測(cè)繪學(xué)院的水下地形抗差匹配算法[9],都為水下定位導(dǎo)航技術(shù)的發(fā)展做出貢獻(xiàn)。

    綜合國內(nèi)外研究進(jìn)展,目前水下地形匹配導(dǎo)航的研究主要集中在地形匹配技術(shù)在不同載體上的應(yīng)用、地形匹配導(dǎo)航系統(tǒng)的濾波器設(shè)計(jì)及匹配定位適配性指標(biāo)分析等;在水下載體的姿態(tài)變化、探頭安裝誤差、深度測(cè)量誤差、聲速誤差、潮位及底圖的時(shí)變性對(duì)定位精度的影響方面尚無系統(tǒng)性的研究成果。文中從載體運(yùn)行狀態(tài)入手,研究水下平臺(tái)距離海底地形高度、姿態(tài)角測(cè)量誤差對(duì)地形匹配定位的影響程度,從而得出各測(cè)量元件為達(dá)到總體精度而需要實(shí)現(xiàn)的精度要求,可為將來實(shí)現(xiàn)更高精度測(cè)量、制定更科學(xué)的技術(shù)指標(biāo)提供參考。

    1 研究方法

    1.1 簡化測(cè)量采樣模型

    水下地形匹配定位是利用潛航器搭載的多波束測(cè)深系統(tǒng)實(shí)時(shí)測(cè)量水下地形數(shù)據(jù),然后與現(xiàn)有底圖進(jìn)行匹配來確定自身位置的技術(shù)。在實(shí)際的多波束測(cè)深系統(tǒng)中,測(cè)量瞬間的結(jié)果為垂直航跡方向的一維線數(shù)據(jù)??紤]到多波束系統(tǒng)采樣頻率較大,即短時(shí)間內(nèi)可發(fā)收多組波束;為簡化模型,近似認(rèn)為多波束獲取的為該時(shí)段平均時(shí)刻采集的一份二維地形數(shù)據(jù)?;驹韴D如圖1所示。

    圖1 潛航器地形匹配原理

    1.2聲線投射點(diǎn)迭代算法

    在與底圖數(shù)據(jù)匹配之前,需要獲取聲線與水下地形的交點(diǎn)位置。由于海底地形曲面具有較強(qiáng)的不規(guī)則性,無法用簡單函數(shù)描述該曲面。參考文獻(xiàn)[10]采用迭代的方式求解投射點(diǎn)坐標(biāo)。同時(shí)顧及計(jì)算效率和準(zhǔn)確率,提出一種改進(jìn)的步長迭代方法,即使用動(dòng)態(tài)步長代替等步長。首先利用大步長迅速找到位于地形面兩側(cè)的兩個(gè)迭代點(diǎn);然后將搜索區(qū)域縮小,以實(shí)現(xiàn)搜索加速。流程如圖2所示。

    圖2 改進(jìn)的步長法流程

    1.3 內(nèi)插算法

    真實(shí)的海底地形是連續(xù)的曲面,而數(shù)字高程地形圖則是離散化的點(diǎn)云數(shù)據(jù)。為得到格網(wǎng)點(diǎn)間的未知點(diǎn)位置,可利用數(shù)字高程中已知格網(wǎng)點(diǎn)的數(shù)據(jù)內(nèi)插獲取。此處,選用內(nèi)插效果較為平滑的雙三次內(nèi)插法:選取與點(diǎn)最近的16個(gè)格網(wǎng)點(diǎn)作為計(jì)算參考值。假設(shè)最近的16個(gè)點(diǎn)的高程值分別為h1,h2,…,h16,則雙三次內(nèi)插法的算式為:

    (1)

    其中,權(quán)系數(shù)Wi通過BiCubic基函數(shù)得到,hi為第i個(gè)格網(wǎng)點(diǎn)的高度;具體基函數(shù)形式如下:

    W(x)=

    (2)

    式中:x為插值點(diǎn)與附近節(jié)點(diǎn)之間格網(wǎng)數(shù)。

    1.4 地形匹配算法

    匹配算法算式為:

    (3)

    式中:σ為當(dāng)前位置地形觀測(cè)高程數(shù)據(jù)的標(biāo)準(zhǔn)差;zi為潛航器聲線到第i個(gè)地形點(diǎn)在縱軸方向的高度,即第i個(gè)觀測(cè)數(shù)據(jù)值;N為聲線數(shù)量;L為該二維地形測(cè)量面與底圖的匹配度,L越大則測(cè)量面與底圖吻合度越大。

    1.5 數(shù)據(jù)仿真與定位精度評(píng)定方法

    由于影響因素眾多,難以直接得到定位誤差與潛航器指標(biāo)的關(guān)系。因此,根據(jù)蒙特·卡羅方法隨機(jī)模擬不同量級(jí)的載體高度、姿態(tài)角等指標(biāo);隨后利用迭代-內(nèi)插的方法,解算出聲線在水下的投射位置;最后再利用匹配算法計(jì)算潛航器的理論位置,并將其與設(shè)定位置相比較,得出定位誤差與各指標(biāo)的關(guān)系。當(dāng)模擬的數(shù)據(jù)量足夠多時(shí),該定位結(jié)果的統(tǒng)計(jì)值便可近似逼近實(shí)際情況。

    在眾多影響水下地形匹配的因素中,文中主要圍繞載體高度、姿態(tài)角,分析其定位效果,并建立數(shù)學(xué)模型。設(shè)定位誤差與各因素的函數(shù)關(guān)系式為:

    ΔP=f(a,d).

    (4)

    式中:d為載體距離海底的高度;a為載體姿態(tài)角;ΔP為定位誤差。

    需要說明的是,載體姿態(tài)角包括俯仰角、橫滾角和艏搖角。假設(shè)姿態(tài)角3個(gè)分量誤差的隨機(jī)分布相同,均符合均值為0、標(biāo)準(zhǔn)差為σ的正態(tài)分布。所以姿態(tài)角轉(zhuǎn)變?yōu)殛P(guān)于單一變量σ的函數(shù),即只利用σ來控制姿態(tài)角的誤差等級(jí)。

    另外,定位誤差ΔP在水平方向分為X和Y兩個(gè)分量。統(tǒng)計(jì)潛航器定位結(jié)果與理論位置,利用貝塞爾公式[11]計(jì)算兩個(gè)分量定位誤差的中誤差,以其平方和作為定位誤差指標(biāo)。

    Δxi=x0-xi,Δyi=y0-yi.

    (5)

    (6)

    (7)

    式中:xi,yi分別是實(shí)驗(yàn)中定位潛航器在(x0,y0)位置的橫縱坐標(biāo)值;Sx,Sy為利用貝塞爾公式計(jì)算得到的X,Y方向的定位誤差中誤差。

    2 數(shù)據(jù)準(zhǔn)備

    2.1 數(shù)字高程底圖

    在我國周邊某海域選取一塊海底地形數(shù)據(jù),利用等比放縮的方式得到600 m×600 m范圍、分辨率0.2 m×0.2 m的實(shí)驗(yàn)數(shù)據(jù)(見圖3)。

    圖3 實(shí)驗(yàn)底圖

    2.2 幾種坐標(biāo)系及其轉(zhuǎn)換

    在匹配仿真的過程中,由于潛航器測(cè)量的坐標(biāo)數(shù)據(jù)與底圖高程數(shù)據(jù)的坐標(biāo)基準(zhǔn)并不統(tǒng)一,還需要統(tǒng)一坐標(biāo)基準(zhǔn)。取潛航器航向方向?yàn)闈摵狡髯鴺?biāo)系的X軸,設(shè)垂直面內(nèi)垂直于X軸向下方向?yàn)閆軸方向,則X,Y,Z軸可構(gòu)成右手坐標(biāo)系(X,Y,Z)q。

    (8)

    式中:R表示兩直角坐標(biāo)系間的旋轉(zhuǎn)矩陣;(ΔX,ΔY,ΔZ)T為兩坐標(biāo)系原點(diǎn)的平移向量。然后可將大地空間直角坐標(biāo)(X,Y,Z)T轉(zhuǎn)換為大地坐標(biāo)(B,L,H)T。具體計(jì)算式可參考文獻(xiàn)[12]。隨后,即可將大地坐標(biāo)轉(zhuǎn)換為海底地形圖坐標(biāo)。

    (9)

    式中:ΔH為大地坐標(biāo)系中大地高與海底地形圖坐標(biāo)系高程之差,B,L為地形點(diǎn)的大地緯度、大地經(jīng)度,D為以理論深度基準(zhǔn)面為起算面且向下為正的垂直分量,即地形點(diǎn)的深度值。

    2.3 水下聲線仿真

    由于潛航器行進(jìn)過程中,不可避免地會(huì)有晃動(dòng),即與參考狀態(tài)間存在一定的姿態(tài)夾角,所以聲吶并不會(huì)向預(yù)設(shè)的角度方向發(fā)射聲波;加之海底地形高低起伏,往往聲線在海底的投射點(diǎn)也相應(yīng)發(fā)生改變。這就需要利用內(nèi)插算法確定格網(wǎng)之間相對(duì)準(zhǔn)確的高程值,而確定投射點(diǎn)位于格網(wǎng)中的平面位置,則要采用迭代的方式精確求解。

    具體仿真過程:在地形格網(wǎng)點(diǎn)中任取一點(diǎn)作為潛航器位置,將該點(diǎn)四周25×25的網(wǎng)格區(qū)域的格網(wǎng)點(diǎn)作為原始聲線的投射點(diǎn)(總計(jì)51×51=2 601個(gè)點(diǎn));而后利用2.2節(jié)步驟進(jìn)行坐標(biāo)轉(zhuǎn)換,以轉(zhuǎn)換后聲線投射點(diǎn)距離潛航器水平面的高差作為觀測(cè)數(shù)據(jù)。

    3 結(jié)果與分析

    3.1 模型構(gòu)建

    分別以載體姿態(tài)、高度為自變量,以定位誤差作為因變量,利用控制變量的方法繪制定位誤差隨姿態(tài)、高度的變化曲線,如圖4、圖5所示。其中圖4為載體距離海底高程1 000 m時(shí)的姿態(tài)角誤差-匹配誤差圖像;圖5為姿態(tài)角誤差σ為0.5°時(shí)的高度-匹配誤差圖像。統(tǒng)計(jì)兩個(gè)自變量與匹配誤差的線性相關(guān)程度,見表1。

    圖4 姿態(tài)角誤差對(duì)定位精度的影響(H=1 000 m)

    圖5 載體高度對(duì)定位精度的影響(σ=0.5°)

    表1 姿態(tài)角和高度關(guān)于定位誤差的線性相關(guān)分析

    從圖4、圖5及表1可知,高度與匹配定位誤差存在顯著的線性關(guān)系,姿態(tài)角誤差與定位誤差的線性程度次之。即文中分析載體高度對(duì)定位精度的影響后,構(gòu)建線性函數(shù)模型。

    在實(shí)驗(yàn)地形數(shù)據(jù)上9個(gè)具有代表性的區(qū)域,各選取一位置作為試驗(yàn)點(diǎn);利用文中仿真方法得到6組不同高程數(shù)據(jù)、7組不同標(biāo)準(zhǔn)差的姿態(tài)角數(shù)據(jù),每組姿態(tài)角生成30個(gè)隨機(jī)姿態(tài)角;共計(jì)11 340條記錄,選取其中4條列于表2。

    表2 不同標(biāo)準(zhǔn)差條件下姿態(tài)角對(duì)于定位的影響(部分結(jié)果)

    根據(jù)前文結(jié)論,采用線性函數(shù)擬合高度-定位誤差曲線??紤]到高程為0時(shí),定位誤差也為0,因此擬合時(shí)線性函數(shù)的零次項(xiàng)設(shè)為0。以姿態(tài)角標(biāo)準(zhǔn)差σ=0.5°為例,繪出各試驗(yàn)點(diǎn)對(duì)應(yīng)的高度-定位誤差曲線見圖6。將9個(gè)試驗(yàn)點(diǎn)的情況取平均,得到該試驗(yàn)底圖中的高度-定位誤差的平均曲線如圖7所示。同理可得其余不同姿態(tài)角標(biāo)準(zhǔn)差擬合得到的對(duì)應(yīng)高度-定位誤差函數(shù)模型一次項(xiàng)系數(shù)見表3,其中p1為一次項(xiàng)系數(shù)。

    圖6 各試驗(yàn)點(diǎn)的高度-定位誤差關(guān)系圖

    圖7 高度-定位誤差的平均線性擬合圖

    根據(jù)表3數(shù)據(jù)繪制一次項(xiàng)系數(shù)p1與σi的關(guān)系如圖8所示,可見兩者構(gòu)成較為明顯的二次函數(shù)關(guān)系。對(duì)p1-σi曲線擬合,可得二次函數(shù)各項(xiàng)系數(shù),見表4,其中,a2,a1,a0分別為p1-σi擬合曲線的二次項(xiàng)、一次項(xiàng)和零次項(xiàng)系數(shù)。

    圖8 p1-σi關(guān)系

    表3 各試驗(yàn)點(diǎn)高度-定位誤差函數(shù)模型參數(shù)

    根據(jù)以上結(jié)論,可推得定位精度與高度、姿態(tài)誤差的關(guān)系模型為:

    ΔP=(a2σ2+a1σ+a0)h.

    (10)

    代入文中實(shí)驗(yàn)數(shù)據(jù),即:

    ΔP=(-0.002 96σ2+0.012 19σ-0.000 35)h.

    (11)

    進(jìn)一步可得,姿態(tài)精度與高度、定位誤差的關(guān)系模型為:

    (12)

    此時(shí),利用潛航器距離海底地形面的高度和潛航器定位精度需求,可推求潛航器姿態(tài)角測(cè)量需滿足的精度。其中σ的負(fù)根被舍棄,式(12)右端其余參數(shù)與表4一致。

    表4 p1-σi模型參數(shù)

    同理有高度與姿態(tài)精度、定位精度的關(guān)系模型:

    (13)

    此時(shí),利用潛航器姿態(tài)角測(cè)量精度指標(biāo)和潛航器對(duì)定位精度的需求,可推求潛航器距離海底地形面需滿足的高度。

    3.2 模型質(zhì)量評(píng)估

    上述3個(gè)關(guān)系式本質(zhì)上是由定位精度與高度、姿態(tài)誤差的關(guān)系模型公式推導(dǎo)得到的,后續(xù)工作僅對(duì)該模型的可靠性進(jìn)行驗(yàn)證。

    實(shí)驗(yàn)在所選取的9個(gè)點(diǎn)圍成的區(qū)域內(nèi)部隨機(jī)選取了1 000個(gè)點(diǎn),每個(gè)點(diǎn)隨機(jī)給定一個(gè)高度值和一個(gè)姿態(tài)角標(biāo)準(zhǔn)差。根據(jù)給定標(biāo)準(zhǔn)差計(jì)算得到30組姿態(tài)角,通過仿真計(jì)算得到30組匹配定位誤差。對(duì)比模型計(jì)算的定位誤差和仿真的參考定位誤差結(jié)果,評(píng)定模型的定位誤差預(yù)測(cè)效果,實(shí)驗(yàn)結(jié)果如圖9所示。因函數(shù)模型計(jì)算的誤差是該區(qū)域的平均水平,所以與參考值相比存在信號(hào)缺失,表現(xiàn)在圖像上顯示為振幅(或能量)弱于參考值。

    圖9 模型計(jì)算的定位誤差與仿真參考誤差

    利用序列相似度公式,可計(jì)算模型計(jì)算誤差與仿真參考誤差之間的相似程度r=81.33%。算式為:

    (14)

    式中:E為仿真得到的匹配定位誤差標(biāo)準(zhǔn)差;E′為模型計(jì)算值;n為檢驗(yàn)點(diǎn)數(shù)。

    另一方面,可根據(jù)文中模型,計(jì)算定位誤差的限差,即2σ的適用性。經(jīng)過計(jì)算,在1 000個(gè)檢驗(yàn)點(diǎn)的匹配誤差中,能夠被該模型計(jì)算的限差值包含在內(nèi)的點(diǎn)有980個(gè),模型適用率為98%,具體結(jié)果如圖10所示。

    圖10 模型計(jì)算定位誤差限差效果

    4 結(jié) 論

    文中圍繞水下定位誤差與載體姿態(tài)、高度的關(guān)系,進(jìn)行仿真實(shí)驗(yàn),得到水下地形匹配定位誤差和載體姿態(tài)、高度的變化關(guān)系,并給出3種參考模型,分別用于不同情況下需求計(jì)算:①已知載體高度及姿態(tài)誤差計(jì)算地形匹配的定位精度;②已知載體高度及定位精度指標(biāo)計(jì)算載體姿態(tài)測(cè)量精度需求;③已知載體姿態(tài)精度指標(biāo)及定位精度指標(biāo)計(jì)算載體需滿足的高度。隨后,針對(duì)構(gòu)建的模型進(jìn)行質(zhì)量評(píng)估,實(shí)驗(yàn)表明模型在隨機(jī)樣本中的預(yù)測(cè)效果優(yōu)于81%,預(yù)測(cè)的定位誤差限差的有效性為98%。

    文中敘述簡化了模型,忽略了潛航器的姿態(tài)漂移量,且認(rèn)為3個(gè)姿態(tài)角的誤差量級(jí)相同;另外,在實(shí)驗(yàn)中未顧及聲線的彎曲效應(yīng)。后續(xù)需要針對(duì)不同量級(jí)且各姿態(tài)角均值不為0的情況進(jìn)一步研究,并完善聲線跟蹤方面的工作。

    猜你喜歡
    底圖定位精度姿態(tài)
    北斗定位精度可達(dá)兩三米
    軍事文摘(2023年4期)2023-04-05 13:57:35
    一種基于實(shí)際GIS底圖的精準(zhǔn)計(jì)費(fèi)方式探究
    攀爬的姿態(tài)
    GPS定位精度研究
    組合導(dǎo)航的AGV定位精度的改善
    全新一代宋的新姿態(tài)
    汽車觀察(2018年9期)2018-10-23 05:46:40
    跑與走的姿態(tài)
    中國自行車(2018年8期)2018-09-26 06:53:44
    底圖檔案管理的科學(xué)化與工作者的綜合素質(zhì)
    企業(yè)底圖檔案的歸檔管理
    陜西檔案(2016年5期)2016-11-26 03:16:46
    如何做好底圖檔案科學(xué)管理
    中文欧美无线码| 久久人妻熟女aⅴ| 国产成人影院久久av| 美女国产高潮福利片在线看| 最近最新中文字幕大全免费视频| 又黄又粗又硬又大视频| 丁香六月欧美| 国产人伦9x9x在线观看| 大香蕉久久成人网| 性高湖久久久久久久久免费观看| 老熟女久久久| 国产精品久久久久成人av| 久久久久视频综合| 亚洲精品久久成人aⅴ小说| 久久久久精品国产欧美久久久 | 黄色 视频免费看| 大型av网站在线播放| av免费在线观看网站| 丝袜人妻中文字幕| 婷婷成人精品国产| av福利片在线| 新久久久久国产一级毛片| 欧美一级毛片孕妇| 久久久精品免费免费高清| 青春草亚洲视频在线观看| 男女床上黄色一级片免费看| 久久国产亚洲av麻豆专区| 亚洲av电影在线进入| 免费在线观看视频国产中文字幕亚洲 | kizo精华| 国产一区二区 视频在线| www.999成人在线观看| 国产成人免费无遮挡视频| 精品久久久久久久毛片微露脸 | 黑人巨大精品欧美一区二区mp4| 成年人黄色毛片网站| 自线自在国产av| 亚洲第一欧美日韩一区二区三区 | 久久热在线av| 亚洲中文av在线| 久久久久网色| 女人久久www免费人成看片| 9191精品国产免费久久| 国产伦人伦偷精品视频| 日韩 亚洲 欧美在线| 免费人妻精品一区二区三区视频| 国产高清视频在线播放一区 | 午夜日韩欧美国产| kizo精华| 丰满人妻熟妇乱又伦精品不卡| 成人av一区二区三区在线看 | 久久人人爽av亚洲精品天堂| 国产欧美日韩精品亚洲av| 亚洲av成人一区二区三| 亚洲av日韩在线播放| 亚洲五月色婷婷综合| 一区二区三区精品91| 18禁观看日本| 十分钟在线观看高清视频www| 久久国产精品男人的天堂亚洲| 亚洲精品国产色婷婷电影| 美女高潮喷水抽搐中文字幕| 中文字幕av电影在线播放| 久久人妻福利社区极品人妻图片| 一边摸一边抽搐一进一出视频| a在线观看视频网站| 在线亚洲精品国产二区图片欧美| 免费不卡黄色视频| 亚洲精品久久成人aⅴ小说| 午夜免费鲁丝| 欧美激情高清一区二区三区| 水蜜桃什么品种好| 免费观看a级毛片全部| 极品少妇高潮喷水抽搐| 免费av中文字幕在线| 亚洲国产日韩一区二区| 一本大道久久a久久精品| 免费人妻精品一区二区三区视频| 国产在线一区二区三区精| 免费久久久久久久精品成人欧美视频| 人人妻人人爽人人添夜夜欢视频| 大片免费播放器 马上看| 欧美精品一区二区大全| 亚洲精品国产区一区二| 亚洲国产欧美网| 在线 av 中文字幕| 免费一级毛片在线播放高清视频 | 亚洲成人国产一区在线观看| 一本综合久久免费| 国产av又大| av欧美777| 99九九在线精品视频| 国产精品偷伦视频观看了| 搡老岳熟女国产| 国产99久久九九免费精品| 永久免费av网站大全| 亚洲 国产 在线| 成人影院久久| 国产成人av激情在线播放| 亚洲国产中文字幕在线视频| 菩萨蛮人人尽说江南好唐韦庄| 91国产中文字幕| 精品福利永久在线观看| 午夜福利免费观看在线| 一级黄色大片毛片| 中文字幕色久视频| 亚洲人成77777在线视频| 亚洲精品在线美女| 美女扒开内裤让男人捅视频| 人人妻人人澡人人爽人人夜夜| 不卡一级毛片| 中文字幕av电影在线播放| 丝袜人妻中文字幕| 欧美97在线视频| 精品高清国产在线一区| 免费人妻精品一区二区三区视频| 欧美97在线视频| 欧美乱码精品一区二区三区| 欧美成狂野欧美在线观看| 80岁老熟妇乱子伦牲交| 一本大道久久a久久精品| 操出白浆在线播放| 黄片小视频在线播放| a级毛片黄视频| 一本—道久久a久久精品蜜桃钙片| 久久久久精品人妻al黑| 亚洲五月色婷婷综合| 老司机午夜福利在线观看视频 | 欧美乱码精品一区二区三区| 亚洲熟女精品中文字幕| 午夜福利一区二区在线看| 久久精品亚洲熟妇少妇任你| 亚洲国产av影院在线观看| 精品少妇黑人巨大在线播放| 午夜激情av网站| 狂野欧美激情性xxxx| 国产精品九九99| 999久久久精品免费观看国产| 女人精品久久久久毛片| 午夜福利免费观看在线| 亚洲欧美激情在线| 国产男女内射视频| 俄罗斯特黄特色一大片| 久久人人97超碰香蕉20202| 淫妇啪啪啪对白视频 | 欧美精品人与动牲交sv欧美| 国产成人系列免费观看| 精品高清国产在线一区| 日本一区二区免费在线视频| 人妻一区二区av| 亚洲精品成人av观看孕妇| 国产精品免费大片| 狠狠狠狠99中文字幕| 国产精品影院久久| 午夜日韩欧美国产| 国产成人精品在线电影| 国产深夜福利视频在线观看| 啦啦啦视频在线资源免费观看| 69精品国产乱码久久久| 亚洲精品美女久久av网站| 国产av精品麻豆| 欧美日韩一级在线毛片| 日韩中文字幕视频在线看片| 日本av免费视频播放| 男女床上黄色一级片免费看| 热re99久久精品国产66热6| 搡老乐熟女国产| 国产黄频视频在线观看| 99精品欧美一区二区三区四区| 动漫黄色视频在线观看| 亚洲情色 制服丝袜| 日本av免费视频播放| 国产在线视频一区二区| 天天躁夜夜躁狠狠躁躁| 大片免费播放器 马上看| 少妇人妻久久综合中文| 老司机影院毛片| 少妇精品久久久久久久| 超碰成人久久| 国产亚洲午夜精品一区二区久久| 亚洲一区二区三区欧美精品| 欧美午夜高清在线| 不卡一级毛片| 99久久综合免费| 亚洲免费av在线视频| 男女下面插进去视频免费观看| 午夜福利乱码中文字幕| 99热全是精品| 国产av又大| 午夜福利免费观看在线| 国产在线观看jvid| 久久国产精品大桥未久av| 精品免费久久久久久久清纯 | 午夜91福利影院| av不卡在线播放| 亚洲欧洲精品一区二区精品久久久| 后天国语完整版免费观看| 午夜福利免费观看在线| 亚洲精品在线美女| 国产精品免费大片| 欧美日韩视频精品一区| 亚洲av成人不卡在线观看播放网 | 一区二区三区激情视频| 真人做人爱边吃奶动态| 欧美一级毛片孕妇| 嫁个100分男人电影在线观看| 亚洲伊人久久精品综合| 久久这里只有精品19| 欧美一级毛片孕妇| 丰满人妻熟妇乱又伦精品不卡| 久久热在线av| 91精品国产国语对白视频| 国产成+人综合+亚洲专区| 亚洲欧美日韩高清在线视频 | 欧美在线一区亚洲| 欧美日本中文国产一区发布| 少妇精品久久久久久久| 日韩中文字幕欧美一区二区| 91精品国产国语对白视频| 熟女少妇亚洲综合色aaa.| 中文欧美无线码| 欧美日韩亚洲综合一区二区三区_| 人妻 亚洲 视频| 少妇粗大呻吟视频| 免费在线观看视频国产中文字幕亚洲 | 国产片内射在线| 久久精品久久久久久噜噜老黄| 亚洲国产精品一区二区三区在线| 亚洲第一青青草原| 精品人妻1区二区| 一边摸一边做爽爽视频免费| 美女福利国产在线| 午夜成年电影在线免费观看| 波多野结衣一区麻豆| 亚洲国产精品一区二区三区在线| 人人妻人人爽人人添夜夜欢视频| 十八禁高潮呻吟视频| 国产精品免费视频内射| 国产成人a∨麻豆精品| 国产激情久久老熟女| 欧美人与性动交α欧美精品济南到| 国产伦人伦偷精品视频| 亚洲久久久国产精品| 国产精品久久久av美女十八| 婷婷成人精品国产| 亚洲自偷自拍图片 自拍| 午夜免费鲁丝| av网站在线播放免费| 欧美日韩福利视频一区二区| 男女下面插进去视频免费观看| 欧美日本中文国产一区发布| 日韩视频在线欧美| 每晚都被弄得嗷嗷叫到高潮| 手机成人av网站| 精品卡一卡二卡四卡免费| 免费av中文字幕在线| 日本撒尿小便嘘嘘汇集6| 国产黄频视频在线观看| 亚洲精品中文字幕一二三四区 | 99国产极品粉嫩在线观看| 国产视频一区二区在线看| 亚洲精品国产av成人精品| 精品国产乱子伦一区二区三区 | 每晚都被弄得嗷嗷叫到高潮| 欧美日韩精品网址| 91九色精品人成在线观看| 欧美大码av| 国产亚洲欧美精品永久| 亚洲欧美成人综合另类久久久| 国产黄色免费在线视频| 欧美xxⅹ黑人| 国产精品一区二区在线观看99| 精品乱码久久久久久99久播| 99国产精品一区二区蜜桃av | 国产av又大| 中国美女看黄片| 中文字幕另类日韩欧美亚洲嫩草| 日韩三级视频一区二区三区| 精品人妻在线不人妻| 极品少妇高潮喷水抽搐| 12—13女人毛片做爰片一| 国产亚洲av高清不卡| 老司机影院成人| 91九色精品人成在线观看| 国产av国产精品国产| 国产精品一区二区在线不卡| 午夜激情久久久久久久| 两人在一起打扑克的视频| 久久香蕉激情| 国产精品成人在线| 一级,二级,三级黄色视频| 亚洲欧美一区二区三区黑人| 老司机影院毛片| 女人久久www免费人成看片| 深夜精品福利| 欧美老熟妇乱子伦牲交| h视频一区二区三区| 丝袜喷水一区| 黑人巨大精品欧美一区二区蜜桃| 不卡一级毛片| 韩国高清视频一区二区三区| 老熟妇乱子伦视频在线观看 | 国产亚洲av高清不卡| 精品福利观看| 岛国在线观看网站| 老司机影院成人| 99国产精品一区二区三区| 中文字幕制服av| 啦啦啦 在线观看视频| 久久天躁狠狠躁夜夜2o2o| 首页视频小说图片口味搜索| 国产99久久九九免费精品| 久久国产精品影院| 日韩视频在线欧美| 91麻豆精品激情在线观看国产 | 日韩一区二区三区影片| 欧美成人午夜精品| 777久久人妻少妇嫩草av网站| 色老头精品视频在线观看| 国产精品熟女久久久久浪| 在线 av 中文字幕| 9热在线视频观看99| 亚洲欧美精品综合一区二区三区| 国产精品熟女久久久久浪| 久久精品亚洲熟妇少妇任你| 水蜜桃什么品种好| 欧美黑人精品巨大| 黄色视频,在线免费观看| 精品熟女少妇八av免费久了| 老司机福利观看| av网站在线播放免费| 夜夜骑夜夜射夜夜干| 亚洲av成人不卡在线观看播放网 | 国产极品粉嫩免费观看在线| 日本撒尿小便嘘嘘汇集6| 日韩中文字幕视频在线看片| 一级片'在线观看视频| 男女高潮啪啪啪动态图| 欧美久久黑人一区二区| 一级毛片精品| 搡老熟女国产l中国老女人| 免费日韩欧美在线观看| 一二三四在线观看免费中文在| 国产日韩欧美视频二区| 成人av一区二区三区在线看 | 久久久久国内视频| 青草久久国产| 别揉我奶头~嗯~啊~动态视频 | 少妇 在线观看| 婷婷成人精品国产| 亚洲伊人色综图| 国产真人三级小视频在线观看| 在线永久观看黄色视频| 亚洲精品美女久久久久99蜜臀| 成人亚洲精品一区在线观看| 久久人妻熟女aⅴ| 亚洲一区二区三区欧美精品| 波多野结衣一区麻豆| 99久久精品国产亚洲精品| 日韩中文字幕视频在线看片| 久久久精品94久久精品| 久久久久久久精品精品| 女人爽到高潮嗷嗷叫在线视频| 手机成人av网站| 欧美日韩黄片免| 国产精品影院久久| 在线观看免费视频网站a站| 黄片小视频在线播放| 午夜福利在线免费观看网站| 欧美精品一区二区大全| 妹子高潮喷水视频| 精品亚洲乱码少妇综合久久| 日韩 欧美 亚洲 中文字幕| 亚洲,欧美精品.| 久久女婷五月综合色啪小说| 人成视频在线观看免费观看| 大香蕉久久成人网| 99久久精品国产亚洲精品| 国产精品欧美亚洲77777| 黄色怎么调成土黄色| 高清欧美精品videossex| 亚洲专区国产一区二区| 老熟妇仑乱视频hdxx| 久久狼人影院| 丝袜喷水一区| 在线观看人妻少妇| 宅男免费午夜| 一本一本久久a久久精品综合妖精| 精品视频人人做人人爽| 欧美在线黄色| 国产又爽黄色视频| 亚洲精品美女久久av网站| 狠狠狠狠99中文字幕| tocl精华| 18禁裸乳无遮挡动漫免费视频| 亚洲国产精品一区三区| 一本久久精品| 国产精品久久久久久精品古装| 十八禁高潮呻吟视频| 宅男免费午夜| 99国产综合亚洲精品| 超色免费av| 久久久久久亚洲精品国产蜜桃av| 成人影院久久| 夜夜骑夜夜射夜夜干| 啦啦啦啦在线视频资源| 亚洲成人免费av在线播放| 99精品欧美一区二区三区四区| 中文字幕另类日韩欧美亚洲嫩草| 中文欧美无线码| 婷婷丁香在线五月| 99国产精品99久久久久| 亚洲欧美日韩另类电影网站| 精品亚洲乱码少妇综合久久| 欧美老熟妇乱子伦牲交| 免费观看av网站的网址| 亚洲色图 男人天堂 中文字幕| 老司机午夜福利在线观看视频 | 天天添夜夜摸| 欧美成狂野欧美在线观看| 国产男人的电影天堂91| 中文字幕av电影在线播放| 97精品久久久久久久久久精品| 国产日韩欧美亚洲二区| 日本精品一区二区三区蜜桃| 交换朋友夫妻互换小说| 午夜老司机福利片| 色视频在线一区二区三区| 亚洲午夜精品一区,二区,三区| 91字幕亚洲| 99国产精品免费福利视频| 日韩人妻精品一区2区三区| 成人国产av品久久久| 亚洲中文av在线| 久久性视频一级片| 欧美中文综合在线视频| 亚洲欧美一区二区三区久久| 悠悠久久av| 欧美午夜高清在线| 纵有疾风起免费观看全集完整版| 日韩一卡2卡3卡4卡2021年| 亚洲 欧美一区二区三区| 一边摸一边做爽爽视频免费| 日韩 欧美 亚洲 中文字幕| 大陆偷拍与自拍| 久久精品久久久久久噜噜老黄| 国产一区二区 视频在线| 97人妻天天添夜夜摸| 国产精品久久久av美女十八| 精品久久久久久久毛片微露脸 | 国产成人a∨麻豆精品| 午夜激情av网站| 欧美97在线视频| 久久人妻熟女aⅴ| 国产区一区二久久| 大片免费播放器 马上看| 午夜精品国产一区二区电影| 亚洲精品美女久久久久99蜜臀| 欧美亚洲 丝袜 人妻 在线| 欧美另类一区| 欧美精品亚洲一区二区| 在线观看免费视频网站a站| 久久99一区二区三区| 日韩中文字幕视频在线看片| 一本久久精品| 亚洲精品久久成人aⅴ小说| 精品欧美一区二区三区在线| 精品人妻1区二区| 亚洲人成电影观看| 亚洲精品国产av成人精品| av视频免费观看在线观看| 免费在线观看黄色视频的| 欧美人与性动交α欧美软件| 99国产极品粉嫩在线观看| 天堂中文最新版在线下载| 亚洲欧美日韩高清在线视频 | 操美女的视频在线观看| 国产成人精品久久二区二区免费| 亚洲伊人色综图| 国产有黄有色有爽视频| 亚洲第一欧美日韩一区二区三区 | 欧美日韩成人在线一区二区| 9色porny在线观看| 欧美亚洲日本最大视频资源| 老汉色av国产亚洲站长工具| 中文字幕人妻熟女乱码| 啦啦啦在线免费观看视频4| 午夜精品国产一区二区电影| 色94色欧美一区二区| av片东京热男人的天堂| 美国免费a级毛片| 丝袜喷水一区| 久久99一区二区三区| 夫妻午夜视频| 国产精品一二三区在线看| 一二三四社区在线视频社区8| 麻豆国产av国片精品| 日日爽夜夜爽网站| 新久久久久国产一级毛片| 黄频高清免费视频| 亚洲精品国产一区二区精华液| a 毛片基地| 精品乱码久久久久久99久播| 亚洲天堂av无毛| 亚洲av片天天在线观看| 国产一区有黄有色的免费视频| 啪啪无遮挡十八禁网站| 久久热在线av| 日韩制服丝袜自拍偷拍| 久久人妻熟女aⅴ| svipshipincom国产片| 久久久久久免费高清国产稀缺| 亚洲中文av在线| 一本色道久久久久久精品综合| 曰老女人黄片| 欧美黑人欧美精品刺激| 日本wwww免费看| 国产老妇伦熟女老妇高清| 欧美午夜高清在线| av福利片在线| 免费在线观看视频国产中文字幕亚洲 | 18禁国产床啪视频网站| 操出白浆在线播放| 免费日韩欧美在线观看| 亚洲精品一二三| 三级毛片av免费| 久久久久网色| 国产精品二区激情视频| 精品一区二区三区四区五区乱码| 色婷婷久久久亚洲欧美| 久久天堂一区二区三区四区| 女警被强在线播放| a级毛片黄视频| 老汉色av国产亚洲站长工具| 最新在线观看一区二区三区| 狂野欧美激情性bbbbbb| 岛国毛片在线播放| 人妻久久中文字幕网| 亚洲精品国产区一区二| 中文字幕另类日韩欧美亚洲嫩草| 正在播放国产对白刺激| 日韩欧美国产一区二区入口| 亚洲av欧美aⅴ国产| 啦啦啦在线免费观看视频4| 在线观看一区二区三区激情| 这个男人来自地球电影免费观看| 精品少妇内射三级| 日韩欧美一区视频在线观看| 精品国产乱码久久久久久小说| av不卡在线播放| 人人妻人人爽人人添夜夜欢视频| 国产成人免费观看mmmm| 亚洲精品成人av观看孕妇| 色94色欧美一区二区| 丰满人妻熟妇乱又伦精品不卡| 亚洲全国av大片| 亚洲精华国产精华精| 中文字幕av电影在线播放| 国产精品香港三级国产av潘金莲| 99精品久久久久人妻精品| 欧美日韩精品网址| 视频区欧美日本亚洲| 一级毛片电影观看| av在线app专区| 青草久久国产| 在线 av 中文字幕| av有码第一页| 欧美精品av麻豆av| 波多野结衣av一区二区av| 欧美日韩视频精品一区| 美女国产高潮福利片在线看| 五月开心婷婷网| av在线播放精品| 精品国产一区二区三区久久久樱花| 欧美乱码精品一区二区三区| 欧美精品人与动牲交sv欧美| 不卡一级毛片| 国产一区二区在线观看av| 亚洲第一av免费看| 日本vs欧美在线观看视频| 少妇 在线观看| 久久人妻熟女aⅴ| 免费一级毛片在线播放高清视频 | 亚洲一码二码三码区别大吗| 午夜福利视频精品| 人妻一区二区av| 夫妻午夜视频| 午夜影院在线不卡| 国产精品秋霞免费鲁丝片| 国产日韩欧美视频二区| 亚洲av日韩在线播放| 亚洲av片天天在线观看| 亚洲精品美女久久久久99蜜臀| 国产精品欧美亚洲77777| 日本91视频免费播放| 18禁黄网站禁片午夜丰满| 性色av一级| 亚洲激情五月婷婷啪啪| 热99国产精品久久久久久7| 欧美日韩中文字幕国产精品一区二区三区 | 大陆偷拍与自拍| 日日夜夜操网爽| 黄片小视频在线播放| 国产精品 欧美亚洲| 啦啦啦免费观看视频1| 色婷婷av一区二区三区视频| www日本在线高清视频| 欧美老熟妇乱子伦牲交| 欧美中文综合在线视频| 国产精品久久久久久精品古装| 99精国产麻豆久久婷婷| 在线观看免费高清a一片| www.999成人在线观看| 国产一区二区 视频在线| 91麻豆精品激情在线观看国产 |