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

    考慮周期性波動因素的中長期空中交通流量預(yù)測

    2015-05-25 03:03:26胡明華張洪海尹嘉男
    西南交通大學(xué)學(xué)報 2015年3期
    關(guān)鍵詞:交通流量周期性季度

    陳 丹, 胡明華, 張洪海, 尹嘉男

    (南京航空航天大學(xué)民航學(xué)院,江蘇南京 211106)

    隨著我國航空運輸業(yè)的飛速發(fā)展,空中交通流量不斷增長,自2006年至今已增長70%左右,年平均增長率達到9%,且未來20年仍將繼續(xù)保持迅猛增長趨勢[1].不斷增長的空中交通需求與有限的空域運行能力之間的不平衡狀況日趨凸顯,空域擁擠與航班延誤等問題頻發(fā),為進一步提升民航服務(wù)品質(zhì),有必要對中長期空中交通流量的總體增長趨勢和周期性波動規(guī)律進行準確預(yù)測,這是實施空中交通流量管理的前提和基礎(chǔ),也是空域管理科學(xué)決策的重要依據(jù).

    中長期空中交通流量預(yù)測的常用方法為時間序列預(yù)測法,主要包括序時平均數(shù)法、指數(shù)平滑法、移動平均法、趨勢預(yù)測法等[2-6].這些方法具有輸入數(shù)據(jù)少、計算量小等特點,但要求時間序列具有穩(wěn)定性,對于非平穩(wěn)時間序列,則需通過差分、取對數(shù)等方法進行適當(dāng)?shù)念A(yù)處理,將其轉(zhuǎn)換為平穩(wěn)時間序列,但預(yù)處理過程可能造成某些重要信息的丟失,降低預(yù)測精度[7].

    為分析惡劣天氣、重大事件等不確定性因素對空中交通流量的影響,將動態(tài)線性模型引入空中交通流量預(yù)測領(lǐng)域,該模型不要求時間序列平穩(wěn)性假設(shè),且可用貝葉斯統(tǒng)計方法求解,不僅可預(yù)測未來發(fā)展趨勢,還能給出相應(yīng)的預(yù)測置信區(qū)間[8-11].

    當(dāng)前預(yù)測方法未考慮氣候、季節(jié)、交通需求等周期性波動因素對空中交通流量的影響[12].本文綜合考慮氣候、季節(jié)、交通需求等周期性波動因素對空中交通流量的影響,通過分析歷史運行數(shù)據(jù),分析空中交通流量中長期總體變化趨勢和周期性波動規(guī)律,建立考慮周期性波動因素的改進動態(tài)線性模型,用于預(yù)測未來某時空范圍內(nèi)的空中交通流量,為戰(zhàn)略流量管理決策提供數(shù)據(jù)支持.

    1 預(yù)測模型

    1.1 一般動態(tài)線性模型

    動態(tài)線性模型是一種特殊的狀態(tài)空間模型,也稱為高斯線性狀態(tài)空間模型.假設(shè)時刻t=0動態(tài)線性模型的初始p維狀態(tài)向量X0服從先驗正態(tài)分布,m0為期望,C0為方差,

    式中:Yt為時刻t空域單元的m維交通流量向量;

    Xt為時刻t交通流的隱含狀態(tài);

    Ft為狀態(tài)向量與觀測向量之間的映射關(guān)系矩陣;

    Gt為相鄰狀態(tài)向量之間的映射關(guān)系矩陣;

    vt(t≥1)和 wt(t≥1)分別為為均值為 0 的方差矩陣Vt(t≥1)和 Wt(t≥1)的獨立的高斯隨機向量,用于描述觀測噪聲和過程噪聲.

    通常將式(2)稱為觀測方程,式(3)稱為狀態(tài)方程,式(1)~(3)循環(huán)遞推,構(gòu)成動態(tài)線性模型[13].

    1.2 線性增長模型

    結(jié)合空中交通流量預(yù)測問題的特點,本文選取一種特殊的動態(tài)線性模型,即線性增長模型[13].此時,m=1,p=2,在標準模型的基礎(chǔ)上,增加了狀態(tài)μt的動態(tài)斜率βt,模型參數(shù)如下:

    對于中長期空中交通流量預(yù)測,μt表示對實際交通流量觀測時序數(shù)據(jù)的估計狀態(tài),βt表示該狀態(tài)的變化率,σμ、σβ為方差,當(dāng) m=1、p=2時,線性增長模型為

    1.3 周期性波動模型

    綜合考慮氣候、季節(jié)、交通需求等周期性波動因素的影響,建立改進模型.假設(shè)時序數(shù)據(jù)的均值為0,即時序數(shù)據(jù)只存在周期性波動.對于均值非0的時序數(shù)據(jù),其變化趨勢項可由1.2節(jié)構(gòu)建的線性增長模型預(yù)測得到.設(shè)波動周期為s,則周期性波動模型(m=1,p=s)的參數(shù)如下:

    Xt為(s-1)維狀態(tài)向量(波動周期為s的周期性波動模型,僅有s-1個自由狀態(tài));

    Ft和 Gt為1×(s-1)維和(s-1)×(s-1)維矩陣;

    W為1×(s-1)維矩陣,表示過程噪聲.

    1.4 預(yù)測精度指標

    為評價模型預(yù)測值與實際值擬合程度的優(yōu)劣,本文選取兩個預(yù)測精度指標,即平均絕對誤差百分比EMAPE和相對均方根誤差百分比ER-RMPSE.這兩個指標廣泛用于空中交通預(yù)測研究[14],

    式中:yt為時刻t的實際觀測值;

    et為時刻t實際觀測值yt與預(yù)測期望值ft之間的偏差.

    由于中長期流量預(yù)測結(jié)果對戰(zhàn)略流量管理具有重要指導(dǎo)意義,因此,要求預(yù)測誤差具有較好的穩(wěn)定性.在統(tǒng)計學(xué)中,常用標準差衡量數(shù)組之間的離散程度,標準差越小,表示數(shù)據(jù)越聚集,穩(wěn)定性越好,反之?dāng)?shù)據(jù)越離散,穩(wěn)定性越差.本文選用預(yù)測誤差的標準差η作為評價預(yù)測性能穩(wěn)定性的指標,

    式中:di為時刻t預(yù)測結(jié)果的相對誤差;

    Emean為各時刻預(yù)測結(jié)果相對誤差的平均值.

    2 預(yù)測流程

    根據(jù)所建模型及求解方法,采用統(tǒng)計分析與計算領(lǐng)域常用的R語言編程工具,對本文提出的預(yù)測方法進行算法實現(xiàn),具體步驟如下:

    (1)分析中長期空中交通流量預(yù)測的各類動靜態(tài)影響因素影響,提煉出具備周期性波動特點的因素;

    (2)建立考慮周期性波動因素的中長期空中交通流量預(yù)測模型,應(yīng)用基于貝葉斯?fàn)顟B(tài)估計與預(yù)測理論的方法求解;

    (3)讀入典型繁忙地區(qū)某空域單元在特定時間范圍內(nèi)的歷史流量觀測時序數(shù)據(jù),并建立觀測向量Yt;

    (4)確定模型參數(shù)Gt、Ft以及波動周期s的值,其中線性增長模型的Gt和Ft見1.2節(jié),改進模型的Gt和Ft見1.3節(jié),波動周期s則取決于用于預(yù)測的觀測數(shù)據(jù)的時間粒度,例如,月交通流量預(yù)測模型中波動周期s=12,季度交通流量預(yù)測模型中波動周期s=4;

    (5)編程設(shè)計R語言程序,對中長期空中交通流量進行預(yù)測,并輸出預(yù)測結(jié)果;

    (6)對預(yù)測結(jié)果進行綜合對比,并分析本文方法的預(yù)測精度.

    3 預(yù)測實例

    選取某地區(qū)2001—2010年交通流量數(shù)據(jù)作為觀測數(shù)據(jù),應(yīng)用線性增長模型和考慮周期性波動因素的改進模型對歷史交通流量時序數(shù)據(jù)建模,并采用貝葉斯?fàn)顟B(tài)估計與預(yù)測方法[15]對模型進行求解,預(yù)測未來兩年該地區(qū)的月交通流量、季度交通流量以及年交通流量,并將預(yù)測結(jié)果與該地區(qū)2011—2012年實際觀測數(shù)據(jù)進行對比分析.

    3.1 月交通流量預(yù)測

    根據(jù)某地區(qū)2001年1月至2010年12月的月交通流量時序數(shù)據(jù),利用1.2節(jié)建立的線性增長模型描述該時序數(shù)據(jù),采用貝葉斯?fàn)顟B(tài)估計和預(yù)測方法對該模型求解,預(yù)測該地區(qū)2011年1月至2012年12月的月交通流量時序數(shù)據(jù)的未來發(fā)展趨勢.預(yù)測結(jié)果如圖1所示.

    圖1 線性增長模型的月交通流量預(yù)測結(jié)果Fig.1 Monthly traffic flow forecasted by the linear growth model

    如圖1所示為根據(jù)2001—2010年的月交通流量觀測數(shù)據(jù)對未來流量變化趨勢的預(yù)測情況.假設(shè)觀測噪聲和過程噪聲均為高斯噪聲,由貝葉斯?fàn)顟B(tài)估計與預(yù)測理論可知,歷史觀測時序數(shù)據(jù)的預(yù)測值同樣服從高斯分布(圖1),圖1中淺色灰線表示20個預(yù)測樣本.

    圖1中,實際時序數(shù)據(jù)具有明顯的周期性波動,但線性增長模型只考慮了時序數(shù)據(jù)的線性發(fā)展趨勢,預(yù)測結(jié)果為一條上升的直線,無法體現(xiàn)時序數(shù)據(jù)的波動性.因此,在線性增長模型的基礎(chǔ)上,考慮周期性波動因素,將線性增長模型作為預(yù)測趨勢項,將周期性波動模型作為波動項,二者相加得到改進的預(yù)測模型,預(yù)測結(jié)果見圖2(s=12).

    圖2 改進模型的月交通流量預(yù)測結(jié)果Fig.2 Monthly traffic flow forecasted by the improved model

    由圖2可見,基于改進模型的月交通流量預(yù)測結(jié)果既能反映時序數(shù)據(jù)的總體發(fā)展趨勢,又能較好體現(xiàn)其周期性波動規(guī)律,與實際觀測時序數(shù)據(jù)更接近.為直觀反映流量數(shù)據(jù)的波動性對預(yù)測結(jié)果的影響,將預(yù)測結(jié)果分解,如圖3所示.

    圖3 改進模型的月交通流量預(yù)測結(jié)果分解Fig.3 Components of the monthly traffic flow forecasted by the improved model

    兩組月交通流量預(yù)測結(jié)果的對比見表1.根據(jù)式(5)和(6)可計算得出線性增長模型預(yù)測結(jié)果的EMAPE和 ER-RMPSE分別為4.00%和5.40%,相對誤差的標準差為3.62%;而改進模型預(yù)測結(jié)果的EMAPE和ER-RMPSE分別為2.62%和3.23%,相對誤差的標準差為1.89%,預(yù)測精度較線性增長模型明顯提高,預(yù)測誤差也更為穩(wěn)定.

    表1 月交通流量預(yù)測結(jié)果對比Tab.1 Comparison of forecast results of monthly traffic flow between models

    3.2 季度交通流量預(yù)測

    根據(jù)2001—2010年的自然季度(春、夏、秋、冬)交通流量時序數(shù)據(jù),運用線性增長模型,預(yù)測2011—2012年的自然季度交通流量,結(jié)果見圖4.

    圖4 線性增長模型的自然季度交通流量預(yù)測結(jié)果Fig.4 Quarterly traffic flow forecasted by the linear growth model

    為準確把握時序數(shù)據(jù)的客觀波動性,采用改進模型對自然季度交通流量進行預(yù)測,波動周期s=4,預(yù)測結(jié)果如圖5所示.為更清晰的顯示改進模型的預(yù)測結(jié)果,圖6分別給出預(yù)測結(jié)果的趨勢項、波動項以及最終預(yù)測結(jié)果.

    圖5 改進模型的自然季度交通流量預(yù)測結(jié)果Fig.5 Quarterly traffic flow forecasted by the improved model

    圖6 改進模型的自然季度交通流量預(yù)測結(jié)果分解Fig.6 Components of the quarterly traffic flow forecasted by the improved model

    自然季度交通流量預(yù)測結(jié)果的對比見表2.根據(jù)1.4節(jié)的預(yù)測精度指標公式可計算得出線性增長模型預(yù)測結(jié)果的EMAPE和ER-RMPSE分別為3.97%和4.94%,相對誤差的標準差為2.93%;而改進模型預(yù)測結(jié)果的 EMAPE和 ER-RMPSE分別為2.91%和3.71%,相對誤差的標準差為2.30%,預(yù)測精度較線性增長模型略有提高,預(yù)測誤差也更為穩(wěn)定.

    利用2001年夏秋季(2001年4月至10月)至2009年冬春季(2009年11月至2010年3月)的航班季度(冬春、夏秋)交通流量時序數(shù)據(jù),用線性增長模型,預(yù)測2010年夏秋季至2012年夏秋季的航班季度交通流量,結(jié)果見圖7.改進模型航班季度交通流量預(yù)測結(jié)果見圖8.

    表2 自然季度交通流量預(yù)測結(jié)果對比Tab.2 Comparison of forecast results of quarterly traffic flow between models

    圖7 線性增長模型的航班季度交通流量預(yù)測結(jié)果Fig.7 Traffic flow of flight season forecasted by the linear growth model

    圖8 改進模型的航班季度交通流量預(yù)測結(jié)果Fig.8 Traffic flow of flight season forecasted by the improved model

    為更清晰地顯示改進模型的預(yù)測結(jié)果,分別給出預(yù)測結(jié)果的趨勢項、波動項以及最終預(yù)測結(jié)果,如圖9所示.

    圖9 改進模型的航班季度交通流量預(yù)測結(jié)果分解Fig.9 Components of traffic flow of flight season forecasted by the improved model

    兩組航班季度交通流量預(yù)測結(jié)果的對比見表3.

    表3 航班季度交通流量預(yù)測結(jié)果對比Tab.3 Comparison of traffic flow of flight season between models

    根據(jù)式(5)和(6)可計算得出線性增長模型預(yù)測結(jié)果的 EMAPE和 ER-RMPSE分別為 21.59%和21.77%,而改進模型預(yù)測結(jié)果的EMAPE和ER-RMPSE分別為6.23%和6.90%,預(yù)測精度較線性增長模型有顯著提高.冬春季航班的執(zhí)行時間為當(dāng)年10月末至翌年3月末,約5個月,夏秋季航班的執(zhí)行時間為當(dāng)年3月末至10月末,約7個月,二者存在2個月的時間差異引起的波動,直接導(dǎo)致不考慮波動性因素的線性增長模型預(yù)測結(jié)果誤差大,此時該模型已不再適用,而改進模型考慮了周期性波動因素,預(yù)測結(jié)果較好.

    3.3 年交通流量預(yù)測

    根據(jù)某地區(qū)2001—2010年的年交通流量時序數(shù)據(jù),用線性增長模型,預(yù)測該地區(qū)2011—2012年的年交通流量,結(jié)果見圖10.

    圖10 線性增長模型的年交通流量預(yù)測結(jié)果Fig.10 Yearly traffic flow forecasted by the linear growth model

    由于年交通流量時序數(shù)據(jù)不存在周期性波動規(guī)律,本文考慮將基于改進模型的月交通流量預(yù)測結(jié)果和自然季度交通流量預(yù)測結(jié)果逐年加和,作為改進模型的年交通流量預(yù)測結(jié)果.3組年交通流量預(yù)測結(jié)果的對比見表4.

    表4 年交通流量預(yù)測結(jié)果對比Tab.4 Comparison of yearly traffic flow between models

    由表4可知,采用改進模型的兩組預(yù)測結(jié)果誤 差要明顯小于線性增長模型的預(yù)測結(jié)果誤差,且前者更穩(wěn)定;由改進模型預(yù)測結(jié)果可見,與季度交通流量時序數(shù)據(jù)的預(yù)測精度結(jié)果相比,月交通流量時序數(shù)據(jù)的預(yù)測精度顯著提升.

    4 結(jié)束語

    考慮到實際空中交通流量受氣候、季節(jié)、交通需求等因素的影響而具有顯著的周期性波動特點,

    本文在一般動態(tài)線性模型的基礎(chǔ)上,建立了考慮周期性波動因素的改進模型用于中長期空中交通流量預(yù)測,并借助貝葉斯?fàn)顟B(tài)估計與預(yù)測方法求解模型.實例分析表明,改進模型不僅能夠反映空中交通流量中長期的總體變化趨勢,而且能夠體現(xiàn)交通流量的周期性波動規(guī)律,因此,得到更為合理的預(yù)測結(jié)果,且相較于現(xiàn)有的未考慮周期性波動時序模型(例如線性增長模型)具有更高的預(yù)測精度和穩(wěn)定性.算例結(jié)果表明,本文方法適用于中長期時間范圍內(nèi)特定空域單元空中交通流量預(yù)測.

    [1]中國民用航空局發(fā)展計劃司.從統(tǒng)計看民航2013[M].北京:中國民航出版社,2013:126-147.

    [2]趙玉環(huán),石新華.基于時間序列的空中交通流量灰預(yù)測模型算法[J].中國民航大學(xué)學(xué)報,2007,25(6):54-57.Z

    HAO Yuhuan,SHI Xinhua.Air traffic flow gray forecast model algorithm based on time series[J].Journal of China Civil Aviation University, 2007,25(6):54-57.

    [3]姜靜逸,韓松臣,王玉婷.新型組合預(yù)測模型在空中交通流量預(yù)測的應(yīng)用[J].中國民航大學(xué)學(xué)報,2009,27(5):4-8.

    JIANG Jingyi,HAN Songchen,WANG Yuting.New combination forecast model in the application of air traffic flow prediction[J]. JournalofChina Civil Aviation University,2009,27(5):4-8.

    [4]?NDER E,KUZU S.Forecasting air traffic volumes using smoothing techniques[J].Journal of Aeronautics and Space Technologies,2014,7(1):65-70.

    [5]趙玉環(huán),郭爽.考慮隨機因素的空中交通流量預(yù)測模型研究[J].中國民航大學(xué)學(xué)報,2008,26(4):59-61.

    ZHAO Yuhuan,GUO Shuang.Research on the forecast method for air traffic flow considering random factors[J].Journal of China Civil Aviation University,2008,26(4):59-61.

    [6]BOUGAS C.Forecasting air passenger traffic flows in canada:an evaluation oftimeseriesmodelsand combinationmethods[D]. [S. l.]:Constantinos Bougas,2013.

    [7]MALLAT S G.A theory for multiresolution signal decomposition:the wavelet representation[J].IEEE Transactions on Pattern Analysis and Machine Intelligence,1989,11(7):674-693.

    [8]MENON P K,SWERIDUK G D,BILIMORIA K D.New approach for modeling,analysis,and control of air traffic flow[J].Journal of Guidance,Control,and Dynamics,2004,27(5):737-744.

    [9]ROY S,SRIDHAR B,VERGHESE G C.An aggregate dynamic stochastic model for an air traffic system[C]∥Proceedings of the 5th Eurocontrol/Federal Aviation Agency Air Traffic Management Research and Development Seminar.Budapest:[s.n.],2003:1-10.

    [10]SRIDHAR B,SONI T,SHETH K,et al.Aggregate flow model for air-traffic management[J].Journal of Guidance,Control,and Dynamics,2006,29(4):992-997.

    [11]SRIDHAR B,CHEN N Y,NG H K.An aggregate sector flow model for air traffic demand forecasting[C]∥ 9th AIAA Aviation Technology,Integration, and Operations Conference(ATIO).[S.l.]:NASA Ames Research Center,2009:1-12.

    [12]張靜,徐肖豪,王飛.天氣季節(jié)性影響的機場到達容量概率分布[J].西南交通大學(xué)學(xué)報,2011,46(1):154-161.

    ZHANG Jing,XU Xiaohao,WANG Fei.Probability distribution for airports'arriving capacity considering seasonal weather effects[J]. Journal of Southwest Jiaotong University,2011,46(1):154-161.

    [13]PETRIS G,PETRONE S,CAMPAGNOLI P.Dynamic linear models with R[M]. New York:Springer-Verlag,2009:41-74.

    [14]COSHALL J T.Combining volatility and smoothing forecasts of UK demand for international tourism[J].Tourism Management,2009,30(4):495-511.

    [15]WEST M,HARRISON P J.Bayesian forecasting and dynamic models[M].2nd ed.New York:Springer-Verlag,1997:20-27.

    猜你喜歡
    交通流量周期性季度
    2021年第4季度航天器發(fā)射統(tǒng)計
    國際太空(2022年2期)2022-03-15 08:03:22
    2021年第3季度航天器發(fā)射統(tǒng)計
    國際太空(2021年11期)2022-01-19 03:27:06
    2021年第2季度航天器發(fā)射統(tǒng)計
    國際太空(2021年8期)2021-11-05 08:32:44
    數(shù)列中的周期性和模周期性
    基于XGBOOST算法的擁堵路段短時交通流量預(yù)測
    基于GA-BP神經(jīng)網(wǎng)絡(luò)的衡大高速公路日交通流量預(yù)測
    一類整數(shù)遞推數(shù)列的周期性
    翻番的1季度與瘋狂的3月
    基于擴頻碼周期性的單通道直擴通信半盲分離抗干擾算法
    基于復(fù)合卡和ETC的交通流量采集研究
    性色avwww在线观看| 亚洲天堂国产精品一区在线| 国产精品无大码| 观看免费一级毛片| 成人无遮挡网站| 午夜影院日韩av| 老女人水多毛片| 我要看日韩黄色一级片| 精品一区二区三区人妻视频| 亚洲久久久久久中文字幕| 嫩草影院新地址| 午夜福利18| 成人高潮视频无遮挡免费网站| 能在线免费观看的黄片| 最近视频中文字幕2019在线8| 一区二区三区高清视频在线| a级一级毛片免费在线观看| 国产麻豆成人av免费视频| 亚洲美女视频黄频| 免费无遮挡裸体视频| 国产精品美女特级片免费视频播放器| 亚洲av美国av| 日本免费一区二区三区高清不卡| 婷婷精品国产亚洲av在线| 欧美成人精品欧美一级黄| 亚洲无线在线观看| 精品人妻偷拍中文字幕| 久久久精品大字幕| 麻豆国产av国片精品| 亚洲自偷自拍三级| 高清午夜精品一区二区三区 | 97人妻精品一区二区三区麻豆| 亚洲中文字幕一区二区三区有码在线看| 亚洲图色成人| 亚洲七黄色美女视频| 一进一出好大好爽视频| 天天一区二区日本电影三级| 亚洲欧美精品综合久久99| 欧美zozozo另类| 久久国产乱子免费精品| av在线亚洲专区| 一本精品99久久精品77| 一个人免费在线观看电影| 亚洲成人中文字幕在线播放| 久久久久久久久久成人| 亚洲精品乱码久久久v下载方式| 春色校园在线视频观看| 床上黄色一级片| 3wmmmm亚洲av在线观看| 熟女电影av网| 中文资源天堂在线| 级片在线观看| 蜜桃亚洲精品一区二区三区| 精品无人区乱码1区二区| 91在线精品国自产拍蜜月| 黑人高潮一二区| 亚洲自拍偷在线| 大香蕉久久网| 成人综合一区亚洲| 一级av片app| 校园人妻丝袜中文字幕| 97超视频在线观看视频| 2021天堂中文幕一二区在线观| 变态另类成人亚洲欧美熟女| 日韩一区二区视频免费看| 亚洲图色成人| 久久精品国产99精品国产亚洲性色| 欧美成人a在线观看| 99久久中文字幕三级久久日本| 人人妻,人人澡人人爽秒播| 国产视频一区二区在线看| 看免费成人av毛片| 国产精品人妻久久久久久| 秋霞在线观看毛片| 欧美一区二区亚洲| 在线播放国产精品三级| 久久久久久久久中文| 丝袜喷水一区| 色综合色国产| 精品日产1卡2卡| 亚洲国产精品久久男人天堂| 亚洲人成网站高清观看| 人妻制服诱惑在线中文字幕| 亚洲三级黄色毛片| 中出人妻视频一区二区| 国产精华一区二区三区| 搡老妇女老女人老熟妇| 深夜精品福利| 国产不卡一卡二| 蜜臀久久99精品久久宅男| 日本一二三区视频观看| 亚洲精华国产精华液的使用体验 | 一级毛片电影观看 | 欧美精品国产亚洲| 2021天堂中文幕一二区在线观| 国产毛片a区久久久久| 天堂av国产一区二区熟女人妻| 国产精品久久电影中文字幕| 精品人妻熟女av久视频| 国产黄色小视频在线观看| 日韩欧美三级三区| 国产成人a∨麻豆精品| 男女边吃奶边做爰视频| 色噜噜av男人的天堂激情| 色5月婷婷丁香| 国产久久久一区二区三区| 日本爱情动作片www.在线观看 | 美女cb高潮喷水在线观看| 精品人妻视频免费看| 欧美日韩国产亚洲二区| 亚洲四区av| 日本爱情动作片www.在线观看 | 嫩草影院精品99| 国产男靠女视频免费网站| 成人亚洲欧美一区二区av| 可以在线观看毛片的网站| 亚洲欧美精品综合久久99| 99在线视频只有这里精品首页| 精品免费久久久久久久清纯| 亚洲国产欧洲综合997久久,| 大型黄色视频在线免费观看| 成人特级黄色片久久久久久久| 欧美高清成人免费视频www| 尤物成人国产欧美一区二区三区| 乱系列少妇在线播放| 精品免费久久久久久久清纯| 国产精品,欧美在线| 黄片wwwwww| 亚洲真实伦在线观看| 精品福利观看| 又黄又爽又免费观看的视频| 亚洲欧美成人精品一区二区| 国产成人a区在线观看| 老师上课跳d突然被开到最大视频| 免费在线观看成人毛片| 国产亚洲精品久久久com| av在线天堂中文字幕| 天堂动漫精品| 亚洲美女视频黄频| 99国产极品粉嫩在线观看| av在线蜜桃| 两个人的视频大全免费| av.在线天堂| 在线观看免费视频日本深夜| .国产精品久久| 国产精品亚洲一级av第二区| 国产视频一区二区在线看| 国产视频一区二区在线看| 大香蕉久久网| 亚洲在线自拍视频| 99riav亚洲国产免费| 高清日韩中文字幕在线| 国产熟女欧美一区二区| 成年av动漫网址| 亚洲精品久久国产高清桃花| 亚州av有码| 婷婷精品国产亚洲av| 欧美极品一区二区三区四区| 黄色视频,在线免费观看| 22中文网久久字幕| 久久久a久久爽久久v久久| 国产精品久久久久久av不卡| 男女下面进入的视频免费午夜| 国产探花极品一区二区| 色视频www国产| 九九热线精品视视频播放| 亚洲欧美精品自产自拍| 精品久久久久久久久久免费视频| 99riav亚洲国产免费| 日韩制服骚丝袜av| 免费av毛片视频| 久久精品夜夜夜夜夜久久蜜豆| www日本黄色视频网| 看非洲黑人一级黄片| 97人妻精品一区二区三区麻豆| 亚洲国产精品sss在线观看| 日本与韩国留学比较| 91av网一区二区| 亚洲成人精品中文字幕电影| 日本与韩国留学比较| 欧美日本亚洲视频在线播放| 最好的美女福利视频网| 久久人人精品亚洲av| 久久精品人妻少妇| 亚洲精品成人久久久久久| 深夜a级毛片| 亚洲真实伦在线观看| 日产精品乱码卡一卡2卡三| 色综合站精品国产| 亚洲美女视频黄频| 亚洲综合色惰| 欧美精品国产亚洲| 男女边吃奶边做爰视频| 亚洲最大成人av| 成人午夜高清在线视频| av视频在线观看入口| 欧美成人一区二区免费高清观看| 久久久久九九精品影院| 日韩欧美精品免费久久| 97热精品久久久久久| 毛片一级片免费看久久久久| 久久天躁狠狠躁夜夜2o2o| 最近在线观看免费完整版| 精品欧美国产一区二区三| 亚洲欧美中文字幕日韩二区| 男插女下体视频免费在线播放| 六月丁香七月| 国产精品久久久久久av不卡| 亚洲高清免费不卡视频| 日韩成人av中文字幕在线观看 | 少妇人妻精品综合一区二区 | 一级a爱片免费观看的视频| 人人妻人人澡欧美一区二区| 日本与韩国留学比较| 国内揄拍国产精品人妻在线| av国产免费在线观看| 久久精品国产亚洲av香蕉五月| 3wmmmm亚洲av在线观看| 色噜噜av男人的天堂激情| 色综合亚洲欧美另类图片| 亚洲美女视频黄频| 九九久久精品国产亚洲av麻豆| 亚洲av.av天堂| 九色成人免费人妻av| 国产精品爽爽va在线观看网站| 黄色视频,在线免费观看| 欧美日韩乱码在线| 日本免费a在线| 九九热线精品视视频播放| 97在线视频观看| 色综合站精品国产| 99九九线精品视频在线观看视频| 国产精品av视频在线免费观看| 男人舔奶头视频| 欧美日韩精品成人综合77777| 男女做爰动态图高潮gif福利片| 亚洲一级一片aⅴ在线观看| 国产一区二区激情短视频| 搡老妇女老女人老熟妇| 国产精品久久久久久av不卡| 亚洲国产色片| 欧美成人免费av一区二区三区| 亚洲美女搞黄在线观看 | 精品熟女少妇av免费看| 午夜福利在线观看免费完整高清在 | av天堂在线播放| 国产精品无大码| 久久婷婷人人爽人人干人人爱| 99热只有精品国产| 变态另类丝袜制服| 我要搜黄色片| 99久久无色码亚洲精品果冻| 久久天躁狠狠躁夜夜2o2o| 亚洲精品影视一区二区三区av| 俄罗斯特黄特色一大片| 岛国在线免费视频观看| 欧美激情在线99| 久久韩国三级中文字幕| 插阴视频在线观看视频| 亚洲成人久久爱视频| 麻豆成人午夜福利视频| 国产视频内射| 又爽又黄无遮挡网站| 欧美另类亚洲清纯唯美| 真实男女啪啪啪动态图| 看免费成人av毛片| 身体一侧抽搐| 少妇人妻一区二区三区视频| 亚洲欧美成人综合另类久久久 | 国产中年淑女户外野战色| 欧美色视频一区免费| 亚洲精品456在线播放app| 人妻丰满熟妇av一区二区三区| 在线国产一区二区在线| 欧美激情久久久久久爽电影| 日本免费一区二区三区高清不卡| 一区福利在线观看| 成人特级av手机在线观看| 亚洲五月天丁香| 国产亚洲欧美98| 成人综合一区亚洲| 国产精品1区2区在线观看.| 国产精品福利在线免费观看| 免费高清视频大片| а√天堂www在线а√下载| 熟女人妻精品中文字幕| 久久久a久久爽久久v久久| 白带黄色成豆腐渣| 午夜福利在线观看免费完整高清在 | 三级毛片av免费| 欧美国产日韩亚洲一区| 亚洲成av人片在线播放无| 午夜福利在线观看免费完整高清在 | 欧美性猛交黑人性爽| 人妻久久中文字幕网| av专区在线播放| 亚洲人成网站在线播| 尾随美女入室| 我的女老师完整版在线观看| 亚洲欧美日韩高清专用| 国产精品嫩草影院av在线观看| 精品国产三级普通话版| 一区二区三区高清视频在线| 国产精品久久久久久亚洲av鲁大| 色吧在线观看| 日本爱情动作片www.在线观看 | 久久久久国产精品人妻aⅴ院| 国产精品综合久久久久久久免费| 久久久成人免费电影| 免费看美女性在线毛片视频| 人妻少妇偷人精品九色| av在线老鸭窝| 亚洲欧美成人精品一区二区| 精品熟女少妇av免费看| 国产高清激情床上av| 欧美性猛交╳xxx乱大交人| 给我免费播放毛片高清在线观看| 老司机福利观看| 国产亚洲精品综合一区在线观看| 你懂的网址亚洲精品在线观看 | 老司机影院成人| 国产v大片淫在线免费观看| 欧美最黄视频在线播放免费| 在线观看免费视频日本深夜| 午夜福利18| 亚洲精品国产成人久久av| 男插女下体视频免费在线播放| 国产91av在线免费观看| 欧美激情国产日韩精品一区| av天堂在线播放| 亚洲精品粉嫩美女一区| 亚洲精品影视一区二区三区av| 久久精品国产亚洲网站| 黄色欧美视频在线观看| 悠悠久久av| 色噜噜av男人的天堂激情| 午夜亚洲福利在线播放| 午夜激情福利司机影院| 看片在线看免费视频| 亚洲中文字幕一区二区三区有码在线看| 亚洲美女视频黄频| av视频在线观看入口| 午夜精品在线福利| 精品人妻一区二区三区麻豆 | 日韩av不卡免费在线播放| 亚洲熟妇中文字幕五十中出| 亚洲精品久久国产高清桃花| 别揉我奶头 嗯啊视频| 天天一区二区日本电影三级| 亚洲精品一区av在线观看| 亚洲国产欧美人成| av女优亚洲男人天堂| 国产高清三级在线| 中文亚洲av片在线观看爽| 菩萨蛮人人尽说江南好唐韦庄 | 中文字幕精品亚洲无线码一区| 精品国内亚洲2022精品成人| 国产一级毛片七仙女欲春2| 国产私拍福利视频在线观看| 欧美成人一区二区免费高清观看| 高清毛片免费看| 一级黄色大片毛片| 日本 av在线| 亚洲精品影视一区二区三区av| 午夜精品国产一区二区电影 | 日本五十路高清| 国国产精品蜜臀av免费| 日韩 亚洲 欧美在线| 中文在线观看免费www的网站| av天堂中文字幕网| 久久精品国产清高在天天线| 国产精品久久视频播放| 久久久久精品国产欧美久久久| 国产老妇女一区| 国产精品久久久久久精品电影| 毛片一级片免费看久久久久| 亚洲久久久久久中文字幕| 欧美不卡视频在线免费观看| 国产高清激情床上av| 亚洲五月天丁香| 最近最新中文字幕大全电影3| 日本-黄色视频高清免费观看| 国产男人的电影天堂91| 国产成人aa在线观看| 天堂av国产一区二区熟女人妻| 91午夜精品亚洲一区二区三区| 欧美又色又爽又黄视频| 国产精品,欧美在线| 国内精品宾馆在线| 黄色欧美视频在线观看| 22中文网久久字幕| 舔av片在线| 久久人人精品亚洲av| 美女xxoo啪啪120秒动态图| 午夜精品国产一区二区电影 | 国产高清视频在线播放一区| 自拍偷自拍亚洲精品老妇| 嫩草影院精品99| 欧美中文日本在线观看视频| 网址你懂的国产日韩在线| 日韩欧美一区二区三区在线观看| 久久久久国内视频| 观看免费一级毛片| 身体一侧抽搐| 九九热线精品视视频播放| 蜜桃久久精品国产亚洲av| 国产亚洲精品久久久久久毛片| 国产成人91sexporn| 18禁裸乳无遮挡免费网站照片| 成人综合一区亚洲| 久久久久久九九精品二区国产| 国产一区二区三区在线臀色熟女| 干丝袜人妻中文字幕| 一个人看的www免费观看视频| 日日撸夜夜添| 男人和女人高潮做爰伦理| 国产激情偷乱视频一区二区| 国产爱豆传媒在线观看| 美女xxoo啪啪120秒动态图| 亚洲图色成人| 亚洲国产高清在线一区二区三| 亚洲国产欧美人成| 久久草成人影院| 亚洲专区国产一区二区| 国产国拍精品亚洲av在线观看| 少妇熟女欧美另类| 久久久久国产网址| 国产不卡一卡二| 最近手机中文字幕大全| 午夜日韩欧美国产| 少妇丰满av| 春色校园在线视频观看| 一级毛片aaaaaa免费看小| 国产男靠女视频免费网站| 久久6这里有精品| a级一级毛片免费在线观看| 国产av在哪里看| 少妇猛男粗大的猛烈进出视频 | 午夜a级毛片| 禁无遮挡网站| 国产亚洲精品综合一区在线观看| 少妇被粗大猛烈的视频| 久久久久国内视频| 国产探花在线观看一区二区| 国产午夜精品久久久久久一区二区三区 | 亚洲欧美成人综合另类久久久 | 国产三级在线视频| 一进一出好大好爽视频| 俄罗斯特黄特色一大片| 国产人妻一区二区三区在| 久久精品国产鲁丝片午夜精品| 国产探花极品一区二区| 久久久午夜欧美精品| 亚洲av熟女| 观看免费一级毛片| 在线观看66精品国产| 美女免费视频网站| a级毛片a级免费在线| 给我免费播放毛片高清在线观看| 别揉我奶头 嗯啊视频| 国产精品国产高清国产av| 2021天堂中文幕一二区在线观| 亚洲精品粉嫩美女一区| 午夜免费激情av| 色5月婷婷丁香| 国产精品国产三级国产av玫瑰| 亚洲成a人片在线一区二区| 欧美日韩综合久久久久久| 亚洲最大成人av| 亚洲av不卡在线观看| 亚洲成人久久性| 一区二区三区四区激情视频 | 蜜桃久久精品国产亚洲av| 三级国产精品欧美在线观看| 天美传媒精品一区二区| 精品久久久噜噜| 精品久久久久久久久亚洲| 国产淫片久久久久久久久| 免费无遮挡裸体视频| a级毛片免费高清观看在线播放| 亚洲图色成人| 国产av在哪里看| 久久九九热精品免费| 天天一区二区日本电影三级| АⅤ资源中文在线天堂| 亚洲人成网站在线播| 日韩人妻高清精品专区| 亚洲高清免费不卡视频| 免费av观看视频| 久久中文看片网| 国产单亲对白刺激| 国产女主播在线喷水免费视频网站 | 婷婷精品国产亚洲av| 国产一区二区激情短视频| 亚洲人成网站在线播| 老熟妇仑乱视频hdxx| 欧美zozozo另类| 成人美女网站在线观看视频| 亚洲成人中文字幕在线播放| 亚洲av免费高清在线观看| 国产伦在线观看视频一区| 1000部很黄的大片| 超碰av人人做人人爽久久| 十八禁国产超污无遮挡网站| 美女大奶头视频| 亚洲人与动物交配视频| av福利片在线观看| 久久久精品欧美日韩精品| 欧美高清性xxxxhd video| 精品人妻偷拍中文字幕| 中国美女看黄片| 看黄色毛片网站| 九九久久精品国产亚洲av麻豆| 在线免费十八禁| 日本成人三级电影网站| 狂野欧美激情性xxxx在线观看| 国产在视频线在精品| 国产精品无大码| 欧美xxxx性猛交bbbb| av福利片在线观看| 欧美性猛交╳xxx乱大交人| 国产欧美日韩精品一区二区| 国产高清视频在线观看网站| 精品99又大又爽又粗少妇毛片| 成人av一区二区三区在线看| 免费人成视频x8x8入口观看| 亚洲中文日韩欧美视频| 香蕉av资源在线| 亚洲av电影不卡..在线观看| 波多野结衣高清无吗| 午夜福利在线观看吧| 欧美日韩国产亚洲二区| a级毛片免费高清观看在线播放| 中文字幕人妻熟人妻熟丝袜美| 精华霜和精华液先用哪个| aaaaa片日本免费| 国产爱豆传媒在线观看| 国产成人影院久久av| 伊人久久精品亚洲午夜| 国产一区二区在线av高清观看| 亚洲天堂国产精品一区在线| 少妇被粗大猛烈的视频| 成人鲁丝片一二三区免费| 欧美+亚洲+日韩+国产| 日韩亚洲欧美综合| 欧美最新免费一区二区三区| 狠狠狠狠99中文字幕| 午夜日韩欧美国产| 天堂√8在线中文| 国内精品久久久久精免费| 香蕉av资源在线| 国产精品久久电影中文字幕| 精华霜和精华液先用哪个| 国产精品一二三区在线看| 欧美+亚洲+日韩+国产| 国产精品久久久久久av不卡| 亚洲在线观看片| 亚洲av中文av极速乱| av视频在线观看入口| 久久久久久久亚洲中文字幕| 此物有八面人人有两片| 99国产精品一区二区蜜桃av| 男女那种视频在线观看| 男女啪啪激烈高潮av片| 日本爱情动作片www.在线观看 | 熟妇人妻久久中文字幕3abv| 久久99热这里只有精品18| 可以在线观看的亚洲视频| 日韩国内少妇激情av| 免费不卡的大黄色大毛片视频在线观看 | 中出人妻视频一区二区| 亚洲欧美精品自产自拍| 深爱激情五月婷婷| 日本一本二区三区精品| av中文乱码字幕在线| 男女视频在线观看网站免费| 成人综合一区亚洲| 精品午夜福利在线看| 欧美一级a爱片免费观看看| 丰满人妻一区二区三区视频av| 午夜久久久久精精品| 久久99热这里只有精品18| 国产精品永久免费网站| 久久精品国产清高在天天线| 久久人妻av系列| 99热只有精品国产| 天堂网av新在线| АⅤ资源中文在线天堂| 成人欧美大片| 欧美一区二区国产精品久久精品| 国产成人a区在线观看| 亚洲av免费在线观看| 国产大屁股一区二区在线视频| 99久久精品热视频| 午夜福利成人在线免费观看| 婷婷精品国产亚洲av在线| 久久6这里有精品| 精品久久久久久久久久免费视频| 美女免费视频网站| 热99re8久久精品国产| 国产熟女欧美一区二区| av在线天堂中文字幕| 成人美女网站在线观看视频| 国产不卡一卡二| 少妇熟女欧美另类| 午夜影院日韩av| 国产69精品久久久久777片| 国产精品无大码| 国产单亲对白刺激| 国产一区二区三区av在线 | 亚洲欧美清纯卡通| 国产老妇女一区| 一本久久中文字幕| 两个人的视频大全免费| 乱系列少妇在线播放|