• <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
    两个人视频免费观看高清| 欧美日韩精品网址| 精品人妻1区二区| 哪里可以看免费的av片| 国产亚洲精品久久久com| 又爽又黄无遮挡网站| 日日摸夜夜添夜夜添小说| 欧美日韩综合久久久久久 | 男女那种视频在线观看| 在线天堂最新版资源| 午夜福利欧美成人| 日韩av在线大香蕉| 特大巨黑吊av在线直播| 日韩精品中文字幕看吧| av中文乱码字幕在线| 国产精品三级大全| 国产成人aa在线观看| 夜夜躁狠狠躁天天躁| 成年版毛片免费区| 村上凉子中文字幕在线| 51午夜福利影视在线观看| 亚洲成a人片在线一区二区| 日日夜夜操网爽| 一区二区三区激情视频| 天堂动漫精品| 国产麻豆成人av免费视频| 精品久久久久久久久久免费视频| 国内精品美女久久久久久| 亚洲激情在线av| 一个人观看的视频www高清免费观看| 99久久精品国产亚洲精品| 天堂av国产一区二区熟女人妻| 动漫黄色视频在线观看| 国产国拍精品亚洲av在线观看 | 日本熟妇午夜| 午夜福利欧美成人| 日韩高清综合在线| 国产成人啪精品午夜网站| 欧美精品啪啪一区二区三区| 色噜噜av男人的天堂激情| 日本免费a在线| 国产视频内射| 香蕉久久夜色| 欧美在线黄色| 狠狠狠狠99中文字幕| 又爽又黄无遮挡网站| 亚洲一区二区三区不卡视频| 国产一区二区在线观看日韩 | 少妇的逼水好多| 亚洲成人久久性| 国产精品一区二区三区四区久久| 可以在线观看的亚洲视频| 成人一区二区视频在线观看| 久久欧美精品欧美久久欧美| 性色av乱码一区二区三区2| 两人在一起打扑克的视频| 午夜免费男女啪啪视频观看 | 99久久综合精品五月天人人| 国产精品国产高清国产av| 国产不卡一卡二| 亚洲欧美一区二区三区黑人| 色尼玛亚洲综合影院| 亚洲av二区三区四区| 国产伦精品一区二区三区视频9 | h日本视频在线播放| 一a级毛片在线观看| 欧美成人性av电影在线观看| 天堂影院成人在线观看| 精品国产亚洲在线| 久久久久久国产a免费观看| 欧美日韩福利视频一区二区| 国产亚洲精品av在线| 国产精品久久久久久亚洲av鲁大| 国产精品精品国产色婷婷| 亚洲欧美激情综合另类| 日韩精品中文字幕看吧| 男女床上黄色一级片免费看| 精品欧美国产一区二区三| 夜夜躁狠狠躁天天躁| 国产中年淑女户外野战色| 脱女人内裤的视频| 久久欧美精品欧美久久欧美| 嫩草影院精品99| 两人在一起打扑克的视频| 欧美绝顶高潮抽搐喷水| 色综合亚洲欧美另类图片| 欧美成人一区二区免费高清观看| 亚洲在线观看片| 亚洲一区二区三区色噜噜| 长腿黑丝高跟| 国产麻豆成人av免费视频| 国产亚洲精品综合一区在线观看| 精品久久久久久,| 91在线精品国自产拍蜜月 | www.999成人在线观看| 国产精华一区二区三区| 99久久精品一区二区三区| 嫩草影视91久久| 免费看美女性在线毛片视频| svipshipincom国产片| 长腿黑丝高跟| 非洲黑人性xxxx精品又粗又长| 国产精品野战在线观看| 亚洲自拍偷在线| 亚洲欧美一区二区三区黑人| 黑人欧美特级aaaaaa片| 在线a可以看的网站| 色噜噜av男人的天堂激情| 亚洲精品粉嫩美女一区| 无遮挡黄片免费观看| 国产亚洲精品一区二区www| 波野结衣二区三区在线 | 久久性视频一级片| netflix在线观看网站| 欧美日本亚洲视频在线播放| 成人无遮挡网站| 国产 一区 欧美 日韩| 欧美日韩黄片免| 有码 亚洲区| 一级黄片播放器| 国产成+人综合+亚洲专区| 欧美bdsm另类| 精品欧美国产一区二区三| 久久久久久久久中文| 亚洲乱码一区二区免费版| а√天堂www在线а√下载| 亚洲国产中文字幕在线视频| 久久香蕉精品热| 十八禁网站免费在线| 无人区码免费观看不卡| 国产乱人伦免费视频| 国产91精品成人一区二区三区| 亚洲成av人片免费观看| 久久精品亚洲精品国产色婷小说| 免费电影在线观看免费观看| 精品无人区乱码1区二区| 男人的好看免费观看在线视频| 成人av在线播放网站| 女人高潮潮喷娇喘18禁视频| 综合色av麻豆| 久久伊人香网站| 一区二区三区激情视频| 欧美日韩中文字幕国产精品一区二区三区| 老司机福利观看| 国产老妇女一区| 亚洲精品粉嫩美女一区| 蜜桃久久精品国产亚洲av| 国产精品久久久久久久电影 | 亚洲无线观看免费| 宅男免费午夜| 夜夜夜夜夜久久久久| 婷婷亚洲欧美| 成年女人永久免费观看视频| 老司机福利观看| 亚洲 欧美 日韩 在线 免费| 国产精品女同一区二区软件 | 欧美3d第一页| 午夜福利在线观看免费完整高清在 | 有码 亚洲区| 夜夜夜夜夜久久久久| 婷婷精品国产亚洲av| 亚洲国产欧洲综合997久久,| 特大巨黑吊av在线直播| 午夜福利高清视频| 老鸭窝网址在线观看| 欧美一级a爱片免费观看看| av天堂在线播放| 国产精品久久久久久人妻精品电影| 久久人妻av系列| 精品电影一区二区在线| 亚洲精品影视一区二区三区av| 啦啦啦韩国在线观看视频| 午夜影院日韩av| 久久久久国内视频| 欧美乱色亚洲激情| 九色成人免费人妻av| 国产69精品久久久久777片| 小说图片视频综合网站| 日本免费a在线| 国产又黄又爽又无遮挡在线| 午夜影院日韩av| 亚洲精品久久国产高清桃花| 国产免费男女视频| 国产欧美日韩精品亚洲av| 两人在一起打扑克的视频| 在线观看66精品国产| 成年女人永久免费观看视频| 中文字幕人妻熟人妻熟丝袜美 | 亚洲久久久久久中文字幕| av福利片在线观看| 在线视频色国产色| 国产探花在线观看一区二区| 搞女人的毛片| 全区人妻精品视频| 最好的美女福利视频网| 波多野结衣高清作品| 亚洲成人久久爱视频| 色av中文字幕| 欧美日韩国产亚洲二区| 亚洲专区中文字幕在线| 精品国产三级普通话版| 男女之事视频高清在线观看| 午夜精品久久久久久毛片777| 99久久成人亚洲精品观看| 成年版毛片免费区| 亚洲精品日韩av片在线观看 | 欧美乱妇无乱码| 亚洲第一欧美日韩一区二区三区| netflix在线观看网站| 久久精品国产自在天天线| 亚洲国产色片| 亚洲专区国产一区二区| 搡老妇女老女人老熟妇| 天堂av国产一区二区熟女人妻| av欧美777| 在线观看免费午夜福利视频| 黄色日韩在线| 一个人免费在线观看的高清视频| 国产高清三级在线| 男女下面进入的视频免费午夜| 在线观看日韩欧美| 亚洲人成网站在线播放欧美日韩| 18禁美女被吸乳视频| 高潮久久久久久久久久久不卡| 无遮挡黄片免费观看| 欧美成狂野欧美在线观看| 搞女人的毛片| 国产一区二区三区在线臀色熟女| 看免费av毛片| 午夜a级毛片| 美女 人体艺术 gogo| 国产精品98久久久久久宅男小说| 少妇裸体淫交视频免费看高清| 国产aⅴ精品一区二区三区波| 亚洲精品一区av在线观看| 嫁个100分男人电影在线观看| 国产精品亚洲美女久久久| 88av欧美| 黄色丝袜av网址大全| 亚洲乱码一区二区免费版| 久久九九热精品免费| 两性午夜刺激爽爽歪歪视频在线观看| 欧美一区二区国产精品久久精品| 亚洲欧美日韩高清在线视频| 国产精品久久久久久精品电影| 国内精品久久久久久久电影| 三级国产精品欧美在线观看| 最后的刺客免费高清国语| 国产亚洲精品av在线| 男女床上黄色一级片免费看| 欧美大码av| 丁香欧美五月| 一级黄片播放器| 精品国产亚洲在线| 最后的刺客免费高清国语| 欧美+日韩+精品| 亚洲精品粉嫩美女一区| 亚洲在线观看片| 国产高清激情床上av| 国产视频内射| 一进一出抽搐gif免费好疼| 一个人观看的视频www高清免费观看| 在线国产一区二区在线| 国产国拍精品亚洲av在线观看 | 在线视频色国产色| av片东京热男人的天堂| 日韩欧美国产一区二区入口| 别揉我奶头~嗯~啊~动态视频| 亚洲欧美日韩东京热| 国产一区二区在线av高清观看| 国产成人啪精品午夜网站| 欧美日本亚洲视频在线播放| 中文字幕精品亚洲无线码一区| 国产精品亚洲一级av第二区| 国产黄色小视频在线观看| 成人欧美大片| 国产欧美日韩精品亚洲av| 国内久久婷婷六月综合欲色啪| 人人妻人人澡欧美一区二区| 非洲黑人性xxxx精品又粗又长| 国产极品精品免费视频能看的| 国产精品久久久久久人妻精品电影| 在线观看舔阴道视频| 中文字幕人妻丝袜一区二区| 国产色爽女视频免费观看| 啦啦啦免费观看视频1| 国内精品久久久久久久电影| 成人国产一区最新在线观看| 午夜精品久久久久久毛片777| 久久香蕉精品热| 亚洲国产精品成人综合色| 国产av麻豆久久久久久久| 午夜免费观看网址| 精品人妻一区二区三区麻豆 | 美女高潮的动态| 日韩免费av在线播放| 午夜老司机福利剧场| 亚洲欧美日韩东京热| 精品久久久久久久毛片微露脸| 夜夜躁狠狠躁天天躁| 久久久成人免费电影| 色视频www国产| 国产黄片美女视频| 国产一区在线观看成人免费| 日韩中文字幕欧美一区二区| 丰满人妻一区二区三区视频av | 少妇人妻精品综合一区二区 | 亚洲av免费在线观看| 久久精品91无色码中文字幕| 色av中文字幕| 一个人看的www免费观看视频| 最新中文字幕久久久久| 国产成人影院久久av| 亚洲中文日韩欧美视频| 亚洲成人中文字幕在线播放| www.www免费av| 一区二区三区国产精品乱码| av国产免费在线观看| 十八禁人妻一区二区| 啦啦啦观看免费观看视频高清| 精品电影一区二区在线| 欧美性感艳星| 国产69精品久久久久777片| 日韩欧美在线二视频| 久久精品91蜜桃| 久久精品91无色码中文字幕| 又紧又爽又黄一区二区| 99国产精品一区二区蜜桃av| 国产精品国产高清国产av| 最近视频中文字幕2019在线8| 男女那种视频在线观看| 热99在线观看视频| 国产麻豆成人av免费视频| 九九在线视频观看精品| 夜夜爽天天搞| 免费人成在线观看视频色| eeuss影院久久| xxxwww97欧美| 亚洲成人精品中文字幕电影| 成年女人毛片免费观看观看9| av中文乱码字幕在线| 无人区码免费观看不卡| 日本五十路高清| 99热6这里只有精品| 亚洲专区中文字幕在线| 一级黄片播放器| 精品久久久久久久久久久久久| 欧美在线黄色| 法律面前人人平等表现在哪些方面| 亚洲精品日韩av片在线观看 | 国产亚洲精品久久久久久毛片| 啪啪无遮挡十八禁网站| 在线a可以看的网站| 久久天躁狠狠躁夜夜2o2o| 少妇高潮的动态图| 男女视频在线观看网站免费| 欧美日韩综合久久久久久 | 婷婷精品国产亚洲av| 国产成年人精品一区二区| 日本黄色片子视频| 婷婷亚洲欧美| 嫩草影院精品99| 最近最新中文字幕大全电影3| 免费看a级黄色片| 国模一区二区三区四区视频| 日日摸夜夜添夜夜添小说| 日韩 欧美 亚洲 中文字幕| 亚洲精品久久国产高清桃花| 国产精品野战在线观看| 久久人人精品亚洲av| 国产亚洲精品久久久com| 欧美日本视频| 国产欧美日韩一区二区三| 欧美日韩一级在线毛片| 在线播放无遮挡| 99精品欧美一区二区三区四区| av片东京热男人的天堂| 亚洲专区国产一区二区| 亚洲av日韩精品久久久久久密| 狂野欧美激情性xxxx| 国产极品精品免费视频能看的| 日日干狠狠操夜夜爽| 国产精品野战在线观看| 精品国产三级普通话版| 热99re8久久精品国产| 波野结衣二区三区在线 | 欧美xxxx黑人xx丫x性爽| 看黄色毛片网站| 久久精品国产自在天天线| 亚洲av电影不卡..在线观看| 色综合亚洲欧美另类图片| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 一二三四社区在线视频社区8| 国产成人a区在线观看| 十八禁网站免费在线| 可以在线观看毛片的网站| 国产一区二区激情短视频| 中文字幕av在线有码专区| 亚洲无线在线观看| 欧美绝顶高潮抽搐喷水| 久久6这里有精品| 老熟妇乱子伦视频在线观看| 午夜免费激情av| 身体一侧抽搐| 午夜激情福利司机影院| 成人三级黄色视频| 欧洲精品卡2卡3卡4卡5卡区| 在线免费观看的www视频| 国产精品一区二区免费欧美| 国产一区二区激情短视频| 免费av观看视频| 少妇裸体淫交视频免费看高清| 在线十欧美十亚洲十日本专区| 欧美成人性av电影在线观看| 女同久久另类99精品国产91| 国产淫片久久久久久久久 | 亚洲片人在线观看| 最近最新中文字幕大全电影3| 国产精品一区二区三区四区久久| 亚洲中文日韩欧美视频| 18禁美女被吸乳视频| 激情在线观看视频在线高清| 日本一二三区视频观看| 日本免费一区二区三区高清不卡| 18禁黄网站禁片午夜丰满| 国产男靠女视频免费网站| 老司机午夜福利在线观看视频| av欧美777| 观看免费一级毛片| 色综合站精品国产| 特级一级黄色大片| 亚洲精华国产精华精| 亚洲电影在线观看av| 成人av在线播放网站| 人人妻,人人澡人人爽秒播| 国产精品亚洲av一区麻豆| 久久国产乱子伦精品免费另类| 一卡2卡三卡四卡精品乱码亚洲| 99久久无色码亚洲精品果冻| 亚洲人成网站在线播| 嫩草影视91久久| 国产精品久久视频播放| 香蕉久久夜色| 天堂√8在线中文| 久久国产精品人妻蜜桃| 精品国产三级普通话版| 久久草成人影院| 久久久久久九九精品二区国产| 亚洲五月婷婷丁香| 亚洲精品一卡2卡三卡4卡5卡| 国产精品自产拍在线观看55亚洲| 午夜精品一区二区三区免费看| 热99在线观看视频| 成年女人毛片免费观看观看9| 91久久精品电影网| 最近最新免费中文字幕在线| 又黄又粗又硬又大视频| 国产一区二区在线av高清观看| 中文在线观看免费www的网站| 亚洲乱码一区二区免费版| 黄色片一级片一级黄色片| 国产视频一区二区在线看| 国产高清有码在线观看视频| 国产亚洲精品综合一区在线观看| 99在线人妻在线中文字幕| 日韩欧美国产在线观看| 午夜a级毛片| 女人被狂操c到高潮| 一进一出抽搐gif免费好疼| 国产精华一区二区三区| 免费av不卡在线播放| 最近最新中文字幕大全免费视频| 亚洲国产欧美网| АⅤ资源中文在线天堂| 男女那种视频在线观看| 97超视频在线观看视频| 精品久久久久久久久久久久久| 男女做爰动态图高潮gif福利片| 国产精品久久久久久久久免 | 日本熟妇午夜| 亚洲性夜色夜夜综合| 国产激情偷乱视频一区二区| 国产高清视频在线观看网站| 精品久久久久久久久久免费视频| 国产免费av片在线观看野外av| 国产亚洲欧美98| 国产精品久久久久久久电影 | 高潮久久久久久久久久久不卡| 午夜影院日韩av| 嫁个100分男人电影在线观看| 成人永久免费在线观看视频| 国产精品三级大全| 99热精品在线国产| 亚洲激情在线av| 国产精品久久久久久久久免 | 日本一本二区三区精品| 99国产极品粉嫩在线观看| 一区二区三区免费毛片| 欧美日本视频| 中文字幕人妻丝袜一区二区| 丰满人妻熟妇乱又伦精品不卡| 久久久久久九九精品二区国产| 偷拍熟女少妇极品色| 亚洲男人的天堂狠狠| 在线天堂最新版资源| 国产成人福利小说| 国产亚洲欧美98| 亚洲一区二区三区色噜噜| 天堂动漫精品| 国内精品久久久久久久电影| 欧美绝顶高潮抽搐喷水| 99久久精品热视频| 日韩亚洲欧美综合| 看黄色毛片网站| 国产真人三级小视频在线观看| 欧美黑人巨大hd| 高潮久久久久久久久久久不卡| 亚洲av成人不卡在线观看播放网| 免费大片18禁| 内地一区二区视频在线| 成年女人看的毛片在线观看| 久久久久久久午夜电影| 老司机午夜福利在线观看视频| 三级男女做爰猛烈吃奶摸视频| 亚洲av成人av| 亚洲一区高清亚洲精品| 欧美不卡视频在线免费观看| 88av欧美| 美女高潮的动态| 国产视频一区二区在线看| 久久久久久国产a免费观看| 中出人妻视频一区二区| 国产视频一区二区在线看| 亚洲国产欧美人成| 色av中文字幕| 国产免费一级a男人的天堂| 国产精品女同一区二区软件 | 久久精品91蜜桃| 国产精品自产拍在线观看55亚洲| 蜜桃久久精品国产亚洲av| 我的老师免费观看完整版| 日本免费a在线| 黑人欧美特级aaaaaa片| 男女午夜视频在线观看| 日韩高清综合在线| 在线播放无遮挡| 高清日韩中文字幕在线| 免费在线观看成人毛片| 手机成人av网站| 国产一区二区亚洲精品在线观看| 91麻豆av在线| 国产高潮美女av| 啪啪无遮挡十八禁网站| 真人做人爱边吃奶动态| 热99在线观看视频| 热99re8久久精品国产| 变态另类成人亚洲欧美熟女| 综合色av麻豆| 一本精品99久久精品77| 亚洲乱码一区二区免费版| 国产亚洲av嫩草精品影院| 熟妇人妻久久中文字幕3abv| 欧美黑人欧美精品刺激| 亚洲第一欧美日韩一区二区三区| av在线蜜桃| 一个人观看的视频www高清免费观看| 丁香六月欧美| 亚洲激情在线av| 成年女人毛片免费观看观看9| 欧美日韩黄片免| 日韩av在线大香蕉| 国产亚洲精品综合一区在线观看| 他把我摸到了高潮在线观看| 欧美丝袜亚洲另类 | 老汉色av国产亚洲站长工具| 亚洲av电影不卡..在线观看| 成人国产一区最新在线观看| 99riav亚洲国产免费| 99国产精品一区二区三区| АⅤ资源中文在线天堂| 丰满人妻熟妇乱又伦精品不卡| 日本与韩国留学比较| av女优亚洲男人天堂| av天堂中文字幕网| 欧美日韩乱码在线| a级一级毛片免费在线观看| 欧美午夜高清在线| 亚洲欧美精品综合久久99| 欧美+亚洲+日韩+国产| 国产亚洲av嫩草精品影院| 老熟妇乱子伦视频在线观看| 深爱激情五月婷婷| 看片在线看免费视频| 亚洲,欧美精品.| 亚洲国产精品久久男人天堂| 色在线成人网| 午夜激情欧美在线| 老汉色∧v一级毛片| 长腿黑丝高跟| 天堂影院成人在线观看| 少妇丰满av| 久久草成人影院| 免费看a级黄色片| 国产麻豆成人av免费视频| 国产一级毛片七仙女欲春2| 国产精品久久视频播放| 日韩精品中文字幕看吧| 少妇裸体淫交视频免费看高清| 久久久久久久久中文| 亚洲内射少妇av| 久久人妻av系列| 国产精品久久久久久精品电影| 欧美日韩亚洲国产一区二区在线观看| 欧美日本亚洲视频在线播放| 美女 人体艺术 gogo|