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

    高精度UT1-UTC差分預報方法研究

    2014-03-06 08:28:19陳略唐歌實胡松杰平勁松許雪晴夏金超
    深空探測學報 2014年3期
    關鍵詞:階數(shù)高精度差分

    陳略,唐歌實,胡松杰,平勁松,許雪晴,夏金超

    (1.北京航天飛行控制中心航天飛行動力學技術重點實驗室;北京100094; 2.中國科學院國家天文臺;北京100012;3.中國科學院上海天文臺;上海200030)

    高精度UT1-UTC差分預報方法研究

    陳略1,唐歌實1,胡松杰1,平勁松2,許雪晴3,夏金超1

    (1.北京航天飛行控制中心航天飛行動力學技術重點實驗室;北京100094; 2.中國科學院國家天文臺;北京100012;3.中國科學院上海天文臺;上海200030)

    針對衛(wèi)星導航所需的高精度地球定向參數(shù)(EOP)中的UT1-UTC預報問題,提出了基于雙差分LS+ AR的UT1-UTC參數(shù)預報方法。對UT1-UTC觀測數(shù)據(jù)進行跳秒檢測、固體地球帶諧潮汐項改正,然后對改正后的UT1-UTC數(shù)據(jù)進行雙差分處理,增強數(shù)據(jù)平穩(wěn)性;采用最小二乘擬合(LS)與自回歸(AR)分析方法對差分處理后的數(shù)據(jù)進行分析與預報;對預報結果進行逆差分處理與潮汐項改正外推、跳秒恢復,獲取高精度的UT1-UTC預報值。通過與國際EOP_PCC預報結果對比表明,UT1-UTC短期預報精度與EOP_PCC較優(yōu)的預報精度相當,其中1天UT1-UTC預報精度優(yōu)于0.03 ms,優(yōu)于EOP_PCC預報結果。介紹了北京航天飛行控制中心的UT1-UTC每日例行預報情況。

    地球定向參數(shù);UT1-UTC預報;差分處理;AR模型;潮汐項改正

    0 引 言

    高精度地球定向參數(shù)(Earth orientation parameters,EOP)是地球參考架與天球參考架實時轉換時的基礎參數(shù)[1]。VLBI、GNSS、DORIS、LLR、SLR等現(xiàn)代空間大地測量技術可精確地測量地球定向參數(shù)[2-3]。高精度EOP在載人航天、深空探測等任務中具有重要價值,尤其是在航天器實時、高精度導航方面不可或缺。目前,由于EOP解算結果在時間上有滯后,因此EOP參數(shù)高精度預報對于航天器導航顯得十分重要,且EOP預報對于與EOP相關的地球物理動力學理論分析具有重要價值。

    EOP預報中的極移預報、UT1-UTC預報是EOP預報中的熱點研究內容。目前,有多種預報方法應用于EOP的預報中,如最小二乘(LS)外推法[4]、最小二乘聯(lián)合自回歸分析法(LS+AR)[5]、神經網(wǎng)絡預報法[6]、譜分析最小二乘外推法[7]、小波分解與自回歸分析方法[8]等。國際地球自轉與參考系服務(IERS)組織于2005年10月發(fā)起了首次全球地球定向參數(shù)預報比賽活動(EOP_PCC),旨在號召全世界EOP預報人員運用不同方法來預報EOP,用以分析現(xiàn)有方法的特點及使用情況。LS+AR預報方法是預報精度最高的方法之一。但運用LS+ AR預報方法時,LS模型參數(shù)的選擇,AR模型階數(shù)的確定,非平穩(wěn)序列AR分析,EOP預報精度評估等是進行EOP預報的關鍵問題。現(xiàn)有LS+AR預報方法在對數(shù)據(jù)平穩(wěn)性要求方面的研究較為欠缺,在一定程度上限制了EOP預報精度。

    本文從分析數(shù)據(jù)的平穩(wěn)特性出發(fā),提出了雙差分LS+AR預報方法,用于實現(xiàn)對EOP中UT1-UTC高精度預報,并獲取了高精度的UT1-UTC短期預報結果。與此同時,同時介紹了與EOP_PCC進行精度比較的結果,此外,介紹了北京航天飛行控制中心EOP參數(shù)高精度預報軟件每日例行情況。

    1 預報理論與方法

    1.1 最小二乘(LS)模型

    對UT1-UTC觀測數(shù)據(jù)進行最小二乘模型擬合時,模型包括線性項與周期項(周年項、半年項、9.3年項、18.6年項等),其最小二乘模型如式(1)所示。

    式中:t為序列的UTC時間,單位為年;A,B,C,D1, D2,E1,E2,…為模型的擬合參數(shù);p1,p2,…為擬合周期,可根據(jù)實際情況確定模型形狀。

    1.2 AR模型

    1.2.1 AR模型定義

    AR模型表示平穩(wěn)隨機序列xt(t=1,2,…,N)與t時刻以前的規(guī)律性變化和t時刻白噪聲的關系,其數(shù)學模型如式(2)所示

    式中:φ1,φ2,…,φp為模型參數(shù);a為白噪聲;p為模型階數(shù),稱上式為p階自回歸模型,簡記為AR(p);為白噪聲的方差。

    1.2.2 AR模型階數(shù)確定

    運用AR模型關鍵在于確定模型階數(shù)p,可采用最終預報誤差(final prediction error criterion, FPE)準則、AIC(akake information criterion)準則、BIC準則、矩陣奇異值分解(singular value decomposition,SVD)準則來確定模型的階數(shù)。

    本文采用解Yule-Walker方程方法來估計模型階數(shù),其基本原理[9]如下:

    根據(jù)自協(xié)方差函數(shù)的定義,一個平穩(wěn)、正態(tài)、零均值的隨機過程{xt}的自協(xié)方差定義為

    其中,E表示集合平均算子,即數(shù)學期望,{xt}的零均值性導致其自協(xié)方差函數(shù)與自相關函數(shù)是一致的,而且作為Rk的特例,當k=0時,得{xt}的方差函數(shù)

    自相關系數(shù)定義為

    顯然,ρ0=1。由許瓦茲不定式,知,的取值范圍為

    對于平穩(wěn)時間序列,當ρk=0時,表示序列前后無關,即為白噪聲;當時,表示序列前后取值的關系不但完全確定,而且前后取值相等,即序列為常值序列。

    對式(2)中左右兩邊乘以xt-k,再取數(shù)學期望并除以R0,得到

    分別令k=1,2,…,p,并注意自協(xié)方差函數(shù)是偶函數(shù)的性質,R-k=Rk,因而ρ-k=ρk,可得線性方程組如下

    分別令k=1,2,…,n,則可得方程的個數(shù)等于未知數(shù)φ1,φ2,…,φn的個數(shù),即上式等效為式(9)。

    此式成為Yule-Walker方程,其中ρ為n×1的列矩陣,T為n階仿陣,稱為Toeplitz矩陣,如式(10)所示。

    T中各元素可按定義估計出。由于T是滿秩的,因此可容易求出模型參數(shù)如式(11)所示。

    本文采用FPE準則來確定AR模型階數(shù)。FPE準則函數(shù)如式(12)所示。

    1.3 精度評估

    為了評估預報精度,采用EOP預報中的通用精度評估準則-平均絕對誤差(mean absulute eroor, MAE)作為預報結果的精度指標,其計算公式如下:

    式中:o為實際觀測值;p為預報值;i為預報跨度;n為預報期數(shù)。

    2 雙差分LS+AR預報流程

    雙差分LS+AR的UT1-UTC預報方法的實現(xiàn)流程如圖1所示。

    圖1 雙差分LS+AR UT1-UTC預報流程圖Fig.1 UT1-UTC prediction process of differential LS+AR method

    3 計算分析

    利用IERS公布的EOP數(shù)據(jù)進行預報與精度驗證。選取EOP 05C04的極移數(shù)據(jù)進行分析,便于與EOP_PCC結果進行比較。時間跨度為1989年1月1日~2009年12月31日,數(shù)據(jù)采樣頻率為1點/天。預報的基礎時間序列為1989年1月1日~2009年10月2日,對2009年10月3日至2009年12月31日的EOP極移進行預報,并將預報值與真實觀測值進行比較,評估預報精度。預報跨度分別從1~30天,預報期數(shù)為90期。具體預報結果如下所述。

    圖2為1989年—2009年的UT1-UTC的時間變化圖,其中圖2(a)為原始的UT1-UTC波形圖;圖2(b)為剔除跳秒后的UT1-UTC圖;圖2(c)為固體地球帶諧潮汐項改正圖,依據(jù)IERS2010規(guī)范進行改正;圖2(d)是經過跳秒剔除與潮汐項改正后的UT1-UTC時間序列。

    圖3為1989年—2009年UT1-UTC雙差分后的LS擬合結果。從圖中可以看出,經過雙差分后的UT1-UTC序列具有明顯的平穩(wěn)性。利用LS所得參數(shù)進行外推,利用AR模型對LS擬合后的殘差進行預報,將LS外推值與AR模型外推值進行相加,即得跨度為1天的預報值。其中,在進行AR預報過程中,AR模型階數(shù)的確定是AR預報最關鍵步驟,本文采用FPE準則確定模型階數(shù),在階數(shù)確定過程中采用迭代搜索方法確定最優(yōu)模型階數(shù),階數(shù)搜索范圍為1~100,采用FPE準則計算出每一項FPE值,其中最小FPE值對應階數(shù)即為最優(yōu)的AR模型階數(shù)。

    按照圖1的處理流程,對預報初值進行逆雙差分處理、固體地球帶諧潮汐項改正外推、跳秒恢復后,即可獲得第一天的UT1-UTC預報值。在得到第1天的預報值之后,將第1天的預報值加入到基礎序列中,重新按照圖1的處理流程進行第2天的UT1-UTC預報。以此類推,可獲取不同預報跨度的UT1-UTC預報結果。

    圖2 1989-2009年UT1-UTC的時間序列分布圖Fig.2 Time sequence of UT1-UTC from 1989 to 2009

    圖3 1989-2009年UT1-UTC雙差分LS擬合結果Fig.3 Differential LS fitting results of UT1-UTC from 1989 to 2009

    基于以上方法進行90期預報,并利用平均絕對誤差(MAE)準則評估UT1-UTC預報精度,其結果如圖4所示,具體數(shù)值如表1所示。

    圖4 BACC的UT1-UTC短期預報結果Fig.4 Short term prediction results of UT1-UTC in Beijing Aerospace Control Center

    表1 BACC的UT1-UTC短期預報MAE結果Table 1 MAE of UT1-UTC short term prediction in Beijing Aerospace Control Center

    全球性的EOP預報比對運動(earth orientation parameters prediction comparison campaign,EOP_PCC),共有8個國家12名時間序列分析領域的頂尖級專家學者參與,涉及20余中預報方法。EOP_PCC運動分為10天超短期預報, 30天短期預報,500天長期預報。

    本文將預報結果與EOP_PCC的結果進行比對分析,分析進行超短期預報比對,與短期預報比對。圖5為EOP_PCC的預報結果,圖6為本文的預報結果。

    從結果可以看出,本文的UT1-UTC短期預報精度與EOP_PCC的預報精度相當。對于1天預報精度,本文UT1-UTC預報精度優(yōu)于EOP_PCC預報精度。本文獲取的1天UT1-UTC預報精度優(yōu)于0.03 ms。EOP_PCC的UT1-UTC一天預報誤差[1]在0.08 ms水平。

    圖5 EOP_PCC的UT1-UTC短期預報結果Fig.5 Short term prediction results of UT1-UTC in EOP_PCC

    圖6 BACC的UT1-UTC短期預報結果Fig.6 Short term prediction results of UT1-UTC in Beijing Aerospace Control Center

    4 UT1-UTC每日例行預報

    北京航天飛行控制中心(BACC)飛行動力學重點實驗室開發(fā)了EOP高精度預報軟件,可實現(xiàn)對UT1-UTC每日例行預報與發(fā)布。該軟件每日從互聯(lián)網(wǎng)上自動化下載UT1-UTC解算數(shù)據(jù),而后進行UT1-UTC預報并生成預報結果,同時輸出BACC/ IERS/美國海軍天文臺(USNO)的UT1-UTC預報的比較結果,用以每日評估BACC的UT1-UTC預報。BACC的EOP預報軟件輸入的UT1-UTC解算數(shù)據(jù)由兩個部分組成:IERS網(wǎng)站公布的EOP 08 C04序列+USNO最近一個月的解算序列。

    圖7 2014-11-03 UT1-UTC 30天預報結果比較Fig.7 UT1-UTC 30 days prediction comparing results on 2014-11-03

    圖8 2014-11-04 UT1-UTC 30天預報結果比較Fig.8 UT1-UTC 30 days prediction comparing results on 2014-11-04

    圖9 2014-11-05 UT1-UTC 30天預報結果比較Fig.9 UT1-UTC 30 days prediction comparing results on 2014-11-05

    圖7~圖9顯示了2014年11月3—5日, BACC的UT1-UTC預報結果與IERS/USNO結果對比圖。在圖7~圖9中,圓圈表示BACC的預報結果,正方形表示USNO的預報結果,梅花形表示IERS的預報結果。從以上圖中可以看出,BACC的UT1-UTC預報結果趨勢與IERS、USNO基本一致,且BACC 30天的UT1-UTC預報結果基本位于IERS與USNO之間,這說明了北京中心UT1-UTC預報結果有效性。

    5 結束語

    本文提出了基于雙差分LS+AR的UT1-UTC預報方法,詳細介紹了理論方法與處理流程,并介紹了BACC的UT1-UTC每日例行預報情況。通過對IERS實際的EOP中UT1-UTC進行預報與比對分析,本文的UT1-UTC預報精度與EOP_PCC較優(yōu)的短期預報精度相當。本文獲取的1天UT1-UTC預報精度優(yōu)于0.03 ms,優(yōu)于EOP_PCC每日預報結果。

    相比于傳統(tǒng)LS+AR的地球定向參數(shù)預報方法,本文重點從AR模型輸入數(shù)據(jù)的平穩(wěn)性要求方面進行考慮,采用了雙差分預處理方法,用以提高UT1-UTC數(shù)據(jù)的平穩(wěn)性,使其更加適應利用AR模型進行預報。此外,AR模型階數(shù)確定采用了迭代搜索方式選擇最優(yōu)模型階數(shù)。通過以上兩項步驟,可有效提高對UT1-UTC的高精度預報。

    后續(xù)將進一步研究本文提出的雙差分LS+AR預報方法在極移、LOD預報上的應用,同時將進一步評估該方法在對EOP參數(shù)長期預報中的精度與適用性。

    [1]Kalarus M,Schuh H,Kosek W,et al.Achievements of the Earth orientation parameters prediction comparison campaign [J].Journal of Geodynamics,2010(62):587-596.

    [2]Tobias Nilsson,Robert Heinkelmann,Maria Karbon,et al. Earth orientation parameters estimated from VLBI during the CONT11 campaign[J].Journal of Geodesy,2014(88): 491-502.

    [3]許雪晴,周永宏.地球定向參數(shù)高精度預報方法研究[J].飛行器測控學報,2010,29(2):70-76.[Xu X Q,Zhou Y H. High precision prediction method of Earth orientation parameters[J].Journal of Spacecraft TT&C Technology, 2010,29(2):70-76.]

    [4]Xu X,Zhou Y,Liao X.Short-term earth orientation parameters predictions by combination of the least-squares, AR model and Kalman filter[J].Journal of Geodynamics, 2010(62):83-86.

    [5]Guo J Y,Li Y B,Dai C L,et al.A technique to improve the accuracy of Earth orientation prediction algorithms based on least squares extrapolation[J].Journal of Geodynamics, 2013(70):36-48.

    [6]Wang Q J,Liao D C,Zhou Y H.Real-time rapid prediction of variations of Earth's rotational rate[J].Chin.Sci.Bull, 2008(53):969-973.

    [7]Akulenko L D,Kumakshev S A,Rykhlova L V.A model for the polar motion of the deformable Earth adequate for astrometric data[J].Astron.Rep.,2002(46):74-82.

    [8]Kosek W,Kalarus M,Johnson T J,et al.A comparison of LOD and UT1-UTC forecasts by different combination prediction techniques[J].Artif.Satell,2005(40):119-125.

    [9]楊叔子.時間序列分析的工程應用[M].武漢:華中科技大學出版社,2007:184-189.[Yang S Z.Time series analysis in engineering application[M].Wuhan:Huazhong University of Science&Technology Press,2007:184-189.]

    通信地址:北京海淀區(qū)5130信箱120號(100094)

    電話:(010)66362814

    E-mail:luechen0912@163.com

    [責任編輯:宋宏]

    High Accuracy Differential Prediction of UT1-UTC

    CHEN Lue1,TANG Geshi1,HU Songjie1,PING Jinsong2,XU Xueqing3,XIA Jinchao1
    (1.National Key Laboratory of Science and Technology on Aerospace Flight Dynamic,Beijing Aerospace Control Center, 100094 Beijing,China;2.National Astronomical Observatory,Chinese Academy of Sciences,100012 Beijing,China; 3.Shanghai Astronomical Observatory,Chinese Academy of Sciences,200030 Shanghai,China)

    This paper proposes a prediction method of UT1-UTC in Earth orientation parameters(EOP)by dual differential least-squares(LS)and autoregressive(AR)model.Firstly,leap seconds are removed in UT1-UTC observations,and Earth zonal harmonic tidal are corrected.Then,the corrected UT1-UTC are processed by dual differential method,the stationarity of polar motion parameters is improved.LS+AR method is utilized to analyze the dual differential UT1-UTC to obtain the preliminary prediction results.Finally,the preliminary prediction results are processed by inverse dual differential method,and tidal correction are extrapolated and leap seconds are recovered to obtain high accuracy UT1-UTC prediction results.The prediction results are compared with EOP prediction comparison campaign(EOP_PCC)results.It shows that the short-term UT1-UTC parameters prediction error is at the same level of EOP_PCC.The one day prediction accuracy of UT1-UTC is less than 0.03ms,which is better than EOP_PCC one day UT1-UTC prediction accuracy.And,the daily routine UT1-UTC prediction in Beijing Aerospace Control Center is introduced.

    Earth orientation parameters;UT1-UTC prediction;dual differential;AR model;earth zonal harmonic tidal correction

    P228

    :A

    :2095-7777(2014)03-0230-06

    10.15982/j.issn.2095-7777.2014.03.012

    陳略(1983—),男,工程師,主要研究方向:地球自轉參數(shù)高精度解算與預報,無線電測量與科學應用。

    2014-07-28

    2014-08-19

    國家自然科學基金(41304026);國家重點基礎研究發(fā)展計劃(973計劃)項目(2015CB857101)

    猜你喜歡
    階數(shù)高精度差分
    關于無窮小階數(shù)的幾點注記
    數(shù)列與差分
    確定有限級數(shù)解的階數(shù)上界的一種n階展開方法
    高抗擾高精度無人機著艦縱向飛行控制
    船載高精度星敏感器安裝角的標定
    基于高精度測角的多面陣航測相機幾何拼接
    高精度免熱處理45鋼的開發(fā)
    山東冶金(2015年5期)2015-12-10 03:27:41
    基于差分隱私的大數(shù)據(jù)隱私保護
    一種新的多址信道有效階數(shù)估計算法*
    電訊技術(2014年1期)2014-09-28 12:25:26
    關于動態(tài)電路階數(shù)的討論
    国产成人欧美在线观看| 村上凉子中文字幕在线| 国产亚洲欧美在线一区二区| 男插女下体视频免费在线播放| 婷婷丁香在线五月| 成人精品一区二区免费| 欧美日韩黄片免| www.色视频.com| 欧美性猛交╳xxx乱大交人| 久9热在线精品视频| 99久久成人亚洲精品观看| 日韩国内少妇激情av| 神马国产精品三级电影在线观看| 丰满的人妻完整版| 色播亚洲综合网| 久久精品国产99精品国产亚洲性色| 亚洲中文日韩欧美视频| 97碰自拍视频| 欧美性猛交╳xxx乱大交人| 俺也久久电影网| 中国美女看黄片| 日本五十路高清| 亚洲精品一区av在线观看| 精品人妻1区二区| 久久国产乱子伦精品免费另类| 一级a爱片免费观看的视频| 国产精品久久久久久精品电影| 亚洲av美国av| 国产av不卡久久| 久久久精品欧美日韩精品| 天天躁日日操中文字幕| 中文资源天堂在线| 女人高潮潮喷娇喘18禁视频| 69人妻影院| 国产精品亚洲av一区麻豆| 日韩欧美在线乱码| 99国产极品粉嫩在线观看| 国产一区二区亚洲精品在线观看| bbb黄色大片| 可以在线观看毛片的网站| 欧美日韩黄片免| 中文字幕人妻丝袜一区二区| 男插女下体视频免费在线播放| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 一区二区三区国产精品乱码| 三级国产精品欧美在线观看| 欧美性猛交黑人性爽| 在线观看日韩欧美| 麻豆一二三区av精品| 99久久久亚洲精品蜜臀av| 免费高清视频大片| 狂野欧美白嫩少妇大欣赏| 国产伦精品一区二区三区四那| 91久久精品电影网| 亚洲狠狠婷婷综合久久图片| 欧美黄色淫秽网站| 老熟妇乱子伦视频在线观看| 一区福利在线观看| 免费搜索国产男女视频| 一本精品99久久精品77| ponron亚洲| 久久香蕉国产精品| 国产亚洲欧美98| 国产熟女xx| 女人高潮潮喷娇喘18禁视频| 国产精品久久久久久人妻精品电影| 精品日产1卡2卡| 免费看光身美女| 男女下面进入的视频免费午夜| 欧美黑人巨大hd| 国产精品 欧美亚洲| 变态另类成人亚洲欧美熟女| 久久精品国产清高在天天线| 国产成人影院久久av| 真实男女啪啪啪动态图| 99久久久亚洲精品蜜臀av| 国产一区二区激情短视频| 日韩有码中文字幕| netflix在线观看网站| 神马国产精品三级电影在线观看| 深爱激情五月婷婷| 国产单亲对白刺激| 国产乱人伦免费视频| avwww免费| 亚洲五月天丁香| 久久久久精品国产欧美久久久| 亚洲中文字幕日韩| 88av欧美| 亚洲成a人片在线一区二区| 欧美乱码精品一区二区三区| 一级a爱片免费观看的视频| 69av精品久久久久久| 色尼玛亚洲综合影院| 熟女电影av网| 嫩草影院精品99| 国产69精品久久久久777片| 国产亚洲精品一区二区www| 成人三级黄色视频| 国产精品 国内视频| 中文在线观看免费www的网站| 一区二区三区激情视频| 国产成人福利小说| 久久草成人影院| 欧美不卡视频在线免费观看| 美女高潮喷水抽搐中文字幕| 国产成人av教育| 搡老熟女国产l中国老女人| 亚洲美女黄片视频| 免费看十八禁软件| 成人欧美大片| 3wmmmm亚洲av在线观看| 给我免费播放毛片高清在线观看| 9191精品国产免费久久| 精华霜和精华液先用哪个| 日韩av在线大香蕉| 亚洲av成人精品一区久久| 12—13女人毛片做爰片一| 国产精品三级大全| 欧美极品一区二区三区四区| 国产一区二区三区视频了| 国产精品三级大全| 欧美中文综合在线视频| 日韩欧美 国产精品| 欧美性猛交黑人性爽| 女人被狂操c到高潮| www日本在线高清视频| 草草在线视频免费看| 男人的好看免费观看在线视频| 免费在线观看日本一区| 好男人电影高清在线观看| 极品教师在线免费播放| 国产久久久一区二区三区| 亚洲成人免费电影在线观看| 国产精品自产拍在线观看55亚洲| 欧美国产日韩亚洲一区| 老汉色∧v一级毛片| 国产一区二区三区在线臀色熟女| 欧美一级毛片孕妇| 日韩欧美 国产精品| 18禁在线播放成人免费| 法律面前人人平等表现在哪些方面| 欧美区成人在线视频| 长腿黑丝高跟| 很黄的视频免费| 欧美日韩瑟瑟在线播放| 久久精品影院6| 免费在线观看亚洲国产| 99视频精品全部免费 在线| 久久香蕉国产精品| 欧美另类亚洲清纯唯美| 午夜福利高清视频| 亚洲精华国产精华精| 高潮久久久久久久久久久不卡| 男人舔奶头视频| 最近视频中文字幕2019在线8| 日韩有码中文字幕| 日本三级黄在线观看| 九色国产91popny在线| 欧美精品啪啪一区二区三区| 国产高潮美女av| 黄片小视频在线播放| 亚洲黑人精品在线| 女人被狂操c到高潮| 在线国产一区二区在线| 国内精品美女久久久久久| 国产精品一及| 夜夜爽天天搞| 亚洲欧美精品综合久久99| 国产成人欧美在线观看| 小蜜桃在线观看免费完整版高清| 国产一级毛片七仙女欲春2| 久久久久久久久久黄片| 日韩欧美精品免费久久 | 性欧美人与动物交配| 波多野结衣巨乳人妻| 国内少妇人妻偷人精品xxx网站| 宅男免费午夜| 久久精品91无色码中文字幕| 欧美极品一区二区三区四区| www日本在线高清视频| 美女大奶头视频| 99国产综合亚洲精品| 在线看三级毛片| 亚洲国产欧美网| 亚洲 国产 在线| 少妇的逼好多水| 欧美乱妇无乱码| 国产在线精品亚洲第一网站| 激情在线观看视频在线高清| 91av网一区二区| 日本a在线网址| 欧美成人免费av一区二区三区| 美女被艹到高潮喷水动态| av在线蜜桃| 国产亚洲欧美在线一区二区| 熟女人妻精品中文字幕| 亚洲av成人精品一区久久| 熟女电影av网| 在线免费观看的www视频| 性色avwww在线观看| 欧美中文综合在线视频| 激情在线观看视频在线高清| 97人妻精品一区二区三区麻豆| 听说在线观看完整版免费高清| 国产美女午夜福利| av专区在线播放| 成人av一区二区三区在线看| 亚洲av成人av| 日韩亚洲欧美综合| 看免费av毛片| 悠悠久久av| 亚洲内射少妇av| 别揉我奶头~嗯~啊~动态视频| 免费人成在线观看视频色| 人妻夜夜爽99麻豆av| 精品国产三级普通话版| 精品欧美国产一区二区三| 欧洲精品卡2卡3卡4卡5卡区| 一级毛片女人18水好多| 色综合欧美亚洲国产小说| 老司机福利观看| 久久久久国产精品人妻aⅴ院| 人人妻人人澡欧美一区二区| 综合色av麻豆| 丰满人妻熟妇乱又伦精品不卡| 大型黄色视频在线免费观看| 午夜精品在线福利| 两性午夜刺激爽爽歪歪视频在线观看| a级毛片a级免费在线| 国产精品永久免费网站| 五月伊人婷婷丁香| bbb黄色大片| 三级男女做爰猛烈吃奶摸视频| 最近视频中文字幕2019在线8| 欧美中文日本在线观看视频| 亚洲一区二区三区不卡视频| 日本与韩国留学比较| 国产精品女同一区二区软件 | av在线蜜桃| 国产一区二区激情短视频| 黄色片一级片一级黄色片| 国产熟女xx| 美女高潮喷水抽搐中文字幕| 99久久综合精品五月天人人| 国产精品影院久久| 一级作爱视频免费观看| 精品福利观看| 超碰av人人做人人爽久久 | 香蕉av资源在线| 首页视频小说图片口味搜索| 伊人久久大香线蕉亚洲五| 色精品久久人妻99蜜桃| 757午夜福利合集在线观看| 国产一区二区亚洲精品在线观看| 久久久久久久亚洲中文字幕 | 女人十人毛片免费观看3o分钟| 国产av麻豆久久久久久久| 久久婷婷人人爽人人干人人爱| 亚洲人成伊人成综合网2020| 国产爱豆传媒在线观看| 国产亚洲av嫩草精品影院| 又黄又爽又免费观看的视频| 欧美日韩中文字幕国产精品一区二区三区| 国产在视频线在精品| 丁香欧美五月| 美女高潮的动态| 成人av一区二区三区在线看| 成年人黄色毛片网站| 网址你懂的国产日韩在线| 淫秽高清视频在线观看| 国产av麻豆久久久久久久| 国产精品一及| av黄色大香蕉| 婷婷六月久久综合丁香| 精品午夜福利视频在线观看一区| 日韩欧美精品v在线| 亚洲人成电影免费在线| 我的老师免费观看完整版| 午夜福利在线观看吧| 国产成+人综合+亚洲专区| 国产成人福利小说| 免费人成在线观看视频色| 精品国产美女av久久久久小说| 欧美3d第一页| 午夜免费成人在线视频| bbb黄色大片| 国产极品精品免费视频能看的| 啪啪无遮挡十八禁网站| 少妇熟女aⅴ在线视频| 国产一区二区三区在线臀色熟女| 小蜜桃在线观看免费完整版高清| 久久久久久国产a免费观看| 男女那种视频在线观看| 久久久久久人人人人人| 免费看美女性在线毛片视频| 狂野欧美激情性xxxx| 内射极品少妇av片p| av欧美777| 国产视频内射| 女人高潮潮喷娇喘18禁视频| 日本撒尿小便嘘嘘汇集6| 欧美一区二区精品小视频在线| 黄片大片在线免费观看| 久久香蕉精品热| 欧美激情在线99| 无限看片的www在线观看| 淫秽高清视频在线观看| 久久久久久久久中文| 国产av不卡久久| 真实男女啪啪啪动态图| 久久久久久久久久黄片| 亚洲人成网站高清观看| 五月玫瑰六月丁香| 国内精品久久久久久久电影| 12—13女人毛片做爰片一| a在线观看视频网站| 一个人看的www免费观看视频| 观看美女的网站| 9191精品国产免费久久| 黄色视频,在线免费观看| 亚洲av中文字字幕乱码综合| 天堂动漫精品| 脱女人内裤的视频| 亚洲精品国产精品久久久不卡| 久久久精品大字幕| 亚洲18禁久久av| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 他把我摸到了高潮在线观看| 99久久成人亚洲精品观看| 国内揄拍国产精品人妻在线| svipshipincom国产片| av国产免费在线观看| e午夜精品久久久久久久| 91字幕亚洲| 在线观看日韩欧美| 伊人久久大香线蕉亚洲五| 国产真实乱freesex| 色哟哟哟哟哟哟| 久久性视频一级片| 三级国产精品欧美在线观看| 一级作爱视频免费观看| 国产乱人视频| 免费在线观看亚洲国产| 亚洲乱码一区二区免费版| 欧美色欧美亚洲另类二区| 精品人妻偷拍中文字幕| 精品无人区乱码1区二区| 91字幕亚洲| 国产美女午夜福利| 午夜免费成人在线视频| 欧美黄色淫秽网站| 欧美av亚洲av综合av国产av| 亚洲av不卡在线观看| 啦啦啦免费观看视频1| 国产v大片淫在线免费观看| 亚洲国产精品999在线| 夜夜夜夜夜久久久久| 最好的美女福利视频网| 久久草成人影院| 少妇人妻一区二区三区视频| 欧美日韩亚洲国产一区二区在线观看| 草草在线视频免费看| 国产成人啪精品午夜网站| 久久精品夜夜夜夜夜久久蜜豆| 午夜福利成人在线免费观看| 男女做爰动态图高潮gif福利片| 欧美另类亚洲清纯唯美| 日韩欧美一区二区三区在线观看| 国产高清videossex| 久久精品国产清高在天天线| 亚洲在线自拍视频| 国产一区二区激情短视频| 免费一级毛片在线播放高清视频| 亚洲av电影在线进入| 在线国产一区二区在线| 免费观看人在逋| 老司机深夜福利视频在线观看| 性色av乱码一区二区三区2| 国产精品99久久久久久久久| 久久6这里有精品| 老熟妇乱子伦视频在线观看| 久久久久国内视频| 精品一区二区三区人妻视频| 日韩有码中文字幕| 欧美又色又爽又黄视频| 日本黄色片子视频| 美女大奶头视频| 色综合欧美亚洲国产小说| 少妇丰满av| 国产久久久一区二区三区| 国产免费av片在线观看野外av| 精品人妻一区二区三区麻豆 | 亚洲国产中文字幕在线视频| 高清日韩中文字幕在线| 九色国产91popny在线| 亚洲18禁久久av| 免费在线观看日本一区| 中文字幕熟女人妻在线| 两性午夜刺激爽爽歪歪视频在线观看| 啪啪无遮挡十八禁网站| av天堂中文字幕网| 一本一本综合久久| 大型黄色视频在线免费观看| 欧美黄色淫秽网站| 久久这里只有精品中国| 最近最新中文字幕大全电影3| 男女床上黄色一级片免费看| 真人做人爱边吃奶动态| 少妇裸体淫交视频免费看高清| 久久精品亚洲精品国产色婷小说| 色播亚洲综合网| 在线a可以看的网站| 禁无遮挡网站| 97碰自拍视频| 精品午夜福利视频在线观看一区| 国产精品久久久久久精品电影| 九色国产91popny在线| 天堂√8在线中文| 小蜜桃在线观看免费完整版高清| 国产黄色小视频在线观看| 成人高潮视频无遮挡免费网站| 老鸭窝网址在线观看| av福利片在线观看| 国产精品嫩草影院av在线观看 | 9191精品国产免费久久| 岛国在线免费视频观看| 中文亚洲av片在线观看爽| 欧美色视频一区免费| 日本免费一区二区三区高清不卡| 亚洲成a人片在线一区二区| 亚洲精品成人久久久久久| 99热这里只有精品一区| 欧美黑人欧美精品刺激| 人人妻人人看人人澡| 叶爱在线成人免费视频播放| 色老头精品视频在线观看| 亚洲精品久久国产高清桃花| 悠悠久久av| 国产日本99.免费观看| 精华霜和精华液先用哪个| www.熟女人妻精品国产| 亚洲av成人av| 国产伦一二天堂av在线观看| 国产精品 国内视频| 精品一区二区三区视频在线观看免费| 日韩有码中文字幕| 97人妻精品一区二区三区麻豆| 日本三级黄在线观看| 亚洲精品色激情综合| 精品久久久久久久人妻蜜臀av| 少妇裸体淫交视频免费看高清| 精品福利观看| 亚洲国产色片| 久久人人精品亚洲av| 1024手机看黄色片| 国产爱豆传媒在线观看| 成年女人毛片免费观看观看9| 99久久精品热视频| 中文字幕久久专区| 欧美日韩乱码在线| 熟女电影av网| 精品无人区乱码1区二区| 五月伊人婷婷丁香| 99热这里只有是精品50| 乱人视频在线观看| 老司机在亚洲福利影院| 久9热在线精品视频| 国产高潮美女av| 亚洲一区高清亚洲精品| 久久99热这里只有精品18| 看免费av毛片| 国产黄色小视频在线观看| 又紧又爽又黄一区二区| 国产成人av教育| 五月伊人婷婷丁香| 天堂网av新在线| av在线天堂中文字幕| 婷婷精品国产亚洲av| 麻豆国产av国片精品| 日本五十路高清| 99热这里只有精品一区| 久久精品影院6| 久久精品亚洲精品国产色婷小说| www.www免费av| 好男人电影高清在线观看| 精品日产1卡2卡| 搡女人真爽免费视频火全软件 | 欧美中文日本在线观看视频| 国产乱人视频| 一本一本综合久久| 国产精品亚洲一级av第二区| 母亲3免费完整高清在线观看| 亚洲最大成人中文| 午夜亚洲福利在线播放| 欧美日韩国产亚洲二区| 久久久精品欧美日韩精品| 日韩欧美在线乱码| 淫秽高清视频在线观看| 看黄色毛片网站| 午夜福利免费观看在线| 天堂av国产一区二区熟女人妻| 国产一区二区激情短视频| 脱女人内裤的视频| 免费高清视频大片| 国产精品电影一区二区三区| 欧美日韩国产亚洲二区| 一边摸一边抽搐一进一小说| 亚洲 国产 在线| 欧美极品一区二区三区四区| 哪里可以看免费的av片| 两性午夜刺激爽爽歪歪视频在线观看| 免费观看的影片在线观看| 国产av在哪里看| svipshipincom国产片| 亚洲av成人av| 成人永久免费在线观看视频| 午夜免费激情av| 国产av在哪里看| 99久久成人亚洲精品观看| 90打野战视频偷拍视频| 欧美最新免费一区二区三区 | 亚洲av电影不卡..在线观看| 日韩欧美三级三区| 搡老熟女国产l中国老女人| www.www免费av| 国产av在哪里看| 一本精品99久久精品77| 国产精品精品国产色婷婷| 国产精品亚洲一级av第二区| 别揉我奶头~嗯~啊~动态视频| av国产免费在线观看| 国内精品久久久久精免费| 日本免费a在线| 亚洲中文字幕一区二区三区有码在线看| 美女高潮的动态| 亚洲av美国av| 伊人久久精品亚洲午夜| 久9热在线精品视频| 3wmmmm亚洲av在线观看| 美女被艹到高潮喷水动态| 丝袜美腿在线中文| 伊人久久大香线蕉亚洲五| svipshipincom国产片| 综合色av麻豆| 久久香蕉精品热| 美女高潮的动态| 99精品欧美一区二区三区四区| 内射极品少妇av片p| tocl精华| 亚洲精品成人久久久久久| 国产探花在线观看一区二区| 亚洲国产精品成人综合色| 中文字幕av成人在线电影| 欧美国产日韩亚洲一区| 天天添夜夜摸| 亚洲成人免费电影在线观看| 女人被狂操c到高潮| 十八禁网站免费在线| 色老头精品视频在线观看| 熟女人妻精品中文字幕| 国产免费av片在线观看野外av| 久久精品国产综合久久久| 男插女下体视频免费在线播放| 国产真实伦视频高清在线观看 | 国产高清videossex| 18禁在线播放成人免费| 不卡一级毛片| 亚洲久久久久久中文字幕| a级毛片a级免费在线| 国产亚洲精品久久久久久毛片| 亚洲欧美激情综合另类| 88av欧美| 国产精品一区二区三区四区免费观看 | 日日摸夜夜添夜夜添小说| 免费搜索国产男女视频| 国产精品影院久久| 国产在线精品亚洲第一网站| 两个人视频免费观看高清| 日韩人妻高清精品专区| 亚洲男人的天堂狠狠| 国产一区二区三区在线臀色熟女| 亚洲电影在线观看av| 欧美成人免费av一区二区三区| 国产亚洲欧美在线一区二区| 国产视频一区二区在线看| 亚洲精品456在线播放app | 亚洲av免费在线观看| 欧美精品啪啪一区二区三区| 国产麻豆成人av免费视频| 国产私拍福利视频在线观看| 国产成+人综合+亚洲专区| 久久久国产成人免费| 黄色女人牲交| 无限看片的www在线观看| 久久久国产成人免费| 久久99热这里只有精品18| 国产成+人综合+亚洲专区| 精品一区二区三区视频在线 | 窝窝影院91人妻| 一边摸一边抽搐一进一小说| 亚洲精品日韩av片在线观看 | 国产真人三级小视频在线观看| 国产aⅴ精品一区二区三区波| 淫秽高清视频在线观看| 婷婷亚洲欧美| 99在线视频只有这里精品首页| 女人高潮潮喷娇喘18禁视频| 成人亚洲精品av一区二区| 宅男免费午夜| a级一级毛片免费在线观看| 18禁美女被吸乳视频| 久久国产精品影院| 国产精品美女特级片免费视频播放器| 美女大奶头视频| 国产aⅴ精品一区二区三区波|