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

    基于解耦消參的動(dòng)對(duì)動(dòng)相對(duì)定位算法

    2015-10-24 02:20:54伍勁實(shí)趙修斌龐春雷余永林
    關(guān)鍵詞:雙差浮點(diǎn)歷元

    伍勁實(shí),趙修斌,龐春雷,余永林,劉 明

    (空軍工程大學(xué)信息與導(dǎo)航學(xué)院,陜西西安710077)

    基于解耦消參的動(dòng)對(duì)動(dòng)相對(duì)定位算法

    伍勁實(shí),趙修斌,龐春雷,余永林,劉 明

    (空軍工程大學(xué)信息與導(dǎo)航學(xué)院,陜西西安710077)

    針對(duì)動(dòng)對(duì)動(dòng)相對(duì)定位基線參量實(shí)時(shí)可變導(dǎo)致模糊度難以準(zhǔn)確求解的問(wèn)題,在比較分析動(dòng)對(duì)動(dòng)與靜態(tài)相對(duì)定位區(qū)別與聯(lián)系的基礎(chǔ)上,提出了基于解耦消參的動(dòng)對(duì)動(dòng)相對(duì)定位算法。該算法將解耦法與最小二乘模糊度降相關(guān)平差(LAMBDA)算法相結(jié)合,采用解耦法變換雙差載波相位觀測(cè)方程,從而消除基線參量并實(shí)時(shí)推算雙差整周模糊度的浮點(diǎn)解及其協(xié)方差矩陣,然后利用LAMBDA算法對(duì)模糊度進(jìn)行固定和確認(rèn)。試驗(yàn)結(jié)果表明,基于解耦消參的動(dòng)對(duì)動(dòng)相對(duì)定位算法初始化模糊度用時(shí)較短,在70 s左右即可正確解算出動(dòng)態(tài)單頻整周模糊度,獲得了厘米級(jí)的相對(duì)定位結(jié)果,適用于短基線高精度動(dòng)態(tài)相對(duì)定位。

    動(dòng)對(duì)動(dòng);相對(duì)定位;整周模糊度;解耦;LAMBDA

    0 引言

    動(dòng)對(duì)動(dòng)高精度相對(duì)定位的關(guān)鍵是在動(dòng)態(tài)情況下快速準(zhǔn)確地確定整周模糊度[1]。由于動(dòng)對(duì)動(dòng)相對(duì)定位基線矢量實(shí)時(shí)可變,屬于變參數(shù)估計(jì)問(wèn)題,且難以找到附加約束條件進(jìn)行輔助[2-3],因此其整周模糊度及其協(xié)方差矩陣的初始化過(guò)程相比于靜態(tài)更為復(fù)雜。目前常用的模糊度浮點(diǎn)解求解方法主要有最小二乘法、遞推最小二乘法、卡爾曼濾波法[4]、雙頻載波相位的直接解法[5]、矩陣解耦法[6]等,其中能夠適用于動(dòng)態(tài)相對(duì)定位的方法主要是最小二乘法、卡爾曼濾波法和矩陣解耦法。最小二乘法屬于批處理算法,難以滿足動(dòng)態(tài)相對(duì)定位的實(shí)時(shí)性要求,且計(jì)算量大;卡爾曼濾波法對(duì)動(dòng)態(tài)模型的建立要求較高,若模型的建立與實(shí)際運(yùn)動(dòng)狀態(tài)不符合,則待估參量的精度難以達(dá)到要求,進(jìn)而難以實(shí)現(xiàn)高精度相對(duì)定位;文獻(xiàn)[6]提出了基于矩陣解耦的模糊度逼近搜索法,該方法能夠較好地用于模糊度浮點(diǎn)解及其協(xié)方差矩陣的實(shí)時(shí)推算,但其對(duì)模糊度浮點(diǎn)解的精度要求較高,需要較長(zhǎng)的時(shí)間才能逼近模糊度真值。本文針對(duì)此問(wèn)題,在比較分析靜態(tài)與動(dòng)態(tài)相對(duì)定位區(qū)別與聯(lián)系的基礎(chǔ)上,提出了基于解耦消參的動(dòng)對(duì)動(dòng)相對(duì)定位算法。

    1 傳統(tǒng)相對(duì)定位模糊度求解

    在短基線條件下,雙差載波相位觀測(cè)方程為:

    其中,φij為雙差載波相位觀測(cè)量,λ為載波波長(zhǎng),l為接收機(jī)至衛(wèi)星的單位矢量,b=bxbyb( )zΤ為基線矢量,Nij為雙差整周模糊度,εij為雙差測(cè)量噪聲。

    若參考站與移動(dòng)站相對(duì)靜止,第k個(gè)歷元觀測(cè)到n顆衛(wèi)星,可構(gòu)造n—1個(gè)雙差觀測(cè)方程,忽略雙差噪聲,將式(1)變?yōu)榉匠探M形式可得

    將k個(gè)歷元雙差觀測(cè)方程進(jìn)行組合,可得到k×(n—1)個(gè)雙差方程:

    由于方程組中含有(n—1)+3=n+2個(gè)未知參量,只需要滿足條件k×(n—1)≥n+2,即可用最小二乘法求解基線矢量和雙差整周模糊度的浮點(diǎn)解。但由于觀測(cè)時(shí)間較短會(huì)造成法矩陣的病態(tài)性[7],因此靜態(tài)相對(duì)定位時(shí),往往通過(guò)增加觀測(cè)時(shí)間來(lái)獲得高精度定位結(jié)果;同時(shí)由于測(cè)站之間相對(duì)靜止,通??删_測(cè)量測(cè)站間基線長(zhǎng)度,將其作為已知約束條件輔助解算整周模糊度[8],提高定位精度和效率。在求解得到模糊度浮點(diǎn)解的前提下,可利用LAMBDA算法對(duì)模糊度進(jìn)行搜索與固定,模糊度一旦正確固定,則可將其帶入雙差方程反解出基線矢量,實(shí)現(xiàn)高精度相對(duì)定位。這種在靜態(tài)情況下求解模糊度的方法一般稱為整周模糊度的靜態(tài)初始化。當(dāng)前,高精度動(dòng)態(tài)相對(duì)定位往往也是采用靜態(tài)初始化的方法,即在靜態(tài)情況下固定住模糊度,然后再開始運(yùn)動(dòng),但這樣不適用于動(dòng)態(tài)情況下的實(shí)時(shí)解算,因此需要尋求模糊度的動(dòng)態(tài)初始化方法。

    2 基于解耦消參的動(dòng)對(duì)動(dòng)相對(duì)定位算法

    2.1 動(dòng)對(duì)動(dòng)模糊度求解特點(diǎn)

    在動(dòng)對(duì)動(dòng)情況下,如果正確解算出模糊度,且不發(fā)生周跳和衛(wèi)星失鎖的情況,則雙差模糊度參量N不隨時(shí)間的變化而變化,但由于運(yùn)載體之間的相對(duì)位置實(shí)時(shí)變化,式(1)中的基線矢量參量b隨時(shí)間實(shí)時(shí)變化。將k個(gè)歷元的雙差觀測(cè)方程進(jìn)行組合,則式(3)變?yōu)槭剑?)。

    由式(4)可知,每一個(gè)歷元均會(huì)增加3個(gè)未知基線矢量參量,則k個(gè)歷元得到的方程組中包含3k+n—1個(gè)未知參量,此時(shí)需要滿足條件k×(n—1)≥3k +n—1;而且隨著觀測(cè)時(shí)間的增加,未知參量個(gè)數(shù)線性增加,觀測(cè)方程系數(shù)矩陣維數(shù)越來(lái)越大,傳統(tǒng)的最小二乘法用于靜態(tài)定位時(shí),模糊度初始化需要至少2 min,如果將其用于動(dòng)對(duì)動(dòng),則會(huì)產(chǎn)生至少360個(gè)未知基線參量,此時(shí)將難以利用最小二乘求解,同時(shí),基線矢量實(shí)時(shí)變化,難以再作為已知約束條件進(jìn)行輔助,這時(shí)應(yīng)尋求適合于動(dòng)態(tài)情況下的模糊度初始化方法。

    由于在實(shí)際求解雙差模糊度過(guò)程中,對(duì)于基線矢量bk并不關(guān)心,因?yàn)橐坏?zhǔn)確求解并固定雙差模糊度,基線矢量就能精確求得。為了滿足動(dòng)對(duì)動(dòng)對(duì)實(shí)時(shí)性的要求,并且充分利用每一歷元的觀測(cè)信息,本文采用基于解耦消參的動(dòng)對(duì)動(dòng)相對(duì)定位算法求解動(dòng)對(duì)動(dòng)模糊度。

    2.2 動(dòng)對(duì)動(dòng)相對(duì)定位算法

    解耦法的基本思想為:將每一歷元的基線矢量參量bk用雙差模糊度參數(shù)N表示,從而消除bk,只剩下N,進(jìn)而僅對(duì)雙差模糊度進(jìn)行推算[9]。

    則第k個(gè)歷元雙差方程為

    基線矢量的最小二乘解為

    根據(jù)式(5),對(duì)前n個(gè)歷元得到的雙差觀測(cè)方程疊加得到

    將式(6)代入式(7)得

    則有

    則模糊度浮點(diǎn)解可表示為:

    根據(jù)誤差傳播定律,可求得模糊度的協(xié)方差矩陣為:

    式(9)和式(10)即為解耦法推算模糊度浮點(diǎn)解及其協(xié)方差矩陣的公式。

    經(jīng)分析可知,采用解耦法實(shí)時(shí)推算模糊度浮點(diǎn)解及其協(xié)方差矩陣,是將當(dāng)前時(shí)刻之前每一歷元的觀測(cè)信息進(jìn)行求和,隨著觀測(cè)時(shí)間的增加,并沒(méi)有造成矩陣維數(shù)的增加。在用解耦法得到雙差模糊度浮點(diǎn)解及其協(xié)方差矩陣的基礎(chǔ)上,利用LAMBDA[10]算法進(jìn)行模糊度的搜索與固定,從而獲得高精度動(dòng)對(duì)動(dòng)相對(duì)定位結(jié)果,滿足動(dòng)對(duì)動(dòng)相對(duì)定位的實(shí)時(shí)性要求。

    從式(9)和式(10)中可以看出,在忽略觀測(cè)噪聲的情況下,雙差模糊度浮點(diǎn)解的求解僅與載波相位觀測(cè)值、接收機(jī)至衛(wèi)星的方向矢量以及載波波長(zhǎng)有關(guān),而接收機(jī)的運(yùn)動(dòng)狀態(tài)雖然對(duì)接收機(jī)的捕獲、跟蹤性能影響較大,但對(duì)接收機(jī)至衛(wèi)星的方向矢量影響較小,理論上,只要接收機(jī)能連續(xù)跟蹤衛(wèi)星信號(hào),提供可靠的載波相位觀測(cè)值,則該算法就具有適用性。

    3 試驗(yàn)結(jié)果與分析

    3.1 試驗(yàn)方法

    本文從兩個(gè)方面驗(yàn)證該算法的可行性:1)如果整個(gè)觀測(cè)過(guò)程不發(fā)生周跳,那么動(dòng)態(tài)整周模糊度應(yīng)與靜態(tài)整周模糊度相同,基于這一思想,先用傳統(tǒng)成熟的方法解算出靜態(tài)時(shí)的模糊度,將其作為參考值,結(jié)果如表1所示,然后用基于解耦消參的動(dòng)對(duì)動(dòng)相對(duì)定位算法對(duì)動(dòng)態(tài)情況下獲得的觀測(cè)數(shù)據(jù)進(jìn)行解算,將求解得到的動(dòng)態(tài)模糊度與參考值作比較,若該算法可行,則二者應(yīng)該一致。2)利用基于解耦消參的動(dòng)對(duì)動(dòng)相對(duì)定位算法求解得到的動(dòng)態(tài)模糊度反解基線長(zhǎng)度,將結(jié)果與其實(shí)際長(zhǎng)度進(jìn)行比較,期望獲得較小的基線誤差。

    表1 靜態(tài)時(shí)解算出的雙差整周模糊度參考值Tab.1 Referencevaluefordoubledifference integerambiguityinstaticcalculation

    若能同時(shí)獲得以上兩個(gè)預(yù)期效果,則說(shuō)明所提算法能較好地適用于動(dòng)態(tài)相對(duì)定位。

    3.2 試驗(yàn)條件

    為驗(yàn)證該算法的準(zhǔn)確性和可行性,采用固定基線情況下動(dòng)態(tài)試驗(yàn)的方法。利用兩個(gè)NovAtel接收機(jī),基準(zhǔn)站和移動(dòng)站分別為SPAN-CPT和OEM628板卡,均連接GPS-703-GGG型號(hào)的雙頻天線,數(shù)據(jù)采樣率為1Hz,兩天線固定于已知長(zhǎng)度的基線兩端,并放置于小車前端,試驗(yàn)前已測(cè)得基線長(zhǎng)為1.91m,取衛(wèi)星截止高度角為15°,觀測(cè)到7顆GPS衛(wèi)星,分別為PRN3,PRN14,PRN20,PRN22,PRN25,PRN31,PRN32,選取PRN3作為參考衛(wèi)星。先靜態(tài)觀測(cè)大約10min,然后開始動(dòng)態(tài)試驗(yàn),運(yùn)動(dòng)速度約為3m/s,繞操場(chǎng)運(yùn)動(dòng)大約5min。

    3.3 試驗(yàn)過(guò)程

    采用解耦法對(duì)動(dòng)態(tài)數(shù)據(jù)推算約1min,這里以雙差模糊度N31為代表(其他雙差模糊度的變化與解算情況與之類似),其浮點(diǎn)解在推算中逼近真值的過(guò)程如圖1所示;然后采用LAMBDA算法固定每一歷元的雙差模糊度,表2和圖2為N31在每一歷元的固定結(jié)果;利用固定得到的雙差模糊度反解每一歷元的基線長(zhǎng)度,所得結(jié)果如圖3所示;將求解得到的基線長(zhǎng)度與實(shí)測(cè)長(zhǎng)度相比較,得到基線誤差曲線圖如圖4所示。

    表2 雙差模糊度N31在每一歷元的固定結(jié)果Tab.2 Thefixedresultfordoubledifference ambiguityN31ineachepoch

    圖1 模糊度浮點(diǎn)解N31逼近真值的過(guò)程Fig.1 FloatsolutionN31 approachingrealvalue

    圖2 模糊度N31在每一歷元的固定結(jié)果Fig.2 FixedambiguityN31 ineachepoch

    圖3 基線長(zhǎng)度Fig.3 Baselinelength

    圖4 基線長(zhǎng)度誤差Fig.4 Baselinelengtherror

    3.4 結(jié)果分析

    從圖1中可以看出,當(dāng)用解耦法推算到大概第70個(gè)歷元時(shí),模糊度N31的浮點(diǎn)解已經(jīng)開始逼近真值,并且之后的推算較為平穩(wěn);結(jié)合表1,由表2和圖2可知,從第69個(gè)歷元開始,在動(dòng)態(tài)情況下的每一歷元均能準(zhǔn)確固定住模糊度N31,該方法首次固定模糊度用時(shí)相比于基于模糊度逼近的搜索法用時(shí)更少的原因,是因?yàn)楹笳邔?duì)模糊度浮點(diǎn)解的精度要求更高,其模糊度解算時(shí)間主要消耗于浮點(diǎn)解逼近真值的過(guò)程;而本文將解耦法與LAMBDA方法相結(jié)合,LAMBDA方法對(duì)浮點(diǎn)解的精度要求不如基于模糊度逼近的搜索法苛刻,因而節(jié)省了浮點(diǎn)解的求解時(shí)間;由圖3和圖4可知,利用解算得到的模糊度反解所得基線長(zhǎng)度與實(shí)驗(yàn)前測(cè)量結(jié)果基本一致,基線長(zhǎng)度誤差為1 cm。

    4 結(jié)論

    本文在分析動(dòng)對(duì)動(dòng)與靜態(tài)相對(duì)定位中整周模糊度初始化的區(qū)別與聯(lián)系的基礎(chǔ)上,提出了基于解耦消參的動(dòng)對(duì)動(dòng)相對(duì)定位算法。該算法將基線參量用雙差模糊度參量代替,僅對(duì)模糊度浮點(diǎn)解及其協(xié)方差矩陣進(jìn)行實(shí)時(shí)推算,然后利用LAMBDA算法進(jìn)行模糊度整數(shù)解的搜索與固定。試驗(yàn)表明,基于解耦消參的動(dòng)對(duì)動(dòng)相對(duì)定位算法初始化模糊度用時(shí)較短,在70 s左右即可正確解算出動(dòng)態(tài)單頻整周模糊度,獲得了厘米級(jí)的相對(duì)定位結(jié)果,能夠較好地適用于短基線高精度動(dòng)態(tài)相對(duì)定位。

    [1]Jacek Paziewski,Pawel Wielgosz.Assessment of GPS+Galileo and multi-frequency Galileo single-epoch precise positioning with network corrections[J].GPS Solution,2014,18(4):571-579.

    [2]王德軍,熊永良,劉寧,等.寬巷組合及軌跡約束下的列車在軌動(dòng)態(tài)單歷元定位算法[J].測(cè)繪學(xué)報(bào),2014,43(2):131-136.

    [3]Wang Bo,Miao Lingiuan,Wang Shunting.A constrained LAMBDA method for GPS attitude determination[J]. GPSSolution,2009(13):97-107.

    [4]段榮,趙修斌,龐春雷,等.改進(jìn)的GPS整周模糊度動(dòng)態(tài)快速解算[J].華中科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,42(9):54-58.

    [5]邱蕾,花向紅,蔡華,等.GPS短基線整周模糊度的直接解法[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版),2009,34(1):97-99.

    [6]俞文伯,高國(guó)江,趙剡.單頻GPS動(dòng)態(tài)相對(duì)定位的模糊度逼近/搜索解法[J].北京航空航天大學(xué)學(xué)報(bào),2002,28(2):242-244.

    [7]龐春雷,趙修斌,盧艷娥,等.一種改進(jìn)型的GPS單頻整周模糊度快速解算方法[J].兵工學(xué)報(bào),2012,33(11):1387-1392.

    [8]Chen Wantong,Qin Honglei.New method for single epoch,single frequency land vehicle attitude determination using low-end GPS receiver[J].GPS Solution,2012,(16):329-338.

    [9]劉志儉.GPS載波相位差分技術(shù)、捷聯(lián)慣性導(dǎo)航系統(tǒng)初始對(duì)準(zhǔn)技術(shù)及其組合技術(shù)研究[D].長(zhǎng)沙:國(guó)防科學(xué)技術(shù)大學(xué),2003.

    [10]Teunissen PJG.The least-squares ambiguity decorrelation adiustment:a method for fast GPSinteger ambiguity estimation[J].Journal of Geodesy,1995(70):65-82.

    Dynamic Relative Positioning Algorithm Based on Decoupling Parameter Elimination

    WU Shaoshi,ZHAO Xiubin,PANG Chunlei,YU Yonglin,LIU Ming
    (Information and Navigation College,Air Force Engineering University,Xi’an 710077,China)

    It is difficult to obtain the integer ambiguity precisely because baseline is real-time changeable in dynamic to dynamic relative positioning,.On the basis of comparison between dynamic to dynamic relative positioning and static relative positioning,an algorithm for dynamic to dynamic relative positioning based on parameter elimination with decoupling method was proposed.The algorithm combined the decoupling method and LAMBDA algorithm.Firstly,double difference carrier phase observation equations were transformed through decoupling method,so that the baseline parameters were removed and the float solution and its covariance matrix of double difference of integer ambiguity were deduced in real-time,then integer ambiguities were fixed and confirmed by LAMBDA algorithm.The test revealed that it took only seventy seconds to calculate the single tone dynamic ambiguity with the proposed algorithm and the baseline error being cm-level,which meant that the proposed algorithm was suitable for short-baseline dynamic scenarios with high-precision.

    dynamic to dynamic;relative positioning;integer ambiguity;decoupling;LAMBDA

    TN967.1

    A

    1008-1194(2015)05-0026-04

    2015-05-25

    國(guó)家自然科學(xué)基金資助項(xiàng)目(61273049)

    伍劭實(shí)(1990—),男,湖北宜昌人,碩士研究生,研究方向:衛(wèi)星導(dǎo)航與定位。E-mail:wushaoshipaper@ 163.com。

    猜你喜歡
    雙差浮點(diǎn)歷元
    虛擬地震臺(tái)陣雙差測(cè)深法及應(yīng)用
    LEO星座增強(qiáng)GNSS PPP模糊度浮點(diǎn)解與固定解性能評(píng)估
    BDS中長(zhǎng)基線三頻RTK算法研究
    BDS參考站間低高度角衛(wèi)星三頻整周模糊度解算方法
    歷元間載波相位差分的GPS/BDS精密單點(diǎn)測(cè)速算法
    基于浮點(diǎn)DSP的鐵路FSK信號(hào)檢測(cè)
    基于雙差的波動(dòng)方程反射波旅行時(shí)反演方法
    Recent advances of TCM treatment of childhood atopic dermatitis
    Clinical observation of Huatan Huoxue Formula in treating coronary heart disease with hyperlipidemia
    Mechanism of sex hormone level in biological clock disorder induced acne and analysis of TCM Pathogenesis
    日本猛色少妇xxxxx猛交久久| 少妇的丰满在线观看| 精品国产乱码久久久久久小说| 男女边吃奶边做爰视频| 精品国产国语对白av| 欧美 日韩 精品 国产| 少妇被粗大的猛进出69影院 | 亚洲熟女精品中文字幕| 校园人妻丝袜中文字幕| 如何舔出高潮| 精品酒店卫生间| 香蕉国产在线看| 精品午夜福利在线看| 一级毛片我不卡| 免费看av在线观看网站| kizo精华| 精品国产一区二区三区久久久樱花| 精品熟女少妇av免费看| av片东京热男人的天堂| 欧美少妇被猛烈插入视频| 成人无遮挡网站| 高清毛片免费看| 黄色一级大片看看| 草草在线视频免费看| 成人影院久久| 最近2019中文字幕mv第一页| 日韩在线高清观看一区二区三区| 精品福利永久在线观看| 视频区图区小说| 日日摸夜夜添夜夜爱| 欧美成人精品欧美一级黄| 高清视频免费观看一区二区| 观看av在线不卡| 国产高清国产精品国产三级| 国产在视频线精品| 亚洲人与动物交配视频| 成年人午夜在线观看视频| 激情视频va一区二区三区| 亚洲图色成人| 亚洲三级黄色毛片| 日本91视频免费播放| 国产精品秋霞免费鲁丝片| 国产黄色免费在线视频| 一区二区三区精品91| 亚洲久久久国产精品| 欧美激情极品国产一区二区三区 | 毛片一级片免费看久久久久| 亚洲av免费高清在线观看| 亚洲五月色婷婷综合| 内地一区二区视频在线| 久久狼人影院| 久久女婷五月综合色啪小说| 亚洲av中文av极速乱| 搡女人真爽免费视频火全软件| 欧美精品国产亚洲| 国产精品免费大片| 国产亚洲最大av| 久久韩国三级中文字幕| 亚洲人与动物交配视频| 熟女人妻精品中文字幕| 欧美亚洲日本最大视频资源| 成人综合一区亚洲| 在线天堂中文资源库| 免费看不卡的av| 日本-黄色视频高清免费观看| 日韩人妻精品一区2区三区| 日本vs欧美在线观看视频| 少妇的逼好多水| 天堂8中文在线网| av在线播放精品| 一级毛片我不卡| 国产精品国产三级国产专区5o| 在线观看人妻少妇| 日韩在线高清观看一区二区三区| 五月开心婷婷网| 美女大奶头黄色视频| 777米奇影视久久| 亚洲国产精品一区三区| www.熟女人妻精品国产 | 男人舔女人的私密视频| 国产精品99久久99久久久不卡 | 国产亚洲欧美精品永久| 美女福利国产在线| 精品99又大又爽又粗少妇毛片| 人妻一区二区av| av在线观看视频网站免费| 亚洲欧洲日产国产| 国产片特级美女逼逼视频| 老女人水多毛片| 91精品三级在线观看| 欧美精品人与动牲交sv欧美| 亚洲av福利一区| 咕卡用的链子| 国产乱来视频区| 久久这里有精品视频免费| 国产av码专区亚洲av| 亚洲欧洲国产日韩| 欧美另类一区| 曰老女人黄片| 亚洲丝袜综合中文字幕| 人体艺术视频欧美日本| 99热国产这里只有精品6| 在线观看美女被高潮喷水网站| 午夜福利影视在线免费观看| 免费黄网站久久成人精品| 乱码一卡2卡4卡精品| 国产一区二区激情短视频 | 一二三四在线观看免费中文在 | 欧美老熟妇乱子伦牲交| 亚洲激情五月婷婷啪啪| 制服丝袜香蕉在线| 国产精品 国内视频| av在线观看视频网站免费| 99热网站在线观看| 王馨瑶露胸无遮挡在线观看| 97精品久久久久久久久久精品| 久久久久久久大尺度免费视频| 免费黄频网站在线观看国产| 国产av一区二区精品久久| 免费大片黄手机在线观看| 人妻一区二区av| 两个人免费观看高清视频| 久久久久精品人妻al黑| av.在线天堂| 日本与韩国留学比较| 精品久久久久久电影网| 亚洲欧洲精品一区二区精品久久久 | 男人添女人高潮全过程视频| 久久毛片免费看一区二区三区| 女性生殖器流出的白浆| 久久精品国产亚洲av涩爱| 日本-黄色视频高清免费观看| 免费高清在线观看日韩| 黄色配什么色好看| 精品国产一区二区三区久久久樱花| 成年人免费黄色播放视频| 亚洲伊人色综图| 在线观看一区二区三区激情| 成年女人在线观看亚洲视频| 尾随美女入室| 国产免费福利视频在线观看| 亚洲精品日韩在线中文字幕| 精品久久久久久电影网| 成人影院久久| 一本久久精品| 丝袜在线中文字幕| 国产av一区二区精品久久| 国产免费现黄频在线看| 2022亚洲国产成人精品| 精品午夜福利在线看| 91精品伊人久久大香线蕉| 咕卡用的链子| 九九爱精品视频在线观看| 欧美变态另类bdsm刘玥| av.在线天堂| 亚洲美女视频黄频| 视频中文字幕在线观看| 中文乱码字字幕精品一区二区三区| 亚洲少妇的诱惑av| 成年人免费黄色播放视频| 交换朋友夫妻互换小说| 亚洲成国产人片在线观看| 亚洲精品国产av成人精品| 丝瓜视频免费看黄片| 伦理电影大哥的女人| 国产极品粉嫩免费观看在线| 日本爱情动作片www.在线观看| 婷婷成人精品国产| 国产一区二区三区综合在线观看 | 99久久精品国产国产毛片| 18禁观看日本| 日本黄色日本黄色录像| 两性夫妻黄色片 | 人妻少妇偷人精品九色| 日本欧美视频一区| 男女下面插进去视频免费观看 | 国产精品麻豆人妻色哟哟久久| 天堂中文最新版在线下载| √禁漫天堂资源中文www| 久久精品夜色国产| 久久婷婷青草| 精品福利永久在线观看| 中国国产av一级| 亚洲精品成人av观看孕妇| 美女福利国产在线| 大香蕉97超碰在线| 天天影视国产精品| 午夜福利网站1000一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 老司机影院成人| 久久久久久久亚洲中文字幕| 大陆偷拍与自拍| 十分钟在线观看高清视频www| 国产精品嫩草影院av在线观看| 精品亚洲成a人片在线观看| 国产av精品麻豆| 十分钟在线观看高清视频www| 高清欧美精品videossex| 涩涩av久久男人的天堂| 美女脱内裤让男人舔精品视频| 国产精品国产三级国产av玫瑰| 最近手机中文字幕大全| 日本欧美视频一区| 成年美女黄网站色视频大全免费| 如日韩欧美国产精品一区二区三区| av片东京热男人的天堂| 日产精品乱码卡一卡2卡三| av又黄又爽大尺度在线免费看| 91精品伊人久久大香线蕉| 黄色 视频免费看| 国产xxxxx性猛交| 日韩人妻精品一区2区三区| 一级片'在线观看视频| 老司机亚洲免费影院| 久久久久久久国产电影| 少妇被粗大的猛进出69影院 | 精品国产一区二区三区久久久樱花| 搡女人真爽免费视频火全软件| 黑人欧美特级aaaaaa片| 亚洲高清免费不卡视频| 秋霞在线观看毛片| videossex国产| 欧美3d第一页| 久久狼人影院| 久久久久久久久久久久大奶| 国产成人精品婷婷| 国产亚洲精品久久久com| 精品亚洲成国产av| 国产成人精品福利久久| av播播在线观看一区| 亚洲国产精品一区三区| 亚洲精品第二区| 另类亚洲欧美激情| 母亲3免费完整高清在线观看 | 国产精品国产三级专区第一集| 免费av中文字幕在线| 欧美xxxx性猛交bbbb| 国产色婷婷99| 一级片免费观看大全| 又粗又硬又长又爽又黄的视频| 我要看黄色一级片免费的| 女的被弄到高潮叫床怎么办| 有码 亚洲区| av天堂久久9| 国产成人aa在线观看| 一个人免费看片子| 赤兔流量卡办理| www.av在线官网国产| 色网站视频免费| 亚洲第一区二区三区不卡| av片东京热男人的天堂| 久久韩国三级中文字幕| 亚洲图色成人| 99国产综合亚洲精品| 爱豆传媒免费全集在线观看| 成人无遮挡网站| 国产一区有黄有色的免费视频| 女性生殖器流出的白浆| 国产精品偷伦视频观看了| 欧美少妇被猛烈插入视频| 蜜桃国产av成人99| 欧美xxⅹ黑人| 大片电影免费在线观看免费| 我要看黄色一级片免费的| 亚洲精品国产av成人精品| 国产精品一区二区在线不卡| 亚洲国产色片| av卡一久久| 啦啦啦啦在线视频资源| 色5月婷婷丁香| 日本与韩国留学比较| 久久99热这里只频精品6学生| 少妇 在线观看| 人妻 亚洲 视频| 黄色配什么色好看| 婷婷色麻豆天堂久久| 欧美人与性动交α欧美精品济南到 | 卡戴珊不雅视频在线播放| 麻豆乱淫一区二区| 久久久欧美国产精品| 如日韩欧美国产精品一区二区三区| 欧美精品av麻豆av| 久久人人爽人人爽人人片va| 亚洲内射少妇av| 九九爱精品视频在线观看| 亚洲人成网站在线观看播放| 丝袜人妻中文字幕| 中文字幕av电影在线播放| av在线app专区| a级毛片在线看网站| 亚洲 欧美一区二区三区| 久久久久久久亚洲中文字幕| 免费观看在线日韩| √禁漫天堂资源中文www| 亚洲五月色婷婷综合| 亚洲,欧美精品.| av播播在线观看一区| 一区二区三区四区激情视频| 欧美成人午夜免费资源| 黑丝袜美女国产一区| 日本vs欧美在线观看视频| 国产精品久久久久成人av| 国产淫语在线视频| 亚洲美女搞黄在线观看| 精品国产露脸久久av麻豆| 免费不卡的大黄色大毛片视频在线观看| 精品国产国语对白av| 伊人久久国产一区二区| 国产黄频视频在线观看| 中文字幕精品免费在线观看视频 | 国产精品一区二区在线不卡| 国产又色又爽无遮挡免| 一级毛片电影观看| 精品人妻在线不人妻| 亚洲一级一片aⅴ在线观看| 国产欧美日韩一区二区三区在线| 18在线观看网站| 丝袜美足系列| 亚洲av福利一区| 热re99久久国产66热| 国产深夜福利视频在线观看| 国产精品一区二区在线观看99| 男女啪啪激烈高潮av片| 国产精品一区二区在线不卡| 在线免费观看不下载黄p国产| 热99久久久久精品小说推荐| 久久久亚洲精品成人影院| 国产成人av激情在线播放| 狠狠精品人妻久久久久久综合| 成人亚洲精品一区在线观看| 日日撸夜夜添| 亚洲国产色片| 亚洲图色成人| 大香蕉久久网| 色吧在线观看| 搡女人真爽免费视频火全软件| 五月玫瑰六月丁香| 亚洲五月色婷婷综合| 亚洲成av片中文字幕在线观看 | 亚洲精品中文字幕在线视频| 人妻一区二区av| 只有这里有精品99| 久久精品国产综合久久久 | 在线观看三级黄色| 性色av一级| 亚洲经典国产精华液单| 日韩av免费高清视频| 欧美亚洲日本最大视频资源| 99视频精品全部免费 在线| 成人18禁高潮啪啪吃奶动态图| 亚洲国产成人一精品久久久| 高清黄色对白视频在线免费看| av一本久久久久| 男女午夜视频在线观看 | 五月天丁香电影| 久久毛片免费看一区二区三区| 成人国产麻豆网| 久久久久精品性色| 99国产精品免费福利视频| 久久久久精品久久久久真实原创| 亚洲国产色片| 色网站视频免费| 亚洲高清免费不卡视频| 看十八女毛片水多多多| 最后的刺客免费高清国语| 90打野战视频偷拍视频| 精品久久久精品久久久| 全区人妻精品视频| 久久久久网色| 91午夜精品亚洲一区二区三区| 9热在线视频观看99| 大话2 男鬼变身卡| 欧美精品亚洲一区二区| 亚洲av电影在线进入| 一级毛片 在线播放| 国产 一区精品| 99久久综合免费| h视频一区二区三区| 免费高清在线观看视频在线观看| 国产精品久久久久久av不卡| 极品少妇高潮喷水抽搐| 午夜日本视频在线| 好男人视频免费观看在线| 亚洲精品久久成人aⅴ小说| 99香蕉大伊视频| 97人妻天天添夜夜摸| 国产一区二区三区综合在线观看 | 日韩一区二区三区影片| 女性生殖器流出的白浆| 黑人猛操日本美女一级片| 天天躁夜夜躁狠狠躁躁| 天天躁夜夜躁狠狠久久av| 一区二区三区四区激情视频| 久久久久久久亚洲中文字幕| 国产 一区精品| 免费日韩欧美在线观看| 九九在线视频观看精品| 亚洲四区av| 亚洲av日韩在线播放| 国产 精品1| 日韩av不卡免费在线播放| 秋霞在线观看毛片| 亚洲熟女精品中文字幕| 免费黄网站久久成人精品| 精品久久久久久电影网| 巨乳人妻的诱惑在线观看| 免费久久久久久久精品成人欧美视频 | 黄色配什么色好看| 午夜91福利影院| 免费大片黄手机在线观看| 国产精品欧美亚洲77777| 国产日韩欧美在线精品| 免费在线观看黄色视频的| 久热久热在线精品观看| 在线观看免费日韩欧美大片| 欧美日韩av久久| 免费看光身美女| 国产精品欧美亚洲77777| 亚洲美女搞黄在线观看| 18在线观看网站| 一本大道久久a久久精品| 在线精品无人区一区二区三| 午夜激情久久久久久久| 亚洲精品国产av蜜桃| 最新的欧美精品一区二区| 欧美国产精品va在线观看不卡| 寂寞人妻少妇视频99o| 亚洲国产色片| 中文欧美无线码| 三级国产精品片| 日本av免费视频播放| 亚洲精品国产av蜜桃| 美女国产视频在线观看| 国产xxxxx性猛交| www.熟女人妻精品国产 | 国产av国产精品国产| 夫妻性生交免费视频一级片| 又黄又粗又硬又大视频| 1024视频免费在线观看| 美女视频免费永久观看网站| 亚洲伊人久久精品综合| 日韩精品免费视频一区二区三区 | 免费黄频网站在线观看国产| 国产精品秋霞免费鲁丝片| 久久久久久人人人人人| 男女啪啪激烈高潮av片| 国精品久久久久久国模美| 乱码一卡2卡4卡精品| 国产精品 国内视频| 成人无遮挡网站| 欧美变态另类bdsm刘玥| 一本—道久久a久久精品蜜桃钙片| 97人妻天天添夜夜摸| 欧美日韩综合久久久久久| 免费观看性生交大片5| 秋霞伦理黄片| av又黄又爽大尺度在线免费看| 精品99又大又爽又粗少妇毛片| 2022亚洲国产成人精品| 少妇人妻精品综合一区二区| 美女主播在线视频| 亚洲少妇的诱惑av| 午夜久久久在线观看| 欧美+日韩+精品| 在线观看免费高清a一片| 国产有黄有色有爽视频| 少妇 在线观看| 欧美人与性动交α欧美精品济南到 | av天堂久久9| 国产精品人妻久久久久久| 欧美少妇被猛烈插入视频| 久久狼人影院| 色吧在线观看| 免费看av在线观看网站| 人成视频在线观看免费观看| 日日摸夜夜添夜夜爱| 日韩成人伦理影院| 国产男人的电影天堂91| 亚洲av在线观看美女高潮| 国产精品成人在线| 男人爽女人下面视频在线观看| 老司机影院毛片| 日日爽夜夜爽网站| 久久鲁丝午夜福利片| 视频中文字幕在线观看| 人妻 亚洲 视频| 精品久久国产蜜桃| 欧美精品高潮呻吟av久久| 亚洲综合色惰| 高清在线视频一区二区三区| 少妇的丰满在线观看| 国产av一区二区精品久久| 综合色丁香网| 美女脱内裤让男人舔精品视频| 下体分泌物呈黄色| 亚洲天堂av无毛| videos熟女内射| 国产视频首页在线观看| 亚洲av电影在线进入| 九色亚洲精品在线播放| 又大又黄又爽视频免费| 午夜福利乱码中文字幕| 一二三四中文在线观看免费高清| 久久人妻熟女aⅴ| 青春草国产在线视频| 少妇人妻久久综合中文| 日韩欧美一区视频在线观看| 国产免费视频播放在线视频| 久久人妻熟女aⅴ| 性色av一级| 十分钟在线观看高清视频www| 国产av一区二区精品久久| 久久国产亚洲av麻豆专区| 永久免费av网站大全| 欧美亚洲日本最大视频资源| 18禁裸乳无遮挡动漫免费视频| 国产成人午夜福利电影在线观看| 国产成人a∨麻豆精品| 母亲3免费完整高清在线观看 | 高清在线视频一区二区三区| 久久久精品94久久精品| 色94色欧美一区二区| 日韩 亚洲 欧美在线| 成人国产麻豆网| 久热这里只有精品99| 啦啦啦中文免费视频观看日本| 日本黄大片高清| 亚洲av综合色区一区| 各种免费的搞黄视频| 午夜福利,免费看| 一二三四中文在线观看免费高清| 久久久久久久久久久久大奶| 高清毛片免费看| 九色成人免费人妻av| 国产欧美亚洲国产| 欧美日韩av久久| 国产国语露脸激情在线看| 一区二区日韩欧美中文字幕 | av视频免费观看在线观看| 日日爽夜夜爽网站| 日韩人妻精品一区2区三区| 欧美成人午夜免费资源| 国产视频首页在线观看| 99视频精品全部免费 在线| 国产精品一二三区在线看| 成人国语在线视频| 色5月婷婷丁香| 亚洲精品国产av成人精品| 街头女战士在线观看网站| 国产亚洲精品第一综合不卡 | 午夜福利乱码中文字幕| 亚洲欧洲精品一区二区精品久久久 | 一级,二级,三级黄色视频| 日韩一区二区三区影片| 中文字幕人妻熟女乱码| 高清在线视频一区二区三区| 99香蕉大伊视频| 性色av一级| 男女边摸边吃奶| 国产精品秋霞免费鲁丝片| 自拍欧美九色日韩亚洲蝌蚪91| 欧美 日韩 精品 国产| 1024视频免费在线观看| 好男人视频免费观看在线| av视频免费观看在线观看| 亚洲av男天堂| 亚洲av综合色区一区| 日韩人妻精品一区2区三区| 久久精品国产综合久久久 | 男女边摸边吃奶| 夜夜骑夜夜射夜夜干| 波野结衣二区三区在线| av国产精品久久久久影院| 在线观看免费视频网站a站| 亚洲欧美成人综合另类久久久| 秋霞在线观看毛片| 色网站视频免费| 精品少妇久久久久久888优播| 青春草视频在线免费观看| 26uuu在线亚洲综合色| a级毛色黄片| 久久久亚洲精品成人影院| 一本大道久久a久久精品| 国产白丝娇喘喷水9色精品| 青春草亚洲视频在线观看| 亚洲三级黄色毛片| 99久久中文字幕三级久久日本| 男女午夜视频在线观看 | 精品视频人人做人人爽| 免费大片18禁| 免费女性裸体啪啪无遮挡网站| 婷婷色av中文字幕| 午夜久久久在线观看| 国产精品一区二区在线观看99| 天堂俺去俺来也www色官网| 少妇的丰满在线观看| 免费大片黄手机在线观看| 51国产日韩欧美| 菩萨蛮人人尽说江南好唐韦庄| 日韩成人伦理影院| 成年人午夜在线观看视频| 久久久久精品久久久久真实原创| 蜜桃在线观看..| 王馨瑶露胸无遮挡在线观看| 国产精品无大码| 在线观看免费视频网站a站| 大片免费播放器 马上看| 91国产中文字幕| 亚洲丝袜综合中文字幕| 性色avwww在线观看| 99香蕉大伊视频| 亚洲精华国产精华液的使用体验| 一本大道久久a久久精品| 国产极品粉嫩免费观看在线| 日韩成人伦理影院|