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

    灰色系統(tǒng)模型在UT1–UTC超短期預報中的應用?

    2016-06-27 08:14:07趙丹寧蔡宏兵
    天文學報 2016年3期
    關鍵詞:實線跨度灰色

    雷 雨 趙丹寧 蔡宏兵

    (1中國科學院國家授時中心西安710600)

    (2中國科學院時間頻率基準重點實驗室西安710600)

    (3中國科學院大學北京100049)

    灰色系統(tǒng)模型在UT1–UTC超短期預報中的應用?

    雷 雨1,2,3?趙丹寧1,3蔡宏兵1,2

    (1中國科學院國家授時中心西安710600)

    (2中國科學院時間頻率基準重點實驗室西安710600)

    (3中國科學院大學北京100049)

    依據(jù)灰色系統(tǒng)理論和UT1–UTC的變化規(guī)律,以較少的觀測樣本建立了預報UT1–UTC的灰色系統(tǒng)模型,并將其與人工神經(jīng)網(wǎng)絡(arti ficial neural network,ANN)、最小二乘(least squares,LS)與自回歸(autoregressive,AR)模型的組合(LS+AR)方法以及地球定向參數(shù)預報比較競賽(Earth Orientation Parameters Prediction Comparison Campaign,EOP PCC)的預報結(jié)果進行對比.結(jié)果表明:灰色系統(tǒng)模型用于UT1–UTC預報是高效可行的,尤其是在1–10 d跨度的超短期預報中預報效果顯著.

    天體測量,時間,方法:其他諸多方面

    1 引言

    地球自轉(zhuǎn)運動表征了地球整體的轉(zhuǎn)動運動狀態(tài),可以用地球定向參數(shù)(Earth orientation parameter,EOP)來描述,包括歲差、章動、極移和UT1–UTC幾個分量. EOP是實現(xiàn)地球參考系與天球參考系相互轉(zhuǎn)換的必需參數(shù),在現(xiàn)代空間導航、深空探測以及天文地球動力學研究等領域都有重要應用.現(xiàn)代空間測地技術如全球衛(wèi)星導航系統(tǒng)(Global Navigation Satellite System,GNSS)、衛(wèi)星激光測距(Satellite Laser Ranging, SLR)和甚長干涉基線測量(Very Long Baseline Interferometry,VLBI)等是獲取EOP的主要手段,但復雜的資料處理過程使得EOP的獲取存在一定的延遲,所以對EOP進行高精度預報是一項非常重要的工作.

    在EOP的幾個分量中,UT1–UTC是變化最快、最難預報的一個分量,因此,UT1–UTC的高精度預報是EOP預報的難點也是當前的研究熱點.目前預報UT1–UTC的方法有多種,如最小二乘外推[1?2]、LS外推和自回歸模型的組合[2?5]、人工神經(jīng)網(wǎng)絡[6?8]、高斯過程(Gaussian processes,GP)[9?10]以及小波分解(wavelet decomposition,WD)和自協(xié)方差(autocovariance,AC)模型的組合(WD+AC)[11]等.為了比較不同方法的預報效果,維也納理工大學大地測量與地球物理研究所從2005年10月1日至2008年2月28日組織了全球性的EOP預報比較競賽,競賽結(jié)果表明,沒有一種方法既適合EOP所有分量的預報又適合所有跨度的預報[12].

    灰色系統(tǒng)是我國鄧聚龍教授于1982年提出的一種研究少數(shù)據(jù)、貧信息不確定性問題的新方法,它是既含已知信息又含未知信息或非確知信息的系統(tǒng),已知的信息為“白色”,未知的信息為“黑色”,不充分、不完全的信息為“灰色”[13].灰色系統(tǒng)通過對原始數(shù)據(jù)實行累加或累減使之成為具有較強規(guī)律的新數(shù)列,然后對生成的新數(shù)列進行建模.它不需要較大的數(shù)據(jù)量,只要原始序列有4個以上數(shù)據(jù)就可以建立灰色系統(tǒng)模型,提供了貧信息情況下解決系統(tǒng)問題的新途徑.

    目前UT1–UTC預報大多是采用國際地球自轉(zhuǎn)與參考系服務(International Earth Rotation and Reference Systems Service,IERS)發(fā)布的EOP C04序列作為基礎序列,在常規(guī)情況下數(shù)據(jù)量是充足的,但在非常時期一旦IERS暫停其EOP服務,就存在數(shù)據(jù)量不足的情形,傳統(tǒng)LS等模型就不再適用,而灰色模型卻為小樣本情況下的UT1–UTC預報提供了一種可行方法;另一方面,近年來國內(nèi)相關單位致力于提供自主的EOP服務,正在建設以數(shù)字天頂筒為主要觀測設備的EOP自主測量與服務系統(tǒng),然而受限于光學測量技術的限制,獲得的EOP測量數(shù)據(jù)是不連續(xù)的,在這種條件下開展EOP預報服務,必須尋求適合于貧信息情況下的時間序列預報方法,而灰色模型為貧信息情況下的UT1–UTC預報提供了一種可行途徑.

    由于受眾多激發(fā)源的影響,UT1–UTC的變化極為復雜,不僅包括年際變化、季節(jié)性變化以及幾十天的短周期變化,同時還包括不規(guī)則的高頻變化,從而很難掌握其復雜細致的變化規(guī)律,但這符合灰色系統(tǒng)理論的特點.因此,本文嘗試將UT1–UTC的變化過程視為灰色系統(tǒng),將灰色系統(tǒng)模型引入UT1–UTC超短期預報研究中,研究表明:灰色系統(tǒng)模型用于UT1–UTC超短期預報是可行的,且預報精度和預報效率較高.

    2 灰色系統(tǒng)模型

    在灰色系統(tǒng)理論中,GM(1,1)模型是外推法的一種新的開拓,是最常用的一種灰色系統(tǒng)模型,是由一個僅包含單變量的1階微分方程構(gòu)成的模型,適合對本身數(shù)據(jù)的預報.

    設有變量X(0)的原始數(shù)據(jù)序列:

    其中,n為原始數(shù)據(jù)個數(shù),與之對應的時間為t(t=1,2,···,n).對原始序列進行累加生成1階累加模塊X(1)

    由模塊X(1)構(gòu)成的微分方程為

    其中,a為發(fā)展系數(shù),b為灰色作用量.根據(jù)導數(shù)離散化形式,(3)式可以寫成矩陣形式:

    利用最小二乘原理,可求得參數(shù)估值:

    求得參數(shù)估值后,可得到微分方程的解:

    離散形式為

    其中,k為參與預報的原始數(shù)據(jù)個數(shù).

    一般形式為

    其中,p為預報點數(shù),則預報后的原始觀測序列為

    在建模前,必須檢查原始數(shù)據(jù)序列符號是否一致,若不一致,則給每個元素均加上一個常數(shù),使序列的符號一致,然后在此基礎上建立GM(1,1)模型并進行預報,最后再從預報值中減去所加常數(shù)即可獲得所需的預報值.所加常數(shù)的取值原則是,符號與原始序列中絕對值最大的數(shù)的符號一致,且其絕對值大于原始序列中絕對值最大數(shù)的絕對值.

    3 實驗分析

    本文實驗數(shù)據(jù)來自IERS發(fā)布的EOP C04序列,數(shù)據(jù)采樣間隔為1 d.在對UT1–UTC建模之前,首先去掉UT1–UTC序列中的閏秒得到UT1–TAI(International Atomic Time,國際原子時)序列,然后根據(jù)IERS協(xié)議扣除UT1–TAI序列中周期為5 d–18.6 yr的固體地球帶諧潮汐項獲得UT1R–TAI序列,其中UT1R表示經(jīng)過固體地球帶諧潮修正后的UT1,最后再利用GM(1,1)模型對UT1R–TAI序列建模預報.圖1(a)–1(d)分別繪出了1990年1月1日至1997年12月31日期間的UT1–UTC原始序列、UT1–TAI序列、UT1R–TAI序列及其累加序列.

    圖1 (a)UT1–UTC原始序列;(b)UT1–TAI原始序列;(c)UT1R–TAI原始序列;(d)UT1R–TAI的累加序列Fig.1 (a)The original series of UT1–UTC;(b)the original series of UT1–TAI;(c)the original series of UT1R–TAI;(d)the accumulated series of UT1R–TAI

    由圖1可以看出,扣除閏秒和潮汐項后的UT1–UTC序列呈不規(guī)則變化,其變化規(guī)律不明顯,但一次累加后的新序列比較光滑,比較符合指數(shù)變化規(guī)律,所以從理論上講,應用GM(1,1)模型對UT1–UTC進行預報是完全可行的.我們首先分析了基礎序列長度對GM(1,1)模型預報效果的影響,結(jié)果發(fā)現(xiàn)基礎序列長度越短,預報效果越好,限于篇幅,圖2僅繪出了基礎序列長度為4、8、12和16時GM(1,1)模型的預報均方根誤差(root mean square error,RMSE),其計算公式為

    圖2中,圓形實線表示基礎序列長度n=4時,GM(1,1)模型的預報RMSE,三角形實線表示基礎序列長度n=8時,GM(1,1)模型的預報RMSE,圓形虛線表示基礎序列長度n=12時,GM(1,1)模型的預報RMSE,三角形虛線表示基礎序列長度n=16時,GM(1,1)模型的預報RMSE.對于每一種建模方案,我們均進行了700期的1–30 d跨度的預報(N=700),預報時間段為1998年1月1日至1999年12月31日,每隔1 d預報1次.從圖2可以看到,基礎數(shù)據(jù)序列長度為4時預報效果最佳,所以在接下來的預報實驗中我們僅采用4個歷元的數(shù)據(jù)建立GM(1,1)模型.

    圖2 UT1-UTC不同基礎序列長度預報精度Fig.2 The comparison of the prediction accuracies for di ff erent data intervals

    由于ANN和LS+AR模型是預報UT1–UTC精度較高的兩種方法,所以為了驗證GM(1,1)模型的預報效果,我們將GM(1,1)模型的預報結(jié)果分別與Schuh等[6]所使用的ANN模型與Xu和Zhou[3]所使用的LS+AR模型的預報結(jié)果進行比較,RMSE比較結(jié)果見圖3–4和表1–2.圖3–4中,圓形實線表示GM(1,1)模型的預報結(jié)果,三角形實線和方形實線分別表示Schuh等[6]與Xu和Zhou[3]的預報結(jié)果.由表1–2可以看出,GM(1,1)模型用于UT1–UTC預報是比較成功的,隨著預報跨度的增大,預報精度有所降低.從表1和圖3的對比中可以看到,對于1–6 d的超短期預報,GM(1,1)模型的預報精度與Schuh等[6]的預報精度大致相當,從第7 d開始,GM(1,1)模型的預報精度低于Schuh等[6]的預報精度并逐漸拉開差距,特別是從第10天開始,前者的預報精度明顯低于后者.從表2和圖4可以看到,對于提前1 d預報,GM(1,1)模型的預報精度明顯優(yōu)于LS+AR模型,對于10 d跨度預報,兩者的預報精度相當,而對于20 d和30 d跨度預報,前者的預報精度則不及后者.

    圖3 GM(1,1)模型預報精度與Schuh等[6]預報精度比較Fig.3 The comparison of the prediction accuracies of GM(1,1)and Schuh et al.[6]

    圖4 GM(1,1)模型預報精度與Xu和Zhou[3]預報精度比較Fig.4 The comparison of the prediction accuracies of GM(1,1)and Xu and Zhou[3]

    表1 GM(1,1)模型預報精度與ANN預報精度(Schuh等[6]的預報結(jié)果)的比較Table 1 The comparison of the prediction accuracies of GM(1,1)and ANN(Schuh et al.[6])

    表2 GM(1,1)模型預報精度與LS+AR預報精度(Xu和Zhou[3]的預報結(jié)果)的比較Table 2 The comparison of the prediction accuracies of GM(1,1)and LS+AR(Xu and Zhou[3])

    為了進一步探討GM(1,1)模型用于UT1–UTC預報的性能,本文也預報了2005年10月1日至2008年2月28日(與EOP PCC時間段相同)不同跨度的UT1–UTC值,統(tǒng)計了預報結(jié)果的平均絕對誤差(mean absolute error,MAE),并同EOP PCC結(jié)果進行了對比,對比結(jié)果見圖5–6,其中MAE的計算公式為

    圖5–6中不同形狀曲線分別為參與EOP PCC不同小組所獲得的預報結(jié)果,其中3位數(shù)字代表各小組編號,EOP PCC中UT1–UTC預報精度較高的兩個小組分別是Gross小組和Gambis小組,這兩個小組均聯(lián)合大氣角動量數(shù)據(jù)進行UT1–UTC預報.圖5–6中圓形點線、方形點線和菱形點線(編號091–093)表示Zotov小組的預報結(jié)果,三角形實線(編號071)表示EOP產(chǎn)品中心的預報結(jié)果,三角形虛線(編號031)表示Gross小組的預報結(jié)果,圓形虛線表示本文預報結(jié)果.圖5中五角星實線、方形實線、菱形實線和圓形實線(編號072–075)表示Gambis小組的預報結(jié)果.圖6中五角星實線(編號053)表示Kosek小組的預報結(jié)果.

    從圖5的比較可以看到,對于1–5 d的超短期預報,GM(1,1)模型的預報精度僅次于排在第1位Gross小組的精度,與排在第2位的Gambis小組的精度相當,從第6天開始, GM(1,1)模型的預報精度低于Gross小組、Gambis小組的精度并逐漸拉開差距,但仍優(yōu)于其他小組的預報精度或與之相當.從圖6可以看出,對于1–30 d的短期預報,當預報跨度在10 d以內(nèi)時,除Gross小組的預報精度超過GM(1,1)模型外,其他小組的預報精度與GM(1,1)模型精度大致相當,當預報跨度超過10 d時,GM(1,1)模型的預報精度則開始變差,不如大部分小組的預報精度.

    圖5 GM(1,1)模型超短期(1–10 d)預報精度與EOP PCC結(jié)果的對比Fig.5The comparison of the prediction accuracies of GM(1,1)and EOP PCC for the ultra short-term (1–10 d)

    圖6 GM(1,1)模型短期(1–30 d)預報精度與EOP PCC結(jié)果的對比Fig.6 The comparison of the prediction accuracies of GM(1,1)and EOP PCC for the short-term (1–30 d)

    4 總結(jié)與展望

    本文通過實例驗證了GM(1,1)模型用于UT1–UTC預報的可行性和有效性.通過與Schuh等[6]所使用的ANN模型以及與Xu和Zhou[3]所使用的LS+AR模型的預報結(jié)果的比較可以發(fā)現(xiàn),雖然在20–30 d跨度的預報中GM(1,1)模型的預報效果遠不及ANN和LS+AR模型,但對于1–10 d跨度預報GM(1,1)模型的預報精度卻與之相當或稍遜.通過與EOP PCC預報結(jié)果的進一步比較,我們同樣發(fā)現(xiàn),GM(1,1)模型在1–10 d超短期跨度預報中效果良好,但當預報跨度超過10 d時效果不甚理想.上述分析表明, UT1–UTC累加序列的短期變化比較符合指數(shù)變化規(guī)律,但長期變化與指數(shù)變化規(guī)律不符,所以利用GM(1,1)模型來描述UT1–UTC的短期變化特性是比較合適的.

    算例表明,GM(1,1)模型的預報精度并非最優(yōu),但與其他方法相比具有可比性,尤其是在1–10 d的超短期預報中效果顯著,這為UT1–UTC超短期預報提供了一種備選方法.此外,對比現(xiàn)有的UT1–UTC預報方法,GM(1,1)模型有2個優(yōu)點:一是它僅僅使用已知的少數(shù)幾個歷元的數(shù)據(jù)便可以建模,減少了要使用的數(shù)據(jù)量,為小樣本情況下的UT1–UTC預報提供了一種有效途徑,這對于恢復我國自主的UT1–UTC服務能力有著重要的現(xiàn)實意義;二是GM(1,1)模型簡單、易于實現(xiàn),無需復雜的調(diào)整參數(shù)就能獲得滿意的預報效果,極大提高了建模效率,這對于UT1–UTC的實時快速預報有著重要的現(xiàn)實意義.由于UT1–UTC與大氣角動量軸向分量有著密切的關系[1?2],下一步的研究我們將會聯(lián)合大氣角動量軸向分量和UT1–UTC序列建立多變量灰色系統(tǒng)模型進行UT1–UTC預報.

    [1]Guo J Y,Li Y B,Dai C L,et al.JGeo,2013,70:36

    [2]Niedzielski T,Kosek W.JGeod,2008,82:83

    [3]Xu X Q,Zhou Y H.AdSpR,2015,56:2248

    [4]Sun Z Z,Xu T H.Geodesy&Geodynamics,2012,3:57

    [5]許雪晴,周永宏.飛行器測控學報,2010,29:70

    [6]Schuh H,Ulrich M,Egger D,et al.JGeod,2002,76:247

    [7]Liao D C,Wang Q J,Zhou Y H,et al.JGeo,2012,62:87

    [8]Lei Y,Zhao D N,Cai H B.Geodesy&Geodynamics,2015,6:151

    [9]雷雨,趙丹寧,高玉平,等.天文學報,2015,56:53

    [10]Lei Y,Zhao D N,Gao Y P,et al.ChA&A,2015,39:368

    [11]Kosek W,Kalarus M,Johnson T J,et al.ArtSa,2005,40:119

    [12]Kalarus M,Schuh H,Kosek W,et al.JGeod,2010,84:587

    [13]鄧聚龍.灰理論基礎.武漢:華中科技大學出版社,2002

    An Application in the Ultra Short-term Prediction of UT1–UTC Based on Grey System Model

    LEI Yu1,2,3ZHAO Dan-ning1,3CAI Hong-bing1,2
    (1 National Time Service Center,Chinese Academy of Sciences,Xi’an 710600)
    (2 Key Laboratory for Time and Frequency Primary Standards,Chinese Academy of Sciences,Xi’an
    710600)
    (3 University of Chinese Academy of Sciences,Beijing 100049)

    This work presents an application of the grey system model in the prediction of UT1–UTC.The short-term prediction of UT1–UTC is studied up to 30 days by means of the grey system model.The EOP(Earth orientation parameter)C04 time series with daily values from the International Earth Rotation and Reference Systems Service(IERS)serve as the data base.The results of the prediction are analyzed and compared with those obtained by the arti ficial neural network(ANN),the combination of least squares(LS)and autoregressive(AR)model(LS+AR),and the Earth Orientation Parameters Prediction Comparison Campaign(EOP PCC).The accuracies of the ultra short-term(1–10 d)prediction are comparable to those obtained by the other prediction methods.The presented method is easy to use.

    astrometry,time,methods:miscellaneous

    P127;

    :A

    10.15940/j.cnki.0001-5245.2016.03.006

    2015-09-28收到原稿,2015-12-02收到修改稿

    ?國家自然科學基金項目(11503031)資助

    ?leiyu@ntsc.ac.cn

    猜你喜歡
    實線跨度灰色
    關于調(diào)整上海道路非必要超長實線及高速監(jiān)控探頭強光燈建議
    緩粘結(jié)預應力技術在大跨度梁中的應用
    大跨度連續(xù)剛構(gòu)橋線形控制分析
    淺灰色的小豬
    秋天來啦
    娃娃畫報(2019年10期)2019-12-17 08:02:09
    戒煙
    詩潮(2019年8期)2019-08-23 05:39:48
    疊疊看 真神奇
    啟蒙(3-7歲)(2019年3期)2019-04-03 01:39:28
    組合鋁合金立柱在超大跨度玻璃幕墻中的應用
    上海建材(2018年4期)2018-11-13 01:08:54
    灰色時代
    Coco薇(2017年2期)2017-04-25 17:59:38
    她、它的灰色時髦觀
    Coco薇(2017年2期)2017-04-25 17:57:49
    免费观看人在逋| 欧美日韩国产亚洲二区| 欧美成人精品欧美一级黄| 国产人妻一区二区三区在| 午夜精品国产一区二区电影 | 黑人高潮一二区| 成人午夜高清在线视频| 99热这里只有是精品在线观看| 亚洲精品一区蜜桃| 精品酒店卫生间| 久久欧美精品欧美久久欧美| 久久精品影院6| 亚洲最大成人av| 免费看光身美女| 久久精品夜色国产| 亚洲人与动物交配视频| 少妇的逼好多水| 26uuu在线亚洲综合色| 婷婷色综合大香蕉| 国产成年人精品一区二区| 噜噜噜噜噜久久久久久91| 亚洲国产精品成人综合色| 久久午夜福利片| 午夜福利在线观看吧| 国产精品一区二区三区四区久久| 国产单亲对白刺激| 好男人在线观看高清免费视频| 一个人观看的视频www高清免费观看| 国语对白做爰xxxⅹ性视频网站| 91久久精品国产一区二区三区| 国产精品久久久久久久久免| 两个人视频免费观看高清| 亚洲欧美成人精品一区二区| 久久精品国产亚洲网站| 99九九线精品视频在线观看视频| 亚洲精品国产成人久久av| 国产亚洲精品久久久com| 亚洲天堂国产精品一区在线| 一级黄色大片毛片| 国产三级在线视频| 纵有疾风起免费观看全集完整版 | 三级国产精品片| 成人av在线播放网站| 男人狂女人下面高潮的视频| 尤物成人国产欧美一区二区三区| 大又大粗又爽又黄少妇毛片口| 成人午夜精彩视频在线观看| 99久久九九国产精品国产免费| av在线蜜桃| 亚洲成人中文字幕在线播放| 长腿黑丝高跟| 丝袜美腿在线中文| 能在线免费观看的黄片| 蜜桃久久精品国产亚洲av| 国产伦精品一区二区三区四那| .国产精品久久| 亚洲精品日韩在线中文字幕| 99久久中文字幕三级久久日本| 亚洲av电影在线观看一区二区三区 | 国产一区二区在线av高清观看| 亚洲av电影不卡..在线观看| 亚洲精品国产av成人精品| 最近中文字幕高清免费大全6| 国产亚洲5aaaaa淫片| 成人性生交大片免费视频hd| 国产亚洲最大av| 国产欧美日韩精品一区二区| 久久人妻av系列| 成人亚洲精品av一区二区| 国产免费又黄又爽又色| 欧美激情国产日韩精品一区| 国产精品乱码一区二三区的特点| 成人毛片60女人毛片免费| 少妇被粗大猛烈的视频| 99在线人妻在线中文字幕| 一区二区三区四区激情视频| 成人av在线播放网站| 亚洲精品乱码久久久v下载方式| av专区在线播放| 欧美最新免费一区二区三区| 国产精品无大码| 亚洲内射少妇av| 国模一区二区三区四区视频| 日韩一区二区三区影片| 亚洲人成网站高清观看| 三级经典国产精品| 亚洲精品自拍成人| 国产免费福利视频在线观看| 在现免费观看毛片| 女人十人毛片免费观看3o分钟| 69av精品久久久久久| 欧美xxxx黑人xx丫x性爽| 91精品伊人久久大香线蕉| 精品一区二区三区视频在线| 亚洲国产色片| 国产精品嫩草影院av在线观看| 青春草亚洲视频在线观看| kizo精华| 久久亚洲国产成人精品v| 精品国产一区二区三区久久久樱花 | 久久这里只有精品中国| 日韩精品有码人妻一区| 国产精品国产三级国产av玫瑰| 国产一区亚洲一区在线观看| 国产成人a区在线观看| 在线播放国产精品三级| 尤物成人国产欧美一区二区三区| 狂野欧美激情性xxxx在线观看| 国内少妇人妻偷人精品xxx网站| 欧美成人午夜免费资源| 熟女人妻精品中文字幕| 干丝袜人妻中文字幕| 真实男女啪啪啪动态图| 日日摸夜夜添夜夜添av毛片| 纵有疾风起免费观看全集完整版 | av天堂中文字幕网| 看十八女毛片水多多多| 久久久久久久亚洲中文字幕| 国产欧美另类精品又又久久亚洲欧美| 亚洲欧洲国产日韩| 亚洲成人中文字幕在线播放| 欧美bdsm另类| 亚洲国产欧美在线一区| 六月丁香七月| 麻豆av噜噜一区二区三区| 老师上课跳d突然被开到最大视频| 亚洲欧美日韩东京热| 日本五十路高清| 天堂√8在线中文| 国产亚洲最大av| 超碰av人人做人人爽久久| 国产黄片视频在线免费观看| 美女高潮的动态| 欧美激情久久久久久爽电影| 亚洲一级一片aⅴ在线观看| 99在线视频只有这里精品首页| 赤兔流量卡办理| 欧美日韩在线观看h| 卡戴珊不雅视频在线播放| 久久人妻av系列| 日本免费一区二区三区高清不卡| 日日摸夜夜添夜夜爱| 精品久久久久久久久久久久久| 天堂中文最新版在线下载 | 国产黄片美女视频| 国产av不卡久久| 韩国高清视频一区二区三区| 美女脱内裤让男人舔精品视频| 亚洲国产欧洲综合997久久,| 欧美xxxx黑人xx丫x性爽| 亚洲精品456在线播放app| videos熟女内射| 日韩在线高清观看一区二区三区| 搞女人的毛片| 18+在线观看网站| 国内精品美女久久久久久| 九色成人免费人妻av| 国产乱来视频区| 国产精品一区二区三区四区免费观看| 亚洲欧洲日产国产| 亚洲国产精品国产精品| 午夜日本视频在线| 国产免费视频播放在线视频 | 三级经典国产精品| 亚洲av免费高清在线观看| 久久久久久久午夜电影| 少妇裸体淫交视频免费看高清| 免费观看的影片在线观看| 亚洲图色成人| 精品午夜福利在线看| 日日撸夜夜添| 黑人高潮一二区| 国产精品久久久久久av不卡| 三级毛片av免费| 亚洲成av人片在线播放无| 两性午夜刺激爽爽歪歪视频在线观看| 欧美三级亚洲精品| 自拍偷自拍亚洲精品老妇| 亚洲国产欧洲综合997久久,| 精品久久久久久久末码| 国产在线一区二区三区精 | 在线观看av片永久免费下载| 亚洲精品色激情综合| 美女cb高潮喷水在线观看| 99热这里只有是精品50| 欧美日韩一区二区视频在线观看视频在线 | 国产精品精品国产色婷婷| 午夜爱爱视频在线播放| 亚洲av中文字字幕乱码综合| 精品国产露脸久久av麻豆 | 亚洲精品乱码久久久久久按摩| h日本视频在线播放| 日韩av不卡免费在线播放| 亚洲欧洲国产日韩| 亚洲高清免费不卡视频| 色尼玛亚洲综合影院| 亚洲成色77777| 亚洲激情五月婷婷啪啪| 日日撸夜夜添| 最近中文字幕高清免费大全6| 91精品国产九色| 特级一级黄色大片| 欧美日韩精品成人综合77777| 狠狠狠狠99中文字幕| 国产精品蜜桃在线观看| 天天一区二区日本电影三级| 亚洲自拍偷在线| 偷拍熟女少妇极品色| 久久久国产成人精品二区| 亚洲av中文字字幕乱码综合| 2022亚洲国产成人精品| or卡值多少钱| 建设人人有责人人尽责人人享有的 | 国产精品嫩草影院av在线观看| 精品国内亚洲2022精品成人| 赤兔流量卡办理| 亚洲一区高清亚洲精品| 色尼玛亚洲综合影院| 联通29元200g的流量卡| 亚洲欧美日韩东京热| 精品人妻偷拍中文字幕| 99热网站在线观看| 国产av一区在线观看免费| 国产大屁股一区二区在线视频| 国产成人freesex在线| 欧美又色又爽又黄视频| 99久久中文字幕三级久久日本| 卡戴珊不雅视频在线播放| 国产精品蜜桃在线观看| 精品人妻一区二区三区麻豆| 国产精品不卡视频一区二区| 国产成人一区二区在线| 亚洲欧洲日产国产| 99热这里只有精品一区| 99久久九九国产精品国产免费| 日韩一本色道免费dvd| 国产精品国产三级国产av玫瑰| 中文字幕熟女人妻在线| 欧美三级亚洲精品| 波多野结衣高清无吗| 国产精品一区二区三区四区久久| 长腿黑丝高跟| 国产白丝娇喘喷水9色精品| 精品国产一区二区三区久久久樱花 | 亚洲精品日韩在线中文字幕| 综合色av麻豆| 亚洲精品影视一区二区三区av| 亚洲四区av| 国产乱来视频区| av线在线观看网站| 九九在线视频观看精品| 国产国拍精品亚洲av在线观看| 精品无人区乱码1区二区| 两性午夜刺激爽爽歪歪视频在线观看| 日韩中字成人| 日日摸夜夜添夜夜爱| 精品欧美国产一区二区三| 国产日韩欧美在线精品| 久久亚洲精品不卡| 亚洲欧美精品综合久久99| 日韩欧美 国产精品| 欧美激情在线99| 18+在线观看网站| 亚州av有码| 成年av动漫网址| 男人舔奶头视频| 18禁裸乳无遮挡免费网站照片| 国产午夜精品论理片| 国产高清视频在线观看网站| 国产 一区 欧美 日韩| 2022亚洲国产成人精品| 欧美性感艳星| 人妻少妇偷人精品九色| 亚洲精品乱码久久久v下载方式| 欧美激情久久久久久爽电影| 久久99热这里只频精品6学生 | 久久热精品热| 纵有疾风起免费观看全集完整版 | 日韩欧美 国产精品| 国产成人福利小说| 一级爰片在线观看| 99热网站在线观看| 国产av一区在线观看免费| 少妇人妻一区二区三区视频| 国产69精品久久久久777片| 成年女人看的毛片在线观看| 欧美成人一区二区免费高清观看| 天天一区二区日本电影三级| 国国产精品蜜臀av免费| 高清在线视频一区二区三区 | 国产中年淑女户外野战色| 日本免费在线观看一区| 成人无遮挡网站| 麻豆成人av视频| 国产私拍福利视频在线观看| 特大巨黑吊av在线直播| 汤姆久久久久久久影院中文字幕 | 亚洲va在线va天堂va国产| 搞女人的毛片| 激情 狠狠 欧美| 国产高清三级在线| 亚洲综合色惰| 国产精品伦人一区二区| 国产精品福利在线免费观看| 午夜激情福利司机影院| 日本黄大片高清| 免费搜索国产男女视频| 久久鲁丝午夜福利片| 嫩草影院精品99| 乱码一卡2卡4卡精品| 欧美一区二区精品小视频在线| 成人欧美大片| 中文字幕av在线有码专区| 成人亚洲欧美一区二区av| 国产不卡一卡二| 亚洲熟妇中文字幕五十中出| 亚洲电影在线观看av| 亚洲欧美成人综合另类久久久 | 在线免费十八禁| 亚洲精品乱久久久久久| 日韩成人av中文字幕在线观看| 免费黄网站久久成人精品| 国语自产精品视频在线第100页| www日本黄色视频网| 亚洲最大成人手机在线| 精品久久久久久久人妻蜜臀av| 精华霜和精华液先用哪个| 一区二区三区乱码不卡18| 狂野欧美白嫩少妇大欣赏| 性色avwww在线观看| 欧美潮喷喷水| 国产探花极品一区二区| 国产大屁股一区二区在线视频| 亚洲国产精品合色在线| 美女大奶头视频| 日日干狠狠操夜夜爽| 特大巨黑吊av在线直播| 99视频精品全部免费 在线| 亚洲av福利一区| 亚洲欧美中文字幕日韩二区| 伊人久久精品亚洲午夜| 精品久久久久久久久亚洲| 一级毛片久久久久久久久女| 亚洲不卡免费看| 婷婷色麻豆天堂久久 | 啦啦啦韩国在线观看视频| 精品一区二区三区人妻视频| 久久韩国三级中文字幕| 中文字幕av在线有码专区| 久久久久久久久大av| 久久99热6这里只有精品| 身体一侧抽搐| 九九热线精品视视频播放| 在现免费观看毛片| 欧美不卡视频在线免费观看| 美女黄网站色视频| 干丝袜人妻中文字幕| 麻豆久久精品国产亚洲av| 欧美成人a在线观看| 亚洲国产日韩欧美精品在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲成av人片在线播放无| 亚洲欧洲国产日韩| 精品无人区乱码1区二区| 在线观看66精品国产| 欧美性感艳星| 国产综合懂色| 日本免费一区二区三区高清不卡| 男女边吃奶边做爰视频| 欧美成人精品欧美一级黄| 亚洲欧美精品自产自拍| АⅤ资源中文在线天堂| 国产精品美女特级片免费视频播放器| 欧美人与善性xxx| 亚洲精品乱久久久久久| 秋霞在线观看毛片| 特级一级黄色大片| 在线观看一区二区三区| 男插女下体视频免费在线播放| 小蜜桃在线观看免费完整版高清| 亚洲成人中文字幕在线播放| 久久久久久久国产电影| 国产 一区精品| 国产淫片久久久久久久久| 99久久九九国产精品国产免费| 淫秽高清视频在线观看| 色尼玛亚洲综合影院| 2022亚洲国产成人精品| 能在线免费看毛片的网站| 啦啦啦韩国在线观看视频| 久久国产乱子免费精品| av国产久精品久网站免费入址| 国产成人a区在线观看| 婷婷色麻豆天堂久久 | 午夜精品一区二区三区免费看| 寂寞人妻少妇视频99o| 欧美zozozo另类| 国产欧美日韩精品一区二区| 日本五十路高清| 国产精品精品国产色婷婷| 中文字幕制服av| av.在线天堂| 久久精品影院6| 中文字幕制服av| 国产一级毛片七仙女欲春2| 欧美精品国产亚洲| 99热这里只有精品一区| 亚洲精品,欧美精品| 1024手机看黄色片| 插逼视频在线观看| 岛国在线免费视频观看| 尾随美女入室| 九色成人免费人妻av| 久久精品夜色国产| 一个人看的www免费观看视频| 天堂√8在线中文| 一级av片app| 两个人的视频大全免费| 卡戴珊不雅视频在线播放| 日韩一区二区视频免费看| 亚洲乱码一区二区免费版| 国产黄片视频在线免费观看| 国产精品电影一区二区三区| 高清av免费在线| 成人午夜精彩视频在线观看| 久久久精品94久久精品| 黄色欧美视频在线观看| 少妇丰满av| 久99久视频精品免费| 欧美最新免费一区二区三区| 精品久久久久久久久av| 亚洲欧美精品综合久久99| 久久欧美精品欧美久久欧美| 国产人妻一区二区三区在| 亚洲电影在线观看av| 偷拍熟女少妇极品色| 亚洲人成网站在线观看播放| 高清日韩中文字幕在线| 最近中文字幕高清免费大全6| 日日摸夜夜添夜夜添av毛片| 免费搜索国产男女视频| 亚洲欧美精品自产自拍| 热99re8久久精品国产| 日本与韩国留学比较| 亚洲av电影不卡..在线观看| 最近最新中文字幕免费大全7| 亚洲在线观看片| 九九热线精品视视频播放| 日日摸夜夜添夜夜添av毛片| 99久久无色码亚洲精品果冻| 精品熟女少妇av免费看| 最近2019中文字幕mv第一页| 免费观看a级毛片全部| 国产在线一区二区三区精 | 久久久久网色| 高清午夜精品一区二区三区| 岛国在线免费视频观看| 秋霞伦理黄片| 久久国产乱子免费精品| 久久久国产成人免费| 人妻制服诱惑在线中文字幕| 少妇的逼好多水| 亚洲人成网站在线观看播放| 嫩草影院新地址| 小说图片视频综合网站| 在线免费十八禁| 日本与韩国留学比较| 亚洲国产成人一精品久久久| 亚洲一级一片aⅴ在线观看| 免费人成在线观看视频色| 欧美一区二区亚洲| 小说图片视频综合网站| 久99久视频精品免费| 两个人的视频大全免费| 国产精华一区二区三区| 国产成人精品一,二区| 又粗又爽又猛毛片免费看| 久久综合国产亚洲精品| 精品99又大又爽又粗少妇毛片| 欧美日韩一区二区视频在线观看视频在线 | 97热精品久久久久久| 国产白丝娇喘喷水9色精品| 岛国在线免费视频观看| 久久精品国产亚洲av涩爱| 最近视频中文字幕2019在线8| 国产色婷婷99| 精华霜和精华液先用哪个| 1024手机看黄色片| 国产一级毛片七仙女欲春2| 国产片特级美女逼逼视频| 在线天堂最新版资源| 久久韩国三级中文字幕| 99热精品在线国产| av在线观看视频网站免费| 99久国产av精品国产电影| av天堂中文字幕网| 在现免费观看毛片| 只有这里有精品99| 国产在视频线在精品| 亚洲国产欧洲综合997久久,| 欧美激情国产日韩精品一区| 成人无遮挡网站| 1000部很黄的大片| 极品教师在线视频| 国产私拍福利视频在线观看| 亚洲精华国产精华液的使用体验| 69av精品久久久久久| 天堂中文最新版在线下载 | 亚洲最大成人av| 国产精品.久久久| 大香蕉97超碰在线| 2021天堂中文幕一二区在线观| 亚洲国产精品久久男人天堂| 亚洲内射少妇av| 久久鲁丝午夜福利片| 午夜免费激情av| 日日啪夜夜撸| 中文字幕免费在线视频6| 国产国拍精品亚洲av在线观看| 久久久久免费精品人妻一区二区| 欧美xxxx黑人xx丫x性爽| 欧美精品一区二区大全| 超碰97精品在线观看| 国产成人免费观看mmmm| a级一级毛片免费在线观看| 69人妻影院| 熟女电影av网| 成年免费大片在线观看| 久久99热这里只有精品18| 天堂av国产一区二区熟女人妻| 少妇裸体淫交视频免费看高清| 国产老妇女一区| 日韩一区二区三区影片| 亚洲av成人精品一二三区| 国产成人a区在线观看| 亚洲国产欧洲综合997久久,| 啦啦啦韩国在线观看视频| 搡女人真爽免费视频火全软件| 国产乱来视频区| 国产 一区精品| 久久国内精品自在自线图片| 国产精品99久久久久久久久| 嫩草影院精品99| 久久久久免费精品人妻一区二区| 免费看美女性在线毛片视频| 免费av观看视频| 国产高清国产精品国产三级 | 91午夜精品亚洲一区二区三区| 久久久久性生活片| 亚洲国产精品sss在线观看| 午夜福利成人在线免费观看| 草草在线视频免费看| 国产成人精品一,二区| 高清av免费在线| 亚洲人成网站在线观看播放| 免费看美女性在线毛片视频| 熟妇人妻久久中文字幕3abv| 国产国拍精品亚洲av在线观看| 久久久a久久爽久久v久久| 国产高清视频在线观看网站| 免费一级毛片在线播放高清视频| 黑人高潮一二区| 中文字幕久久专区| 色哟哟·www| 国产又黄又爽又无遮挡在线| 亚洲欧美成人综合另类久久久 | 91久久精品电影网| 最后的刺客免费高清国语| 一个人看视频在线观看www免费| 97在线视频观看| 高清视频免费观看一区二区 | 午夜精品一区二区三区免费看| 日本av手机在线免费观看| 国产免费福利视频在线观看| 2021少妇久久久久久久久久久| 亚洲自拍偷在线| 欧美一级a爱片免费观看看| 搡老妇女老女人老熟妇| 在线天堂最新版资源| 青春草国产在线视频| 亚洲一级一片aⅴ在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲欧美日韩高清专用| 久久精品夜夜夜夜夜久久蜜豆| 99热精品在线国产| 精品久久久久久久末码| av在线观看视频网站免费| 久久久精品大字幕| 亚洲aⅴ乱码一区二区在线播放| 精品人妻偷拍中文字幕| 精品人妻一区二区三区麻豆| 18禁在线播放成人免费| av在线蜜桃| 成年av动漫网址| 最近视频中文字幕2019在线8| 国产午夜福利久久久久久| 久久草成人影院| 国产男人的电影天堂91| 赤兔流量卡办理| 午夜亚洲福利在线播放| 九草在线视频观看| 丝袜美腿在线中文| 国产精品.久久久| 亚洲va在线va天堂va国产| 亚洲成人av在线免费| 能在线免费看毛片的网站| 国产一级毛片七仙女欲春2| 久久久国产成人免费| 国产极品天堂在线| 亚洲av电影不卡..在线观看| 亚洲国产欧洲综合997久久,| 亚洲丝袜综合中文字幕| 69人妻影院| 夜夜看夜夜爽夜夜摸| 中文字幕精品亚洲无线码一区|