• <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的交通流量采集研究
    人人妻人人爽人人添夜夜欢视频| 一级a爱视频在线免费观看| 成年av动漫网址| a级片在线免费高清观看视频| 国产福利在线免费观看视频| 国产精品人妻久久久影院| 亚洲婷婷狠狠爱综合网| 国产日韩欧美亚洲二区| 成人亚洲精品一区在线观看| 免费在线观看黄色视频的| 黑人猛操日本美女一级片| 又黄又粗又硬又大视频| av在线老鸭窝| 国产有黄有色有爽视频| 日韩视频在线欧美| 国产免费福利视频在线观看| 午夜久久久在线观看| 久久午夜福利片| 欧美 日韩 精品 国产| 尾随美女入室| 青草久久国产| 综合色丁香网| 中文字幕亚洲精品专区| 久久久久国产网址| 亚洲欧美成人综合另类久久久| xxxhd国产人妻xxx| 欧美 日韩 精品 国产| 亚洲,欧美,日韩| 大香蕉久久成人网| 精品少妇久久久久久888优播| 黄网站色视频无遮挡免费观看| 青春草视频在线免费观看| 久久人妻熟女aⅴ| 最近最新中文字幕免费大全7| 亚洲国产欧美网| 亚洲精品国产av蜜桃| 在线精品无人区一区二区三| 中文欧美无线码| av天堂久久9| 亚洲五月色婷婷综合| 日韩欧美一区视频在线观看| 丝袜人妻中文字幕| 丰满迷人的少妇在线观看| 狂野欧美激情性bbbbbb| 伦精品一区二区三区| 国产黄频视频在线观看| 亚洲情色 制服丝袜| 麻豆精品久久久久久蜜桃| 久久毛片免费看一区二区三区| 欧美精品一区二区大全| 国产视频首页在线观看| 侵犯人妻中文字幕一二三四区| 精品视频人人做人人爽| 婷婷色麻豆天堂久久| 男女国产视频网站| 国产成人欧美| 在线天堂中文资源库| av免费观看日本| 午夜免费鲁丝| 久久精品国产a三级三级三级| 成人亚洲精品一区在线观看| 色吧在线观看| 久久国内精品自在自线图片| 久久这里有精品视频免费| 伊人亚洲综合成人网| 国产男女内射视频| 精品国产超薄肉色丝袜足j| 亚洲,欧美精品.| 街头女战士在线观看网站| 成人18禁高潮啪啪吃奶动态图| 亚洲精品第二区| 少妇的逼水好多| 亚洲精品视频女| 国产人伦9x9x在线观看 | 熟女电影av网| 亚洲精品国产一区二区精华液| 9热在线视频观看99| 久久这里有精品视频免费| 国产精品 国内视频| 精品一区在线观看国产| 晚上一个人看的免费电影| 黄色配什么色好看| 亚洲成国产人片在线观看| 2021少妇久久久久久久久久久| 久久久国产欧美日韩av| 精品人妻偷拍中文字幕| 两个人看的免费小视频| 久久久久久人人人人人| 国产精品.久久久| 国产精品秋霞免费鲁丝片| 嫩草影院入口| 欧美97在线视频| 天天躁夜夜躁狠狠久久av| 成人影院久久| 另类精品久久| videosex国产| 2018国产大陆天天弄谢| 国产精品久久久久久精品古装| 久久精品国产亚洲av涩爱| 波多野结衣一区麻豆| 免费观看无遮挡的男女| 亚洲第一av免费看| 国产精品 国内视频| 搡老乐熟女国产| 韩国av在线不卡| 欧美黄色片欧美黄色片| 亚洲伊人久久精品综合| 欧美日韩精品网址| 久久av网站| av免费观看日本| 色播在线永久视频| 精品一品国产午夜福利视频| 在线观看免费日韩欧美大片| 国产免费一区二区三区四区乱码| 精品视频人人做人人爽| 欧美人与性动交α欧美精品济南到 | 日本免费在线观看一区| 日产精品乱码卡一卡2卡三| 免费人妻精品一区二区三区视频| 天天躁夜夜躁狠狠久久av| 少妇人妻 视频| 国产精品一区二区在线不卡| 亚洲综合精品二区| 老汉色∧v一级毛片| 日韩电影二区| 亚洲人成77777在线视频| 中文欧美无线码| 国产精品免费大片| 日韩欧美精品免费久久| 另类精品久久| 久久久欧美国产精品| 高清黄色对白视频在线免费看| 国产成人91sexporn| 国产精品嫩草影院av在线观看| 赤兔流量卡办理| 街头女战士在线观看网站| 卡戴珊不雅视频在线播放| 亚洲国产精品国产精品| 日韩av在线免费看完整版不卡| 久久热在线av| videos熟女内射| 电影成人av| 久久热在线av| 精品亚洲乱码少妇综合久久| 2022亚洲国产成人精品| 青春草亚洲视频在线观看| 麻豆精品久久久久久蜜桃| 黄色一级大片看看| 免费久久久久久久精品成人欧美视频| 精品国产超薄肉色丝袜足j| 人人妻人人澡人人爽人人夜夜| 中文精品一卡2卡3卡4更新| 97精品久久久久久久久久精品| 蜜桃在线观看..| 五月伊人婷婷丁香| av女优亚洲男人天堂| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 毛片一级片免费看久久久久| 亚洲三级黄色毛片| 久久精品亚洲av国产电影网| 黄色怎么调成土黄色| 亚洲少妇的诱惑av| 国产精品偷伦视频观看了| 精品少妇内射三级| 精品人妻在线不人妻| 制服诱惑二区| av女优亚洲男人天堂| 人人妻人人爽人人添夜夜欢视频| 亚洲av男天堂| 另类精品久久| 交换朋友夫妻互换小说| 国产综合精华液| 女人高潮潮喷娇喘18禁视频| 国产精品偷伦视频观看了| 大片免费播放器 马上看| 精品视频人人做人人爽| 精品久久久久久电影网| 在线观看www视频免费| 99精国产麻豆久久婷婷| 久久综合国产亚洲精品| 青春草视频在线免费观看| 亚洲天堂av无毛| 美国免费a级毛片| 少妇被粗大的猛进出69影院| 色吧在线观看| 亚洲国产成人一精品久久久| 电影成人av| 久久 成人 亚洲| 日韩电影二区| 亚洲av.av天堂| 侵犯人妻中文字幕一二三四区| 亚洲中文av在线| 又粗又硬又长又爽又黄的视频| 麻豆精品久久久久久蜜桃| 亚洲成人av在线免费| 国产精品欧美亚洲77777| 午夜精品国产一区二区电影| freevideosex欧美| 欧美国产精品va在线观看不卡| 亚洲久久久国产精品| 亚洲成国产人片在线观看| 一区二区三区激情视频| 国产成人精品无人区| 久久久久久久精品精品| 国产av一区二区精品久久| 国产极品粉嫩免费观看在线| 国产成人aa在线观看| 女性生殖器流出的白浆| 久久久久久久久免费视频了| 亚洲精品久久成人aⅴ小说| 少妇被粗大猛烈的视频| 一边摸一边做爽爽视频免费| 国产av精品麻豆| 欧美精品av麻豆av| 国产免费又黄又爽又色| 亚洲天堂av无毛| 国产片特级美女逼逼视频| 岛国毛片在线播放| 国产精品女同一区二区软件| 秋霞伦理黄片| 国产97色在线日韩免费| 久久久久久久大尺度免费视频| 肉色欧美久久久久久久蜜桃| 啦啦啦视频在线资源免费观看| 亚洲国产看品久久| 性少妇av在线| 欧美亚洲日本最大视频资源| 日韩av不卡免费在线播放| 午夜日韩欧美国产| 亚洲视频免费观看视频| 十分钟在线观看高清视频www| 另类精品久久| 最新的欧美精品一区二区| 性色avwww在线观看| 亚洲成av片中文字幕在线观看 | 亚洲美女搞黄在线观看| 国产精品无大码| 99久国产av精品国产电影| 尾随美女入室| 欧美少妇被猛烈插入视频| av线在线观看网站| 观看av在线不卡| 国产 精品1| 欧美最新免费一区二区三区| 在线观看国产h片| 欧美精品高潮呻吟av久久| 哪个播放器可以免费观看大片| 色婷婷久久久亚洲欧美| 亚洲成色77777| 亚洲欧洲精品一区二区精品久久久 | 最近2019中文字幕mv第一页| 久久女婷五月综合色啪小说| 日日摸夜夜添夜夜爱| 国产精品国产av在线观看| 国产97色在线日韩免费| 女的被弄到高潮叫床怎么办| 免费人妻精品一区二区三区视频| 最近的中文字幕免费完整| 国产有黄有色有爽视频| 王馨瑶露胸无遮挡在线观看| 看免费成人av毛片| 综合色丁香网| 欧美黄色片欧美黄色片| 国产高清国产精品国产三级| 午夜福利一区二区在线看| 色婷婷av一区二区三区视频| 精品99又大又爽又粗少妇毛片| 国产极品天堂在线| 国产 精品1| 一级毛片黄色毛片免费观看视频| 好男人视频免费观看在线| 欧美+日韩+精品| 国产精品欧美亚洲77777| 免费黄频网站在线观看国产| 99热网站在线观看| 欧美亚洲日本最大视频资源| 黄片无遮挡物在线观看| 色播在线永久视频| 老熟女久久久| 天天躁日日躁夜夜躁夜夜| 超碰97精品在线观看| 国产又爽黄色视频| 一级,二级,三级黄色视频| 国产亚洲欧美精品永久| 王馨瑶露胸无遮挡在线观看| 欧美 日韩 精品 国产| 午夜福利乱码中文字幕| 黄色一级大片看看| 最新中文字幕久久久久| 午夜老司机福利剧场| 成年人免费黄色播放视频| 中文字幕最新亚洲高清| 少妇被粗大的猛进出69影院| 国产精品欧美亚洲77777| 久久久a久久爽久久v久久| 最黄视频免费看| 青春草国产在线视频| 老鸭窝网址在线观看| 一本大道久久a久久精品| 久久久久久久久久久免费av| 91精品伊人久久大香线蕉| 精品少妇黑人巨大在线播放| 男女高潮啪啪啪动态图| 少妇被粗大的猛进出69影院| 激情五月婷婷亚洲| 成年动漫av网址| av片东京热男人的天堂| 老汉色∧v一级毛片| 超色免费av| 最近2019中文字幕mv第一页| 91午夜精品亚洲一区二区三区| 色婷婷久久久亚洲欧美| 国产一区二区激情短视频 | 成年动漫av网址| 成年女人毛片免费观看观看9 | 亚洲 欧美一区二区三区| 亚洲视频免费观看视频| 制服人妻中文乱码| 蜜桃在线观看..| 国产精品欧美亚洲77777| 美女午夜性视频免费| 一级毛片 在线播放| 成年av动漫网址| 五月伊人婷婷丁香| 中文字幕制服av| 成人手机av| 最近2019中文字幕mv第一页| 亚洲精品自拍成人| 丝袜脚勾引网站| 大话2 男鬼变身卡| www日本在线高清视频| 你懂的网址亚洲精品在线观看| 99久久综合免费| 岛国毛片在线播放| 久久韩国三级中文字幕| 久久女婷五月综合色啪小说| 精品国产国语对白av| 五月开心婷婷网| 久久97久久精品| 精品人妻熟女毛片av久久网站| 成人影院久久| 高清视频免费观看一区二区| 色视频在线一区二区三区| 亚洲男人天堂网一区| 亚洲精品国产色婷婷电影| 日韩欧美一区视频在线观看| 伦理电影免费视频| 亚洲国产精品999| 国产高清国产精品国产三级| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 黑人巨大精品欧美一区二区蜜桃| a级毛片在线看网站| 在线观看美女被高潮喷水网站| 飞空精品影院首页| 成人午夜精彩视频在线观看| 1024视频免费在线观看| 国产人伦9x9x在线观看 | 一级毛片黄色毛片免费观看视频| 99久久中文字幕三级久久日本| 久久国内精品自在自线图片| 黄片小视频在线播放| 男女啪啪激烈高潮av片| 亚洲欧洲精品一区二区精品久久久 | 欧美激情极品国产一区二区三区| 超色免费av| 亚洲色图综合在线观看| 天天影视国产精品| 91aial.com中文字幕在线观看| 老司机亚洲免费影院| 亚洲国产日韩一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品国产三级国产专区5o| 捣出白浆h1v1| 亚洲综合色网址| 狂野欧美激情性bbbbbb| 欧美激情高清一区二区三区 | 欧美日韩视频高清一区二区三区二| 丝袜在线中文字幕| 日本-黄色视频高清免费观看| 国产精品女同一区二区软件| 熟妇人妻不卡中文字幕| 久久久久精品人妻al黑| 成人国产麻豆网| 亚洲精品乱久久久久久| 欧美国产精品va在线观看不卡| 日韩,欧美,国产一区二区三区| 各种免费的搞黄视频| 亚洲三区欧美一区| 黄频高清免费视频| 久久97久久精品| 亚洲人成网站在线观看播放| av片东京热男人的天堂| 精品国产超薄肉色丝袜足j| 寂寞人妻少妇视频99o| av在线播放精品| 亚洲视频免费观看视频| 精品久久蜜臀av无| 久久久久久人妻| 亚洲av国产av综合av卡| 最近手机中文字幕大全| 亚洲成av片中文字幕在线观看 | 最近2019中文字幕mv第一页| 人体艺术视频欧美日本| 亚洲 欧美一区二区三区| 国产精品国产av在线观看| 在线观看三级黄色| 99久久人妻综合| 免费观看a级毛片全部| 999久久久国产精品视频| 亚洲国产av新网站| 少妇精品久久久久久久| 黑人猛操日本美女一级片| 尾随美女入室| 久久国产亚洲av麻豆专区| 日韩熟女老妇一区二区性免费视频| 好男人视频免费观看在线| 欧美日韩一区二区视频在线观看视频在线| 久久久久精品性色| 亚洲精品,欧美精品| 香蕉精品网在线| 又大又黄又爽视频免费| 精品亚洲成国产av| 男女边吃奶边做爰视频| 亚洲久久久国产精品| 啦啦啦视频在线资源免费观看| 久久国产精品大桥未久av| 婷婷色综合www| 欧美精品一区二区大全| 久久精品久久久久久久性| av福利片在线| 波多野结衣av一区二区av| 中文精品一卡2卡3卡4更新| 少妇的逼水好多| 我的亚洲天堂| 亚洲成色77777| 欧美 亚洲 国产 日韩一| 老鸭窝网址在线观看| 日韩中文字幕欧美一区二区 | 2022亚洲国产成人精品| 三上悠亚av全集在线观看| 一级毛片 在线播放| 久久精品熟女亚洲av麻豆精品| 国产精品一区二区在线不卡| 亚洲中文av在线| 最近手机中文字幕大全| 麻豆av在线久日| 在线观看www视频免费| 国产极品粉嫩免费观看在线| 三上悠亚av全集在线观看| 成人黄色视频免费在线看| 麻豆乱淫一区二区| 老司机影院毛片| 精品少妇久久久久久888优播| 亚洲视频免费观看视频| 一本—道久久a久久精品蜜桃钙片| 一级a爱视频在线免费观看| 国产精品一区二区在线观看99| 欧美国产精品一级二级三级| 国产爽快片一区二区三区| 国产成人精品无人区| 日韩一本色道免费dvd| 日本爱情动作片www.在线观看| 亚洲视频免费观看视频| 久久精品久久精品一区二区三区| 90打野战视频偷拍视频| 在现免费观看毛片| 成人手机av| 乱人伦中国视频| 2018国产大陆天天弄谢| 久久精品国产亚洲av涩爱| 午夜91福利影院| 久久99一区二区三区| 一边摸一边做爽爽视频免费| 又大又黄又爽视频免费| 精品国产一区二区三区四区第35| 久久人人爽av亚洲精品天堂| 亚洲成人av在线免费| 亚洲欧美中文字幕日韩二区| 国产 一区精品| 国产亚洲精品第一综合不卡| 91久久精品国产一区二区三区| 亚洲精品美女久久久久99蜜臀 | 欧美日韩视频高清一区二区三区二| 亚洲av电影在线进入| 国产精品久久久久久精品电影小说| 久久99蜜桃精品久久| 伊人久久大香线蕉亚洲五| 欧美精品国产亚洲| 亚洲在久久综合| 色播在线永久视频| 日韩人妻精品一区2区三区| 国产成人欧美| 精品午夜福利在线看| 国产伦理片在线播放av一区| 天天躁日日躁夜夜躁夜夜| 一级,二级,三级黄色视频| 中文字幕制服av| 亚洲国产欧美在线一区| 高清视频免费观看一区二区| 最近手机中文字幕大全| 只有这里有精品99| 亚洲国产精品国产精品| 亚洲精品国产色婷婷电影| 女性生殖器流出的白浆| 下体分泌物呈黄色| 国产欧美日韩综合在线一区二区| 国产免费又黄又爽又色| 亚洲精品久久午夜乱码| 永久网站在线| 老汉色∧v一级毛片| 日韩三级伦理在线观看| 久久av网站| 一本久久精品| 国产精品国产三级专区第一集| 老熟女久久久| 国产视频首页在线观看| 国产高清不卡午夜福利| 一二三四在线观看免费中文在| 视频区图区小说| 成年动漫av网址| 高清不卡的av网站| 久久久久久久久久久久大奶| 久久久久久人妻| 国产亚洲午夜精品一区二区久久| 国产一区二区激情短视频 | 欧美日韩一级在线毛片| 日韩精品免费视频一区二区三区| 免费观看无遮挡的男女| 美女大奶头黄色视频| 欧美成人精品欧美一级黄| 精品久久蜜臀av无| 亚洲精品美女久久av网站| 国产精品嫩草影院av在线观看| 波多野结衣一区麻豆| 欧美在线黄色| av在线播放精品| 天天操日日干夜夜撸| 国产一级毛片在线| 熟女少妇亚洲综合色aaa.| 大香蕉久久网| a级毛片黄视频| 免费在线观看黄色视频的| 午夜福利乱码中文字幕| 嫩草影院入口| 汤姆久久久久久久影院中文字幕| 91成人精品电影| 国产综合精华液| 成年美女黄网站色视频大全免费| 国产精品久久久久久精品古装| 久久久久国产一级毛片高清牌| 不卡视频在线观看欧美| 老司机影院成人| 26uuu在线亚洲综合色| 午夜福利乱码中文字幕| 精品卡一卡二卡四卡免费| 91午夜精品亚洲一区二区三区| 亚洲精品一区蜜桃| 美女xxoo啪啪120秒动态图| 日本免费在线观看一区| 欧美亚洲日本最大视频资源| 熟女少妇亚洲综合色aaa.| 伦精品一区二区三区| 男男h啪啪无遮挡| 国产成人精品久久久久久| 免费在线观看视频国产中文字幕亚洲 | 国产在线一区二区三区精| 国产精品秋霞免费鲁丝片| 日韩制服丝袜自拍偷拍| 国产极品粉嫩免费观看在线| 999精品在线视频| 汤姆久久久久久久影院中文字幕| 麻豆精品久久久久久蜜桃| 日本av免费视频播放| 国产欧美日韩一区二区三区在线| 只有这里有精品99| 在线观看www视频免费| 久久99精品国语久久久| 少妇熟女欧美另类| 一级片免费观看大全| 少妇人妻精品综合一区二区| 日韩一本色道免费dvd| 亚洲欧美一区二区三区国产| 人成视频在线观看免费观看| 97人妻天天添夜夜摸| 91国产中文字幕| 亚洲美女搞黄在线观看| 精品午夜福利在线看| 亚洲人成网站在线观看播放| 黄色怎么调成土黄色| 亚洲成人av在线免费| 97人妻天天添夜夜摸| 制服诱惑二区| av在线app专区| av一本久久久久| 国产一区二区三区综合在线观看| 精品久久久精品久久久| 亚洲av欧美aⅴ国产| 97精品久久久久久久久久精品| 麻豆av在线久日| 一二三四在线观看免费中文在| 日本午夜av视频| 欧美av亚洲av综合av国产av | 人妻人人澡人人爽人人| 日韩一本色道免费dvd| 香蕉精品网在线| av免费在线看不卡| 亚洲精品久久成人aⅴ小说| 欧美日韩视频高清一区二区三区二| 久久女婷五月综合色啪小说| av女优亚洲男人天堂| 国产精品一国产av| 丁香六月天网| 日本猛色少妇xxxxx猛交久久| 欧美+日韩+精品| 人妻人人澡人人爽人人|