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

    基于GPS測站坐標(biāo)殘差序列的ARMA建模方法研究*

    2012-11-14 13:48:24孫付平李楚陽
    大地測量與地球動力學(xué) 2012年2期
    關(guān)鍵詞:階數(shù)測站殘差

    田 亮 孫付平 李楚陽

    (1)解放軍61363部隊,西安 710054 2)解放軍信息工程大學(xué)測繪學(xué)院,鄭州450000)

    基于GPS測站坐標(biāo)殘差序列的ARMA建模方法研究*

    田 亮1)孫付平2)李楚陽1)

    (1)解放軍61363部隊,西安 710054 2)解放軍信息工程大學(xué)測繪學(xué)院,鄭州450000)

    以GPS測站坐標(biāo)殘差序列為實驗對象,采用ARMA模型兩種經(jīng)典建模方法(Box-Jenkins方法和Pandit-Wu方法)進(jìn)行了建模擬合的比較與分析。實驗比較發(fā)現(xiàn),當(dāng)殘差數(shù)據(jù)比較豐富時兩種建模方法擬合效果理想,無明顯差異,當(dāng)實驗數(shù)據(jù)較少時,Pandit-Wu方法擬合效果更好一些,剩余殘差振幅較Box-Jenkins方法減少約,且分布更加平穩(wěn)。

    GPS;坐標(biāo)殘差;ARMA模型;時間序列分析;ITRF2008

    1 引言

    隨著各種空間定位技術(shù)的迅速發(fā)展與成熟,國際地球參考框架的精度和穩(wěn)定性取得了很大提高[1]。但是為了更加精確地描述地球表面各種毫米量級地球動力學(xué)變化,國際大地測量學(xué)界又提出了地心坐標(biāo)精度接近或達(dá)到1 mm的新目標(biāo)。要達(dá)到如此高的精度目標(biāo),勢必對各種空間技術(shù)的觀測精度以及數(shù)據(jù)處理方法和模型提出了更高的要求。在這樣的背景下國際地球參考框架ITRF2008于2010年成功推出。該地球參考框架無論精度還是穩(wěn)定性方面都大大優(yōu)于以往的ITRF系列[2]。但是通過分析ITRF2008測站坐標(biāo)剩余殘差序列,發(fā)現(xiàn)依然包含許多周期性的誤差規(guī)律沒有進(jìn)行模型改正。本文曾對這些尚未模型改正的誤差規(guī)律進(jìn)行了建模擬合研究[3],發(fā)現(xiàn)ARMA建模擬合效果明顯優(yōu)于三角函數(shù)法以及小波分析等擬合方法。因此本文對ARMA模型的不同建模方法做了進(jìn)一步比較研究,對常用的Box-Jenkins方法和Pandit-Wu方法進(jìn)行實驗比較與分析。

    2 計算方法

    ARMA模型定義為[4,5]

    其中{εt}~WN(0,σ2),p、q≥0為整數(shù),(p,q)為模型的階數(shù),記為{Xt}~ARMA(p,q)。

    值得注意的是ARMA模型建模時要求時間序列必須為平穩(wěn)序列,但在實際測量中,時間序列往往包含一些非平穩(wěn)的趨勢項和周期項。為此在使用模型擬合前要先提取趨勢項和周期部分或者采用差分方法對時間序列進(jìn)行預(yù)處理,使其滿足平穩(wěn)性的要求。基于差分后的時間序列所建立的模型又稱為ARIMA模型,記為{Xt}~ARIMA(p,d,q)。

    2.1 Box-Jenkins方法建模

    經(jīng)數(shù)據(jù)預(yù)處理后殘差序列采用Box-Jenkins方法建模擬合,主要建模步驟如下:

    1)殘差序列平穩(wěn)性(ADF)檢驗。確定殘差序列是否平穩(wěn),如果殘差序列為非平穩(wěn)序列,則進(jìn)行差分處理。

    2)模型識別。判斷殘差時間序列的自相關(guān)和偏自相關(guān)函數(shù)的形式,如果模型自相關(guān)函數(shù)拖尾,偏自相關(guān)函數(shù)截尾,采用AR模型;如果模型自相關(guān)函數(shù)截尾,偏自相關(guān)函數(shù)拖尾,采用MA模型;如果兩者皆拖尾,采用ARMA模型。

    3)模型定階采用AIC最小信息準(zhǔn)則。

    4)模型參數(shù)估計采用非線性最小二乘估計方法。

    5)模型的適應(yīng)性檢驗,即檢驗剩余殘差是否為白噪聲序列。

    2.2 Pandit-Wu方法建模

    具體建模步驟如下:

    1)ADF檢驗。確定殘差序列是否平穩(wěn),如果檢驗殘差序列為非平穩(wěn)序列,將進(jìn)行差分處理。

    2)從n=1開始,逐漸增加模型階數(shù),擬合ARMA(2n,2n-1)模型,依次為ARMA(2,1)、ARMA (4,3)、ARMA(6,5)…即以兩階為單位逐漸增加模型階數(shù),直到F檢驗表明當(dāng)再增加模型階數(shù)時,剩余殘差平方和不再顯著減小為止。如果系統(tǒng)的確是階數(shù)為奇數(shù)的模型,那么擬合的較高階數(shù)為偶數(shù)的模型第2n個自回歸系數(shù)的絕對值必然接近于零。那么我們可以刪除這樣的小參數(shù)得到較低奇數(shù)階的模型。

    3)模型適應(yīng)性檢驗。

    4)求最優(yōu)模型。采用F檢驗判斷模型是否包含小參數(shù),若有參數(shù)可以刪除,擬合較低階的模型,從而得到最優(yōu)模型。

    3 算例分析

    3.1 數(shù)據(jù)選取

    選取的實驗數(shù)據(jù)為ITRF2008解的衍生產(chǎn)品GPS測站坐標(biāo)殘差序列[6-8],絕大部分GPS坐標(biāo)殘差序列包含了1997—2009年約12年的殘差數(shù)據(jù),采樣間隔為7天,部分殘差序列存在間斷點或粗差。因此首先對選取的殘差序列進(jìn)行數(shù)據(jù)預(yù)處理,對于含有間斷點的殘差序列進(jìn)行插值擬合,存在粗差的情況,設(shè)定閾值(通過殘差文件中給出的formal error來判斷)進(jìn)行剔除。

    3.2 實驗結(jié)果比較與分析

    由于全球GPS測站坐標(biāo)殘差序列數(shù)據(jù)量非常大,受篇幅限制本文僅選取部分測站作為實驗對象,HOFN測站(64°N,164°W)數(shù)據(jù)累積豐富且觀測數(shù)據(jù)質(zhì)量較高,沒有明顯的外部因素突發(fā)性影響,因此以HOFN測站坐標(biāo)殘差序列為實驗對象非常理想。

    首先利用Box-Jenkins方法對HOFN測站坐標(biāo)殘差序列進(jìn)行建模擬合,定階準(zhǔn)則采用AIC信息準(zhǔn)則,階數(shù)選取如表1所示,實際建模效果如圖1~3所示。

    表1 ARMA模型(Box-Jenkins方法)階數(shù)選取Tab.1 Selection of the order number with ARMA model (Box-Jenkins)

    利用Pandit-Wu方法建模階數(shù)選取如表2所示,實際建模效果如圖4~6所示。

    通過觀察擬合圖和剩余殘差圖可以發(fā)現(xiàn),兩種建模方法的擬合效果都比較理想,沒有顯著差異。細(xì)微來看Pandit-Wu方法擬合后剩余殘差振幅比Box-Jenkins方法稍小,但是效果不明顯。為了進(jìn)一步區(qū)分比較兩種方法,對上述實驗進(jìn)行改動,實驗數(shù)據(jù)截取一半即選取HOFN測站殘差序列5年的殘差數(shù)據(jù)(1997—2001年的時間序列)作為實驗對象,采用上述兩種方法進(jìn)行建模擬合比較。擬合效果如圖7~9所示。

    表2 ARMA模型(Pandit-Wu方法)階數(shù)選取Tab.2 Selection of the order number of ARMA model (Pandit-Wu)

    圖1 垂向擬合結(jié)果Fig.1 Vertical fitting results with Box-Jenkins method

    圖2 東向擬合結(jié)果Fig.2 Eastward fitting results with Box-Jenkins method

    圖3 北向擬合結(jié)果Fig.3 Northward fitting results with Box-Jenkins method

    圖4 垂向擬合結(jié)果Fig.4 Vertical fitting results with Pandit-Wu method

    圖5 東向擬合結(jié)果Fig.5 Eastward fitting results with Pandit-Wu method

    圖6 北向擬合結(jié)果Fig.6 Northward fitting results with Pandit-Wu method

    圖7 垂向擬合剩余殘差比較Fig.7 Comparison between the residuals after vertical fitting with the two methods

    圖8 北向擬合剩余殘差比較Fig.8 Comparison between residuals after northward fitting with the two methods

    通過3個方向的剩余殘差圖比較可以發(fā)現(xiàn),兩種建模方法在實驗數(shù)據(jù)減少一半的情況下擬合結(jié)果出現(xiàn)了較大差別:Pandit-Wu方法擬合后垂向剩余殘差振幅在 6 mm左右,Box-Jenkins方法達(dá)到了10mm;Pandit-Wu方法擬合后殘差分布仍然無明顯周期性振動規(guī)律,而Box-Jenkins方法擬合剩余殘差明顯存在一定的周期性振動規(guī)律。東向和北向剩余殘差經(jīng)Pandit-Wu方法擬合后振幅在2mm左右,Box-Jenkins方法達(dá)到了3 mm。由此可見Pandit-Wu方法的擬合效果優(yōu)于Box-Jenkins方法。

    圖9 東向擬合剩余殘差比較Fig.9 Comparison between residuals after eastward fitting with the two methods

    之所以Pandit-Wu方法擬合效果更加理想,從原理上分析主要是因為Box-Jenkins建模方法是以時間序列的自相關(guān)函數(shù)和偏自相關(guān)函數(shù)的統(tǒng)計特性為依據(jù),但是在建模之前我們并不知道時間序列的理論自相關(guān)函數(shù),而是用樣本自相關(guān)函數(shù)和偏自相關(guān)函數(shù)來近似,因此不可避免地會產(chǎn)生誤差。尤其是當(dāng)采樣數(shù)據(jù)較少時,樣本無法完全代表總體殘差分布,誤差將不可避免。而Pandit-Wu方法無需了解殘差序列特性,只需要逐漸增加模型階數(shù),直到剩余殘差平方和不再顯著減小為止??梢娫跓o法預(yù)知樣本內(nèi)部信息的前提下,Pandit-Wu方法更加方便可靠。

    5 結(jié)束語

    通過對GPS測站坐標(biāo)殘差序列采用ARMA兩種經(jīng)典建模方法進(jìn)行建模比較發(fā)現(xiàn),Pandit-Wu方法的建模效果無論是在數(shù)據(jù)豐富還是缺乏時都能保持得很好,采用該法對剩余殘差序列進(jìn)行建模擬合對于進(jìn)一步提高測站地心坐標(biāo)精度,提高國際地球參考框架的精度和穩(wěn)定性有重要意義。

    1 陳俊勇.大地坐標(biāo)框架理論和實踐的進(jìn)展[J].大地測量與地球動力學(xué),2007,(1):1-6.(Chen Junyong.Progress in theory and practice for geodetic reference coordinate frame[J].Journal of Geodesy and Geodynamics,2007,(1):1-6)

    2 黃立人.地殼運動的參考框架[J].大地測量與地球動力學(xué),2002,(3):102-108.(Huang Liren.The reference frame of earth deformation[J].Journal of Geodesy and Geodynamics,2002,(3):102-108)

    3 田亮,等.GPS測站坐標(biāo)非線性變化研究方法的比較與分析[J].測繪工程,2011,20(2):27-29.(Tian Liang,et al.Comparative analysis of studying on GPS station non-linear variation[J].Engineering of Surveying and Mapping,2011,20(2):27-29)

    4 王振龍,胡永宏.應(yīng)用時間序列分析[M].北京:科學(xué)出版社,2007.(Wang Zhenlong and Hu Yonghong.The application of time series analysis[M].Beijing:Science Press, 2007)

    5 陳杰.MATLAB寶典[M].北京:電子工業(yè)出版社,2006.(Chen Jie.MATLAB canon[M].Beijing:Electronics Industry Press,2006)

    6 http://itrf.ensg.ign.fr/ITRF_solutions/2008/ITRF2008.php,2010.

    7 http://itrf.ensg.ign.fr/ITRF_solutions/2008/doc/ ITRF2008_GPS.SSC.txt,2010.

    8 Dennis D and McCarthy D D(eds.).IERS Conventions (2003),IERS Technical Note No.32,2003.

    STUDY ON ARMA MODEL BASED ON GPS STATION COORDINATE RESIDUAL SERIES

    Tian Liang1),Sun Fuping2)and Li Chuyang1)

    (1)61363 Troops of PLA,Xi’an 710054 2)PLA Information Engineering University,Zhengzhou 450000)

    As to GPS station coordinate residuals,we used two classical ARMA modeling methods(Box-Jenkins method and Pandit-Wu method)to fit the GPS station coordinate residuals,and made comparative study.According to the experiment,we found out that the two methods have a little difference while there is abundant data.If the data is poor,with the Pandit-Wu method the results will be better.The amplitude of residuals reduces aboutand the distribution is more uniform.

    GPS;coordinate residuals;ARMA model;time series analysis;ITRF2008

    1671-5942(2012)02-0124-04

    2011-11-30

    國家自然科學(xué)基金(41074011)

    田亮,1987年生,工學(xué)碩士,助理工程師,主要研究方向為空間大地測量.E-mail:tl200431610237@126.com

    P207

    A

    猜你喜歡
    階數(shù)測站殘差
    GNSS鐘差估計中的兩種測站選取策略分析
    基于雙向GRU與殘差擬合的車輛跟馳建模
    關(guān)于無窮小階數(shù)的幾點注記
    確定有限級數(shù)解的階數(shù)上界的一種n階展開方法
    基于殘差學(xué)習(xí)的自適應(yīng)無人機目標(biāo)跟蹤算法
    基于遞歸殘差網(wǎng)絡(luò)的圖像超分辨率重建
    全球GPS測站垂向周年變化統(tǒng)計改正模型的建立
    測站分布對GPS解算ERP的影響分析
    城市勘測(2018年1期)2018-03-15 03:37:02
    平穩(wěn)自相關(guān)過程的殘差累積和控制圖
    河南科技(2015年8期)2015-03-11 16:23:52
    一種新的多址信道有效階數(shù)估計算法*
    少妇裸体淫交视频免费看高清 | 91麻豆精品激情在线观看国产 | 欧美成人免费av一区二区三区 | 久久亚洲精品不卡| 女性生殖器流出的白浆| 欧美黑人欧美精品刺激| 伦理电影免费视频| 女人被狂操c到高潮| 国产色视频综合| 欧美日韩视频精品一区| 大香蕉久久成人网| 日本一区二区免费在线视频| 一进一出抽搐动态| 69精品国产乱码久久久| 国产国语露脸激情在线看| 国产成人免费观看mmmm| 免费av中文字幕在线| 日日摸夜夜添夜夜添小说| 国产精品久久久久成人av| 亚洲专区国产一区二区| 91字幕亚洲| 国产成人系列免费观看| 大陆偷拍与自拍| 国产精品美女特级片免费视频播放器 | 一级a爱片免费观看的视频| 老司机影院毛片| bbb黄色大片| 中文字幕av电影在线播放| 国产一区二区三区视频了| 91麻豆精品激情在线观看国产 | 精品第一国产精品| 亚洲成人免费电影在线观看| 人人妻人人爽人人添夜夜欢视频| 99riav亚洲国产免费| 高清在线国产一区| 精品国产超薄肉色丝袜足j| 国产精品国产av在线观看| 久久久久久久久久久久大奶| 亚洲片人在线观看| 老鸭窝网址在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 中文字幕精品免费在线观看视频| 一级a爱片免费观看的视频| 午夜亚洲福利在线播放| 欧美另类亚洲清纯唯美| 国产视频一区二区在线看| 国产精品久久久人人做人人爽| 色播亚洲综合网| 精品久久久久久成人av| 在线观看免费视频日本深夜| www.999成人在线观看| 亚洲av免费高清在线观看| 九色国产91popny在线| 久久香蕉精品热| 国产精品,欧美在线| 在线观看午夜福利视频| 免费观看精品视频网站| 亚洲精品在线观看二区| 69人妻影院| 久久草成人影院| 熟妇人妻久久中文字幕3abv| 久久精品国产亚洲av香蕉五月| 国产精品 国内视频| 亚洲不卡免费看| 欧美日韩精品网址| 久久婷婷人人爽人人干人人爱| 一个人观看的视频www高清免费观看| 色av中文字幕| 精品一区二区三区视频在线观看免费| АⅤ资源中文在线天堂| 97超级碰碰碰精品色视频在线观看| 99久久99久久久精品蜜桃| 最好的美女福利视频网| 免费人成视频x8x8入口观看| 男女那种视频在线观看| 亚洲专区中文字幕在线| 观看美女的网站| 久久草成人影院| 国产毛片a区久久久久| 国产不卡一卡二| 亚洲国产精品成人综合色| 成熟少妇高潮喷水视频| 欧美区成人在线视频| 91麻豆精品激情在线观看国产| 午夜福利视频1000在线观看| 中文字幕高清在线视频| 人人妻人人澡欧美一区二区| 亚洲黑人精品在线| 国产色婷婷99| 日本免费一区二区三区高清不卡| 嫩草影院精品99| 夜夜夜夜夜久久久久| 久久精品国产综合久久久| 国产精品一区二区三区四区免费观看 | 在线观看av片永久免费下载| 国产亚洲精品久久久久久毛片| 国产一级毛片七仙女欲春2| 色哟哟哟哟哟哟| 国产三级在线视频| 亚洲成人中文字幕在线播放| 国产欧美日韩一区二区三| 男女下面进入的视频免费午夜| 男人舔奶头视频| 伊人久久大香线蕉亚洲五| ponron亚洲| 搡女人真爽免费视频火全软件 | 天堂动漫精品| 日本免费a在线| 久久久久九九精品影院| 97超级碰碰碰精品色视频在线观看| 国产aⅴ精品一区二区三区波| 无限看片的www在线观看| 精品一区二区三区视频在线观看免费| 中文字幕人成人乱码亚洲影| 综合色av麻豆| 欧美成人a在线观看| 国产69精品久久久久777片| 啦啦啦观看免费观看视频高清| 国语自产精品视频在线第100页| 最好的美女福利视频网| 欧美性猛交黑人性爽| 午夜福利视频1000在线观看| 亚洲精品乱码久久久v下载方式 | 叶爱在线成人免费视频播放| 毛片女人毛片| 国产高清videossex| 国产黄片美女视频| 欧美在线黄色| 女生性感内裤真人,穿戴方法视频| 天天躁日日操中文字幕| 久久草成人影院| 99久久99久久久精品蜜桃| 每晚都被弄得嗷嗷叫到高潮| 国产主播在线观看一区二区| 日本精品一区二区三区蜜桃| 亚洲av日韩精品久久久久久密| 久久精品91无色码中文字幕| 国产精品 国内视频| 麻豆成人午夜福利视频| av片东京热男人的天堂| 亚洲欧美精品综合久久99| 亚洲在线观看片| 91在线观看av| 亚洲精品在线美女| 伊人久久大香线蕉亚洲五| 国产精品国产高清国产av| 亚洲最大成人中文| 黄色女人牲交| 久久久久国内视频| 香蕉久久夜色| 午夜久久久久精精品| 内地一区二区视频在线| av福利片在线观看| 免费在线观看日本一区| 操出白浆在线播放| 18禁裸乳无遮挡免费网站照片| 日本免费a在线| 91在线精品国自产拍蜜月 | 国产又黄又爽又无遮挡在线| 欧美在线一区亚洲| 国产精品三级大全| 国产精品久久久人人做人人爽| 国产不卡一卡二| 91麻豆精品激情在线观看国产| 老鸭窝网址在线观看| 国产一区二区在线观看日韩 | 最新在线观看一区二区三区| 九色国产91popny在线| 亚洲avbb在线观看| 女人被狂操c到高潮| 成人鲁丝片一二三区免费| 欧美国产日韩亚洲一区| 亚洲在线观看片| 香蕉久久夜色| 成人精品一区二区免费| 可以在线观看毛片的网站| 亚洲av五月六月丁香网| 中文亚洲av片在线观看爽| 欧美日韩黄片免| 亚洲av成人不卡在线观看播放网| 少妇裸体淫交视频免费看高清| 成人午夜高清在线视频| 90打野战视频偷拍视频| 久久久久久国产a免费观看| 亚洲国产欧洲综合997久久,| 亚洲成人中文字幕在线播放| 男女做爰动态图高潮gif福利片| 午夜影院日韩av| 午夜免费成人在线视频| 757午夜福利合集在线观看| 精品一区二区三区av网在线观看| 欧美色欧美亚洲另类二区| а√天堂www在线а√下载| 国产精品,欧美在线| 久久人人精品亚洲av| 88av欧美| 国产一区二区三区在线臀色熟女| 黄片大片在线免费观看| 欧美国产日韩亚洲一区| 国产精品98久久久久久宅男小说| 日韩欧美国产一区二区入口| 变态另类成人亚洲欧美熟女| 全区人妻精品视频| 日韩国内少妇激情av| 手机成人av网站| 亚洲人成网站在线播放欧美日韩| 在线播放无遮挡| 夜夜看夜夜爽夜夜摸| 99久久精品热视频| 亚洲欧美日韩卡通动漫| 91在线精品国自产拍蜜月 | 最近最新中文字幕大全电影3| 色综合婷婷激情| 两个人的视频大全免费| 国产伦精品一区二区三区视频9 | 最近最新中文字幕大全免费视频| 91在线精品国自产拍蜜月 | 麻豆一二三区av精品| 亚洲在线观看片| 啦啦啦观看免费观看视频高清| 午夜久久久久精精品| 欧美精品啪啪一区二区三区| 国产主播在线观看一区二区| 最近最新中文字幕大全免费视频| 午夜福利成人在线免费观看| 岛国在线观看网站| 亚洲美女视频黄频| 一区二区三区免费毛片| 久久人妻av系列| 日韩中文字幕欧美一区二区| 老汉色av国产亚洲站长工具| 美女 人体艺术 gogo| 国产欧美日韩精品一区二区| 国产伦精品一区二区三区四那| 亚洲成av人片在线播放无| av在线蜜桃| 中文字幕人妻丝袜一区二区| 国产精品自产拍在线观看55亚洲| 欧美乱妇无乱码| 国产69精品久久久久777片| 一卡2卡三卡四卡精品乱码亚洲| 最近最新中文字幕大全电影3| 无限看片的www在线观看| 亚洲中文字幕一区二区三区有码在线看| 亚洲 国产 在线| 搡女人真爽免费视频火全软件 | 黄色成人免费大全| 精品福利观看| 好男人电影高清在线观看| 女生性感内裤真人,穿戴方法视频| 亚洲av二区三区四区| 成人国产综合亚洲| 此物有八面人人有两片| 日本 av在线| 18禁在线播放成人免费| 一区二区三区免费毛片| 国产av一区在线观看免费| 国产精品一区二区三区四区免费观看 | 青草久久国产| www.熟女人妻精品国产| 国内精品一区二区在线观看| 亚洲精品成人久久久久久| 99国产综合亚洲精品| 夜夜夜夜夜久久久久| 国产精品亚洲美女久久久| 久久欧美精品欧美久久欧美| 51午夜福利影视在线观看| 夜夜躁狠狠躁天天躁| 久久精品夜夜夜夜夜久久蜜豆| 操出白浆在线播放| av在线蜜桃| 久久天躁狠狠躁夜夜2o2o| 国产激情偷乱视频一区二区| 午夜影院日韩av| 成人午夜高清在线视频| 亚洲最大成人手机在线| 人妻丰满熟妇av一区二区三区| 欧美不卡视频在线免费观看| 色吧在线观看| 在线观看日韩欧美| 免费人成在线观看视频色| 一个人看的www免费观看视频| 国产爱豆传媒在线观看| 亚洲av美国av| 免费大片18禁| 毛片女人毛片| 俄罗斯特黄特色一大片| 国产欧美日韩一区二区精品| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 成人无遮挡网站| 法律面前人人平等表现在哪些方面| 午夜精品一区二区三区免费看| 国产99白浆流出| 亚洲,欧美精品.| 久久午夜亚洲精品久久| 国产又黄又爽又无遮挡在线| 中文亚洲av片在线观看爽| 18禁国产床啪视频网站| 国产精品久久视频播放| 高清在线国产一区| 黄色日韩在线| 日本 av在线| 成人鲁丝片一二三区免费| 麻豆久久精品国产亚洲av| 亚洲精品久久国产高清桃花| 亚洲色图av天堂| 一个人免费在线观看的高清视频| 美女被艹到高潮喷水动态| 欧美午夜高清在线| 亚洲精品乱码久久久v下载方式 | 夜夜躁狠狠躁天天躁| 他把我摸到了高潮在线观看| 国产高清videossex| 亚洲精品成人久久久久久| 亚洲成a人片在线一区二区| 91久久精品电影网| 欧美色视频一区免费| 亚洲国产高清在线一区二区三| 欧美激情在线99| 欧美乱色亚洲激情| 深爱激情五月婷婷| 国产亚洲精品综合一区在线观看| 国产一区二区在线av高清观看| 国产成人a区在线观看| 日本 av在线| 亚洲av不卡在线观看| 亚洲成人精品中文字幕电影| 久久久国产成人精品二区| 日本五十路高清| 欧美日韩一级在线毛片| 悠悠久久av| 嫁个100分男人电影在线观看| 欧美最新免费一区二区三区 | 欧美乱妇无乱码| 久久香蕉精品热| 天天躁日日操中文字幕| 在线十欧美十亚洲十日本专区| 露出奶头的视频| 真人做人爱边吃奶动态| 久久久久久国产a免费观看| 99精品在免费线老司机午夜| 天堂√8在线中文| 中文字幕人妻熟人妻熟丝袜美 | 一级毛片女人18水好多| 免费电影在线观看免费观看| 婷婷丁香在线五月| 国产成人系列免费观看| 狂野欧美激情性xxxx| 欧美又色又爽又黄视频| 美女 人体艺术 gogo| 久久精品人妻少妇| 午夜激情福利司机影院| 最好的美女福利视频网| 亚洲精品456在线播放app | 国产爱豆传媒在线观看| 在线免费观看不下载黄p国产 | 国内精品美女久久久久久| 亚洲精品久久国产高清桃花| 欧美三级亚洲精品| 一本一本综合久久| av在线蜜桃| 国产淫片久久久久久久久 | www日本黄色视频网| 亚洲一区高清亚洲精品| 啦啦啦免费观看视频1| 女人被狂操c到高潮| 无遮挡黄片免费观看| 性欧美人与动物交配| 日韩成人在线观看一区二区三区| 亚洲av不卡在线观看| 欧美bdsm另类| 99久久成人亚洲精品观看| 国产主播在线观看一区二区| 欧美大码av| 欧美黄色片欧美黄色片| 国产精品电影一区二区三区| 精品久久久久久久毛片微露脸| 亚洲国产欧美网| 国产精品自产拍在线观看55亚洲| 老司机深夜福利视频在线观看| 一夜夜www| 岛国在线观看网站| 国模一区二区三区四区视频| 麻豆成人av在线观看| 久久久久九九精品影院| 不卡一级毛片| 欧美日韩福利视频一区二区| 97碰自拍视频| 国产亚洲欧美98| 亚洲欧美日韩无卡精品| 精品日产1卡2卡| 精华霜和精华液先用哪个| 成人av在线播放网站| 黄色成人免费大全| 精品国产美女av久久久久小说| 亚洲avbb在线观看| 老司机深夜福利视频在线观看| 精品乱码久久久久久99久播| a在线观看视频网站| 免费观看精品视频网站| 国产成人系列免费观看| 国产精品久久久久久人妻精品电影| a级一级毛片免费在线观看| 精品国产三级普通话版| 岛国在线免费视频观看| 久久精品国产99精品国产亚洲性色| 香蕉久久夜色| 中文字幕人妻丝袜一区二区| 最新在线观看一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 国产高清视频在线观看网站| 久久草成人影院| 久久久久久人人人人人| 69av精品久久久久久| 99热精品在线国产| 一进一出好大好爽视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产真人三级小视频在线观看| 日韩欧美在线二视频| 亚洲国产精品999在线| 99久国产av精品| 久久这里只有精品中国| 久久国产精品人妻蜜桃| 十八禁网站免费在线| 小蜜桃在线观看免费完整版高清| 国产精品女同一区二区软件 | 国产亚洲精品综合一区在线观看| 国产高清有码在线观看视频| 免费av不卡在线播放| www.www免费av| 动漫黄色视频在线观看| 操出白浆在线播放| 91av网一区二区| 日本在线视频免费播放| 亚洲av免费在线观看| 尤物成人国产欧美一区二区三区| 十八禁网站免费在线| 中文字幕高清在线视频| 午夜亚洲福利在线播放| 一进一出好大好爽视频| 伊人久久大香线蕉亚洲五| 亚洲av免费高清在线观看| 一区二区三区激情视频| 性色avwww在线观看| 小说图片视频综合网站| 欧美zozozo另类| 国产精品一区二区三区四区免费观看 | 国产精品久久久久久精品电影| 在线观看日韩欧美| 又黄又爽又免费观看的视频| 99在线人妻在线中文字幕| 黑人欧美特级aaaaaa片| ponron亚洲| 欧美又色又爽又黄视频| 国产亚洲精品久久久com| 亚洲av二区三区四区| 成人午夜高清在线视频| 人人妻人人看人人澡| 久久久久久大精品| 成人欧美大片| 亚洲精品粉嫩美女一区| 高清在线国产一区| 婷婷六月久久综合丁香| 熟女少妇亚洲综合色aaa.| 亚洲乱码一区二区免费版| 国产成人aa在线观看| 国产精品综合久久久久久久免费| 久久久久国内视频| 少妇熟女aⅴ在线视频| 18禁在线播放成人免费| 久久精品国产综合久久久| 欧美性猛交黑人性爽| 国产伦精品一区二区三区四那| 99在线人妻在线中文字幕| 夜夜躁狠狠躁天天躁| 国产成年人精品一区二区| 国产高清三级在线| 精品国产三级普通话版| 精品99又大又爽又粗少妇毛片 | xxx96com| 日本 av在线| 国产精品久久电影中文字幕| 欧美bdsm另类| 精品99又大又爽又粗少妇毛片 | 老司机深夜福利视频在线观看| 我要搜黄色片| www日本在线高清视频| 一个人观看的视频www高清免费观看| 搞女人的毛片| 日本a在线网址| 国产精品乱码一区二三区的特点| 国产一区二区在线av高清观看| 成人一区二区视频在线观看| 国产成年人精品一区二区| 精品久久久久久久久久免费视频| 亚洲av免费在线观看| 久久久久久久亚洲中文字幕 | 免费大片18禁| 欧美乱色亚洲激情| 国产伦人伦偷精品视频| 国产成+人综合+亚洲专区| 久久久久久大精品| 在线播放无遮挡| 免费看美女性在线毛片视频| 欧美乱码精品一区二区三区| 亚洲精品国产精品久久久不卡| 国产伦精品一区二区三区四那| 国产精品日韩av在线免费观看| 美女被艹到高潮喷水动态| 亚洲国产欧美人成| 小蜜桃在线观看免费完整版高清| 一个人免费在线观看的高清视频| 波多野结衣巨乳人妻| 久久久精品欧美日韩精品| 人人妻人人看人人澡| 亚洲性夜色夜夜综合| 又黄又粗又硬又大视频| 我要搜黄色片| 欧美日韩瑟瑟在线播放| 亚洲欧美日韩卡通动漫| h日本视频在线播放| 国产精品 国内视频| 精品不卡国产一区二区三区| 午夜老司机福利剧场| 亚洲电影在线观看av| 最近在线观看免费完整版| 两个人的视频大全免费| bbb黄色大片| 人人妻人人澡欧美一区二区| 日本黄色片子视频| www日本黄色视频网| 国产久久久一区二区三区| 宅男免费午夜| 一级毛片女人18水好多| 丰满乱子伦码专区| 看黄色毛片网站| 琪琪午夜伦伦电影理论片6080| 亚洲天堂国产精品一区在线| 搡女人真爽免费视频火全软件 | 校园春色视频在线观看| 国产三级在线视频| 熟女人妻精品中文字幕| av在线蜜桃| 99精品久久久久人妻精品| 不卡一级毛片| 好男人在线观看高清免费视频| 国产精品久久久久久久电影 | 国产99白浆流出| 亚洲av免费高清在线观看| 亚洲专区国产一区二区| av视频在线观看入口| 色综合婷婷激情| 免费看美女性在线毛片视频| 色综合婷婷激情| 国产成人福利小说| 可以在线观看毛片的网站| 色综合站精品国产| 一级毛片高清免费大全| 色综合站精品国产| 欧美乱色亚洲激情| 女人十人毛片免费观看3o分钟| 亚洲av五月六月丁香网| 一个人观看的视频www高清免费观看| 日韩国内少妇激情av| 亚洲aⅴ乱码一区二区在线播放| 久久久久国产精品人妻aⅴ院| 狠狠狠狠99中文字幕| 国产亚洲精品一区二区www| 天美传媒精品一区二区| 欧美日本视频| 乱人视频在线观看| 亚洲午夜理论影院| a级毛片a级免费在线| 久久精品国产亚洲av涩爱 | 美女黄网站色视频| 精品国内亚洲2022精品成人| 舔av片在线| 国产精品国产高清国产av| 久久久久九九精品影院| 国产爱豆传媒在线观看| 成人鲁丝片一二三区免费| 亚洲美女视频黄频| 国产成人av激情在线播放| 18禁美女被吸乳视频| 亚洲一区二区三区色噜噜| 国产精品久久电影中文字幕| 国产精品自产拍在线观看55亚洲| 精品一区二区三区视频在线观看免费| 成人亚洲精品av一区二区| a级毛片a级免费在线| 精品国产亚洲在线| 精品久久久久久成人av| 丰满乱子伦码专区| 精品国产亚洲在线| 亚洲久久久久久中文字幕| 成人特级av手机在线观看| 欧美日韩综合久久久久久 | 亚洲av第一区精品v没综合| 日韩人妻高清精品专区| 国产 一区 欧美 日韩| 一个人免费在线观看电影| 精品电影一区二区在线| 在线播放国产精品三级| 国产黄色小视频在线观看| 国产乱人视频| 国产高清三级在线| 免费搜索国产男女视频| 国产野战对白在线观看| 国产激情欧美一区二区| 高清毛片免费观看视频网站| 日韩 欧美 亚洲 中文字幕| 在线播放无遮挡| 亚洲国产色片| 亚洲成人中文字幕在线播放| 午夜免费男女啪啪视频观看 |