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

    一種基于SARIMA的經(jīng)驗(yàn)對流層延遲模型

    2018-04-08 01:47:52劉中流甘祥前梁春麗
    測繪通報(bào) 2018年3期
    關(guān)鍵詞:對流層電離層差分

    劉中流,任 超,甘祥前,梁春麗

    (1. 桂林理工大學(xué)測繪地理信息學(xué)院,廣西 桂林 541004; 2. 廣西空間信息與測繪重點(diǎn)實(shí)驗(yàn)室,廣西 桂林 541004)

    GNSS精密定位技術(shù)已廣泛應(yīng)用于測量改正的各個(gè)領(lǐng)域。然而GNSS信號穿過大氣層時(shí)會引起大氣延遲誤差,這往往對定位精度會產(chǎn)生很大的影響。大氣延遲誤差主要分為電離層延遲誤差和對流層延遲誤差。由于電離層對GNSS信號表現(xiàn)出色散介質(zhì)特性,故產(chǎn)生的電離層延遲與GNSS信號頻率相關(guān),因而目前在GNSS常規(guī)測量中常采用雙頻觀測方法來消除電離層延遲的影響,且在一定范圍內(nèi),其對電離層延遲的改正可達(dá)90%以上[1]。然而對流層對GNSS信號表現(xiàn)出了非色散介質(zhì)特性,因此雙頻觀測方法不適用于對流層延遲改正。對對流層延遲改正的常用方法有外部修正法、參數(shù)估計(jì)法和模型改正法[2-3],它們可以通過采集或運(yùn)算的方法估計(jì)出ZTD值,再代入GNSS觀測方程中便可進(jìn)行延遲誤差改正[4]。

    改進(jìn)對流層模型一直是減弱對流層延遲影響的有效方法[5-6],而近年來,部分學(xué)者在預(yù)測估計(jì)ZTD的領(lǐng)域進(jìn)行了研究,并獲得了一定的成效:王勇等使用BP神經(jīng)網(wǎng)絡(luò)的方法,在試驗(yàn)中取得了76%的預(yù)測值與真實(shí)值偏差(bias)在3 cm以內(nèi)的成果[7];李劍鋒等在試驗(yàn)中使用組合模型的方法,使得試驗(yàn)中的平均預(yù)測值精度(均方根誤差,下同)達(dá)到了14.88 mm[8];呂慧珠等使用譜分析和AR補(bǔ)償結(jié)合的方法,使得試驗(yàn)中的平均預(yù)測值精度達(dá)到了19.36 mm[9]。此類方法也可以較好地計(jì)算出對流層延遲,并代入觀測方程中進(jìn)行誤差改正。許多學(xué)者將自回歸移動平均模型(ARIMA)應(yīng)用到電離層電子含量的預(yù)測中,并取得了較好的成果[10-12],而在對流層預(yù)測中對其的應(yīng)用卻很少有。ARIMA模型對非平穩(wěn)序列具有較好的預(yù)測效果,且適用范圍廣,但卻容易在捕捉序列線性與非線性復(fù)合性特征時(shí)出現(xiàn)精度下降[13]。針對上述情況,本文將ARIMA模型應(yīng)用于ZTD預(yù)測,并引入一種基于季節(jié)性自回歸移動平均模型(SARIMA)的對流層延遲估計(jì)新方法,該方法可分別針對不同的原始ZTD序列進(jìn)行參數(shù)構(gòu)建,然后進(jìn)行預(yù)測以得到最終結(jié)果;針對我國不同地區(qū)的ZTD值進(jìn)行預(yù)測,并統(tǒng)計(jì)分析它們的預(yù)測精度。

    1 模型介紹

    1.1 自回歸移動平均模型(ARIMA)理論

    ARIMA是一個(gè)可以衡量內(nèi)生變量與其滯后性與改變量關(guān)系的系統(tǒng)模型,是對移動平均模型(ARMA)的擴(kuò)展[14]。

    ARIMA(p,d,q)的實(shí)質(zhì)是ARMA(p,q)的d階單整(d次差分),其效果是將一個(gè)非平穩(wěn)的時(shí)間序列轉(zhuǎn)化為平穩(wěn)時(shí)間序列[15]。設(shè)Xt為等時(shí)間間隔的ZTD序列,則可將ARMA(p,q)模型表達(dá)為

    φp(B)Xt=θq(B)εt

    (1)

    式中,B為滯后算子;εt為白噪聲序列;φp(B)為p階自回歸系數(shù)多項(xiàng)式,定義其相應(yīng)特征方程為φp(x);θq(B)為q階滑動平均系數(shù)多項(xiàng)式,定義其相應(yīng)特征方程為θq(x)。

    當(dāng)原ZTD序列{Xt}不平穩(wěn)時(shí),可考慮通過式(10)進(jìn)行差分

    (2)

    式中,u為差分算子;d為差分階數(shù)。{Xt}在經(jīng)過d次差分后變?yōu)槠椒€(wěn)序列,即可建立ARIMA(p,d,q)模型的一般形式

    (3)

    1.2 季節(jié)性自回歸移動平均模型(SARIMA)理論

    SARIMA模型又是對ARIMA模型的拓展。雖然ARIMA(p,d,q)模型已經(jīng)可以解決許多時(shí)間序列的擬合預(yù)測問題,但在實(shí)際操作中,又另有許多時(shí)間序列的變化具有明顯的周期性特征,我們稱之為季節(jié)性時(shí)間序列。由于時(shí)間序列的周期性變化,它們在不同周期內(nèi)的同時(shí)期數(shù)據(jù)將處于一個(gè)相當(dāng)?shù)乃剑藭r(shí)可以恰當(dāng)引入乘積季節(jié)性SARIMA(p,d,q)×(P,D,Q)模型[16]。當(dāng)序列間的長期趨勢、季節(jié)效應(yīng)和隨機(jī)波動間不是簡單的相加關(guān)系時(shí),可以使用ARIMA模型提取短期相關(guān)性,用周期步長為S的SARIMA模型提取季節(jié)周期性,再采用二者的乘積組合模型對ZTD序列進(jìn)行預(yù)測,從而得到乘積季節(jié)模型的完整結(jié)構(gòu)為

    (4)

    圖1 SARIMA模型進(jìn)行ZTD預(yù)測的基本流程

    2 試驗(yàn)分析

    2.1 數(shù)據(jù)準(zhǔn)備

    本文分別選取IGS中心發(fā)布的長春站(chan)和上海佘山站(shao)的2016年年積日1—30 d和92—121 d的ZTD數(shù)據(jù)作為本次分析的原始ZTD數(shù)據(jù),數(shù)據(jù)時(shí)間分辨率變?yōu)? h。采用ARIMA和SARIMA模型分別對它們進(jìn)行預(yù)測分析,并評定其精度。預(yù)測分析原則為,選取這兩個(gè)時(shí)段前29 d逐小時(shí)觀測的ZTD數(shù)據(jù)作為訓(xùn)練數(shù)據(jù),并對未來24個(gè)ZTD值逐小時(shí)地進(jìn)行預(yù)測。

    2.2 精度指標(biāo)

    實(shí)現(xiàn)對各個(gè)ZTD分量序列的SARIMA模型構(gòu)建后,選取前29 d的數(shù)據(jù)進(jìn)行擬合分析,并對后24 h的數(shù)據(jù)進(jìn)行預(yù)測。本文引入精度指標(biāo):均方根誤差(root mean square error,RMSE)、平均絕對誤差(mean absolute error,MAE)、平均相對誤差(mean relative error,MRE),它們的數(shù)學(xué)表達(dá)式分別為

    (5)

    (6)

    (7)

    2.3 建模分析

    在進(jìn)行擬合預(yù)測前,應(yīng)對SARIMA模型(或ARIMA模型)進(jìn)行定階。此處以上海佘山站1—29 d 原始ZTD數(shù)據(jù)為例進(jìn)行模型構(gòu)建分析,此時(shí)段數(shù)據(jù)如圖2所示。

    圖2 上海佘山站年積日1—29 d ZTD數(shù)據(jù)

    通過觀察易于發(fā)現(xiàn),該序列波幅在趨勢上下波動且無明顯增大或減小傾向,應(yīng)屬于趨勢平穩(wěn)(trend stationary),此時(shí)再引入ADF單位根檢驗(yàn)法,判斷序列是否嚴(yán)平穩(wěn)。ADF檢驗(yàn)結(jié)果見表1。

    表1 ZTD序列ADF檢驗(yàn)結(jié)果

    結(jié)果顯示,原序列即使在10%的置信水平上也不能拒絕單位根假設(shè)。此時(shí)需要對原序列進(jìn)行一階差分,并再次進(jìn)行ADF檢驗(yàn),結(jié)果見表2。

    表2 一階差分ZTD序列ADF檢驗(yàn)結(jié)果

    此時(shí)結(jié)果顯示,即使在1%的置信水平下,也可以拒絕單位根假設(shè),代表序列平穩(wěn),故d值取1。下一步應(yīng)進(jìn)行相關(guān)性分析,作出該一階差分序列滯后36階的ACF、PACF函數(shù)如圖3所示。

    由圖3可知,其ACF、PACF函數(shù)均表現(xiàn)出拖尾性,滿足AR(p)、MA(q)過程。此處PACF值三步截尾,一階滯后處PACF值較為顯著,且24階滯后處再次上升至一個(gè)相對峰值,隨后更高階滯后的自相關(guān)系數(shù)明顯減小,此類型表明AR(3)×SAR(24)基本符合該序列(此處SAR(24)僅為計(jì)算機(jī)的取向概念,告知計(jì)算機(jī)季節(jié)組中的每24 h為一個(gè)周期,而該模型應(yīng)正確描述為p=3,P=1),同理可得q=8,P=1。故可基本選定SARIMA(3,1,1)×(8,0,1)模型。但此時(shí)定階的模型僅代表它能很大程度上迫近原序列,并不代表就是最優(yōu)模型,此時(shí)可選取該模型已確認(rèn)參數(shù)的一定鄰域內(nèi)的正整數(shù)替換原參數(shù),并一起進(jìn)行AIC準(zhǔn)則運(yùn)算,選擇使得AIC最小的參數(shù)組合確認(rèn)為最優(yōu)模型參數(shù)。本文通過AIC驗(yàn)證最終確認(rèn)模型為SARIMA(3,1,1)×(8,0,1),并以此進(jìn)行回歸,再進(jìn)行殘差白噪聲檢驗(yàn),殘差A(yù)CF、PACF函數(shù)圖如圖4所示。

    圖3 一階差分ZTD序列ACF、PACF圖

    顯而易見,殘差為白噪聲序列,故可認(rèn)為該模型可靠,可進(jìn)一步用于ZTD數(shù)據(jù)預(yù)測。同理,其他序列也以此過程進(jìn)行模型定階,并作預(yù)測分析。

    2.4 預(yù)測分析

    對所準(zhǔn)備的原始序列進(jìn)行預(yù)測分析,最終得到兩模型預(yù)測精度對比分別見表3、表4、表5。

    圖4 殘差序列ACF、PACF圖

    表4 兩模型預(yù)測成果MAE對比 mm

    表5 兩模型預(yù)測成果MRE (%)

    本文對上述各表的預(yù)測結(jié)果進(jìn)行分析,并總結(jié)得到以下結(jié)論:

    (1) 加入季節(jié)性分析后,對ZTD序列的預(yù)測精度有較好的提高,對年積日1—30 d的數(shù)據(jù)預(yù)測結(jié)果顯示,平均RMSE降低了27.7%,對年積日92—121 d的數(shù)據(jù)預(yù)測結(jié)果總體精度較前者雖有所下降,但同期平均RMSE也降低了25.7%。

    (2) SARIMA模型預(yù)測結(jié)果顯示,平均偏差不大,且占原值比很小,年積日92—121 d的數(shù)據(jù)預(yù)測結(jié)果的平均MRE僅1%左右,年積日1—30 d的數(shù)據(jù)預(yù)測結(jié)果平均MRE僅0.15%。

    (3) SARIMA模型對ZTD數(shù)據(jù)的總體估計(jì)精度較領(lǐng)域內(nèi)當(dāng)前的其他模型(前文提及)的估計(jì)精度要高,總體平均RMSE降低了19.9%。

    (4) 在部分時(shí)期、地區(qū)出現(xiàn)了預(yù)測精度降低的情況,由于對流層延遲受濕度影響較大,故本文猜測與同期氣象變化存在關(guān)系,但不于此處進(jìn)行研究。

    3 結(jié) 論

    本文綜合分析了中國長春、上海地區(qū)的ZTD變化,并分別以ARIMA和SARIMA模型,對兩個(gè)地區(qū)不同時(shí)期的ZTD序列進(jìn)行了短期(24 h)預(yù)測。綜合不同地區(qū)、時(shí)段的預(yù)測成果,結(jié)論如下:

    (1) ARIMA和SARIMA模型對ZTD序列均有較好的預(yù)測效果,而由于ZTD序列存在一定的季節(jié)性特征,故總體而言加入季節(jié)性分析后預(yù)測精度均有較大提高。

    (2) 本文在定階試驗(yàn)時(shí)發(fā)現(xiàn),存在個(gè)別ZTD序列不具有較強(qiáng)的季節(jié)性特征的情況,此時(shí)不能盲目加入季節(jié)性影響,否則將降低模型精度。

    (3) SARIMA模型對ZTD序列的預(yù)測精度較高,可滿足許多ZTD數(shù)值估計(jì)的需求,具有一定的實(shí)用意義。

    (4) 關(guān)于SARIMA模型在不同地區(qū)、時(shí)期對ZTD序列的預(yù)測的應(yīng)用效果需要進(jìn)一步研究。

    參考文獻(xiàn):

    [1] 李征航,黃勁松.GPS測量與數(shù)據(jù)處理[M].武漢:武漢大學(xué)出版社,2005.

    [2] 趙章明,馮徑,洪亮.衛(wèi)星定位中對流層延遲模型對比分析[J].測繪通報(bào),2016(11):18-21.

    [3] 李昭,邱衛(wèi)寧,邱蕾,等.幾種對流層延遲改正模型的分析與比較[J].測繪通報(bào),2009(7):16-18.

    [4] 丁曉光.對流層延遲改正在GPS數(shù)據(jù)處理中的應(yīng)用與研究[D].西安:長安大學(xué),2009.

    [5] 姚宜斌,胡羽豐,余琛.一種改進(jìn)的全球?qū)α鲗犹祉斞舆t模型[J].測繪學(xué)報(bào),2015,44(3):242-249.

    [6] 姚宜斌,張豹,嚴(yán)鳳,等.兩種精化的對流層延遲改正模型[J].地球物理學(xué)報(bào),2015,58(5):1492-1501.

    [7] 王勇,張立輝,楊晶.基于BP神經(jīng)網(wǎng)絡(luò)的對流層延遲預(yù)測研究[J].大地測量與地球動力學(xué),2011,31(3):134-137.

    [8] 李劍鋒,王永前,郭俊元.預(yù)測模型在對流層延遲計(jì)算中的應(yīng)用研究[J].測繪科學(xué)技術(shù)學(xué)報(bào),2015,32(5):450-454.

    [9] 呂慧珠,黃文德,聞德保.一種基于頻譜分析和AR補(bǔ)償?shù)膶α鲗友舆t預(yù)報(bào)模型[J].大地測量與地球動力學(xué),2015,35(2):283-286.

    [10]劉軍,柴洪洲,劉先冬,等.基于ARIMA(p,1,1)的電離層預(yù)報(bào)模型[J].大地測量與地球動力學(xué),2010,30(3):79-82.

    [11]劉立龍,陳軍,黃良珂,等.基于小波-ARIMA電離層短期總電子含量預(yù)報(bào)[J].桂林理工大學(xué)學(xué)報(bào),2016,36(2):294-299.

    [12]張小紅,任曉東,吳風(fēng)波,等.自回歸移動平均模型的電離層總電子含量短期預(yù)報(bào)[J].測繪學(xué)報(bào),2014,43(2):118-124.

    [13]朱幫助,林健.基于ARIMA和LSSVM的非線性集成預(yù)測模型[J].數(shù)學(xué)的實(shí)踐與認(rèn)識,2009,39(12):34-40.

    [14]ZHANG G P.Time Series Forecasting Using a Hybrid ARIMA and Neural Network Model[J].Neurocomputing,2003,50(1):159-175.

    [15]李子奈.計(jì)量經(jīng)濟(jì)學(xué)[M].北京:高等教育出版社,2000.

    [16]張婷.CPI的 SARIMA 模型與X-12季節(jié)調(diào)整模型對比預(yù)測分析[J].經(jīng)濟(jì)問題,2014(12):37-41.

    猜你喜歡
    對流層電離層差分
    一種電離層TEC格點(diǎn)預(yù)測模型
    郴州地區(qū)對流層頂氣候概況
    數(shù)列與差分
    Kalman濾波估算電離層延遲的一種優(yōu)化方法
    電離層對中高軌SAR影響機(jī)理研究
    實(shí)時(shí)干涉測量中對流層延遲與鐘差精修正建模
    載人航天(2016年4期)2016-12-01 06:56:24
    成都地區(qū)2005~2015年對流層NO2柱濃度趨勢與時(shí)空分布
    Linux Shell語言在電離層解算中的應(yīng)用
    基于差分隱私的大數(shù)據(jù)隱私保護(hù)
    相對差分單項(xiàng)測距△DOR
    太空探索(2014年1期)2014-07-10 13:41:50
    国产亚洲精品久久久久久毛片| 最新中文字幕久久久久| 免费电影在线观看免费观看| 午夜老司机福利剧场| 欧美成人免费av一区二区三区| 国产伦精品一区二区三区四那| 欧美最新免费一区二区三区 | 国语自产精品视频在线第100页| 日本与韩国留学比较| 国产一区二区激情短视频| 亚洲成av人片免费观看| 真实男女啪啪啪动态图| 国产99白浆流出| 国产欧美日韩一区二区三| 免费看美女性在线毛片视频| 国产精品久久久久久久久免 | 亚洲av五月六月丁香网| 国产主播在线观看一区二区| 69av精品久久久久久| 亚洲人成伊人成综合网2020| 国产精品久久电影中文字幕| 中出人妻视频一区二区| 午夜福利在线在线| 国产老妇女一区| 在线观看免费视频日本深夜| 中文资源天堂在线| 激情在线观看视频在线高清| 久久精品国产自在天天线| 精品久久久久久久毛片微露脸| 偷拍熟女少妇极品色| 最近视频中文字幕2019在线8| 国产黄a三级三级三级人| 丰满乱子伦码专区| 制服丝袜大香蕉在线| 亚洲av熟女| 久久久色成人| 成年免费大片在线观看| 黄色片一级片一级黄色片| 欧美成狂野欧美在线观看| 69av精品久久久久久| 搡老妇女老女人老熟妇| 88av欧美| 国产一区二区亚洲精品在线观看| 亚洲精品国产精品久久久不卡| 成年女人永久免费观看视频| 亚洲av日韩精品久久久久久密| 欧美一区二区亚洲| 亚洲黑人精品在线| 69av精品久久久久久| 午夜老司机福利剧场| 精品久久久久久久毛片微露脸| 久久亚洲精品不卡| 3wmmmm亚洲av在线观看| 亚洲18禁久久av| 叶爱在线成人免费视频播放| 精品免费久久久久久久清纯| 中亚洲国语对白在线视频| 99热6这里只有精品| 欧美成人免费av一区二区三区| 午夜视频国产福利| 一本精品99久久精品77| 搡老熟女国产l中国老女人| 日本黄大片高清| 国产精品av视频在线免费观看| or卡值多少钱| 成人无遮挡网站| 日本黄色片子视频| 国产精品久久电影中文字幕| 欧美绝顶高潮抽搐喷水| 两个人的视频大全免费| 午夜福利欧美成人| 哪里可以看免费的av片| 亚洲熟妇中文字幕五十中出| 最近最新中文字幕大全免费视频| 久久精品91无色码中文字幕| 51国产日韩欧美| 亚洲精品456在线播放app | 国产午夜福利久久久久久| 动漫黄色视频在线观看| 国产男靠女视频免费网站| 久久精品国产清高在天天线| 欧美性猛交╳xxx乱大交人| 亚洲精华国产精华精| 三级毛片av免费| 蜜桃久久精品国产亚洲av| 国产成人a区在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 动漫黄色视频在线观看| 亚洲精品日韩av片在线观看 | 日韩亚洲欧美综合| 色综合婷婷激情| 男女午夜视频在线观看| 欧美zozozo另类| 丰满人妻熟妇乱又伦精品不卡| 欧美区成人在线视频| 免费电影在线观看免费观看| 丝袜美腿在线中文| 午夜精品久久久久久毛片777| 成年人黄色毛片网站| 狂野欧美激情性xxxx| av专区在线播放| 动漫黄色视频在线观看| 97人妻精品一区二区三区麻豆| 女生性感内裤真人,穿戴方法视频| 51午夜福利影视在线观看| 久久九九热精品免费| 少妇人妻精品综合一区二区 | 亚洲国产欧洲综合997久久,| 天堂√8在线中文| 久久久久久久午夜电影| 女人十人毛片免费观看3o分钟| 久久九九热精品免费| 国产精品嫩草影院av在线观看 | 两个人的视频大全免费| 亚洲美女视频黄频| 欧美一区二区亚洲| 国产欧美日韩一区二区精品| 在线播放国产精品三级| 亚洲无线观看免费| 欧美日韩黄片免| 亚洲国产精品sss在线观看| 变态另类丝袜制服| 在线观看午夜福利视频| 精品午夜福利视频在线观看一区| 少妇的逼水好多| 在线天堂最新版资源| 一二三四社区在线视频社区8| 18禁黄网站禁片午夜丰满| 中文资源天堂在线| 日韩欧美一区二区三区在线观看| 51国产日韩欧美| 特级一级黄色大片| 黄片小视频在线播放| 一级作爱视频免费观看| 成人高潮视频无遮挡免费网站| 九九久久精品国产亚洲av麻豆| 久久久久九九精品影院| 18+在线观看网站| 日韩 欧美 亚洲 中文字幕| 99riav亚洲国产免费| 成年女人毛片免费观看观看9| 在线观看舔阴道视频| 人人妻人人澡欧美一区二区| 综合色av麻豆| 十八禁人妻一区二区| 国产69精品久久久久777片| 成人一区二区视频在线观看| 国产黄片美女视频| 亚洲美女黄片视频| 天堂√8在线中文| 18禁裸乳无遮挡免费网站照片| 99久国产av精品| 亚洲国产欧美网| 国内精品美女久久久久久| 久久亚洲精品不卡| 日日干狠狠操夜夜爽| 欧美色欧美亚洲另类二区| 丁香六月欧美| 嫩草影视91久久| 18禁美女被吸乳视频| 尤物成人国产欧美一区二区三区| 又粗又爽又猛毛片免费看| 岛国在线免费视频观看| 琪琪午夜伦伦电影理论片6080| 夜夜夜夜夜久久久久| 亚洲色图av天堂| 国产高清视频在线观看网站| 色av中文字幕| 久久欧美精品欧美久久欧美| 欧美激情在线99| 久久久久久久亚洲中文字幕 | 久久久国产精品麻豆| 蜜桃亚洲精品一区二区三区| 熟女电影av网| 日本免费a在线| 黄色丝袜av网址大全| 亚洲 国产 在线| 欧美zozozo另类| 国内精品久久久久精免费| 男女午夜视频在线观看| 一个人观看的视频www高清免费观看| 天堂影院成人在线观看| 在线视频色国产色| 久久久国产成人免费| 在线观看舔阴道视频| 丰满的人妻完整版| av片东京热男人的天堂| 悠悠久久av| 亚洲av第一区精品v没综合| 亚洲不卡免费看| 亚洲成人久久性| 少妇裸体淫交视频免费看高清| 精品久久久久久久久久久久久| 国产高潮美女av| 男女之事视频高清在线观看| 丁香六月欧美| 国产精品,欧美在线| 国内精品久久久久久久电影| 99久久九九国产精品国产免费| 99在线人妻在线中文字幕| 狂野欧美激情性xxxx| 国产午夜福利久久久久久| 在线观看美女被高潮喷水网站 | 看黄色毛片网站| 精品国产三级普通话版| aaaaa片日本免费| 在线观看免费视频日本深夜| 两个人视频免费观看高清| 啦啦啦观看免费观看视频高清| 丁香欧美五月| 一级黄片播放器| 又紧又爽又黄一区二区| 90打野战视频偷拍视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 日韩欧美精品v在线| 少妇裸体淫交视频免费看高清| 色视频www国产| 人人妻人人看人人澡| 不卡一级毛片| 男女午夜视频在线观看| 一二三四社区在线视频社区8| 给我免费播放毛片高清在线观看| 日韩欧美国产在线观看| 麻豆成人av在线观看| avwww免费| 青草久久国产| 18+在线观看网站| 午夜免费成人在线视频| 成人特级av手机在线观看| 国产精品美女特级片免费视频播放器| 脱女人内裤的视频| 国产一区二区三区视频了| 久久中文看片网| www.www免费av| 淫妇啪啪啪对白视频| 老熟妇仑乱视频hdxx| 69人妻影院| 99热精品在线国产| 亚洲激情在线av| 国产视频内射| 伊人久久大香线蕉亚洲五| 欧美成人一区二区免费高清观看| 欧美日本视频| АⅤ资源中文在线天堂| 亚洲国产中文字幕在线视频| 少妇人妻一区二区三区视频| av女优亚洲男人天堂| 老熟妇乱子伦视频在线观看| 亚洲av一区综合| 99热精品在线国产| 美女高潮喷水抽搐中文字幕| 亚洲不卡免费看| 免费看a级黄色片| 久久九九热精品免费| 色老头精品视频在线观看| 丰满乱子伦码专区| 99精品在免费线老司机午夜| 黄片小视频在线播放| 特大巨黑吊av在线直播| 日本与韩国留学比较| 欧美乱妇无乱码| 国内精品一区二区在线观看| 免费观看人在逋| 窝窝影院91人妻| 久久婷婷人人爽人人干人人爱| 老司机午夜十八禁免费视频| 欧美一区二区精品小视频在线| 一边摸一边抽搐一进一小说| 欧美绝顶高潮抽搐喷水| 综合色av麻豆| 欧美成人免费av一区二区三区| 又紧又爽又黄一区二区| 老汉色av国产亚洲站长工具| 一本久久中文字幕| 亚洲天堂国产精品一区在线| 精品国产超薄肉色丝袜足j| 亚洲精品粉嫩美女一区| 久久久久性生活片| 日本熟妇午夜| 国产免费av片在线观看野外av| 午夜免费激情av| 怎么达到女性高潮| 波多野结衣高清无吗| 国产高清有码在线观看视频| 国产久久久一区二区三区| 全区人妻精品视频| 国产精品日韩av在线免费观看| 国产97色在线日韩免费| av福利片在线观看| 麻豆国产97在线/欧美| 成年女人毛片免费观看观看9| 国产亚洲精品久久久久久毛片| 欧美成人性av电影在线观看| 国产亚洲精品av在线| 久久精品影院6| 国产精品永久免费网站| 少妇丰满av| 国产视频内射| av中文乱码字幕在线| 别揉我奶头~嗯~啊~动态视频| 精品不卡国产一区二区三区| 天天躁日日操中文字幕| 国产精品99久久99久久久不卡| 午夜亚洲福利在线播放| 国产伦在线观看视频一区| 有码 亚洲区| 观看免费一级毛片| 一进一出好大好爽视频| 国产高清videossex| 欧美在线黄色| 成熟少妇高潮喷水视频| 久久精品国产自在天天线| 欧美丝袜亚洲另类 | 日韩国内少妇激情av| 天天添夜夜摸| 国产亚洲精品av在线| 久久精品综合一区二区三区| 国产精品乱码一区二三区的特点| 最近视频中文字幕2019在线8| 窝窝影院91人妻| 五月玫瑰六月丁香| 久久久久久久亚洲中文字幕 | 久久久久九九精品影院| 成人特级黄色片久久久久久久| 老鸭窝网址在线观看| 亚洲最大成人中文| 在线观看美女被高潮喷水网站 | 国产乱人伦免费视频| 搡老妇女老女人老熟妇| 国内精品一区二区在线观看| 91久久精品国产一区二区成人 | 久久久国产精品麻豆| 精品国内亚洲2022精品成人| 国产欧美日韩精品亚洲av| 天堂影院成人在线观看| 99久久精品热视频| 国产在视频线在精品| 99久久精品热视频| 国产aⅴ精品一区二区三区波| 啦啦啦免费观看视频1| 欧美中文日本在线观看视频| 国产男靠女视频免费网站| 在线免费观看不下载黄p国产 | 久久香蕉精品热| 最近最新免费中文字幕在线| 真人做人爱边吃奶动态| 久久精品国产自在天天线| 一本久久中文字幕| 国产成人av激情在线播放| 大型黄色视频在线免费观看| 欧美性猛交╳xxx乱大交人| 日韩欧美在线二视频| 我的老师免费观看完整版| 叶爱在线成人免费视频播放| 在线播放国产精品三级| 久久久久亚洲av毛片大全| 99riav亚洲国产免费| 91在线观看av| 丰满人妻一区二区三区视频av | 欧美+日韩+精品| 丰满人妻一区二区三区视频av | 757午夜福利合集在线观看| 国产高清视频在线观看网站| 91麻豆av在线| 欧美一区二区精品小视频在线| 免费搜索国产男女视频| 国产亚洲精品综合一区在线观看| 97超视频在线观看视频| 校园春色视频在线观看| 国产高清激情床上av| 欧美极品一区二区三区四区| 草草在线视频免费看| 日本在线视频免费播放| 美女黄网站色视频| 99热精品在线国产| or卡值多少钱| 日韩国内少妇激情av| 亚洲七黄色美女视频| www日本在线高清视频| 狠狠狠狠99中文字幕| 日韩国内少妇激情av| 欧美成人性av电影在线观看| 婷婷精品国产亚洲av| 18美女黄网站色大片免费观看| 免费av不卡在线播放| 一区二区三区免费毛片| 精品一区二区三区视频在线观看免费| 成人av在线播放网站| 中文字幕av成人在线电影| 在线免费观看不下载黄p国产 | 老汉色∧v一级毛片| 国产淫片久久久久久久久 | 3wmmmm亚洲av在线观看| 色尼玛亚洲综合影院| 国产精品久久久人人做人人爽| 尤物成人国产欧美一区二区三区| 99精品欧美一区二区三区四区| 中文字幕久久专区| 国产一区二区亚洲精品在线观看| 国产亚洲精品一区二区www| 国产精品女同一区二区软件 | 久久久久久久久中文| 国产精品爽爽va在线观看网站| 亚洲国产日韩欧美精品在线观看 | 国产视频一区二区在线看| 桃红色精品国产亚洲av| ponron亚洲| 51国产日韩欧美| 国产男靠女视频免费网站| 欧美中文日本在线观看视频| 国产高清激情床上av| 精品人妻一区二区三区麻豆 | 日韩欧美在线二视频| 精品国产美女av久久久久小说| 国产色爽女视频免费观看| www.999成人在线观看| 99久久精品一区二区三区| 琪琪午夜伦伦电影理论片6080| netflix在线观看网站| 国产黄色小视频在线观看| 婷婷亚洲欧美| 最后的刺客免费高清国语| 夜夜看夜夜爽夜夜摸| 国产淫片久久久久久久久 | 国产主播在线观看一区二区| xxx96com| 久久久久久大精品| 欧美精品啪啪一区二区三区| 18禁裸乳无遮挡免费网站照片| 亚洲一区二区三区色噜噜| 别揉我奶头~嗯~啊~动态视频| 国产精品久久久久久亚洲av鲁大| 99久国产av精品| 九九久久精品国产亚洲av麻豆| 日本一本二区三区精品| 欧美性猛交黑人性爽| 中文字幕久久专区| 美女高潮喷水抽搐中文字幕| 99久久精品一区二区三区| 最近视频中文字幕2019在线8| 伊人久久精品亚洲午夜| 久久精品91无色码中文字幕| 成人av一区二区三区在线看| 国产日本99.免费观看| 少妇人妻一区二区三区视频| 尤物成人国产欧美一区二区三区| 九色成人免费人妻av| 亚洲国产欧美网| 欧美成人性av电影在线观看| 亚洲 欧美 日韩 在线 免费| 午夜激情欧美在线| 久久精品国产自在天天线| 中文在线观看免费www的网站| 五月玫瑰六月丁香| 在线观看舔阴道视频| 99国产精品一区二区蜜桃av| 999久久久精品免费观看国产| 美女 人体艺术 gogo| 日韩欧美精品v在线| 日本在线视频免费播放| 国产黄a三级三级三级人| 香蕉久久夜色| 啪啪无遮挡十八禁网站| 国产一区二区三区视频了| 中文字幕精品亚洲无线码一区| 欧美av亚洲av综合av国产av| 中文在线观看免费www的网站| 九色成人免费人妻av| 欧美日韩国产亚洲二区| 一区二区三区免费毛片| 在线观看av片永久免费下载| 国产v大片淫在线免费观看| 欧美午夜高清在线| 少妇的逼水好多| 最近最新中文字幕大全免费视频| 亚洲精品在线美女| 夜夜夜夜夜久久久久| 搡女人真爽免费视频火全软件 | www.色视频.com| 国内精品久久久久精免费| 亚洲性夜色夜夜综合| 网址你懂的国产日韩在线| 久久九九热精品免费| 欧美又色又爽又黄视频| 日韩欧美在线二视频| 国产欧美日韩精品亚洲av| 欧美日韩国产亚洲二区| 性色avwww在线观看| 好男人电影高清在线观看| 精品日产1卡2卡| 免费人成视频x8x8入口观看| 老汉色∧v一级毛片| 美女高潮的动态| 色综合亚洲欧美另类图片| 国产精品野战在线观看| 69av精品久久久久久| 欧美日韩黄片免| 日本黄色片子视频| 91麻豆精品激情在线观看国产| 波多野结衣巨乳人妻| 99国产精品一区二区蜜桃av| 国产免费一级a男人的天堂| 欧美色欧美亚洲另类二区| 九九久久精品国产亚洲av麻豆| 可以在线观看毛片的网站| 91在线精品国自产拍蜜月 | 成年免费大片在线观看| 香蕉av资源在线| 小蜜桃在线观看免费完整版高清| 欧美在线一区亚洲| 亚洲精品在线美女| 女人被狂操c到高潮| 他把我摸到了高潮在线观看| 一区二区三区国产精品乱码| 男女视频在线观看网站免费| 一进一出抽搐gif免费好疼| 欧美一级毛片孕妇| 天堂√8在线中文| 亚洲最大成人中文| 国产欧美日韩一区二区精品| 啦啦啦免费观看视频1| 午夜精品一区二区三区免费看| 我要搜黄色片| 999久久久精品免费观看国产| 国产精品亚洲av一区麻豆| 香蕉av资源在线| 9191精品国产免费久久| 99riav亚洲国产免费| 日本五十路高清| 尤物成人国产欧美一区二区三区| 中文字幕熟女人妻在线| 欧美区成人在线视频| 亚洲欧美激情综合另类| 高清日韩中文字幕在线| 亚洲国产日韩欧美精品在线观看 | 男人舔奶头视频| 日本三级黄在线观看| 国产亚洲欧美在线一区二区| 老司机在亚洲福利影院| 可以在线观看的亚洲视频| 国产不卡一卡二| 精品国内亚洲2022精品成人| 亚洲国产精品成人综合色| 国产高清视频在线播放一区| 亚洲久久久久久中文字幕| 国产一区二区激情短视频| 欧美日韩精品网址| 中文字幕精品亚洲无线码一区| 99久久精品热视频| 欧美日本视频| 成年人黄色毛片网站| 搡老岳熟女国产| 免费人成在线观看视频色| 亚洲国产精品成人综合色| 大型黄色视频在线免费观看| 久久久久精品国产欧美久久久| 亚洲性夜色夜夜综合| 久久精品91蜜桃| 蜜桃亚洲精品一区二区三区| 99在线人妻在线中文字幕| 白带黄色成豆腐渣| 亚洲人成伊人成综合网2020| 国产私拍福利视频在线观看| 又黄又粗又硬又大视频| 久久午夜亚洲精品久久| 少妇丰满av| 中文字幕av在线有码专区| 18禁裸乳无遮挡免费网站照片| 亚洲自拍偷在线| 亚洲av成人精品一区久久| 日日干狠狠操夜夜爽| 亚洲专区中文字幕在线| 久久久久精品国产欧美久久久| 99国产综合亚洲精品| 999久久久精品免费观看国产| 亚洲av不卡在线观看| 成年女人永久免费观看视频| 小说图片视频综合网站| 日本黄色片子视频| АⅤ资源中文在线天堂| 网址你懂的国产日韩在线| 免费在线观看成人毛片| 亚洲av美国av| 欧美不卡视频在线免费观看| 很黄的视频免费| 99精品欧美一区二区三区四区| 国产精品国产高清国产av| 国产中年淑女户外野战色| 黄片小视频在线播放| 欧美精品啪啪一区二区三区| 亚洲人与动物交配视频| 丰满人妻一区二区三区视频av | 亚洲av电影在线进入| 国产视频一区二区在线看| 久久精品国产99精品国产亚洲性色| 少妇的逼好多水| 又紧又爽又黄一区二区| www国产在线视频色| 国产精品三级大全| 小说图片视频综合网站| 国产精品久久久久久久久免 | 91久久精品国产一区二区成人 | 亚洲国产色片| 国产精品久久久久久精品电影| 欧美+日韩+精品| 女警被强在线播放| 精品乱码久久久久久99久播| 日本熟妇午夜| 亚洲人成伊人成综合网2020| 日本精品一区二区三区蜜桃| 网址你懂的国产日韩在线| 亚洲一区高清亚洲精品| 九色成人免费人妻av| 欧美区成人在线视频| 91麻豆av在线|