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

    結(jié)合IGS分析中心產(chǎn)品的軌道綜合算法及其精度分析

    2018-06-04 03:06:02陸軼材高成發(fā)
    測繪通報 2018年5期
    關(guān)鍵詞:框架軌道衛(wèi)星

    陸軼材,高成發(fā),郭 奇

    (東南大學(xué)交通學(xué)院,江蘇 南京 210096)

    IGS軌道和鐘差產(chǎn)品是精密定位等重要的基礎(chǔ)數(shù)據(jù),也是維持衛(wèi)星導(dǎo)航系統(tǒng)時空基準(zhǔn)的重要產(chǎn)品。目前IGS運行的數(shù)據(jù)處理分析中心不少于12家,主要包括歐洲定軌中心(CODE)、加拿大自然資源部(EMR)、歐洲空間局(ESA)、德國波茲坦地學(xué)研究中心(GFZ)、法國空間研究中心(GRG)、美國噴氣動力實驗室(JPL)、美國麻省理工學(xué)院(MIT)、美國國家大地測量局(NGS)、美國斯克里普斯海洋學(xué)研究所(SIO)、美國海軍天文臺(USN)、捷克大地實驗室(GOP)和武漢大學(xué)(WHU)均能獨立解算,生成精度優(yōu)于5 cm的GNSS精密軌道[1-3]。由于各家分析中心使用的基準(zhǔn)站坐標(biāo)、數(shù)據(jù)處理軟件及定軌策略等不盡相同,定軌結(jié)果會存在明顯的差異[4]。IGS對多個分析中心產(chǎn)品進(jìn)行質(zhì)量分析,然后通過綜合解算生成參考框架和時間基準(zhǔn)穩(wěn)定、統(tǒng)一且1 d精度優(yōu)于2.5 cm的綜合產(chǎn)品[5-7]。

    隨著IGS分析中心的發(fā)展,IGS產(chǎn)品綜合的方法不斷得到改進(jìn)和完善。文獻(xiàn)[4]表明分析中心的軌道坐標(biāo)存在系統(tǒng)偏差,且該偏差與極移參數(shù)之差相關(guān);文獻(xiàn)[8]提出了一種利用L1G范數(shù)法求取轉(zhuǎn)換參數(shù)進(jìn)行軌道產(chǎn)品加權(quán)平均的綜合算法。文獻(xiàn)[9—11]探討和分析了IGS產(chǎn)品的一致性問題,提出利用旋轉(zhuǎn)變換和附加改正等方法維持IGS最終軌道/鐘差和周解SINEX等產(chǎn)品的一致性。文獻(xiàn)[12]總結(jié)了IGS軌道產(chǎn)品綜合的3種常用的加權(quán)平均定權(quán)方法:等權(quán)法、基于Helmert轉(zhuǎn)換的驗后方差定權(quán)法及方差因子定權(quán)法。文獻(xiàn)[13]通過選擇最優(yōu)的分析中心組合和框架差異消除方法進(jìn)行精密軌道綜合,極大提升了合成軌道的精度?,F(xiàn)有文獻(xiàn)中常用的定權(quán)方法均是以各分析中心等權(quán)平均軌道為參考軌道,進(jìn)行參數(shù)估計,很少基于分析中心框架、軌道RMS趨于穩(wěn)定的方向?qū)壍谰C合的研究。IGS最終軌道產(chǎn)品在每周四發(fā)布且具有12~18 d的時延性,相比之下分析中心的軌道產(chǎn)品則提前2~3 d,即IGS的最終綜合產(chǎn)品通常需要在分析中心產(chǎn)品全部提交后的2~3 d后才能得到,因此對分析中心的軌道產(chǎn)品進(jìn)行綜合生成精度更優(yōu)的產(chǎn)品是極具工程意義的。故本文提出了基于IGS分析中心產(chǎn)品的軌道綜合算法,并自行編寫程序?qū)PS衛(wèi)星軌道進(jìn)行綜合,為北斗系統(tǒng)中不同類型衛(wèi)星(如GEO/IGSO/MEO)軌道綜合方法提供參考,為iGMAS中心的建設(shè)提供建議[14-16]。

    1 軌道綜合原理及數(shù)學(xué)模型

    目前IGS僅發(fā)布GPS、GLONASS綜合軌道產(chǎn)品,產(chǎn)品綜合方法有兩種:動力學(xué)法和加權(quán)平均法。第一種方法是將各分析中心產(chǎn)品作為偽觀測數(shù)據(jù)進(jìn)行軌道確定,估計衛(wèi)星初始狀態(tài)和其他動力模型參數(shù),保持了軌道動力學(xué)特性,可靠性較高,但增加了算法復(fù)雜度,且軌道和鐘差綜合無法采用同一算法。第二種方法是將分析中心產(chǎn)品首先轉(zhuǎn)換到同一參考框架,然后進(jìn)行加權(quán)平均得到綜合產(chǎn)品,是IGS目前采用的方法。該方法簡單靈活,軌道和鐘差綜合可采用類似算法,軌道綜合產(chǎn)品在特定條件下滿足動力學(xué)特性,但喪失了軌道動力學(xué)的高階項特性。本文將使用加權(quán)平均法進(jìn)行軌道綜合研究。

    雖然分析中心定軌策略與解算軟件不同,導(dǎo)致精度存在明顯差異,但從GPS 1506周后,分析中心軌道精度趨于穩(wěn)定,WRMS曲線趨于平緩(如圖1所示),分析中心的權(quán)值也在整體上趨于穩(wěn)定。

    圖1 分析中心軌道相對于IGS綜合軌道的WRMS

    1.1 IGS軌道產(chǎn)品的綜合

    基于文獻(xiàn)[2,8,17—20],歸納IGS軌道產(chǎn)品綜合的主要步驟如下:①統(tǒng)一分析中心軌道坐標(biāo)的參考框架,通常選取IERS的極移解(Bulletin A/B)或IGS分析中心的ERP均值作為“真值”,分析中心的軌道坐標(biāo)根據(jù)其極移解與“真值”的差值分別繞X軸和Y軸作相應(yīng)旋轉(zhuǎn);②進(jìn)行Helmert轉(zhuǎn)換,進(jìn)一步統(tǒng)一分析中心的參考框架并確定分析中心的初始權(quán);③基于加權(quán)平均得到組合軌道。

    1.1.1 參考框架統(tǒng)一

    由于不同分析中心的軌道坐標(biāo)參考框架之間存在旋轉(zhuǎn)、平移等系統(tǒng)性偏差,在軌道綜合之前需先通過旋轉(zhuǎn)變換消除坐標(biāo)參考框架不一致的影響。本文最終軌道綜合選取IGS發(fā)布的周極移解為“真值”,以當(dāng)天中間時刻ERP參數(shù)為參考,計算分析中心極移解與“真值”的偏差,然后對軌道坐標(biāo)作旋轉(zhuǎn)變換。由于極移偏差均為小角度,則旋轉(zhuǎn)變換的公式可按下式簡化

    (1)

    轉(zhuǎn)之后的分析中心衛(wèi)星坐標(biāo);dXp、dYp分別為各家分析中心極移解相較于IGS發(fā)布的周極移在X、Y方向分量的差值。

    1.1.2 Helmert轉(zhuǎn)換定初始權(quán)

    在進(jìn)行極移旋轉(zhuǎn)變換后,分析中心衛(wèi)星坐標(biāo)系統(tǒng)間仍存在系統(tǒng)偏差。以IGS發(fā)布的最終軌道為“真值”,將分析中心軌道坐標(biāo)通過Helmert轉(zhuǎn)換統(tǒng)一到相同的坐標(biāo)參考框架消除偏差,并依據(jù)最小二乘準(zhǔn)則,用間接平差法求解參數(shù),誤差方程如下

    V3×N,1=B3×N,7·X7,1-L3×N,1

    (2)

    (3)

    基于Helmert旋轉(zhuǎn)后的驗后方差對各分析中心進(jìn)行定權(quán)

    (4)

    (5)

    1.1.3 綜合軌道合成

    (6)

    (7)

    1.2 軌道綜合精度評定

    根據(jù)各分析中心變換處理后的坐標(biāo)與綜合軌道間X、Y、Z3個方向坐標(biāo)分量的關(guān)系,可得各分析中心與參考軌道間的RMS及各顆衛(wèi)星的RMS,公式如下

    (8)

    (9)

    式中,mAC為分析中心AC的衛(wèi)星個數(shù)。

    2 試驗及結(jié)果分析

    按照圖2流程,利用自編的軌道綜合算法,利用2017年2月26日—2017年4月8日(GPS 1938—1940周)的9家分析中心事后精密軌道為試驗數(shù)據(jù),進(jìn)行精密軌道綜合試驗。其中分別以第K周的數(shù)據(jù)作為先驗觀測數(shù)據(jù),進(jìn)行轉(zhuǎn)換參數(shù)的求解與定權(quán)處理,對第K+1周的數(shù)據(jù)進(jìn)行軌道合成。

    2.1 IGS軌道產(chǎn)品的綜合

    定權(quán)結(jié)果見表1—表2,為了易于比較和分析,將分析中心ESA的權(quán)值定為單位權(quán)(1.000)。經(jīng)分析得到,分析中心NGS、CODE、GFZ權(quán)值基本大于1,即相對精度較好,剩余分析中心的相對精度則稍差,其中分析中心NGS在所有分析中心的相對權(quán)值最大,即相對精度最高。該結(jié)果與IGS發(fā)布的Weighted Average圖具有一致性,驗證了此算法的可行性(DOY:年積日)。

    2.2 軌道合成結(jié)果及分析

    對GPS 1939周各個分析中心的軌道產(chǎn)品進(jìn)行軌道綜合,圖3為PRN08號衛(wèi)星7天內(nèi)所有歷元的綜合軌道與IGS最終軌道在X、Y、Z方向的差值,其余衛(wèi)星的效果與圖3一致。圖4為綜合軌道、分析中心與IGS最終軌道的單天三維RMS(GPS周1939)。圖5為分析中心綜合軌道相對于IGS最終軌道各顆衛(wèi)星單天平均三維RMS(GPS周1939),其中PRN 04為非健康狀態(tài)。表3給出了合成軌道的單天平均三維RMS(GPS周1939—1943)。

    圖2 軌道合成程序流程

    ACDOY057058059060061062063064065066067068069070CODE1.0631.0941.5281.5681.4501.4561.7151.8151.5841.1931.4991.5120.9951.231EMR0.3770.7310.6960.6210.6110.4820.7930.6940.5330.3120.5280.5540.4070.529ESA1.0001.0001.0001.0001.0001.0001.0001.0001.0001.0001.0001.0001.0001.000GFZ1.2611.3541.3451.3481.0151.1611.3471.2201.5341.3020.9881.1950.9561.218GRG0.5530.6090.8560.6840.6660.3860.7230.6220.3180.7090.6250.7870.4830.723JPL0.4061.0450.8050.9620.9140.8670.9250.7990.4390.6100.8290.8220.7650.950MIT0.5000.5550.5650.5100.5300.6400.5800.5510.6690.4820.6040.6840.4860.647NGS1.5401.5021.4221.3281.4331.7631.6821.2590.4601.4131.7531.5681.1861.466SIO0.5360.5420.5670.6420.5690.5490.7270.5560.6170.6050.6480.8780.6940.907

    表2 分析中心軌道產(chǎn)品定權(quán)(1940—1941周)

    圖3 綜合軌道與IGS最終軌道的差值

    圖4 綜合軌道、分析中心與IGS最終軌道的單天3D RMS

    分析以上數(shù)據(jù)處理結(jié)果,圖4、圖5表明該算法得到的綜合軌道相比于各分析中心的軌道精度有顯著的提高。從圖4可以發(fā)現(xiàn),分析中心軌道與IGS最終軌道比較的平均三維RMS均大于15 mm,即表明分析中心與IGS發(fā)布的綜合軌道間存在框架差異問題。圖5中表明綜合軌道與IGS最終軌道相比,各GPS衛(wèi)星的軌道精度在3~7 mm之間,絕大部分優(yōu)于5 mm。

    圖5 綜合軌道各顆衛(wèi)星單天平均3D RMS

    表3以周為單位,對綜合軌道的整體精度進(jìn)行統(tǒng)計計算,結(jié)果表明綜合軌道單周精度優(yōu)于4.5 mm,證明了對分析中心軌道產(chǎn)品進(jìn)行加權(quán)綜合的有效性。

    表3 合成軌道精度統(tǒng)計(1939—1944)

    3 結(jié) 論

    本文研究了IGS分析中心軌道產(chǎn)品的綜合方法,提出了基于IGS分析中心產(chǎn)品的軌道綜合算法。利用自編算法對IGS各個分析中心2017年2月26日—2017年4月8日精密軌道進(jìn)行綜合,獲得GPS綜合軌道。同時以IGS最終軌道為參考,比較分析了綜合軌道的精度。結(jié)果表明,IGS進(jìn)行軌道綜合的9個分析中心中,NGS、GFZ、CODE和ESA 4個分析中心的軌道產(chǎn)品精度相對較好,剩余5個分析中心的軌道產(chǎn)品相對較差;分析中心軌道與IGS發(fā)布的綜合軌道間存在框架差異性;使用本文算法計算得到的合成軌道結(jié)果與IGS綜合軌道作比較,二者三維差異小于5 mm,優(yōu)于快速星歷,證明了本文算法的可行性。目前,我國將建成多個北斗數(shù)據(jù)分析中心并提供軌道等各類產(chǎn)品,本文研究將為分析中心軌道產(chǎn)品的綜合提供有益的參考[21]。

    參考文獻(xiàn):

    [1] DOW J M,NEILAN R E,GENDT G.The International GPS Service:Celebrating the 10th Anniversary and Looking to the Next Decade[J].Advances in Space Research,2005,36(3):320-326.

    [2] BEUTLER G,KOUBA J,SPRINGER T A.Combining the Orbits of the IGS Analysis Centers[J].Bulletin Géodésique,1995,69(4):200-222.

    [3] LI X,GE M,DAI X,et al.Accuracy and Reliability of Multi-GNSS Real-time Precise Positioning:GPS,GLONASS,BeiDou,and Galileo[J].Journal of Geodesy,2015,89(6):607-635.

    [4] SPRINGER T A,BEUTLER G.Towards an Official IGS Orbit by Combining the Results of All IGS Processing Centers[C]∥Proceedings of the 1993 IGS Workshop.[S.l.]:IGS,1993:24-26.

    [5] GRIFFITHS J.Status of IGS Core Products [C]∥2013 AGU Fall Meeting.[S.l.]:IGS,2013.

    [6] 李征航,黃勁松.GPS測量與數(shù)據(jù)處理[M].3版.武漢:武漢大學(xué)出版社,2016.

    [7] 高成發(fā),胡伍生.衛(wèi)星導(dǎo)航定位原理與應(yīng)用[M].北京:人民交通出版社,2011.

    [8] KOUBA J,MIREAULT Y,LAHAYE F.IGS Orbit/Clock Combination and Evaluation:Appendix I of the Analysis Coordinator Report,International GPS Service for Geo-dynamics (IGS) 1994 Annual Report[M].[S.l.]:Jet Propulsion Laboratory Publication,1994:70-94.

    [9] SPRINGER T A,ZUMBERGE J F,KOUBA J.The IGS Analysis Products and the Consistency of the Combined Solutions[C]∥Proceedings of the 1998 IGS Analysis Center Workshop.[S.l.]:IGS,1998.

    [10] 姚宜斌.GPS精密定位定軌后處理算法與實現(xiàn)[D].武漢:武漢大學(xué),2004.

    [11] 施闖,鄒蓉,姚宜斌,等.基于SINEX解的數(shù)據(jù)組合及系統(tǒng)誤差分析[J].武漢大學(xué)學(xué)報(信息科學(xué)版),2008,33(6):608-611.

    [12] 陳康慷.IGS分析中心軌道/鐘差產(chǎn)品綜合及網(wǎng)解模式動態(tài)精密單點定位研究[D].西安:長安大學(xué),2014.

    [13] 耿濤,徐夏炎.IGS 分析中心軌道綜合算法實現(xiàn)及精度分析[J].大地測量與地球動力學(xué),2017,37(4):369-373.

    [14] STEIGENBERGER P,HUGENTOBLER U,HAUSCHILD A,et al.Orbit and Clock Analysis of Compass GEO and IGSO Satellites[J].Journal of Geodesy,2013,87(6):515-525.

    [15] 譚暢,陳國,張強,等.iGMAS軌道產(chǎn)品綜合及精度初步分析[J].武漢大學(xué)學(xué)報(信息科學(xué)版),2016,41(11):1469-1475.

    [16] CHEN Kangkang,XU Tianhe,Chen G,et al.The Orbit and Clock Combination of iGMAS Analysis Centers and the Analysis of Their Precision[C]∥中國衛(wèi)星導(dǎo)航學(xué)術(shù)年會.西安:[s.n.],2015.

    [17] GRIFFITHS J.A Review of the IGS Final Orbit Combination Procedures (2012),Report from the Analysis Coordinator[EB/OL].[2017-08-07].http:∥acc.igs.org.

    [18] WEBER R,F(xiàn)RANGER E.Combination GLONASS Orbits [C]∥IGS International GLONASS Experiment IGEX-98 Workshop.Tennessee:IGS,1999.

    [19] 劉偉平,郝金明,于合理,等.導(dǎo)航衛(wèi)星精密軌道與鐘差確定方法研究及精度分析[J].測繪通報,2014(5):5-7.

    [20] FERLAND R,PIRASZEWSKI M.The IGS-combined Station Coordinates,Earth Rotation Parameters and Apparent Geo-center[J].Journal of Geodesy,2009,83(3):385-392.

    [21] 于素梅.iGMAS產(chǎn)品精度評定和一致性分析[D].西安:長安大學(xué),2016.

    猜你喜歡
    框架軌道衛(wèi)星
    框架
    miniSAR遙感衛(wèi)星
    廣義框架的不相交性
    基于單純形法的TLE軌道確定
    CryoSat提升軌道高度與ICESat-2同步運行
    朝美重回“相互羞辱軌道”?
    靜止衛(wèi)星派
    科學(xué)家(2019年3期)2019-08-18 09:47:43
    WTO框架下
    法大研究生(2017年1期)2017-04-10 08:55:06
    Puma" suede shoes with a focus on the Product variables
    一種基于OpenStack的云應(yīng)用開發(fā)框架
    日韩一区二区视频免费看| 爱豆传媒免费全集在线观看| 亚洲精品一二三| 国产一级毛片在线| a级毛片黄视频| 秋霞在线观看毛片| 久久久久久免费高清国产稀缺| 国产乱来视频区| 久久久久久久国产电影| 老女人水多毛片| av在线老鸭窝| 精品一区二区三卡| 日韩欧美一区视频在线观看| 国产精品二区激情视频| 精品亚洲成国产av| 捣出白浆h1v1| 国产熟女午夜一区二区三区| 卡戴珊不雅视频在线播放| 国产精品久久久久久久久免| 亚洲精品美女久久av网站| 免费人妻精品一区二区三区视频| 我的亚洲天堂| 亚洲精品国产av蜜桃| 亚洲精品美女久久久久99蜜臀 | 97在线人人人人妻| 亚洲成色77777| 一区在线观看完整版| 9色porny在线观看| 超碰97精品在线观看| 在线观看一区二区三区激情| 天堂中文最新版在线下载| 亚洲av男天堂| 女性被躁到高潮视频| 欧美av亚洲av综合av国产av | www.av在线官网国产| 老鸭窝网址在线观看| 女的被弄到高潮叫床怎么办| 少妇 在线观看| av视频免费观看在线观看| 亚洲在久久综合| 90打野战视频偷拍视频| a级毛片在线看网站| 国产伦理片在线播放av一区| 曰老女人黄片| 一区福利在线观看| 老司机影院成人| 亚洲av免费高清在线观看| 最近的中文字幕免费完整| 欧美日韩亚洲高清精品| 久久久久久久亚洲中文字幕| 亚洲精品一二三| 亚洲经典国产精华液单| 中文天堂在线官网| 69精品国产乱码久久久| 精品午夜福利在线看| 国产成人av激情在线播放| 国产精品欧美亚洲77777| 中国三级夫妇交换| 亚洲欧美清纯卡通| 久久久精品94久久精品| 黑人欧美特级aaaaaa片| 亚洲综合精品二区| 欧美成人午夜精品| 国产日韩欧美亚洲二区| 18禁国产床啪视频网站| 久久久国产精品麻豆| 成人毛片a级毛片在线播放| 99香蕉大伊视频| 又粗又硬又长又爽又黄的视频| 精品少妇久久久久久888优播| 不卡av一区二区三区| 一区二区三区激情视频| 欧美激情 高清一区二区三区| 亚洲精品久久久久久婷婷小说| 欧美精品av麻豆av| 午夜精品国产一区二区电影| av.在线天堂| 精品久久蜜臀av无| 国产精品国产三级国产专区5o| 欧美老熟妇乱子伦牲交| 少妇被粗大猛烈的视频| 午夜91福利影院| 丁香六月天网| 天天躁夜夜躁狠狠久久av| 欧美国产精品一级二级三级| 热re99久久精品国产66热6| 精品少妇黑人巨大在线播放| 一区在线观看完整版| 中文字幕av电影在线播放| 三上悠亚av全集在线观看| 九草在线视频观看| 久久久久国产一级毛片高清牌| 在线 av 中文字幕| 一个人免费看片子| videosex国产| 成人国产av品久久久| 亚洲一级一片aⅴ在线观看| 欧美日韩国产mv在线观看视频| 国产亚洲精品第一综合不卡| 高清黄色对白视频在线免费看| 少妇的丰满在线观看| 精品久久久久久电影网| 9热在线视频观看99| 天美传媒精品一区二区| 欧美人与性动交α欧美软件| 久久久久久免费高清国产稀缺| 久久久a久久爽久久v久久| 尾随美女入室| 久久人妻熟女aⅴ| 啦啦啦中文免费视频观看日本| 亚洲欧美清纯卡通| 中文字幕制服av| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲在久久综合| 人人妻人人澡人人看| 亚洲国产最新在线播放| 少妇的逼水好多| 国产精品成人在线| 亚洲成人手机| 少妇的逼水好多| 制服诱惑二区| 不卡av一区二区三区| 久久精品aⅴ一区二区三区四区 | 中文字幕av电影在线播放| 久久人人97超碰香蕉20202| 久久久国产精品麻豆| 亚洲精品乱久久久久久| 人妻系列 视频| 成年人免费黄色播放视频| 精品一区二区免费观看| 精品少妇一区二区三区视频日本电影 | 久久精品夜色国产| av有码第一页| 一级片'在线观看视频| 免费黄色在线免费观看| 老司机影院毛片| 亚洲综合色网址| 99热全是精品| 叶爱在线成人免费视频播放| 欧美中文综合在线视频| 精品国产超薄肉色丝袜足j| 岛国毛片在线播放| 亚洲内射少妇av| 国产亚洲一区二区精品| 亚洲成色77777| 啦啦啦在线观看免费高清www| 亚洲欧美日韩另类电影网站| 99精国产麻豆久久婷婷| 中国国产av一级| 午夜精品国产一区二区电影| a级毛片在线看网站| 日本爱情动作片www.在线观看| 久久精品国产亚洲av天美| 寂寞人妻少妇视频99o| 亚洲欧洲精品一区二区精品久久久 | 免费在线观看视频国产中文字幕亚洲 | 日本vs欧美在线观看视频| 亚洲一区中文字幕在线| 国产在线免费精品| 欧美国产精品一级二级三级| 亚洲国产精品国产精品| 国产极品天堂在线| 91精品伊人久久大香线蕉| 亚洲欧美色中文字幕在线| 亚洲精品自拍成人| 精品久久久精品久久久| 亚洲精品一二三| 日本wwww免费看| 9热在线视频观看99| 亚洲在久久综合| 十分钟在线观看高清视频www| 一二三四在线观看免费中文在| 成年人免费黄色播放视频| 一本色道久久久久久精品综合| 亚洲欧美精品综合一区二区三区 | 国产精品一区二区在线不卡| av在线老鸭窝| 欧美日韩一级在线毛片| 一二三四在线观看免费中文在| 欧美激情极品国产一区二区三区| 如日韩欧美国产精品一区二区三区| 26uuu在线亚洲综合色| 另类精品久久| 久久精品国产亚洲av天美| 男女无遮挡免费网站观看| 美女脱内裤让男人舔精品视频| 国产极品粉嫩免费观看在线| 人人妻人人添人人爽欧美一区卜| 人妻系列 视频| 免费看不卡的av| 久久免费观看电影| 亚洲内射少妇av| 国产 精品1| 国产淫语在线视频| 午夜福利网站1000一区二区三区| 亚洲av福利一区| 亚洲国产精品成人久久小说| 国产av一区二区精品久久| 日韩av在线免费看完整版不卡| xxxhd国产人妻xxx| 国产av一区二区精品久久| 国产国语露脸激情在线看| 欧美日韩av久久| 大香蕉久久网| 亚洲国产精品一区三区| 国产精品香港三级国产av潘金莲 | 美女高潮到喷水免费观看| 国产男人的电影天堂91| 精品亚洲乱码少妇综合久久| 日本色播在线视频| 在线免费观看不下载黄p国产| 人人妻人人澡人人爽人人夜夜| 一边亲一边摸免费视频| 新久久久久国产一级毛片| 亚洲精品自拍成人| 老司机亚洲免费影院| 不卡视频在线观看欧美| 亚洲精品国产色婷婷电影| 亚洲av在线观看美女高潮| 日韩一本色道免费dvd| 嫩草影院入口| 日本黄色日本黄色录像| www.熟女人妻精品国产| 免费少妇av软件| 99香蕉大伊视频| 午夜激情av网站| 国产精品一区二区在线观看99| 免费久久久久久久精品成人欧美视频| 日本91视频免费播放| 只有这里有精品99| 国产精品久久久久成人av| 国产一区二区三区综合在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 国产成人欧美| 美女脱内裤让男人舔精品视频| 亚洲色图 男人天堂 中文字幕| 自线自在国产av| 成年女人在线观看亚洲视频| 国产熟女欧美一区二区| 亚洲一级一片aⅴ在线观看| 亚洲四区av| 青春草视频在线免费观看| 亚洲一码二码三码区别大吗| 一级片'在线观看视频| 欧美精品亚洲一区二区| 午夜激情av网站| 青春草视频在线免费观看| 麻豆av在线久日| 青青草视频在线视频观看| 国产精品久久久久久久久免| 久久国产精品大桥未久av| 久久这里只有精品19| 日日爽夜夜爽网站| 丁香六月天网| 搡女人真爽免费视频火全软件| 三级国产精品片| 欧美xxⅹ黑人| 日韩免费高清中文字幕av| 欧美日本中文国产一区发布| 亚洲国产精品一区三区| 少妇的逼水好多| 亚洲精品乱久久久久久| av在线app专区| 亚洲色图 男人天堂 中文字幕| 欧美日韩av久久| 国产精品av久久久久免费| 日韩成人av中文字幕在线观看| 久久精品熟女亚洲av麻豆精品| 日韩不卡一区二区三区视频在线| 飞空精品影院首页| 国产精品无大码| 久久久久精品性色| 寂寞人妻少妇视频99o| 91精品伊人久久大香线蕉| 久久av网站| 午夜激情久久久久久久| 亚洲第一av免费看| 亚洲成国产人片在线观看| 女的被弄到高潮叫床怎么办| 久久99精品国语久久久| 侵犯人妻中文字幕一二三四区| 中文字幕人妻丝袜制服| 日韩欧美精品免费久久| 亚洲av电影在线进入| 色视频在线一区二区三区| 久久久精品免费免费高清| 日韩三级伦理在线观看| 丁香六月天网| 色吧在线观看| 新久久久久国产一级毛片| 日本色播在线视频| 三级国产精品片| 母亲3免费完整高清在线观看 | 久久婷婷青草| 一级片'在线观看视频| 新久久久久国产一级毛片| 少妇被粗大猛烈的视频| 久久久久久久大尺度免费视频| 免费人妻精品一区二区三区视频| 一边摸一边做爽爽视频免费| 亚洲激情五月婷婷啪啪| 一二三四在线观看免费中文在| 男女免费视频国产| 欧美国产精品va在线观看不卡| 国产男女超爽视频在线观看| 大片免费播放器 马上看| 亚洲精品美女久久av网站| 国产日韩一区二区三区精品不卡| a级片在线免费高清观看视频| 一边亲一边摸免费视频| 久久婷婷青草| 免费久久久久久久精品成人欧美视频| 久久婷婷青草| 国产亚洲午夜精品一区二区久久| 在线观看免费视频网站a站| 中文字幕色久视频| 日本欧美视频一区| 亚洲国产欧美在线一区| 亚洲欧美色中文字幕在线| 99热网站在线观看| av在线app专区| 亚洲美女视频黄频| 中文字幕制服av| 制服诱惑二区| 女人被躁到高潮嗷嗷叫费观| 欧美日韩一级在线毛片| 999精品在线视频| 汤姆久久久久久久影院中文字幕| a 毛片基地| 国产免费福利视频在线观看| 久久国内精品自在自线图片| 大话2 男鬼变身卡| 曰老女人黄片| 人成视频在线观看免费观看| 国产成人aa在线观看| 亚洲精品aⅴ在线观看| 久久国内精品自在自线图片| 一级毛片黄色毛片免费观看视频| 性色avwww在线观看| 免费高清在线观看日韩| 国产免费一区二区三区四区乱码| 午夜免费男女啪啪视频观看| 欧美日韩成人在线一区二区| 久热这里只有精品99| 最近中文字幕2019免费版| 十分钟在线观看高清视频www| 欧美国产精品va在线观看不卡| 久热这里只有精品99| 日韩中字成人| 免费在线观看视频国产中文字幕亚洲 | 亚洲成人一二三区av| 欧美激情高清一区二区三区 | 国产精品久久久久久av不卡| 在线观看免费视频网站a站| 永久网站在线| a级片在线免费高清观看视频| 日本wwww免费看| 久久人人爽人人片av| 欧美激情高清一区二区三区 | 有码 亚洲区| 欧美 亚洲 国产 日韩一| 欧美97在线视频| 天天操日日干夜夜撸| 午夜激情久久久久久久| 成年人午夜在线观看视频| 亚洲在久久综合| 亚洲精品乱久久久久久| 老汉色∧v一级毛片| 黄频高清免费视频| 人人妻人人添人人爽欧美一区卜| 国产亚洲av片在线观看秒播厂| 久久久亚洲精品成人影院| 天堂8中文在线网| 午夜福利乱码中文字幕| 街头女战士在线观看网站| 高清不卡的av网站| 午夜福利视频在线观看免费| 中文欧美无线码| 久久国产亚洲av麻豆专区| 在线亚洲精品国产二区图片欧美| 亚洲精品视频女| 国产精品久久久久成人av| 18+在线观看网站| 色婷婷av一区二区三区视频| 国产成人aa在线观看| 性高湖久久久久久久久免费观看| 国产视频首页在线观看| 日韩视频在线欧美| 日本av免费视频播放| 午夜福利影视在线免费观看| 午夜免费观看性视频| 最新的欧美精品一区二区| 美国免费a级毛片| 成年人免费黄色播放视频| 国产日韩一区二区三区精品不卡| 久久久a久久爽久久v久久| 少妇人妻 视频| 美女视频免费永久观看网站| 国产成人精品久久久久久| 中文精品一卡2卡3卡4更新| 亚洲av国产av综合av卡| 黄片播放在线免费| 亚洲av免费高清在线观看| 一个人免费看片子| 国产在视频线精品| 久久久欧美国产精品| av片东京热男人的天堂| 亚洲欧美一区二区三区国产| 不卡视频在线观看欧美| 2018国产大陆天天弄谢| 精品亚洲成a人片在线观看| 少妇人妻久久综合中文| 一级黄片播放器| 国产av精品麻豆| 免费观看性生交大片5| 人妻 亚洲 视频| 久久久精品国产亚洲av高清涩受| 久久久久久人人人人人| 老鸭窝网址在线观看| 咕卡用的链子| 国产探花极品一区二区| 欧美精品国产亚洲| 在线亚洲精品国产二区图片欧美| 欧美日韩精品成人综合77777| 国产精品久久久久久av不卡| 高清视频免费观看一区二区| 精品久久蜜臀av无| 亚洲伊人色综图| 午夜日本视频在线| 一级毛片我不卡| 大片免费播放器 马上看| 老司机影院毛片| 男女高潮啪啪啪动态图| 精品人妻一区二区三区麻豆| 日本av免费视频播放| 国产精品99久久99久久久不卡 | 亚洲,欧美,日韩| 97在线视频观看| 亚洲第一av免费看| 亚洲av.av天堂| 久久久久久人人人人人| 亚洲精品久久久久久婷婷小说| 啦啦啦啦在线视频资源| 国产极品粉嫩免费观看在线| 欧美精品高潮呻吟av久久| 亚洲av福利一区| 1024香蕉在线观看| 一二三四中文在线观看免费高清| 香蕉丝袜av| 欧美日韩精品网址| 欧美精品一区二区免费开放| 久久精品国产鲁丝片午夜精品| 看免费av毛片| videosex国产| a级毛片黄视频| 免费人妻精品一区二区三区视频| 九九爱精品视频在线观看| 一本—道久久a久久精品蜜桃钙片| 青春草视频在线免费观看| 久久影院123| av有码第一页| 在线 av 中文字幕| 免费人妻精品一区二区三区视频| 久久精品国产亚洲av天美| 亚洲av福利一区| 宅男免费午夜| 热re99久久国产66热| 天堂中文最新版在线下载| 少妇 在线观看| 中文字幕人妻熟女乱码| 男女午夜视频在线观看| 男女高潮啪啪啪动态图| 精品酒店卫生间| 夜夜骑夜夜射夜夜干| 一级,二级,三级黄色视频| 如何舔出高潮| 日韩中文字幕视频在线看片| 欧美在线黄色| 国产一区亚洲一区在线观看| 女人高潮潮喷娇喘18禁视频| 免费黄频网站在线观看国产| 国产 一区精品| 久久久久久免费高清国产稀缺| 一级,二级,三级黄色视频| 日本欧美国产在线视频| 叶爱在线成人免费视频播放| 久久韩国三级中文字幕| 国产一区二区三区综合在线观看| 亚洲国产精品一区三区| 少妇的逼水好多| 免费久久久久久久精品成人欧美视频| 亚洲图色成人| 最黄视频免费看| 丝袜脚勾引网站| 免费女性裸体啪啪无遮挡网站| 久久 成人 亚洲| 岛国毛片在线播放| 中文字幕av电影在线播放| 国产一区二区激情短视频 | 超碰成人久久| 99国产综合亚洲精品| 国产成人免费无遮挡视频| 中文字幕人妻丝袜制服| 一区二区三区激情视频| 99热网站在线观看| 1024香蕉在线观看| 咕卡用的链子| 亚洲精品日韩在线中文字幕| 久久精品国产a三级三级三级| 一边亲一边摸免费视频| 99国产精品免费福利视频| 最近手机中文字幕大全| 国产精品成人在线| 欧美av亚洲av综合av国产av | 97精品久久久久久久久久精品| 欧美人与性动交α欧美精品济南到 | 久久久精品免费免费高清| av在线播放精品| 国产福利在线免费观看视频| 女人被躁到高潮嗷嗷叫费观| 久久久国产一区二区| 国产一区二区三区综合在线观看| 老司机亚洲免费影院| 亚洲熟女精品中文字幕| 国产免费又黄又爽又色| 久久99热这里只频精品6学生| 中文字幕另类日韩欧美亚洲嫩草| 日韩中文字幕欧美一区二区 | 老鸭窝网址在线观看| 纯流量卡能插随身wifi吗| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲,欧美精品.| 午夜91福利影院| 2022亚洲国产成人精品| 国产亚洲av片在线观看秒播厂| 久久热在线av| 老司机影院毛片| 日韩一区二区三区影片| 国产精品亚洲av一区麻豆 | 欧美日韩视频精品一区| 国产色婷婷99| 日日摸夜夜添夜夜爱| 亚洲av欧美aⅴ国产| 国产爽快片一区二区三区| 国产精品 欧美亚洲| 色94色欧美一区二区| 少妇人妻精品综合一区二区| 亚洲综合色惰| 深夜精品福利| 在线 av 中文字幕| 国产极品粉嫩免费观看在线| 少妇被粗大的猛进出69影院| 精品国产国语对白av| 汤姆久久久久久久影院中文字幕| 在线观看人妻少妇| 校园人妻丝袜中文字幕| 边亲边吃奶的免费视频| 日韩制服丝袜自拍偷拍| 久久ye,这里只有精品| 最近中文字幕2019免费版| 性色av一级| 国产一区亚洲一区在线观看| 在线天堂中文资源库| 欧美日韩av久久| 亚洲国产av影院在线观看| 熟女av电影| 99久国产av精品国产电影| 亚洲天堂av无毛| 日韩电影二区| 精品卡一卡二卡四卡免费| 国产极品天堂在线| 国产精品麻豆人妻色哟哟久久| 免费观看性生交大片5| 美女大奶头黄色视频| 九九爱精品视频在线观看| 搡老乐熟女国产| 黄频高清免费视频| 午夜精品国产一区二区电影| 伦理电影免费视频| 免费观看在线日韩| 三上悠亚av全集在线观看| 久久久精品94久久精品| 香蕉丝袜av| 欧美日韩亚洲高清精品| av一本久久久久| 狠狠婷婷综合久久久久久88av| 美女视频免费永久观看网站| 亚洲欧洲日产国产| 国产av码专区亚洲av| 999精品在线视频| 夫妻午夜视频| 叶爱在线成人免费视频播放| 国产精品久久久久久久久免| 久久久精品94久久精品| av.在线天堂| 狠狠精品人妻久久久久久综合| 国产人伦9x9x在线观看 | 久久久久久免费高清国产稀缺| 一级片'在线观看视频| 久久女婷五月综合色啪小说| 精品国产超薄肉色丝袜足j| 国产男女内射视频| 大陆偷拍与自拍| 欧美亚洲日本最大视频资源| 国产精品欧美亚洲77777| 欧美激情高清一区二区三区 | 欧美最新免费一区二区三区| 啦啦啦啦在线视频资源| 亚洲欧洲日产国产| 在线天堂最新版资源| 久久久久久久国产电影| 成人漫画全彩无遮挡| 天天操日日干夜夜撸| 国产亚洲精品第一综合不卡| 国精品久久久久久国模美|