• <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ù)的討論
    超碰成人久久| 亚洲av成人av| 国内久久婷婷六月综合欲色啪| 日本三级黄在线观看| 在线播放国产精品三级| netflix在线观看网站| 成人一区二区视频在线观看| 欧美最黄视频在线播放免费| 亚洲国产中文字幕在线视频| 亚洲欧美精品综合久久99| 国产精品永久免费网站| 真人做人爱边吃奶动态| 久久久久国产一级毛片高清牌| 亚洲无线在线观看| 欧美成人一区二区免费高清观看 | 又黄又粗又硬又大视频| av福利片在线观看| 久久欧美精品欧美久久欧美| 中文字幕人成人乱码亚洲影| 在线免费观看的www视频| 舔av片在线| 老鸭窝网址在线观看| 性色av乱码一区二区三区2| 1024视频免费在线观看| 国产熟女午夜一区二区三区| 欧美色视频一区免费| 1024视频免费在线观看| x7x7x7水蜜桃| 亚洲第一欧美日韩一区二区三区| 级片在线观看| 在线国产一区二区在线| 男女那种视频在线观看| 久久国产精品人妻蜜桃| 在线国产一区二区在线| 久久国产精品影院| 亚洲全国av大片| 成人三级黄色视频| 欧美午夜高清在线| 琪琪午夜伦伦电影理论片6080| 免费看美女性在线毛片视频| 桃红色精品国产亚洲av| 午夜老司机福利片| 国产精品永久免费网站| 一级毛片女人18水好多| 国产精品亚洲一级av第二区| 一二三四社区在线视频社区8| 欧美成狂野欧美在线观看| 观看免费一级毛片| 午夜福利高清视频| 叶爱在线成人免费视频播放| 很黄的视频免费| 在线国产一区二区在线| 国产一区在线观看成人免费| a级毛片在线看网站| 亚洲男人的天堂狠狠| 日韩精品中文字幕看吧| 国产av一区二区精品久久| 亚洲五月天丁香| 国内久久婷婷六月综合欲色啪| av免费在线观看网站| 国内揄拍国产精品人妻在线| 久久九九热精品免费| 亚洲在线自拍视频| 天天躁夜夜躁狠狠躁躁| 听说在线观看完整版免费高清| xxxwww97欧美| 国产精品乱码一区二三区的特点| 色综合亚洲欧美另类图片| 色精品久久人妻99蜜桃| 欧美国产日韩亚洲一区| 精品一区二区三区视频在线观看免费| 19禁男女啪啪无遮挡网站| 久久久国产成人免费| 一级毛片女人18水好多| 精品人妻1区二区| 嫩草影院精品99| 欧美色欧美亚洲另类二区| 国产午夜精品论理片| 色在线成人网| 亚洲精品久久国产高清桃花| 啦啦啦韩国在线观看视频| 国产精品综合久久久久久久免费| 99国产综合亚洲精品| 国产97色在线日韩免费| 一本一本综合久久| 91老司机精品| 国产爱豆传媒在线观看 | 国产午夜精品论理片| 天天躁狠狠躁夜夜躁狠狠躁| 老鸭窝网址在线观看| 高清在线国产一区| 亚洲精品中文字幕在线视频| 黄色视频,在线免费观看| av片东京热男人的天堂| 国产精品国产高清国产av| 男人舔女人下体高潮全视频| 桃红色精品国产亚洲av| 一级毛片高清免费大全| 成人三级黄色视频| www.www免费av| 两个人看的免费小视频| 欧美+亚洲+日韩+国产| 最近视频中文字幕2019在线8| 桃色一区二区三区在线观看| 色老头精品视频在线观看| 一区福利在线观看| 在线观看免费日韩欧美大片| 麻豆一二三区av精品| 日本五十路高清| 91av网站免费观看| 黄色丝袜av网址大全| 亚洲熟女毛片儿| 天堂√8在线中文| 在线观看免费午夜福利视频| 亚洲av成人一区二区三| 亚洲国产精品sss在线观看| 国产成人精品无人区| 欧洲精品卡2卡3卡4卡5卡区| 午夜两性在线视频| 床上黄色一级片| aaaaa片日本免费| 免费在线观看视频国产中文字幕亚洲| 久久婷婷人人爽人人干人人爱| 99riav亚洲国产免费| 91老司机精品| 亚洲美女黄片视频| 一级片免费观看大全| 久久国产精品人妻蜜桃| 欧美日本亚洲视频在线播放| 最近最新免费中文字幕在线| 高潮久久久久久久久久久不卡| 天天躁夜夜躁狠狠躁躁| 99国产精品一区二区蜜桃av| 久久久久久久久免费视频了| 亚洲美女黄片视频| 午夜福利在线观看吧| 久久久久久久久免费视频了| 国产精品久久久久久亚洲av鲁大| 午夜福利在线观看吧| 中文字幕久久专区| 97人妻精品一区二区三区麻豆| a级毛片在线看网站| 久久久精品大字幕| 免费看a级黄色片| 欧美成人性av电影在线观看| 两个人看的免费小视频| 久久人人精品亚洲av| 大型av网站在线播放| 久久人妻av系列| 午夜福利高清视频| 亚洲七黄色美女视频| 黄色丝袜av网址大全| 久久久久久九九精品二区国产 | 国产精品免费视频内射| 国产三级黄色录像| 欧美另类亚洲清纯唯美| 亚洲全国av大片| 日日干狠狠操夜夜爽| 最新美女视频免费是黄的| 精品电影一区二区在线| 麻豆av在线久日| 日韩大尺度精品在线看网址| 国产精品久久久av美女十八| a在线观看视频网站| 成人欧美大片| 男人舔女人下体高潮全视频| 国产麻豆成人av免费视频| 成人永久免费在线观看视频| 免费电影在线观看免费观看| 又爽又黄无遮挡网站| av免费在线观看网站| 99在线人妻在线中文字幕| 成年免费大片在线观看| 夜夜看夜夜爽夜夜摸| 欧美日韩中文字幕国产精品一区二区三区| 久久久久精品国产欧美久久久| 老司机靠b影院| 久久中文字幕人妻熟女| 亚洲在线自拍视频| 中文字幕人成人乱码亚洲影| 99在线视频只有这里精品首页| 欧美乱色亚洲激情| 香蕉av资源在线| 国产麻豆成人av免费视频| 欧美激情久久久久久爽电影| 麻豆国产97在线/欧美 | 午夜福利18| 又大又爽又粗| 欧美乱色亚洲激情| 亚洲精品国产精品久久久不卡| 国产成人欧美在线观看| 日本黄大片高清| av在线天堂中文字幕| www.熟女人妻精品国产| av中文乱码字幕在线| 啪啪无遮挡十八禁网站| 久久久久久久久中文| av欧美777| 1024手机看黄色片| 黄色成人免费大全| 听说在线观看完整版免费高清| 欧美国产日韩亚洲一区| 欧美av亚洲av综合av国产av| 一级作爱视频免费观看| 757午夜福利合集在线观看| 亚洲熟妇中文字幕五十中出| 欧美3d第一页| 午夜福利高清视频| 男插女下体视频免费在线播放| 国产熟女xx| 久久婷婷人人爽人人干人人爱| 麻豆成人午夜福利视频| 成人国产综合亚洲| 九九热线精品视视频播放| 熟女少妇亚洲综合色aaa.| 12—13女人毛片做爰片一| 精品免费久久久久久久清纯| 天堂av国产一区二区熟女人妻 | xxxwww97欧美| 免费在线观看影片大全网站| 中亚洲国语对白在线视频| 国产精品野战在线观看| 色综合亚洲欧美另类图片| x7x7x7水蜜桃| 老汉色∧v一级毛片| 午夜亚洲福利在线播放| 一级毛片精品| 亚洲,欧美精品.| 桃色一区二区三区在线观看| 欧美三级亚洲精品| 亚洲国产精品成人综合色| 亚洲av中文字字幕乱码综合| 在线永久观看黄色视频| 99热6这里只有精品| 国内精品一区二区在线观看| 好男人在线观看高清免费视频| 岛国在线观看网站| 亚洲成人国产一区在线观看| 俺也久久电影网| 中文字幕高清在线视频| 成年版毛片免费区| 人妻久久中文字幕网| 欧美丝袜亚洲另类 | 丁香欧美五月| 哪里可以看免费的av片| 禁无遮挡网站| 免费观看精品视频网站| 一进一出抽搐gif免费好疼| 日韩三级视频一区二区三区| 99热6这里只有精品| 又黄又粗又硬又大视频| 国产一区二区在线观看日韩 | 好看av亚洲va欧美ⅴa在| 97人妻精品一区二区三区麻豆| 欧美zozozo另类| 国产精品 欧美亚洲| 午夜两性在线视频| 床上黄色一级片| 国产亚洲av嫩草精品影院| 久久热在线av| 亚洲最大成人中文| 国产一区在线观看成人免费| 国产精品久久视频播放| 欧美大码av| 国产精品98久久久久久宅男小说| 亚洲精品色激情综合| 女警被强在线播放| 中文字幕av在线有码专区| 99国产精品一区二区三区| 国产成人系列免费观看| 夜夜爽天天搞| 三级毛片av免费| 日韩国内少妇激情av| 国产精品影院久久| 亚洲一码二码三码区别大吗| 欧美国产日韩亚洲一区| 人妻夜夜爽99麻豆av| 国产亚洲欧美在线一区二区| 欧美3d第一页| 又紧又爽又黄一区二区| 丰满人妻一区二区三区视频av | 欧美日韩精品网址| 中文字幕人妻丝袜一区二区| 日韩有码中文字幕| 不卡av一区二区三区| 国产私拍福利视频在线观看| 99热这里只有精品一区 | 精品国产美女av久久久久小说| 久久精品国产99精品国产亚洲性色| 男人舔女人的私密视频| 午夜精品在线福利| 女人爽到高潮嗷嗷叫在线视频| 一a级毛片在线观看| 国产精品爽爽va在线观看网站| 亚洲色图 男人天堂 中文字幕| 9191精品国产免费久久| 97人妻精品一区二区三区麻豆| 国产成人精品久久二区二区91| 级片在线观看| 欧美日本亚洲视频在线播放| 黄频高清免费视频| 精品午夜福利视频在线观看一区| 男男h啪啪无遮挡| 久久欧美精品欧美久久欧美| 一二三四社区在线视频社区8| 国内精品久久久久久久电影| 欧美精品啪啪一区二区三区| 日韩欧美一区二区三区在线观看| 91国产中文字幕| 真人做人爱边吃奶动态| av在线播放免费不卡| 91国产中文字幕| netflix在线观看网站| 亚洲美女视频黄频| 久久精品成人免费网站| 日日夜夜操网爽| 国产成人系列免费观看| 1024视频免费在线观看| 国产精品电影一区二区三区| 麻豆久久精品国产亚洲av| 欧美日韩亚洲综合一区二区三区_| 国产精品电影一区二区三区| 91麻豆精品激情在线观看国产| 成人三级做爰电影| 18禁黄网站禁片免费观看直播| 国产视频内射| 色在线成人网| 少妇的丰满在线观看| 波多野结衣高清无吗| 久久久国产成人免费| 国产午夜福利久久久久久| 国模一区二区三区四区视频 | 五月伊人婷婷丁香| 性欧美人与动物交配| 国产精品久久久久久亚洲av鲁大| 精品熟女少妇八av免费久了| 日韩免费av在线播放| 青草久久国产| 两个人的视频大全免费| 麻豆成人午夜福利视频| 天天躁夜夜躁狠狠躁躁| 免费在线观看亚洲国产| 亚洲精品一区av在线观看| 国产真实乱freesex| 麻豆久久精品国产亚洲av| 两性夫妻黄色片| 亚洲 欧美 日韩 在线 免费| 国产亚洲av嫩草精品影院| 夜夜看夜夜爽夜夜摸| 在线观看午夜福利视频| 精品国产亚洲在线| 高清毛片免费观看视频网站| 全区人妻精品视频| 后天国语完整版免费观看| 亚洲电影在线观看av| 桃红色精品国产亚洲av| 夜夜看夜夜爽夜夜摸| www.www免费av| 小说图片视频综合网站| 给我免费播放毛片高清在线观看| 黄色 视频免费看| tocl精华| 亚洲全国av大片| 女警被强在线播放| 亚洲精品一区av在线观看| 在线看三级毛片| 亚洲真实伦在线观看| 国产成人一区二区三区免费视频网站| av天堂在线播放| 在线观看日韩欧美| 18禁观看日本| 99久久精品国产亚洲精品| 女人高潮潮喷娇喘18禁视频| 国产69精品久久久久777片 | 国产成人精品久久二区二区91| 亚洲最大成人中文| 69av精品久久久久久| 亚洲最大成人中文| 亚洲精华国产精华精| 国内精品一区二区在线观看| 中文字幕人妻丝袜一区二区| 伦理电影免费视频| 91麻豆精品激情在线观看国产| 欧美大码av| 久久久久久久久久黄片| 日韩大尺度精品在线看网址| 视频区欧美日本亚洲| 成人欧美大片| 黄色视频,在线免费观看| 成人永久免费在线观看视频| 69av精品久久久久久| 精品福利观看| 日本一二三区视频观看| 久久婷婷成人综合色麻豆| 好男人在线观看高清免费视频| 12—13女人毛片做爰片一| 欧美日韩福利视频一区二区| 人人妻人人澡欧美一区二区| 可以免费在线观看a视频的电影网站| 搡老岳熟女国产| 久久香蕉精品热| 午夜福利视频1000在线观看| 麻豆久久精品国产亚洲av| 亚洲人成电影免费在线| 国产三级在线视频| 精品欧美一区二区三区在线| 日韩欧美一区二区三区在线观看| 国产成人av激情在线播放| 亚洲熟妇中文字幕五十中出| 中文字幕高清在线视频| 精品免费久久久久久久清纯| av福利片在线观看| 国产精华一区二区三区| 午夜视频精品福利| 国产熟女午夜一区二区三区| 老熟妇仑乱视频hdxx| 91国产中文字幕| 亚洲精品中文字幕在线视频| 给我免费播放毛片高清在线观看| 黄色成人免费大全| 国产成人影院久久av| 亚洲国产高清在线一区二区三| 51午夜福利影视在线观看| 日韩欧美在线二视频| 免费电影在线观看免费观看| 99热这里只有是精品50| 天堂√8在线中文| 一区二区三区国产精品乱码| 国产黄色小视频在线观看| 曰老女人黄片| 成年人黄色毛片网站| 日韩免费av在线播放| xxxwww97欧美| 欧美黑人精品巨大| 久久久久久久久免费视频了| 视频区欧美日本亚洲| 亚洲专区字幕在线| 母亲3免费完整高清在线观看| 国产精品亚洲一级av第二区| 在线观看免费午夜福利视频| 欧美av亚洲av综合av国产av| 亚洲欧洲精品一区二区精品久久久| 黄片小视频在线播放| 国产99久久九九免费精品| 中文字幕熟女人妻在线| 久久午夜综合久久蜜桃| 精品高清国产在线一区| 欧美三级亚洲精品| 欧美日韩亚洲综合一区二区三区_| 日日夜夜操网爽| 黄色视频不卡| 观看免费一级毛片| 国产精品av视频在线免费观看| 亚洲成人久久爱视频| 国产视频一区二区在线看| 亚洲七黄色美女视频| 久久久久免费精品人妻一区二区| 99在线人妻在线中文字幕| 一本精品99久久精品77| 一a级毛片在线观看| 亚洲第一欧美日韩一区二区三区| 18禁观看日本| 欧美人与性动交α欧美精品济南到| 午夜激情av网站| 一个人观看的视频www高清免费观看 | 国产成人啪精品午夜网站| 国产亚洲精品第一综合不卡| 国产成人影院久久av| 亚洲av片天天在线观看| 亚洲,欧美精品.| 中出人妻视频一区二区| 久久国产乱子伦精品免费另类| 中文字幕精品亚洲无线码一区| 又黄又爽又免费观看的视频| 天天添夜夜摸| 成人18禁在线播放| 伊人久久大香线蕉亚洲五| 久久午夜亚洲精品久久| 脱女人内裤的视频| 国产欧美日韩一区二区三| e午夜精品久久久久久久| 我的老师免费观看完整版| 欧美黑人巨大hd| 超碰成人久久| 久久中文字幕人妻熟女| 亚洲av片天天在线观看| 久久香蕉国产精品| 亚洲欧洲精品一区二区精品久久久| 国产又色又爽无遮挡免费看| 97超级碰碰碰精品色视频在线观看| 国产高清视频在线播放一区| 国产成人影院久久av| 亚洲国产精品合色在线| 午夜福利免费观看在线| 精品一区二区三区av网在线观看| 五月伊人婷婷丁香| 国产97色在线日韩免费| 女人高潮潮喷娇喘18禁视频| 中文字幕精品亚洲无线码一区| 老熟妇仑乱视频hdxx| 好男人在线观看高清免费视频| 777久久人妻少妇嫩草av网站| 搡老妇女老女人老熟妇| 国产伦一二天堂av在线观看| a级毛片a级免费在线| 成人欧美大片| 午夜a级毛片| 成人亚洲精品av一区二区| 久久国产乱子伦精品免费另类| 亚洲人成伊人成综合网2020| 男人舔女人的私密视频| 日韩有码中文字幕| 身体一侧抽搐| 最近在线观看免费完整版| 亚洲精品中文字幕在线视频| 男插女下体视频免费在线播放| 黑人巨大精品欧美一区二区mp4| 1024视频免费在线观看| 欧美精品亚洲一区二区| 国产精品爽爽va在线观看网站| 脱女人内裤的视频| 香蕉av资源在线| 国产激情久久老熟女| 亚洲自拍偷在线| 久久久精品国产亚洲av高清涩受| av中文乱码字幕在线| 国产真实乱freesex| 成人国产一区最新在线观看| 真人一进一出gif抽搐免费| 老司机靠b影院| 国产真实乱freesex| 男人舔女人的私密视频| 在线观看免费午夜福利视频| 日本免费一区二区三区高清不卡| 成人手机av| 国产精品免费视频内射| 免费人成视频x8x8入口观看| 香蕉久久夜色| 香蕉丝袜av| 少妇熟女aⅴ在线视频| 搡老妇女老女人老熟妇| 两性午夜刺激爽爽歪歪视频在线观看 | 制服丝袜大香蕉在线| 中文字幕av在线有码专区| 国产欧美日韩精品亚洲av| 国产精品电影一区二区三区| 亚洲精品国产精品久久久不卡| 50天的宝宝边吃奶边哭怎么回事| www日本黄色视频网| 亚洲av电影不卡..在线观看| 男人舔女人下体高潮全视频| 成人亚洲精品av一区二区| 国内揄拍国产精品人妻在线| 99热这里只有精品一区 | 亚洲精品中文字幕在线视频| 成人av在线播放网站| 男人舔奶头视频| 欧美3d第一页| 亚洲片人在线观看| 国产av一区在线观看免费| 亚洲男人天堂网一区| 夜夜躁狠狠躁天天躁| 亚洲精品粉嫩美女一区| 国产成人精品久久二区二区免费| www日本在线高清视频| www.999成人在线观看| 午夜精品一区二区三区免费看| 亚洲中文av在线| 好看av亚洲va欧美ⅴa在| 国产成人精品久久二区二区91| 久久国产精品影院| 亚洲精华国产精华精| 精品一区二区三区视频在线观看免费| 麻豆成人午夜福利视频| 神马国产精品三级电影在线观看 | 婷婷精品国产亚洲av在线| 女人被狂操c到高潮| 日韩精品免费视频一区二区三区| 亚洲精品国产一区二区精华液| 日韩欧美一区二区三区在线观看| 国产aⅴ精品一区二区三区波| 黄色视频不卡| 亚洲午夜理论影院| 国产黄色小视频在线观看| 国产精品99久久99久久久不卡| 亚洲色图 男人天堂 中文字幕| 91字幕亚洲| 国产v大片淫在线免费观看| 日韩欧美一区二区三区在线观看| 成人一区二区视频在线观看| 久久中文字幕一级| 国语自产精品视频在线第100页| a在线观看视频网站| 天天躁夜夜躁狠狠躁躁| 国产成人欧美在线观看| 巨乳人妻的诱惑在线观看| 黄色丝袜av网址大全| 19禁男女啪啪无遮挡网站| 久久精品亚洲精品国产色婷小说| 久久亚洲真实| 黄色 视频免费看| 夜夜夜夜夜久久久久| 午夜福利欧美成人| 九色国产91popny在线| 欧美日韩中文字幕国产精品一区二区三区| 亚洲专区字幕在线| 最近视频中文字幕2019在线8| 国内精品久久久久精免费| 天天添夜夜摸| 777久久人妻少妇嫩草av网站| 国产亚洲av嫩草精品影院| 97碰自拍视频| 亚洲av电影不卡..在线观看| 亚洲男人天堂网一区| 天堂av国产一区二区熟女人妻 |