• <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
    女的被弄到高潮叫床怎么办| 黄频高清免费视频| 丝袜脚勾引网站| 丰满乱子伦码专区| 亚洲国产av影院在线观看| 男女免费视频国产| 另类亚洲欧美激情| 97在线人人人人妻| 韩国高清视频一区二区三区| 日韩欧美精品免费久久| 国产激情久久老熟女| 搡女人真爽免费视频火全软件| 亚洲成人av在线免费| 99re6热这里在线精品视频| 国产又爽黄色视频| 欧美成人午夜免费资源| 少妇人妻精品综合一区二区| 中文字幕av电影在线播放| 免费不卡的大黄色大毛片视频在线观看| 亚洲熟女精品中文字幕| 亚洲欧美一区二区三区黑人 | 久久国产精品大桥未久av| 中文字幕人妻丝袜一区二区 | 亚洲国产成人一精品久久久| 亚洲成人手机| 美女国产视频在线观看| 汤姆久久久久久久影院中文字幕| 欧美最新免费一区二区三区| 亚洲男人天堂网一区| 国产精品香港三级国产av潘金莲 | 久久免费观看电影| 午夜福利网站1000一区二区三区| 黄色配什么色好看| 日韩在线高清观看一区二区三区| 亚洲三级黄色毛片| 女人被躁到高潮嗷嗷叫费观| 青春草国产在线视频| 我要看黄色一级片免费的| 亚洲国产欧美网| 久久精品久久精品一区二区三区| 少妇的丰满在线观看| 视频区图区小说| 日韩av不卡免费在线播放| 国产精品av久久久久免费| 色94色欧美一区二区| 成年动漫av网址| videos熟女内射| 亚洲国产欧美日韩在线播放| 人妻系列 视频| av不卡在线播放| 男男h啪啪无遮挡| 久久久久人妻精品一区果冻| 成人漫画全彩无遮挡| 久久久久视频综合| 久久久精品国产亚洲av高清涩受| 日韩一本色道免费dvd| 老汉色av国产亚洲站长工具| 最近2019中文字幕mv第一页| 精品卡一卡二卡四卡免费| 日韩成人av中文字幕在线观看| 一二三四中文在线观看免费高清| 人人妻人人添人人爽欧美一区卜| 久久久久久久精品精品| 久久午夜综合久久蜜桃| 成人漫画全彩无遮挡| 亚洲情色 制服丝袜| 91久久精品国产一区二区三区| 久久久久人妻精品一区果冻| 男女高潮啪啪啪动态图| 一级片'在线观看视频| 日韩制服骚丝袜av| 大香蕉久久成人网| 香蕉国产在线看| 亚洲精品日韩在线中文字幕| 边亲边吃奶的免费视频| 伊人亚洲综合成人网| 精品福利永久在线观看| 久久人妻熟女aⅴ| 香蕉丝袜av| 99精国产麻豆久久婷婷| 亚洲中文av在线| 亚洲在久久综合| 国产淫语在线视频| 午夜老司机福利剧场| 国产熟女午夜一区二区三区| 大话2 男鬼变身卡| 国产片特级美女逼逼视频| 久久午夜综合久久蜜桃| 久久久精品国产亚洲av高清涩受| 在线看a的网站| 久热这里只有精品99| 人妻少妇偷人精品九色| 波多野结衣一区麻豆| 高清黄色对白视频在线免费看| 欧美人与善性xxx| 男女午夜视频在线观看| 国产欧美日韩一区二区三区在线| 国产1区2区3区精品| 国产一区二区三区av在线| 午夜福利在线免费观看网站| 91成人精品电影| 欧美日韩国产mv在线观看视频| 亚洲成国产人片在线观看| 中文字幕另类日韩欧美亚洲嫩草| 91精品三级在线观看| 欧美国产精品va在线观看不卡| 观看av在线不卡| 久久人妻熟女aⅴ| 热99久久久久精品小说推荐| 天堂俺去俺来也www色官网| 国产成人av激情在线播放| 男女下面插进去视频免费观看| 久久精品国产亚洲av涩爱| 十八禁高潮呻吟视频| 欧美日韩av久久| av在线播放精品| 一级毛片 在线播放| 母亲3免费完整高清在线观看 | 亚洲精品中文字幕在线视频| 啦啦啦视频在线资源免费观看| 色婷婷av一区二区三区视频| 日韩人妻精品一区2区三区| videosex国产| 久久狼人影院| 丝瓜视频免费看黄片| 麻豆精品久久久久久蜜桃| 亚洲综合精品二区| 久久久国产一区二区| 精品国产乱码久久久久久小说| 黄频高清免费视频| 不卡视频在线观看欧美| 国产日韩欧美在线精品| 日韩一区二区三区影片| 韩国av在线不卡| 久久国产亚洲av麻豆专区| 亚洲av成人精品一二三区| 免费在线观看黄色视频的| 亚洲熟女精品中文字幕| 一本—道久久a久久精品蜜桃钙片| 久久久久国产一级毛片高清牌| 777米奇影视久久| 美女xxoo啪啪120秒动态图| 好男人视频免费观看在线| 狂野欧美激情性bbbbbb| 国产爽快片一区二区三区| 日本色播在线视频| 女性生殖器流出的白浆| 成人午夜精彩视频在线观看| 午夜福利乱码中文字幕| 精品少妇黑人巨大在线播放| 大陆偷拍与自拍| 欧美激情高清一区二区三区 | 麻豆av在线久日| 国产成人精品久久二区二区91 | 观看av在线不卡| 精品久久久久久电影网| 一级毛片 在线播放| 国产亚洲欧美精品永久| 欧美精品av麻豆av| 国产av国产精品国产| 久久99热这里只频精品6学生| 午夜福利乱码中文字幕| 欧美精品国产亚洲| 激情视频va一区二区三区| 精品一区二区三卡| 汤姆久久久久久久影院中文字幕| 亚洲国产欧美网| 精品一品国产午夜福利视频| 最黄视频免费看| 午夜91福利影院| 乱人伦中国视频| 欧美+日韩+精品| 美女国产视频在线观看| 一级毛片我不卡| 成人二区视频| 久久久久国产网址| 欧美日韩精品成人综合77777| 大片电影免费在线观看免费| 亚洲精品在线美女| 老汉色av国产亚洲站长工具| 美女主播在线视频| 精品少妇久久久久久888优播| 久久久久久免费高清国产稀缺| 考比视频在线观看| 汤姆久久久久久久影院中文字幕| 色播在线永久视频| 国产精品久久久av美女十八| 亚洲欧美精品自产自拍| 欧美日韩精品成人综合77777| 国产成人精品婷婷| 国产精品久久久久成人av| 欧美精品高潮呻吟av久久| 国产精品av久久久久免费| 亚洲av成人精品一二三区| 男男h啪啪无遮挡| 亚洲欧美日韩另类电影网站| 国产精品一区二区在线观看99| 欧美bdsm另类| 高清不卡的av网站| 日韩制服丝袜自拍偷拍| www.av在线官网国产| 国产熟女午夜一区二区三区| 看免费av毛片| 久久人人爽人人片av| 在线观看人妻少妇| 永久网站在线| 午夜激情久久久久久久| 午夜激情av网站| 黄色配什么色好看| 搡女人真爽免费视频火全软件| 毛片一级片免费看久久久久| 国产亚洲av片在线观看秒播厂| 久久国产精品大桥未久av| 国产一区二区在线观看av| 日日摸夜夜添夜夜爱| 午夜精品国产一区二区电影| 久久鲁丝午夜福利片| 亚洲国产欧美在线一区| 天美传媒精品一区二区| 国产乱来视频区| 欧美激情 高清一区二区三区| 老鸭窝网址在线观看| 久久99热这里只频精品6学生| 欧美在线黄色| 91国产中文字幕| 国产一区二区 视频在线| 超碰97精品在线观看| 日本黄色日本黄色录像| 久久久久精品久久久久真实原创| 一二三四在线观看免费中文在| 免费看av在线观看网站| 久久久国产精品麻豆| 欧美少妇被猛烈插入视频| 亚洲欧美一区二区三区久久| 亚洲欧美色中文字幕在线| 免费人妻精品一区二区三区视频| 日本av手机在线免费观看| 天天躁夜夜躁狠狠躁躁| 国产熟女欧美一区二区| 一本色道久久久久久精品综合| 午夜福利视频在线观看免费| 永久网站在线| 中文字幕最新亚洲高清| 亚洲欧美色中文字幕在线| 午夜福利在线观看免费完整高清在| 综合色丁香网| 欧美97在线视频| 亚洲激情五月婷婷啪啪| 国语对白做爰xxxⅹ性视频网站| 97在线人人人人妻| 一级,二级,三级黄色视频| 人妻人人澡人人爽人人| 久久精品aⅴ一区二区三区四区 | 大码成人一级视频| 黄色视频在线播放观看不卡| 亚洲第一青青草原| 中文字幕色久视频| 久久久久国产精品人妻一区二区| 99re6热这里在线精品视频| 高清av免费在线| 一二三四在线观看免费中文在| 亚洲男人天堂网一区| 亚洲熟女精品中文字幕| 国产精品免费视频内射| 九草在线视频观看| 日韩一卡2卡3卡4卡2021年| 亚洲精品国产av成人精品| 日本欧美视频一区| 亚洲熟女精品中文字幕| 欧美日韩国产mv在线观看视频| 人妻一区二区av| 少妇被粗大猛烈的视频| 免费在线观看视频国产中文字幕亚洲 | 亚洲欧美一区二区三区国产| 久久亚洲国产成人精品v| 超色免费av| 亚洲国产精品一区三区| www.自偷自拍.com| 在线看a的网站| 中文字幕最新亚洲高清| 成年av动漫网址| 99国产综合亚洲精品| 国产精品成人在线| 国产成人免费无遮挡视频| 国产精品久久久av美女十八| 亚洲精品久久久久久婷婷小说| 黄片播放在线免费| 日韩制服丝袜自拍偷拍| 亚洲精品视频女| 亚洲欧洲精品一区二区精品久久久 | 国产人伦9x9x在线观看 | 天天躁夜夜躁狠狠久久av| 免费av中文字幕在线| 一级a爱视频在线免费观看| 久久这里有精品视频免费| 五月伊人婷婷丁香| 咕卡用的链子| 另类亚洲欧美激情| 亚洲av成人精品一二三区| 99久久精品国产国产毛片| 久久狼人影院| 精品国产乱码久久久久久小说| 亚洲欧洲精品一区二区精品久久久 | 男女边摸边吃奶| 日韩成人av中文字幕在线观看| 亚洲三区欧美一区| 久久久久久久久久久免费av| 国产成人免费无遮挡视频| 丝袜喷水一区| 国产成人精品无人区| 在线观看美女被高潮喷水网站| 老熟女久久久| 伊人亚洲综合成人网| 亚洲精品久久午夜乱码| 美女午夜性视频免费| 久久99蜜桃精品久久| 丝袜美足系列| 99久久精品国产国产毛片| 亚洲,欧美精品.| 91aial.com中文字幕在线观看| 国产免费视频播放在线视频| 免费大片黄手机在线观看| 女人精品久久久久毛片| 久久久久国产精品人妻一区二区| 国产成人精品福利久久| 亚洲精品国产一区二区精华液| 秋霞在线观看毛片| 国产精品一区二区在线观看99| 日日爽夜夜爽网站| 精品国产一区二区三区四区第35| 国产精品一国产av| 一区二区三区乱码不卡18| 久热这里只有精品99| 五月伊人婷婷丁香| 亚洲av免费高清在线观看| 欧美+日韩+精品| 99久久中文字幕三级久久日本| 国产黄色视频一区二区在线观看| 五月开心婷婷网| xxx大片免费视频| av国产久精品久网站免费入址| 亚洲三级黄色毛片| 亚洲美女黄色视频免费看| 国产免费一区二区三区四区乱码| 亚洲视频免费观看视频| 精品一区在线观看国产| 精品亚洲成国产av| 边亲边吃奶的免费视频| 一级毛片 在线播放| 男人舔女人的私密视频| 亚洲国产精品一区二区三区在线| 黄色怎么调成土黄色| 欧美av亚洲av综合av国产av | 国产精品一国产av| 免费日韩欧美在线观看| 日本黄色日本黄色录像| 视频区图区小说| 天天躁夜夜躁狠狠躁躁| a级毛片在线看网站| 中文字幕av电影在线播放| 亚洲国产日韩一区二区| 黄色怎么调成土黄色| 久久久久久久亚洲中文字幕| av网站免费在线观看视频| 男人爽女人下面视频在线观看| 亚洲精品美女久久久久99蜜臀 | 男女啪啪激烈高潮av片| 中文字幕av电影在线播放| 一区福利在线观看| 色婷婷av一区二区三区视频| 亚洲精品第二区| 边亲边吃奶的免费视频| 性色avwww在线观看| 熟女av电影| 99久久中文字幕三级久久日本| 黄网站色视频无遮挡免费观看| 国产精品二区激情视频| 国产男人的电影天堂91| 18禁观看日本| 亚洲成人av在线免费| 日本猛色少妇xxxxx猛交久久| 午夜福利乱码中文字幕| 菩萨蛮人人尽说江南好唐韦庄| 少妇 在线观看| 亚洲成人av在线免费| 精品一区在线观看国产| 99国产精品免费福利视频| 久久久久国产一级毛片高清牌| 香蕉精品网在线| 交换朋友夫妻互换小说| 天天躁日日躁夜夜躁夜夜| 人妻 亚洲 视频| 男人操女人黄网站| 亚洲欧美精品综合一区二区三区 | 美女xxoo啪啪120秒动态图| 女人高潮潮喷娇喘18禁视频| 91国产中文字幕| av电影中文网址| 人妻一区二区av| 在线观看三级黄色| 国产成人aa在线观看| 精品人妻偷拍中文字幕| 丝袜美足系列| 91精品三级在线观看| 大片电影免费在线观看免费| 欧美成人午夜免费资源| 成人国语在线视频| 999精品在线视频| 男的添女的下面高潮视频| 精品国产国语对白av| 亚洲男人天堂网一区| 亚洲成国产人片在线观看| 欧美97在线视频| 国产麻豆69| 久久精品久久久久久久性| 男女免费视频国产| 99九九在线精品视频| av视频免费观看在线观看| 91久久精品国产一区二区三区| 亚洲国产精品国产精品| 亚洲精品久久久久久婷婷小说| 热re99久久精品国产66热6| 午夜福利网站1000一区二区三区| 少妇 在线观看| 成年美女黄网站色视频大全免费| 91精品伊人久久大香线蕉| 女性生殖器流出的白浆| 在线观看美女被高潮喷水网站| 九草在线视频观看| 国产日韩一区二区三区精品不卡| 日本欧美视频一区| www.精华液| 国产亚洲av片在线观看秒播厂| 欧美精品国产亚洲| 涩涩av久久男人的天堂| 美女主播在线视频| 久久久国产欧美日韩av| 精品一区二区三区四区五区乱码 | 麻豆av在线久日| 高清欧美精品videossex| 91精品三级在线观看| 精品少妇内射三级| 黄色配什么色好看| 男女国产视频网站| 美女午夜性视频免费| 国产精品秋霞免费鲁丝片| 亚洲国产欧美日韩在线播放| 纵有疾风起免费观看全集完整版| 一本久久精品| 国产白丝娇喘喷水9色精品| 夫妻午夜视频| 老汉色∧v一级毛片| 日韩精品有码人妻一区| 久久久精品区二区三区| 日韩av在线免费看完整版不卡| 日本欧美国产在线视频| 欧美最新免费一区二区三区| 日韩成人av中文字幕在线观看| 国产不卡av网站在线观看| 高清av免费在线| 涩涩av久久男人的天堂| 人人妻人人澡人人看| 另类精品久久| 久久久久久久久久久免费av| 亚洲av免费高清在线观看| 亚洲成色77777| 精品少妇久久久久久888优播| 日韩欧美一区视频在线观看| 啦啦啦中文免费视频观看日本| a级片在线免费高清观看视频| 亚洲综合色网址| 国产又色又爽无遮挡免| 午夜福利视频精品| 亚洲精品aⅴ在线观看| 久久国内精品自在自线图片| 成人午夜精彩视频在线观看| 亚洲精品av麻豆狂野| 久久亚洲国产成人精品v| 欧美中文综合在线视频| 午夜影院在线不卡| 各种免费的搞黄视频| 久久久久久人妻| 久久精品亚洲av国产电影网| 少妇熟女欧美另类| 国产精品 国内视频| 亚洲精品日本国产第一区| 久久久国产精品麻豆| 国产视频首页在线观看| 亚洲成色77777| 国产免费视频播放在线视频| 亚洲精品美女久久av网站| 在线观看美女被高潮喷水网站| 精品一品国产午夜福利视频| av在线播放精品| 91精品三级在线观看| 我的亚洲天堂| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 热99国产精品久久久久久7| 亚洲经典国产精华液单| 免费观看性生交大片5| 欧美激情 高清一区二区三区| 黑人猛操日本美女一级片| 美女午夜性视频免费| 大香蕉久久成人网| 极品少妇高潮喷水抽搐| 人人妻人人澡人人看| 在线亚洲精品国产二区图片欧美| 哪个播放器可以免费观看大片| 久久人人爽人人片av| 国产精品秋霞免费鲁丝片| 久久人人97超碰香蕉20202| 午夜av观看不卡| 久久人人爽人人片av| 这个男人来自地球电影免费观看 | 成人18禁高潮啪啪吃奶动态图| 久久精品aⅴ一区二区三区四区 | 亚洲国产精品成人久久小说| 欧美xxⅹ黑人| 又黄又粗又硬又大视频| 免费观看在线日韩| 精品国产国语对白av| 亚洲精品自拍成人| 最近2019中文字幕mv第一页| 电影成人av| av女优亚洲男人天堂| 午夜福利在线观看免费完整高清在| 日本av手机在线免费观看| 亚洲第一区二区三区不卡| 人妻人人澡人人爽人人| 肉色欧美久久久久久久蜜桃| 亚洲四区av| 波多野结衣一区麻豆| 香蕉精品网在线| 国产精品无大码| 在线观看一区二区三区激情| 国产一区亚洲一区在线观看| 精品人妻一区二区三区麻豆| 中文精品一卡2卡3卡4更新| 侵犯人妻中文字幕一二三四区| 亚洲欧洲国产日韩| 高清在线视频一区二区三区| 国产麻豆69| kizo精华| 久久热在线av| 亚洲精品一二三| 国产精品女同一区二区软件| 伦精品一区二区三区| 人人妻人人添人人爽欧美一区卜| 看免费成人av毛片| 天堂8中文在线网| 亚洲第一区二区三区不卡| www日本在线高清视频| 欧美日韩成人在线一区二区| 韩国高清视频一区二区三区| 日本av免费视频播放| 成年av动漫网址| 男女国产视频网站| 久久亚洲国产成人精品v| av天堂久久9| videos熟女内射| 在线观看美女被高潮喷水网站| 亚洲国产av新网站| 国产亚洲精品第一综合不卡| 我的亚洲天堂| 美女高潮到喷水免费观看| 日日摸夜夜添夜夜爱| 国产成人精品在线电影| 男男h啪啪无遮挡| 18在线观看网站| 国产乱来视频区| 国产免费福利视频在线观看| 成年人午夜在线观看视频| 男的添女的下面高潮视频| 一级a爱视频在线免费观看| 日韩av不卡免费在线播放| 亚洲 欧美一区二区三区| 国产成人精品在线电影| 久久午夜福利片| 国产一区二区激情短视频 | av女优亚洲男人天堂| 26uuu在线亚洲综合色| 啦啦啦啦在线视频资源| 成人黄色视频免费在线看| 婷婷色av中文字幕| 91国产中文字幕| 99热全是精品| 午夜av观看不卡| 亚洲精品在线美女| 中文字幕人妻丝袜制服| 中文字幕制服av| 美女国产视频在线观看| 亚洲第一av免费看| 777米奇影视久久| 亚洲经典国产精华液单| 午夜老司机福利剧场| 久久99一区二区三区| 国产深夜福利视频在线观看| 亚洲精品乱久久久久久| 精品一区二区三卡| 国产视频首页在线观看| 亚洲欧美精品自产自拍| 最近的中文字幕免费完整| 国产黄频视频在线观看| 欧美亚洲日本最大视频资源| 亚洲精品国产一区二区精华液| 国产又色又爽无遮挡免| 亚洲精品久久午夜乱码| 亚洲经典国产精华液单| 亚洲成人一二三区av| 久久精品久久久久久噜噜老黄| videos熟女内射| 国产视频首页在线观看| 不卡视频在线观看欧美| 99九九在线精品视频| av免费观看日本|