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

    兩種偽距定位精度分析及計(jì)算程序的實(shí)現(xiàn)

    2021-04-22 02:53:36楊久東龔櫟澎
    導(dǎo)航定位學(xué)報(bào) 2021年2期
    關(guān)鍵詞:星間坐標(biāo)值歷元

    李 韌,楊久東,龔櫟澎,梁 鵬

    兩種偽距定位精度分析及計(jì)算程序的實(shí)現(xiàn)

    李 韌,楊久東,龔櫟澎,梁 鵬

    (華北理工大學(xué) 礦業(yè)工程學(xué)院,河北 唐山 063210)

    針對卡爾曼濾波偽距定位方法難以建立準(zhǔn)確的定位模型,容易造成濾波發(fā)散,導(dǎo)致定位精度降低的問題,比較分析了卡爾曼濾波偽距定位和星頻雙差偽距定位兩種偽距定位模型,利用三個(gè)點(diǎn)的靜態(tài)觀測數(shù)據(jù),采用全球定位系統(tǒng)(GPS)的衛(wèi)星偽距觀測值和廣播星歷文件計(jì)算得到測站點(diǎn)坐標(biāo),將兩種偽距定位模型計(jì)算的坐標(biāo)與赫戈(HGO)軟件靜態(tài)解算的坐標(biāo)進(jìn)行精度對比與分析。結(jié)果顯示,卡爾曼濾波偽距定位方法的誤差在6 m內(nèi),星頻雙差偽距定位方法的誤差在3 m內(nèi),對總體均方根誤差(RMSE)進(jìn)行分析,星頻雙差偽距定位結(jié)果低于卡爾曼濾波偽距定位,證明星頻雙差偽距定位精度高于卡爾曼濾波偽距定位精度,星頻雙差偽距定位數(shù)學(xué)模型是可行的。

    偽距定位;卡爾曼濾波偽距定位;星間單差;星頻雙差;編程實(shí)現(xiàn)

    0 引言

    隨著科技的發(fā)展,衛(wèi)星導(dǎo)航定位為社會(huì)做出了重大貢獻(xiàn),偽距定位是其中不可或缺的一部分。它利用接收機(jī)的偽距觀測數(shù)據(jù)及廣播星歷文件計(jì)算測站坐標(biāo),是導(dǎo)航定位中最常用的定位方法,主要用于導(dǎo)航定位中初始概略位置的定位[1]。隨著對精密單點(diǎn)定位(precise point positioning, PPP)的研究深入,作為其基礎(chǔ)方法及定位理論的偽距單點(diǎn)定位方法及理論也是極其重要的[2-3]。偽距定位已經(jīng)不僅僅局限于傳統(tǒng)的偽距定位,通過不斷深入研究,已經(jīng)衍生出多種偽距定位解算模型[4-5]。研究發(fā)現(xiàn),采用最小二乘迭代法對偽距進(jìn)行解算,其迭代次數(shù)少但定位精度較低[6-7];采用卡爾曼(Kalman)濾波法解算,可使精度有所提高,但迭代次數(shù)過多[8];很多學(xué)者研究使用最小二乘法和Kalman濾波法組合使用進(jìn)行偽距解算,以減少迭代次數(shù),提高定位精度[9]。偽距單點(diǎn)定位方法的精度一般為米級(jí),組合偽距觀測值定位精度可高于單一偽距觀測值精度[10],全球衛(wèi)星導(dǎo)航系統(tǒng)(global navigation satellite system, GNSS)靜態(tài)精密單點(diǎn)定位模型,點(diǎn)位精度可達(dá)毫米級(jí)[11-13]。

    偽距單點(diǎn)定位方法及理論的研究雖然很多,但關(guān)于偽距定位的數(shù)學(xué)模型的表述卻不夠詳細(xì),筆者通過對偽距定位理論的進(jìn)一步研究,整理出Kalman濾波偽距定位和星間頻間偽距定位模型,以派森(Python)3.8為工具進(jìn)行編程,并利用實(shí)測數(shù)據(jù)驗(yàn)證該數(shù)學(xué)模型的可行性、正確性及程序的穩(wěn)定性,同時(shí)將兩種偽距定位方法的定位結(jié)果進(jìn)行比較,分析兩種偽距定位方法的定位效果及精度,證明星頻雙差偽距定位精度高于卡爾曼濾波偽距定位精度,星頻雙差偽距定位數(shù)學(xué)模型是可行的。

    1 偽距定位模型

    1.1 Kalman濾波偽距定位模型

    Kalman濾波方法是一種利用線性系統(tǒng)狀態(tài)方程,通過系統(tǒng)輸入輸出觀測數(shù)據(jù),對系統(tǒng)狀態(tài)進(jìn)行最優(yōu)估計(jì)的算法,在GNSS定位中應(yīng)用較為廣泛,適用于高精度衛(wèi)星導(dǎo)航定位領(lǐng)域。首先建立Kalman濾波模型,即

    濾波器的遞推公式為

    通過Kalman濾波計(jì)算出測站的位置分量后,偽距定位坐標(biāo)計(jì)算公式為

    由于單獨(dú)每個(gè)歷元的數(shù)據(jù)都可解算一個(gè)測站點(diǎn)坐標(biāo),各歷元數(shù)據(jù)可視作是獨(dú)立觀測值,可通過均方根誤差(root mean square error, RMSE)指標(biāo),判別各歷元解算坐標(biāo)與平均值的變化幅度,坐標(biāo)各分量的RMSE值和坐標(biāo)的RMSE值的計(jì)算方法為

    1.2 星間頻間偽距定位模型

    根據(jù)偽距定位原理,可將偽距公式進(jìn)行變形整理后得

    參考衛(wèi)星的觀測方程與其余衛(wèi)星的觀測方程進(jìn)行星間作差,可得

    GPS衛(wèi)星信號(hào)播發(fā)頻率L1是由調(diào)制基本頻率154×10.23 MHz獲得的頻率,L2是由調(diào)制基本頻率120×10.23 MHz獲得的頻率,所以有

    因?yàn)樾情g頻間偽距定位方法是依據(jù)L1、L2頻率間倍率關(guān)系作差進(jìn)行計(jì)算定位的,但衛(wèi)星的頻率發(fā)射不穩(wěn)定,存在一定的頻率波動(dòng),在后期解算中,需要對系數(shù)進(jìn)行一系列調(diào)整,才能利用最小二乘計(jì)算出較為精確的測站的偽距定位坐標(biāo)。最小二乘解算中的權(quán)陣采用同一歷元下衛(wèi)星高度角定權(quán)模型,權(quán)陣為方差陣的逆,其方差計(jì)算方法為

    觀測誤差方程為

    由最小二乘求解得

    由式(17)可以求得坐標(biāo)改正值,將坐標(biāo)改正值代入式(18)可以求得接收機(jī)坐標(biāo)為

    2 定位坐標(biāo)解算及精度分析

    在利用Python語言進(jìn)行兩種偽距定位模型的程序設(shè)計(jì)時(shí),建立了不同的模塊來實(shí)現(xiàn)偽距定位坐標(biāo)的計(jì)算,主要包括文件讀取、坐標(biāo)轉(zhuǎn)換、時(shí)間統(tǒng)一、衛(wèi)星位置計(jì)算、偽距定位模型計(jì)算、精度評(píng)價(jià)等6大模塊,定位程序主要處理流程如圖1所示。

    圖1 程序設(shè)計(jì)流程圖

    在計(jì)算偽距定位時(shí),需要注意坐標(biāo)系轉(zhuǎn)換。由于GPS定位求得的各種參數(shù)為1984世界大地坐標(biāo)系(world geodetic coordinate system 1984, WGS84)下的結(jié)果,在測站坐標(biāo)計(jì)算的過程中,需要根據(jù)實(shí)際測繪工作需要,將WGS84坐標(biāo)轉(zhuǎn)換為站心坐標(biāo)系。

    運(yùn)用本文所述的兩種偽距定位模型,選擇2019年11月河北省唐山市曹妃甸生態(tài)城控制網(wǎng)中C401、C410和GPS1三個(gè)控制點(diǎn)的靜態(tài)觀測數(shù)據(jù)進(jìn)行試驗(yàn),C401與C410的距離為5493 m,C401與GPS1的距離為1950 m,C410與GPS1距離為4348 m。同步觀測時(shí)間為100 min,靜態(tài)接收機(jī)型號(hào)為海星達(dá)IRTK2型雙頻接收機(jī),靜態(tài)數(shù)據(jù)文件為RINEX3.02文件,采樣間隔1 s,高度截止角為10°,因觀測時(shí)間較長,故選取觀測歷元中的300個(gè)歷元參與模型計(jì)算。程序自動(dòng)計(jì)算1~300個(gè)歷元、兩個(gè)偽距定位模型中的定位坐標(biāo),同時(shí)使用中海達(dá)赫戈(HGO)基線解算軟件進(jìn)行控制網(wǎng)解算,經(jīng)嚴(yán)密二維約束平差后,最弱點(diǎn)中誤差的最大值為1.19 mm,符合精度要求。由于相對偽距定位模型的精度為米級(jí),可以將HGO解算的二維約束平差結(jié)果作為真值使用。圖2至圖10為對兩種偽距定位模型進(jìn)行比較的結(jié)果。

    圖2至圖4為測站C401的坐標(biāo)值(、、),圖5至圖7為測站C410的坐標(biāo)值(、、),圖8至圖10為測站GPS1的坐標(biāo)值(、、)。從圖2至圖10可知,三個(gè)測站的Kalman濾波偽距定位模型和星間頻間偽距定位模型的每個(gè)歷元解算的定位坐標(biāo),都出現(xiàn)了不同程度的波動(dòng),符合偽距定位的解算規(guī)律。在兩種定位方法模型解算的300個(gè)歷元坐標(biāo)中:Kalman濾波偽距定位方法在前10個(gè)歷元處波動(dòng)較大,原因在于,在動(dòng)態(tài)濾波過程中,Kalman濾波需要一定的歷元數(shù)進(jìn)行收斂,而收斂后的圖線波動(dòng)較小,圖形表現(xiàn)較為平穩(wěn);星間頻間偽距定位模型解算的300歷元定位坐標(biāo),波動(dòng)較大,大部分圖線波動(dòng)在精密坐標(biāo)圖線的上方與下方,其平均值與精密坐標(biāo)圖線更為接近。

    圖2 測站C401的X坐標(biāo)值

    圖3 測站C401的Y坐標(biāo)值

    圖4 測站C401的Z坐標(biāo)值

    圖5 測站C410的X坐標(biāo)值

    圖6 測站C410的Y坐標(biāo)值

    圖7 測站C410的Z坐標(biāo)值

    圖8 測站GPS1的X坐標(biāo)值

    圖9 測站GPS1的Y坐標(biāo)值

    圖10 測站GPS1的Z坐標(biāo)值

    表1 測站定位坐標(biāo)誤差及坐標(biāo)變化范圍

    為了更直觀地反映觀測值偏離真實(shí)有效數(shù)據(jù)值的波動(dòng)情況,將其結(jié)果列于表2,RMSE值越小,表明誤差就越小,精度越高,偽距定位越準(zhǔn)確。通過對表2進(jìn)行分析可知,用兩種偽距定位方法在三個(gè)測站進(jìn)行300個(gè)歷元的解算,三個(gè)測站的RMSE值表現(xiàn)一致,均是卡爾曼濾波偽距定位的坐標(biāo)解算精度高于星間頻間偽距定位,而星間頻間偽距定位在三個(gè)測站上的RMSE、RMSE、RMSE的值均低于卡爾曼濾波偽距定位,其中測站C401、測站C410和測站GPS1的星間頻間偽距定位的總體均方根誤差RMSE比Kalman濾波偽距定位分別提升了48%、68%和24%。綜合分析,星間頻間偽距定位方法的定位精度要高于Kalman濾波偽距定位,但是坐標(biāo)解算精度低于Kalman濾波偽距定位精度,主要原因可能在于選擇的解算模型不夠精確,在定位計(jì)算中發(fā)生了偏移,導(dǎo)致定位精度不高。

    表2 測站RMSE值

    3 結(jié)束語

    本文對兩種偽距定位的數(shù)學(xué)模型進(jìn)行分析和研究,用Python程序?qū)崿F(xiàn)了Kalman濾波偽距定位模型、星間頻間偽距定位模型的計(jì)算,并利用靜態(tài)實(shí)測數(shù)據(jù)進(jìn)行分析,測站點(diǎn)的偽距定位精度在米級(jí),Kalman濾波偽距定位方法計(jì)算的、、坐標(biāo)誤差在6 m內(nèi),星間頻間偽距定位方法計(jì)算的、、坐標(biāo)誤差在3 m內(nèi),星間頻間偽距定位的精度高于Kalman濾波偽距定位。在經(jīng)過大量實(shí)測數(shù)據(jù)的計(jì)算檢驗(yàn)后,表明筆者編制的偽距定位程序可以滿足要求。由于實(shí)驗(yàn)數(shù)據(jù)區(qū)域較小,屬于短基線解算,對流層及電離層誤差影響較小,而且僅僅采用GPS單系統(tǒng)的偽距觀測值進(jìn)行了偽距計(jì)算,未加入其它衛(wèi)星系統(tǒng)的數(shù)據(jù)進(jìn)行聯(lián)合解算。若在后期能加入其它衛(wèi)星系統(tǒng)的數(shù)據(jù)進(jìn)行有效組合偽距定位,則能進(jìn)一步提高偽距定位精度及穩(wěn)定性。

    [1]衛(wèi)柳艷. GPS網(wǎng)絡(luò)RTK系統(tǒng)的算法及定位精度研究[D]. 南京: 河海大學(xué), 2006: 2-6.

    [2]CAI Changsheng, GAO Yang. Modeling and assessment of combined GPS/GLONASS precise point positioning[J]. GPS Solutions. 2013, 17(3): 223-236.

    [3]CAI C S, KUANG C L, DAI W J. Performance analysis of combined GPS/GLONASS precise point positioning[J]. Science of Surveying & Mapping, 2012, 37(3): 183-185.

    [4]原洪峰, 楊龍, 蔣德. 不同實(shí)驗(yàn)條件下GPS偽距定位精度分析[J]. 全球定位系統(tǒng), 2005, 17(5): 23-26.

    [5]DING He, SUN Fuping, LI Yaping, et al. Modeling and performance analysis of combined BDS/GPS/GLONASS precise point positioning[J]. Journal of Geodesy & Geodynamics, 2016, 36(4): 303-307.

    [6]王祖蔭. 卡爾曼濾波與經(jīng)典最小二乘法[J]. 巖礦測試, 1993, 12(1): 23-26.

    [7]LI P, ZHANG X H. Modeling and performance analysis of GPS/GLONASS/BDS precise point positioning[C]//China Satellite Navigation Conference (CSNC). Proceedings of Lecture Notes in Electrical Engineerin: Volume III. Berlin, Heidelberg: Springer, 2014: 251-26. DOI: https: //doi. org/10.1007/978-3-642-54740-9_23[2020-05-18].

    [8]王小冬. 卡爾曼濾波在衛(wèi)星導(dǎo)航中的應(yīng)用研究[D]. 大連: 大連海事大學(xué), 2008: 5-9.

    [9] 劉春, 馬穎. 改進(jìn)卡爾曼濾波在北斗偽距定位中的研究[J]. 電子測量與儀器學(xué)報(bào), 2016, 23(5): 779-785.

    [10] 李旺, 喬靈娜. 約束條件下GPS偽距定位精度分析[J]. 測繪與空間地理信息, 2019, 42(4): 76-78

    [11] 王趁香. GPS/BDS精密單點(diǎn)定位算法研究與實(shí)現(xiàn)[D]. 阜新: 遼寧工程技術(shù)大學(xué), 2017: 5-9

    [12] 崔立魯, 杜石, 余寒, 等. BDS/GPS組合偽距單點(diǎn)定位性能分析與評(píng)價(jià)[J]. 測繪工程, 2019, 28(2): 17-20, 25.

    [13] 楊徐, 徐愛功, 秦小茜, 等. 高度角定權(quán)模型的BDS/GPS偽距單點(diǎn)定位分析[J]. 導(dǎo)航定位學(xué)報(bào), 2017, 5(2): 72-78, 85.

    [14] 李征航, 張小紅. 衛(wèi)星導(dǎo)航定位新技術(shù)及高精度數(shù)據(jù)處理方法[M]. 武漢: 武漢大學(xué)出版社, 2009: 65-68.

    [15] 符華年, 張旭東, 胡玉坤. GPS靜態(tài)精密單點(diǎn)定位算法精度分析[J]. 測繪通報(bào), 2019(2): 95-98. DOI: 10. 13474 /j. cnki. 11-2246. 2019. 0051.

    [16] 郭平, 周適, 段太生, 等. 一種偽距單點(diǎn)定位的數(shù)學(xué)模型研究及程序?qū)崿F(xiàn)[J]. 全球定位系統(tǒng), 2019, 44(4): 53-60.

    Accuracy analysis and program implementation of two kinds of pseudo range positioning

    LI Ren,YANG Jiudong,GONG Lipeng,LIANG Peng

    ( Mining Engineering College,North China University of Technology,Tangshan,Hebei 063210,China)

    The Kalman filter pseudo-range positioning method finds it difficult to establish an exact positioning model and filtering divergence could be easily caused which reduce the accuracy of position. In point of these problems, we have compared and analyzed two pseudo-range positioning model,namely Kalman filter pseudo-range positioning model and star frequency double difference pseudo-distance positioning, taken advantage of static observation data measured by engineering at 3 points, gotten the satellite pseudo-distance measured data of the Global Positioning System (GPS) and measured position coordinates computed by broadcast ephemeris.Comparing and analyzing the accuracy of coordinates computed by the pseudo-distance positioning model and by the static HGO software, we can find that the error range of Kalman filter pseudo-range positioning model is within 6 m and the star frequency double difference pseudo-distance positioning is within 3 m. Analyzing the value of General Root Mean Square Error (RMSE), we can see that the positioning results of star frequency double difference pseudo-distance are all lower than the Kalman filter pseudodistance positioning. It has demonstrated that the accuracy of star frequency double difference pseudo-distance positioning is higher than the Kalman filter pseudo-distance positioning and its feasibility for star frequency double difference pseudo-distance positioning to locate math models.

    pseudo range positioning;Kalman filter pseudo range positioning;single difference between stars;double difference between stars;programming implementation

    P228

    A

    2095-4999(2021)02-0083-07

    李韌,楊久東,龔櫟澎,等. 兩種偽距定位精度分析及計(jì)算程序的實(shí)現(xiàn)[J]. 導(dǎo)航定位學(xué)報(bào), 2021, 9(2): 83-89.(LI Ren, YANG Jiudong, GONG Lipeng, et al. Accuracy analysis and program implementation of two kinds of pseudo range positioning[J]. Journal of Navigation and Positioning, 2021, 9(2): 83-89.)

    10.16547/j.cnki.10-1096.20210213.

    2020-06-17

    國家自然科學(xué)基金項(xiàng)目(51904105)。

    李韌(1996—),男,河北廊坊人,碩士研究生,研究方向?yàn)榭臻g大地測量。

    龔櫟澎(1995—),男,河北唐山人,碩士研究生,研究方向?yàn)榭臻g大地測量。

    猜你喜歡
    星間坐標(biāo)值歷元
    萌寵家族之星間飛行
    麥弗遜懸架主銷軸線對半軸滑移的影響
    北京汽車(2023年1期)2023-03-03 00:50:38
    基于星間鏈路的導(dǎo)航衛(wèi)星時(shí)間自主恢復(fù)策略
    歷元間載波相位差分的GPS/BDS精密單點(diǎn)測速算法
    基于二分法迭代的凸模數(shù)控銑削加工編程*
    星地星間聯(lián)合時(shí)間比對與衛(wèi)星鐘預(yù)報(bào)
    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
    基于星間測距的導(dǎo)航衛(wèi)星自主時(shí)間維持
    国内揄拍国产精品人妻在线 | 狂野欧美激情性xxxx| 国产高清videossex| 不卡av一区二区三区| 免费在线观看影片大全网站| 啦啦啦 在线观看视频| 啪啪无遮挡十八禁网站| 欧美日韩亚洲国产一区二区在线观看| 在线观看舔阴道视频| 亚洲精品一卡2卡三卡4卡5卡| 亚洲欧美日韩无卡精品| 男女之事视频高清在线观看| 中出人妻视频一区二区| 亚洲午夜精品一区,二区,三区| 999久久久国产精品视频| 久久天堂一区二区三区四区| 国产午夜精品久久久久久| 成人国产一区最新在线观看| 精品一区二区三区av网在线观看| xxx96com| 国产熟女xx| 白带黄色成豆腐渣| 亚洲一区中文字幕在线| 欧美不卡视频在线免费观看 | 看黄色毛片网站| 欧美激情久久久久久爽电影| 精品国产超薄肉色丝袜足j| 2021天堂中文幕一二区在线观 | 91麻豆av在线| 精品久久久久久久久久久久久 | 97碰自拍视频| 国产高清激情床上av| 国产成年人精品一区二区| 国产精品一区二区三区四区久久 | 成人亚洲精品av一区二区| 色播在线永久视频| 一级黄色大片毛片| 国产亚洲精品一区二区www| 悠悠久久av| 少妇被粗大的猛进出69影院| 叶爱在线成人免费视频播放| 精品久久久久久久末码| 欧美色欧美亚洲另类二区| 免费在线观看日本一区| 熟女电影av网| 亚洲自拍偷在线| 国产精品久久久久久精品电影 | 人人妻人人看人人澡| 欧美 亚洲 国产 日韩一| 亚洲 欧美 日韩 在线 免费| 一边摸一边做爽爽视频免费| 搡老岳熟女国产| 18禁国产床啪视频网站| 欧美日韩乱码在线| 大香蕉久久成人网| 狂野欧美激情性xxxx| 久久久久久久久免费视频了| 久久精品亚洲精品国产色婷小说| 日韩欧美 国产精品| 丁香欧美五月| 免费高清在线观看日韩| 一级作爱视频免费观看| 午夜福利成人在线免费观看| 国产久久久一区二区三区| 怎么达到女性高潮| 国产aⅴ精品一区二区三区波| 亚洲中文日韩欧美视频| 又紧又爽又黄一区二区| 97碰自拍视频| 国产免费av片在线观看野外av| 国产精品,欧美在线| 岛国视频午夜一区免费看| 叶爱在线成人免费视频播放| 1024视频免费在线观看| 一进一出抽搐动态| 丰满人妻熟妇乱又伦精品不卡| 日日摸夜夜添夜夜添小说| 亚洲最大成人中文| 夜夜看夜夜爽夜夜摸| 亚洲欧美日韩无卡精品| 亚洲熟妇中文字幕五十中出| 国产成人精品久久二区二区91| 久久久国产成人免费| 成人国语在线视频| 黄色成人免费大全| 国产精品九九99| 美女 人体艺术 gogo| 久久久久久久久中文| 2021天堂中文幕一二区在线观 | 欧美精品亚洲一区二区| 久久久水蜜桃国产精品网| 亚洲色图 男人天堂 中文字幕| 高潮久久久久久久久久久不卡| 身体一侧抽搐| 免费观看人在逋| 精品久久蜜臀av无| 欧美激情久久久久久爽电影| 99精品欧美一区二区三区四区| 午夜精品久久久久久毛片777| 啦啦啦免费观看视频1| 成人午夜高清在线视频 | 亚洲激情在线av| 大型av网站在线播放| 国产成人啪精品午夜网站| 国产视频一区二区在线看| 亚洲,欧美精品.| 国产av又大| 可以在线观看毛片的网站| 午夜福利在线在线| 久久亚洲真实| 两个人看的免费小视频| 国产三级在线视频| 一级片免费观看大全| 国产高清videossex| 国产99白浆流出| 欧美激情极品国产一区二区三区| 国产成人欧美在线观看| 午夜两性在线视频| 国产精品爽爽va在线观看网站 | 精品电影一区二区在线| 亚洲av熟女| 国产久久久一区二区三区| 亚洲 欧美 日韩 在线 免费| 免费观看人在逋| www.精华液| 夜夜看夜夜爽夜夜摸| 变态另类丝袜制服| 久久久国产成人免费| 在线观看舔阴道视频| 午夜福利欧美成人| 别揉我奶头~嗯~啊~动态视频| 国产亚洲精品av在线| 午夜福利在线在线| 日韩大码丰满熟妇| 精品国产国语对白av| 国产一区在线观看成人免费| 亚洲色图 男人天堂 中文字幕| 成人精品一区二区免费| 婷婷六月久久综合丁香| 99在线视频只有这里精品首页| 亚洲国产欧美网| 精华霜和精华液先用哪个| 1024视频免费在线观看| 韩国av一区二区三区四区| www日本黄色视频网| 最近最新免费中文字幕在线| 麻豆一二三区av精品| 村上凉子中文字幕在线| 一区二区三区精品91| 日日爽夜夜爽网站| 一区福利在线观看| 1024手机看黄色片| www日本在线高清视频| 一进一出好大好爽视频| 99久久精品国产亚洲精品| 757午夜福利合集在线观看| 在线观看一区二区三区| 无人区码免费观看不卡| 一级a爱视频在线免费观看| 国产精品一区二区三区四区久久 | 国内少妇人妻偷人精品xxx网站 | 免费观看人在逋| 丰满人妻熟妇乱又伦精品不卡| 亚洲国产欧美日韩在线播放| 99国产精品一区二区蜜桃av| 精品国产乱子伦一区二区三区| 欧美又色又爽又黄视频| www.www免费av| 国产精品野战在线观看| 一进一出好大好爽视频| 一进一出抽搐动态| 亚洲人成网站在线播放欧美日韩| 俄罗斯特黄特色一大片| 夜夜躁狠狠躁天天躁| 午夜福利在线在线| 两性夫妻黄色片| netflix在线观看网站| av在线播放免费不卡| 国产精华一区二区三区| 成人精品一区二区免费| 亚洲精品国产区一区二| 欧美 亚洲 国产 日韩一| 国产精品野战在线观看| 日韩中文字幕欧美一区二区| 黄片小视频在线播放| 久久性视频一级片| 亚洲欧美日韩高清在线视频| 黄网站色视频无遮挡免费观看| 黄片小视频在线播放| 午夜福利在线观看吧| 精品无人区乱码1区二区| 亚洲 欧美一区二区三区| 亚洲精品国产区一区二| 亚洲久久久国产精品| 久久久久九九精品影院| 高清在线国产一区| 欧美精品亚洲一区二区| av在线播放免费不卡| 亚洲一码二码三码区别大吗| 夜夜夜夜夜久久久久| 黑人操中国人逼视频| 欧美又色又爽又黄视频| 亚洲无线在线观看| 午夜激情福利司机影院| 人人妻人人看人人澡| 国产男靠女视频免费网站| 淫秽高清视频在线观看| xxx96com| 日韩欧美一区视频在线观看| 国产真实乱freesex| 欧美性长视频在线观看| 亚洲无线在线观看| 久久久久精品国产欧美久久久| 97碰自拍视频| 免费看十八禁软件| 男男h啪啪无遮挡| 午夜免费成人在线视频| 国产精品一区二区精品视频观看| 国产乱人伦免费视频| 国产精品香港三级国产av潘金莲| 美女高潮到喷水免费观看| 一本久久中文字幕| 免费观看精品视频网站| 亚洲一区二区三区色噜噜| 国产精品久久久久久精品电影 | 看免费av毛片| 亚洲中文字幕一区二区三区有码在线看 | 久久天堂一区二区三区四区| 免费看十八禁软件| 精品第一国产精品| 精品一区二区三区视频在线观看免费| 精品乱码久久久久久99久播| 在线观看免费日韩欧美大片| 国产成人av激情在线播放| 青草久久国产| 久久午夜综合久久蜜桃| 国产成人欧美| 中亚洲国语对白在线视频| 精品国产国语对白av| 精品乱码久久久久久99久播| 日本在线视频免费播放| 波多野结衣av一区二区av| 精品欧美国产一区二区三| 美女国产高潮福利片在线看| 亚洲九九香蕉| 麻豆一二三区av精品| 日韩欧美一区二区三区在线观看| 成熟少妇高潮喷水视频| 国产野战对白在线观看| 久久午夜亚洲精品久久| 欧美成人午夜精品| 夜夜爽天天搞| 日韩精品中文字幕看吧| 欧美性猛交黑人性爽| 国产aⅴ精品一区二区三区波| 国产精品久久久av美女十八| 久久这里只有精品19| 国产亚洲av高清不卡| 香蕉丝袜av| 怎么达到女性高潮| 欧美黄色片欧美黄色片| 日韩欧美国产在线观看| 一区福利在线观看| 精品熟女少妇八av免费久了| 亚洲国产日韩欧美精品在线观看 | 黄片小视频在线播放| 香蕉国产在线看| 欧美最黄视频在线播放免费| av中文乱码字幕在线| 青草久久国产| 国产亚洲精品久久久久久毛片| 亚洲av电影在线进入| 成人三级做爰电影| 亚洲欧美精品综合一区二区三区| 亚洲激情在线av| 一区二区三区精品91| 热re99久久国产66热| 免费av毛片视频| 婷婷丁香在线五月| 亚洲美女黄片视频| 桃色一区二区三区在线观看| 久久久久亚洲av毛片大全| 大香蕉久久成人网| 欧美成人午夜精品| 久久精品成人免费网站| 老熟妇乱子伦视频在线观看| 精品免费久久久久久久清纯| 国产亚洲精品av在线| 亚洲精品美女久久av网站| 午夜久久久久精精品| 亚洲国产看品久久| 老司机午夜福利在线观看视频| 韩国精品一区二区三区| 三级毛片av免费| 精品国产亚洲在线| 国产色视频综合| 一二三四在线观看免费中文在| 免费看美女性在线毛片视频| 亚洲午夜理论影院| 成在线人永久免费视频| 看黄色毛片网站| 哪里可以看免费的av片| 国产成人av激情在线播放| 国产免费av片在线观看野外av| 欧美一级毛片孕妇| 精品久久久久久,| 99热6这里只有精品| 亚洲精品中文字幕在线视频| 两性夫妻黄色片| 女同久久另类99精品国产91| 搡老妇女老女人老熟妇| 欧美久久黑人一区二区| 欧美一区二区精品小视频在线| 90打野战视频偷拍视频| 日本三级黄在线观看| 嫩草影视91久久| av有码第一页| 国产国语露脸激情在线看| 看免费av毛片| 久久久国产欧美日韩av| 成年免费大片在线观看| 午夜影院日韩av| 精品一区二区三区四区五区乱码| 韩国精品一区二区三区| 很黄的视频免费| 国产精品综合久久久久久久免费| 久久香蕉激情| 精品不卡国产一区二区三区| 美女高潮到喷水免费观看| 成熟少妇高潮喷水视频| 久久婷婷成人综合色麻豆| 久久精品91无色码中文字幕| 少妇 在线观看| 久久精品国产亚洲av高清一级| 亚洲专区字幕在线| 国产黄色小视频在线观看| 国产亚洲精品综合一区在线观看 | 1024手机看黄色片| 欧美一区二区精品小视频在线| 后天国语完整版免费观看| 亚洲电影在线观看av| 18禁裸乳无遮挡免费网站照片 | 18禁美女被吸乳视频| 黑人操中国人逼视频| 90打野战视频偷拍视频| 精品国产国语对白av| 久久久国产精品麻豆| 成人一区二区视频在线观看| 久久久久久久午夜电影| 热re99久久国产66热| 琪琪午夜伦伦电影理论片6080| 精品国产乱子伦一区二区三区| 亚洲精品色激情综合| xxx96com| 亚洲熟妇熟女久久| 久久 成人 亚洲| 久久香蕉国产精品| 亚洲av片天天在线观看| 香蕉丝袜av| 欧美乱色亚洲激情| 99国产精品一区二区蜜桃av| 亚洲精品一卡2卡三卡4卡5卡| 男女视频在线观看网站免费 | 国产蜜桃级精品一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 国产欧美日韩一区二区三| 色尼玛亚洲综合影院| 欧美在线黄色| 丰满人妻熟妇乱又伦精品不卡| 国产av一区二区精品久久| 国产精品 欧美亚洲| 97人妻精品一区二区三区麻豆 | 天天躁夜夜躁狠狠躁躁| 中文字幕精品亚洲无线码一区 | 91国产中文字幕| 国产午夜福利久久久久久| 桃红色精品国产亚洲av| 精品国产亚洲在线| 欧美日韩福利视频一区二区| www.自偷自拍.com| xxxwww97欧美| 国产成人系列免费观看| 国产成人啪精品午夜网站| 两个人免费观看高清视频| 国产成人啪精品午夜网站| 久久婷婷人人爽人人干人人爱| 精品欧美一区二区三区在线| 伊人久久大香线蕉亚洲五| 夜夜爽天天搞| 啦啦啦韩国在线观看视频| 亚洲国产精品999在线| 99在线视频只有这里精品首页| 欧美激情久久久久久爽电影| 视频区欧美日本亚洲| 正在播放国产对白刺激| 午夜免费成人在线视频| 午夜精品在线福利| 免费无遮挡裸体视频| 午夜福利免费观看在线| 亚洲天堂国产精品一区在线| 精品久久久久久,| 午夜精品久久久久久毛片777| 黄色成人免费大全| 国产一区二区三区视频了| 男女之事视频高清在线观看| 99国产精品99久久久久| 国产精品免费一区二区三区在线| 欧美成人免费av一区二区三区| 好男人电影高清在线观看| 一本久久中文字幕| 一级a爱视频在线免费观看| 18禁裸乳无遮挡免费网站照片 | 大型av网站在线播放| 黄色成人免费大全| 真人一进一出gif抽搐免费| 麻豆久久精品国产亚洲av| 午夜精品久久久久久毛片777| 欧美+亚洲+日韩+国产| 桃色一区二区三区在线观看| avwww免费| 正在播放国产对白刺激| 天堂√8在线中文| 热re99久久国产66热| 久久人妻福利社区极品人妻图片| 色老头精品视频在线观看| 亚洲精华国产精华精| 在线观看免费视频日本深夜| 国产熟女午夜一区二区三区| 最好的美女福利视频网| www.精华液| 欧美激情 高清一区二区三区| 淫妇啪啪啪对白视频| 成人免费观看视频高清| 不卡av一区二区三区| 日本 欧美在线| 动漫黄色视频在线观看| 可以在线观看的亚洲视频| 国产男靠女视频免费网站| 特大巨黑吊av在线直播 | 真人做人爱边吃奶动态| 精品国产乱码久久久久久男人| 91九色精品人成在线观看| 国产成人精品久久二区二区91| 精品不卡国产一区二区三区| 精品熟女少妇八av免费久了| 欧美av亚洲av综合av国产av| 日本在线视频免费播放| 亚洲熟妇熟女久久| 中文字幕精品免费在线观看视频| 麻豆av在线久日| 欧美亚洲日本最大视频资源| 别揉我奶头~嗯~啊~动态视频| 久久久久免费精品人妻一区二区 | 啪啪无遮挡十八禁网站| 中文字幕精品免费在线观看视频| 天天添夜夜摸| 一级毛片高清免费大全| 老汉色av国产亚洲站长工具| 91老司机精品| 亚洲av美国av| 人人澡人人妻人| 99热只有精品国产| 国产人伦9x9x在线观看| 国产三级在线视频| 男女床上黄色一级片免费看| 国产av一区在线观看免费| 怎么达到女性高潮| 国产一区二区在线av高清观看| 成人精品一区二区免费| 亚洲无线在线观看| 中文字幕高清在线视频| 国产精品日韩av在线免费观看| 热re99久久国产66热| 99国产综合亚洲精品| 两性午夜刺激爽爽歪歪视频在线观看 | 午夜视频精品福利| 青草久久国产| 久久精品91蜜桃| 嫁个100分男人电影在线观看| 变态另类成人亚洲欧美熟女| 大型av网站在线播放| 欧美不卡视频在线免费观看 | 在线看三级毛片| 亚洲国产毛片av蜜桃av| 精品久久蜜臀av无| 一区福利在线观看| 高清在线国产一区| 久久久久久久久久黄片| 淫妇啪啪啪对白视频| 在线天堂中文资源库| 国产精品电影一区二区三区| 视频区欧美日本亚洲| 欧美乱码精品一区二区三区| 免费看美女性在线毛片视频| 国产精品乱码一区二三区的特点| 欧美成人性av电影在线观看| 欧美精品亚洲一区二区| a在线观看视频网站| 老汉色∧v一级毛片| 在线观看日韩欧美| 久久久久久九九精品二区国产 | 国产亚洲欧美精品永久| 一级a爱片免费观看的视频| 欧美一级毛片孕妇| 精华霜和精华液先用哪个| 特大巨黑吊av在线直播 | 人妻久久中文字幕网| 国产精品香港三级国产av潘金莲| 黄色丝袜av网址大全| 88av欧美| 50天的宝宝边吃奶边哭怎么回事| 国产欧美日韩一区二区三| 观看免费一级毛片| 久久精品国产亚洲av香蕉五月| 精品欧美一区二区三区在线| 午夜免费成人在线视频| 国产成人av教育| 日韩欧美国产一区二区入口| aaaaa片日本免费| 欧美在线黄色| 精品久久久久久成人av| av欧美777| 欧洲精品卡2卡3卡4卡5卡区| 久久亚洲精品不卡| 成人欧美大片| 一个人观看的视频www高清免费观看 | 在线观看日韩欧美| 午夜a级毛片| 丁香六月欧美| 极品教师在线免费播放| 国产精品亚洲美女久久久| 在线观看www视频免费| 成在线人永久免费视频| 成人永久免费在线观看视频| 性色av乱码一区二区三区2| 一本大道久久a久久精品| 亚洲七黄色美女视频| 777久久人妻少妇嫩草av网站| 亚洲国产欧美一区二区综合| 欧美绝顶高潮抽搐喷水| 免费高清在线观看日韩| 不卡av一区二区三区| 人成视频在线观看免费观看| 女性被躁到高潮视频| 欧美成人性av电影在线观看| 亚洲一区中文字幕在线| 国产成人精品久久二区二区91| 国产精品国产高清国产av| 久9热在线精品视频| 中文字幕人成人乱码亚洲影| 久久青草综合色| 亚洲欧美一区二区三区黑人| 12—13女人毛片做爰片一| 美女扒开内裤让男人捅视频| 女生性感内裤真人,穿戴方法视频| 欧美不卡视频在线免费观看 | www.精华液| 露出奶头的视频| 午夜两性在线视频| 少妇的丰满在线观看| 久久精品影院6| 白带黄色成豆腐渣| 日本一区二区免费在线视频| 在线永久观看黄色视频| 精品国产国语对白av| 午夜福利视频1000在线观看| 亚洲第一av免费看| 免费一级毛片在线播放高清视频| 久久精品91无色码中文字幕| 男女午夜视频在线观看| 我的亚洲天堂| 久久久久久久精品吃奶| 亚洲av成人不卡在线观看播放网| 最新美女视频免费是黄的| 欧美日韩亚洲综合一区二区三区_| 国产99白浆流出| 校园春色视频在线观看| 琪琪午夜伦伦电影理论片6080| 在线十欧美十亚洲十日本专区| 亚洲一卡2卡3卡4卡5卡精品中文| 国产又色又爽无遮挡免费看| 99国产精品99久久久久| 韩国av一区二区三区四区| 俺也久久电影网| 在线视频色国产色| 在线播放国产精品三级| 国产精品二区激情视频| 成人亚洲精品av一区二区| 精品卡一卡二卡四卡免费| 制服丝袜大香蕉在线| 成人国产综合亚洲| 午夜福利在线观看吧| 色哟哟哟哟哟哟| x7x7x7水蜜桃| 美女大奶头视频| e午夜精品久久久久久久| 国产真人三级小视频在线观看| 成人18禁在线播放| 99精品欧美一区二区三区四区| 黄色a级毛片大全视频| 精品久久久久久久毛片微露脸| www日本在线高清视频| 桃红色精品国产亚洲av| 日韩精品青青久久久久久| 首页视频小说图片口味搜索| 国产欧美日韩一区二区三| 丁香六月欧美| 亚洲精品中文字幕一二三四区| 国产亚洲精品久久久久久毛片| 人人妻人人澡人人看| 精品国产一区二区三区四区第35| 亚洲欧洲精品一区二区精品久久久| 欧美激情极品国产一区二区三区| 日本熟妇午夜| 国产亚洲av嫩草精品影院|