• <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ā)框架
    国产精品偷伦视频观看了| 精品久久久久久电影网| 国产精品福利在线免费观看| 又大又黄又爽视频免费| 精品久久久久久电影网| 一边亲一边摸免费视频| 多毛熟女@视频| 免费大片18禁| 国产精品久久久久久av不卡| av不卡在线播放| 最近的中文字幕免费完整| 亚洲av综合色区一区| 亚洲不卡免费看| 青春草视频在线免费观看| 亚洲无线观看免费| 国产成人精品福利久久| 精品国产一区二区三区久久久樱花| 91久久精品电影网| 亚洲天堂av无毛| 国内少妇人妻偷人精品xxx网站| 午夜福利网站1000一区二区三区| 国产免费一级a男人的天堂| 久久99精品国语久久久| 3wmmmm亚洲av在线观看| 亚洲av男天堂| 热re99久久国产66热| 下体分泌物呈黄色| 少妇人妻精品综合一区二区| 韩国av在线不卡| 国产在线一区二区三区精| 久久狼人影院| 97超碰精品成人国产| 黄色欧美视频在线观看| 精品一品国产午夜福利视频| 91成人精品电影| 亚洲av在线观看美女高潮| 51国产日韩欧美| 中文字幕制服av| 免费少妇av软件| kizo精华| 亚洲成色77777| 国产永久视频网站| 日本色播在线视频| 中国三级夫妇交换| 国产亚洲最大av| 亚洲欧美日韩卡通动漫| 一本一本综合久久| 欧美丝袜亚洲另类| av在线app专区| 国产精品人妻久久久久久| 欧美最新免费一区二区三区| 美女大奶头黄色视频| 欧美日韩精品成人综合77777| 97超视频在线观看视频| 午夜免费鲁丝| 在线 av 中文字幕| 国产精品伦人一区二区| 国产一区有黄有色的免费视频| 在线播放无遮挡| www.av在线官网国产| 十八禁网站网址无遮挡 | 日本黄大片高清| 伊人亚洲综合成人网| 中文欧美无线码| 人人妻人人澡人人看| 有码 亚洲区| 亚洲欧美日韩卡通动漫| a 毛片基地| 如何舔出高潮| 国产成人精品一,二区| 亚洲精华国产精华液的使用体验| 色网站视频免费| 久久99精品国语久久久| 日韩大片免费观看网站| 国产日韩一区二区三区精品不卡 | 汤姆久久久久久久影院中文字幕| 日韩不卡一区二区三区视频在线| √禁漫天堂资源中文www| 欧美另类一区| 男女边吃奶边做爰视频| 久久久久久伊人网av| 大话2 男鬼变身卡| 80岁老熟妇乱子伦牲交| 欧美日韩av久久| av天堂久久9| 狂野欧美激情性xxxx在线观看| 美女内射精品一级片tv| 99热全是精品| 久久久久久久亚洲中文字幕| 亚洲国产精品999| 久久国产精品大桥未久av | 精品久久久噜噜| 免费看av在线观看网站| 国产视频内射| av黄色大香蕉| 欧美+日韩+精品| 日本黄大片高清| 最新中文字幕久久久久| 亚洲av电影在线观看一区二区三区| 亚洲va在线va天堂va国产| 成人特级av手机在线观看| 午夜av观看不卡| 日本-黄色视频高清免费观看| 插逼视频在线观看| 69精品国产乱码久久久| 国产亚洲午夜精品一区二区久久| 午夜老司机福利剧场| 欧美丝袜亚洲另类| 日韩欧美一区视频在线观看 | 国产欧美日韩一区二区三区在线 | 男女边摸边吃奶| 久久精品国产自在天天线| √禁漫天堂资源中文www| 国产免费一区二区三区四区乱码| 丰满饥渴人妻一区二区三| 在线观看一区二区三区激情| 国产在线免费精品| 亚洲精品色激情综合| 国产高清国产精品国产三级| 美女内射精品一级片tv| 国产精品国产三级国产专区5o| 国产一区有黄有色的免费视频| 亚洲国产日韩一区二区| 黑丝袜美女国产一区| 国产熟女欧美一区二区| 婷婷色综合大香蕉| 亚洲天堂av无毛| 国产精品欧美亚洲77777| 午夜老司机福利剧场| 亚洲精品中文字幕在线视频 | 亚洲无线观看免费| 交换朋友夫妻互换小说| 一级毛片 在线播放| av国产精品久久久久影院| 一级毛片久久久久久久久女| 亚洲国产色片| 极品少妇高潮喷水抽搐| 制服丝袜香蕉在线| 亚洲美女视频黄频| 亚洲伊人久久精品综合| 大香蕉97超碰在线| 亚洲图色成人| 久久精品国产亚洲av涩爱| 午夜福利在线观看免费完整高清在| 只有这里有精品99| 日韩强制内射视频| 国产精品一区www在线观看| 国产高清三级在线| 久久免费观看电影| av又黄又爽大尺度在线免费看| 在线观看免费高清a一片| 色视频www国产| 免费看光身美女| 自拍偷自拍亚洲精品老妇| 又粗又硬又长又爽又黄的视频| av卡一久久| 免费黄色在线免费观看| 91精品一卡2卡3卡4卡| 男人舔奶头视频| 国产亚洲91精品色在线| 一级毛片我不卡| 国产亚洲一区二区精品| 久久久久久久久久久丰满| 人妻制服诱惑在线中文字幕| 在线观看免费日韩欧美大片 | 男女国产视频网站| 欧美 亚洲 国产 日韩一| 天堂俺去俺来也www色官网| 午夜福利,免费看| a级毛片在线看网站| 亚洲国产欧美在线一区| 亚洲成人手机| 国产视频内射| 成人亚洲欧美一区二区av| 精品久久久久久久久亚洲| 亚洲av电影在线观看一区二区三区| 亚洲四区av| 18+在线观看网站| 久久久午夜欧美精品| 国国产精品蜜臀av免费| 久久久久网色| 久久久国产欧美日韩av| 蜜臀久久99精品久久宅男| 亚洲内射少妇av| 大香蕉97超碰在线| 国产日韩欧美亚洲二区| a级毛片在线看网站| 啦啦啦在线观看免费高清www| 热re99久久国产66热| 亚洲精品第二区| 精品一区二区三卡| 全区人妻精品视频| 美女国产视频在线观看| 午夜激情福利司机影院| 欧美高清成人免费视频www| 人妻夜夜爽99麻豆av| 精品人妻偷拍中文字幕| 综合色丁香网| 欧美3d第一页| 亚洲精品一二三| 成人国产麻豆网| 黄色怎么调成土黄色| 三级国产精品片| 18+在线观看网站| av.在线天堂| 日韩三级伦理在线观看| 一本一本综合久久| 国产乱人偷精品视频| 久久久国产精品麻豆| 精品久久久久久久久亚洲| 国产黄片视频在线免费观看| 特大巨黑吊av在线直播| 免费看光身美女| 欧美日韩精品成人综合77777| 国产高清有码在线观看视频| 亚洲精品中文字幕在线视频 | 久久综合国产亚洲精品| 国产精品久久久久成人av| 日本91视频免费播放| 精品国产露脸久久av麻豆| 日本vs欧美在线观看视频 | 亚洲,欧美,日韩| 久久久午夜欧美精品| 我要看日韩黄色一级片| 国产69精品久久久久777片| 少妇被粗大的猛进出69影院 | 午夜福利影视在线免费观看| 老熟女久久久| 久久女婷五月综合色啪小说| 女人精品久久久久毛片| 一级毛片我不卡| 国产精品麻豆人妻色哟哟久久| 麻豆成人午夜福利视频| 国产欧美日韩精品一区二区| 亚洲精品自拍成人| 搡老乐熟女国产| 亚洲四区av| 青春草视频在线免费观看| 亚洲,欧美,日韩| 亚洲美女黄色视频免费看| 美女xxoo啪啪120秒动态图| 天堂俺去俺来也www色官网| 熟女av电影| 亚洲真实伦在线观看| 亚洲av电影在线观看一区二区三区| 丰满饥渴人妻一区二区三| 伊人久久精品亚洲午夜| 26uuu在线亚洲综合色| 国产免费视频播放在线视频| 国产精品伦人一区二区| 国产探花极品一区二区| 一区在线观看完整版| 久久97久久精品| 麻豆乱淫一区二区| 日韩欧美一区视频在线观看 | 国产极品粉嫩免费观看在线 | 国国产精品蜜臀av免费| 日日爽夜夜爽网站| 熟妇人妻不卡中文字幕| 欧美日韩综合久久久久久| 不卡视频在线观看欧美| 99热国产这里只有精品6| 日日摸夜夜添夜夜爱| 狠狠精品人妻久久久久久综合| 国产高清不卡午夜福利| 国产亚洲91精品色在线| 久久人人爽人人片av| 国产探花极品一区二区| 99久久精品热视频| 午夜老司机福利剧场| 久久人人爽人人片av| 99re6热这里在线精品视频| 91aial.com中文字幕在线观看| 国精品久久久久久国模美| 黑人巨大精品欧美一区二区蜜桃 | 久久精品国产自在天天线| 男男h啪啪无遮挡| 亚洲国产毛片av蜜桃av| 纵有疾风起免费观看全集完整版| 久久久精品94久久精品| 超碰97精品在线观看| 国产在线视频一区二区| 精品国产乱码久久久久久小说| 久久国产亚洲av麻豆专区| 一级爰片在线观看| 亚洲精品日韩在线中文字幕| 一级二级三级毛片免费看| 久热久热在线精品观看| 夫妻性生交免费视频一级片| 日韩av免费高清视频| 六月丁香七月| 日韩视频在线欧美| 精品人妻偷拍中文字幕| 嫩草影院新地址| 黑人巨大精品欧美一区二区蜜桃 | 最近2019中文字幕mv第一页| 久久久久久久国产电影| 色网站视频免费| 欧美日韩视频精品一区| 亚洲欧美日韩卡通动漫| 免费黄频网站在线观看国产| 国产永久视频网站| 亚洲精品久久午夜乱码| 新久久久久国产一级毛片| 久久99一区二区三区| 成人18禁高潮啪啪吃奶动态图 | 亚州av有码| 水蜜桃什么品种好| 成年人免费黄色播放视频 | 国产在视频线精品| 国产视频内射| 国产午夜精品一二区理论片| 精品久久久久久久久亚洲| 久久97久久精品| 涩涩av久久男人的天堂| 亚洲精品aⅴ在线观看| av一本久久久久| 另类精品久久| 女人久久www免费人成看片| 国产成人91sexporn| 9色porny在线观看| 国产乱来视频区| 亚洲av二区三区四区| 精品国产一区二区久久| 日日撸夜夜添| 丝袜脚勾引网站| 人人妻人人爽人人添夜夜欢视频 | 亚洲av男天堂| 人妻少妇偷人精品九色| 一个人看视频在线观看www免费| 国产精品一区二区性色av| 我要看日韩黄色一级片| 在现免费观看毛片| 色5月婷婷丁香| 久久久久久久久大av| 国产淫片久久久久久久久| 久久久精品免费免费高清| 在线观看免费日韩欧美大片 | 一级二级三级毛片免费看| 在线亚洲精品国产二区图片欧美 | 午夜免费观看性视频| 亚洲美女黄色视频免费看| 99久久综合免费| 精品久久久噜噜| 国产黄频视频在线观看| 国产免费视频播放在线视频| 国产免费又黄又爽又色| 只有这里有精品99| 99九九在线精品视频 | 日本黄色日本黄色录像| 99久久中文字幕三级久久日本| 男人爽女人下面视频在线观看| 亚洲无线观看免费| 久久久久久久久大av| 亚洲人与动物交配视频| av有码第一页| 国产一区亚洲一区在线观看| 亚洲精品乱久久久久久| 最近的中文字幕免费完整| 亚洲在久久综合| 亚洲精品日韩av片在线观看| 91久久精品国产一区二区三区| 国产亚洲一区二区精品| 精品国产露脸久久av麻豆| 一区二区av电影网| 国产中年淑女户外野战色| 精品少妇黑人巨大在线播放| 国产在线视频一区二区| 国产亚洲欧美精品永久| videossex国产| 少妇被粗大的猛进出69影院 | 日本猛色少妇xxxxx猛交久久| 性色av一级| .国产精品久久| 97在线视频观看| 国产精品蜜桃在线观看| 日本-黄色视频高清免费观看| 亚洲av电影在线观看一区二区三区| 国产在线一区二区三区精| 天天躁夜夜躁狠狠久久av| 日日摸夜夜添夜夜爱| 美女中出高潮动态图| 国产一区二区在线观看日韩| 男人爽女人下面视频在线观看| 99久久精品一区二区三区| 国产欧美日韩综合在线一区二区 | 亚洲精品乱码久久久v下载方式| 国产精品一区二区性色av| 久久久久久久精品精品| 国产欧美日韩精品一区二区| 欧美bdsm另类| 秋霞在线观看毛片| 九九在线视频观看精品| 国产日韩一区二区三区精品不卡 | xxx大片免费视频| 亚洲精品中文字幕在线视频 | 高清不卡的av网站| 午夜免费男女啪啪视频观看| 一区二区av电影网| 精品少妇内射三级| 免费人成在线观看视频色| 精品久久国产蜜桃| av在线app专区| 美女主播在线视频| 国产美女午夜福利| 高清午夜精品一区二区三区| 啦啦啦中文免费视频观看日本| 欧美97在线视频| 久久精品国产自在天天线| 国产黄色视频一区二区在线观看| 国产男女超爽视频在线观看| 欧美另类一区| 纯流量卡能插随身wifi吗| 色网站视频免费| 18禁在线播放成人免费| 欧美一级a爱片免费观看看| 69精品国产乱码久久久| 亚洲,一卡二卡三卡| 国产免费又黄又爽又色| 又大又黄又爽视频免费| 欧美xxxx性猛交bbbb| 国产黄片美女视频| 狂野欧美激情性xxxx在线观看| 亚洲国产毛片av蜜桃av| 亚洲欧洲日产国产| 欧美高清成人免费视频www| 青青草视频在线视频观看| 国产欧美亚洲国产| 美女内射精品一级片tv| 制服丝袜香蕉在线| 97精品久久久久久久久久精品| 国产精品国产三级专区第一集| 国产视频内射| 国产欧美日韩综合在线一区二区 | 交换朋友夫妻互换小说| 国产精品国产三级国产专区5o| 亚洲熟女精品中文字幕| 亚洲欧美一区二区三区国产| 天堂中文最新版在线下载| 在线观看免费高清a一片| 少妇的逼水好多| 草草在线视频免费看| 欧美xxxx性猛交bbbb| 久久久久国产网址| 黑人高潮一二区| 久久人人爽av亚洲精品天堂| 精品久久久噜噜| 极品人妻少妇av视频| 国产精品久久久久久久久免| 亚洲真实伦在线观看| 久久久久精品性色| 精品国产露脸久久av麻豆| 日本爱情动作片www.在线观看| 美女中出高潮动态图| 一级毛片 在线播放| 国产91av在线免费观看| 中文欧美无线码| 亚洲国产精品成人久久小说| 免费观看av网站的网址| 王馨瑶露胸无遮挡在线观看| 日本黄色日本黄色录像| 丰满乱子伦码专区| 黄色欧美视频在线观看| 精品国产一区二区久久| 丝袜在线中文字幕| 亚洲精品自拍成人| 好男人视频免费观看在线| 亚洲国产精品999| 九色成人免费人妻av| 亚洲精品乱码久久久v下载方式| 国产精品偷伦视频观看了| 精品久久久精品久久久| 国产成人精品福利久久| 天堂8中文在线网| 成人午夜精彩视频在线观看| 亚洲一级一片aⅴ在线观看| 中文字幕免费在线视频6| 自拍偷自拍亚洲精品老妇| av一本久久久久| 精品一区二区三卡| 五月玫瑰六月丁香| 久久久国产欧美日韩av| 日本wwww免费看| 三上悠亚av全集在线观看 | videossex国产| 亚洲欧美日韩东京热| 亚洲情色 制服丝袜| 两个人免费观看高清视频 | 久久久久人妻精品一区果冻| 蜜桃久久精品国产亚洲av| 亚洲真实伦在线观看| 97精品久久久久久久久久精品| 边亲边吃奶的免费视频| 2021少妇久久久久久久久久久| 欧美日韩综合久久久久久| 国精品久久久久久国模美| 18+在线观看网站| 国产亚洲欧美精品永久| 国产免费又黄又爽又色| 亚州av有码| 乱码一卡2卡4卡精品| 极品人妻少妇av视频| 成年av动漫网址| 国产真实伦视频高清在线观看| 丝袜喷水一区| 日韩制服骚丝袜av| 成年av动漫网址| 久久99蜜桃精品久久| 免费黄色在线免费观看| 大话2 男鬼变身卡| 亚洲色图综合在线观看| 欧美+日韩+精品| 国产精品女同一区二区软件| 国产成人一区二区在线| 欧美性感艳星| 亚洲性久久影院| 国产精品一区二区在线不卡| 赤兔流量卡办理| 97在线视频观看| 人妻 亚洲 视频| 国产精品三级大全| 久久久久久久国产电影| 久久久久精品久久久久真实原创| 嘟嘟电影网在线观看| 国产精品久久久久久久电影| 久久午夜综合久久蜜桃| 日本91视频免费播放| 精品人妻偷拍中文字幕| 欧美xxxx性猛交bbbb| 成人毛片a级毛片在线播放| 一区二区三区乱码不卡18| 人妻制服诱惑在线中文字幕| 一区二区三区四区激情视频| 伊人久久精品亚洲午夜| av在线播放精品| 亚洲精华国产精华液的使用体验| 有码 亚洲区| 99热这里只有精品一区| 午夜日本视频在线| 亚洲婷婷狠狠爱综合网| 夫妻性生交免费视频一级片| 精品卡一卡二卡四卡免费| www.av在线官网国产| 久久精品久久精品一区二区三区| 一个人看视频在线观看www免费| 精品少妇久久久久久888优播| 久久av网站| 国产精品国产三级国产av玫瑰| 亚洲成人一二三区av| 国产午夜精品久久久久久一区二区三区| 啦啦啦啦在线视频资源| 亚洲精品久久久久久婷婷小说| a级毛色黄片| 亚洲欧美精品专区久久| 丝袜在线中文字幕| 一区二区av电影网| 老熟女久久久| 亚洲一级一片aⅴ在线观看| 国产成人精品一,二区| 五月玫瑰六月丁香| 成年女人在线观看亚洲视频| 日韩av免费高清视频| 80岁老熟妇乱子伦牲交| 亚洲精品国产成人久久av| 国产精品国产av在线观看| 亚洲精品国产色婷婷电影| a级一级毛片免费在线观看| 精华霜和精华液先用哪个| 欧美丝袜亚洲另类| 欧美日韩在线观看h| 日韩一区二区三区影片| 久久99精品国语久久久| 嘟嘟电影网在线观看| 国产极品粉嫩免费观看在线 | 男女国产视频网站| 免费看不卡的av| 欧美日韩视频高清一区二区三区二| 欧美+日韩+精品| 久久99热这里只频精品6学生| 国产91av在线免费观看| 亚洲欧洲国产日韩| 99热这里只有精品一区| 日本wwww免费看| 精品99又大又爽又粗少妇毛片| 亚洲欧美成人综合另类久久久| 赤兔流量卡办理| 精品久久久精品久久久| 噜噜噜噜噜久久久久久91| 久久久久久久久久成人| 国产淫片久久久久久久久| 国内揄拍国产精品人妻在线| 亚洲精华国产精华液的使用体验| 久久国内精品自在自线图片| av福利片在线| kizo精华| 在现免费观看毛片| 少妇人妻 视频| 久久久欧美国产精品| 一区二区三区乱码不卡18| 亚洲天堂av无毛| 午夜视频国产福利| 亚洲国产最新在线播放| 免费人妻精品一区二区三区视频| 亚洲第一av免费看| 91精品伊人久久大香线蕉| 国产黄色免费在线视频| 精品亚洲成国产av| 国精品久久久久久国模美| 国产黄色免费在线视频| 国产成人精品一,二区| 国产成人a∨麻豆精品| 在线观看免费日韩欧美大片 | .国产精品久久| 一级av片app| 一级a做视频免费观看|