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

    基于ADPSO算法優(yōu)化LSSVM的高速公路交通量預(yù)測(cè)方法

    2017-06-05 09:34:14司文靜封喜波耿立艷張占福北華航天工業(yè)學(xué)院建筑工程系河北廊坊06000河北省高速公路廊坊北三縣管理處河北廊坊06000石家莊鐵道大學(xué)經(jīng)濟(jì)管理學(xué)院河北石家莊00043曼徹斯特城市大學(xué)商學(xué)院曼徹斯特英國(guó)6BH石家莊鐵道大學(xué)四方學(xué)院河北石家莊03
    關(guān)鍵詞:交通量適應(yīng)度粒子

    司文靜,封喜波,耿立艷,張占福(.北華航天工業(yè)學(xué)院 建筑工程系,河北 廊坊 06000;.河北省高速公路廊坊北三縣管理處,河北 廊坊 06000;3.石家莊鐵道大學(xué) 經(jīng)濟(jì)管理學(xué)院,河北 石家莊 00043;4.曼徹斯特城市大學(xué) 商學(xué)院,曼徹斯特 英國(guó) M 6BH;.石家莊鐵道大學(xué) 四方學(xué)院,河北 石家莊 03)

    ?

    基于ADPSO算法優(yōu)化LSSVM的高速公路交通量預(yù)測(cè)方法

    司文靜1,封喜波2,耿立艷3,4,張占福5
    (1.北華航天工業(yè)學(xué)院 建筑工程系,河北 廊坊 065000;2.河北省高速公路廊坊北三縣管理處,河北 廊坊 065000;3.石家莊鐵道大學(xué) 經(jīng)濟(jì)管理學(xué)院,河北 石家莊 050043;4.曼徹斯特城市大學(xué) 商學(xué)院,曼徹斯特 英國(guó) M15 6BH;5.石家莊鐵道大學(xué) 四方學(xué)院,河北 石家莊 051132)

    針對(duì)高速公路交通量與其經(jīng)濟(jì)影響因素之間的復(fù)雜非線性關(guān)系,將最小二乘支持向量機(jī)(least squares support vector machines,LSSVM)與自適應(yīng)動(dòng)態(tài)粒子群優(yōu)化(adaptive dynamic particle swarm optimization,ADPSO)算法相結(jié)合,提出一種ADPSO算法優(yōu)化LSSVM的高速公路交通量新型預(yù)測(cè)方法.將建模簡(jiǎn)單、精度高的LSSVM作為預(yù)測(cè)模型,通過尋優(yōu)能力優(yōu)異的ADPSO算法選擇LSSVM最優(yōu)參數(shù).以某市高速公路交通量為例驗(yàn)證模型的有效性.結(jié)果表明,所提方法的預(yù)測(cè)性能較好,適合于高速公路交通量的短期預(yù)測(cè).

    高速公路;交通量預(yù)測(cè);自適應(yīng)動(dòng)態(tài)粒子群優(yōu)化算法;最小二乘支持向量機(jī)

    高速公路交通量的準(zhǔn)確預(yù)測(cè)是高速公路規(guī)劃與設(shè)計(jì)的重要基礎(chǔ)工作.由于高速公路交通量往往與某一區(qū)域的社會(huì)經(jīng)濟(jì)發(fā)展水平聯(lián)系緊密,經(jīng)濟(jì)因素是影響某地區(qū)高速交通量的主要因素,而經(jīng)濟(jì)影響因素對(duì)高速公路交通量的作用形式不同,這就使得高速公路交通量與其經(jīng)濟(jì)影響因素之間存在復(fù)雜的非線性關(guān)系.高速公路交通量預(yù)測(cè)方法中,較早使用的方法主要有四階段法[1]、時(shí)間序列預(yù)測(cè)法[2]、灰色預(yù)測(cè)法[3]等.四階段法需要大量的調(diào)查數(shù)據(jù),不適合大規(guī)模應(yīng)用.時(shí)間序列預(yù)測(cè)法和灰色預(yù)測(cè)法均為單變量預(yù)測(cè)模型,計(jì)算簡(jiǎn)單且易于實(shí)現(xiàn),但難以揭示高速公路交通量與經(jīng)濟(jì)影響因素之間的復(fù)雜非線性關(guān)系.

    近年來,神經(jīng)網(wǎng)絡(luò)被引入到高速公路交通量預(yù)測(cè)中.神經(jīng)網(wǎng)絡(luò)是一種數(shù)據(jù)驅(qū)動(dòng)的非參數(shù)預(yù)測(cè)方法,可靈活反映高速公路交通量與經(jīng)濟(jì)影響因素之間的復(fù)雜非線性關(guān)系,取得了較好的預(yù)測(cè)效果[4].神經(jīng)網(wǎng)絡(luò)以經(jīng)驗(yàn)風(fēng)險(xiǎn)最小化為基礎(chǔ),在實(shí)際應(yīng)用中容易遇到局部極優(yōu)值、過擬合等問題,而且需要大量的數(shù)據(jù)樣本進(jìn)行網(wǎng)絡(luò)訓(xùn)練.作為一種基于統(tǒng)計(jì)學(xué)習(xí)理論和結(jié)構(gòu)風(fēng)險(xiǎn)最小化原則的學(xué)習(xí)算法[5],支持向量機(jī)(support vector machines,SVM)較好地解決了神經(jīng)網(wǎng)絡(luò)的大樣本訓(xùn)練問題,具有良好的非線性擬合能力和泛化能力,在數(shù)據(jù)樣本較少的情況下依然能獲得較高的交通量預(yù)測(cè)精度[6-7].由于SVM需求解二次規(guī)劃方程,使用的樣本數(shù)量越多,SVM的計(jì)算復(fù)雜程度越高.最小二乘支持向量機(jī)[8](least squares support vector machines,LSSVM)是一種以SVM為基礎(chǔ)的學(xué)習(xí)算法,它將最小二乘線性系統(tǒng)作為損失函數(shù),簡(jiǎn)化了計(jì)算過程,提高了求解速度[9-10],更適合于高速公路交通量預(yù)測(cè)研究.LSSVM的性能依賴于自身參數(shù)的選擇,傳統(tǒng)參數(shù)選擇方法大多具有主觀性,難以達(dá)到滿意的預(yù)測(cè)效果.在眾多智能優(yōu)化算法中,粒子群優(yōu)化(particle swarm optimization,PSO)算法由于具有實(shí)現(xiàn)簡(jiǎn)單、尋優(yōu)效果好等特點(diǎn)[11],在LSSVM參數(shù)選擇方面的有效性已得到證明.傳統(tǒng)PSO算法中,慣性權(quán)重、加速系數(shù)的取值對(duì)算法的優(yōu)化性能起到非常重要的作用.國(guó)內(nèi)外學(xué)者對(duì)這2個(gè)參數(shù)進(jìn)行了不同的修正,形成了不同的改進(jìn)PSO算法,其中,有代表性的有基于慣性權(quán)重自適應(yīng)調(diào)整的自適應(yīng)慣性權(quán)重粒子群優(yōu)化(adaptive inertia weight particle swarm optimization,AIWPSO)算法[12]和基于加速系數(shù)動(dòng)態(tài)調(diào)整的動(dòng)態(tài)加速系數(shù)粒子群優(yōu)化(dynamic acceleration coefficients particle swarm optimization,DACPSO)算法[13],這2種改進(jìn)PSO算法均在一定程度上改善了傳統(tǒng)PSO算法的尋優(yōu)效率.為進(jìn)一步提高PSO算法的尋優(yōu)能力,本文基于這2種改進(jìn)PSO算法,提出自適應(yīng)動(dòng)態(tài)粒子群優(yōu)化(adaptive dynamic particle swarm optimization,ADPSO)算法,并利用ADPSO算法選擇LSSVM最優(yōu)參數(shù).通過對(duì)高速公路交通量的實(shí)例計(jì)算,驗(yàn)證該方法的有效性.

    1 最小二乘支持向量機(jī)

    LSSVM是SVM的一種改進(jìn)形式,將標(biāo)準(zhǔn)SVM中的不等式約束條件轉(zhuǎn)換為等式約束條件.設(shè)有N組訓(xùn)練樣本{(xk,yk)|k=1,2,…,N},其中,xk∈Rd為d維輸入向量,yk∈R為對(duì)應(yīng)的一維輸出.LSSVM函數(shù)估計(jì)的優(yōu)化問題如下:

    (1)

    其中,φ(x)為一非線性映射函數(shù),用于將原始輸入數(shù)據(jù)映射到高維特征空間;ω、b分別為權(quán)重向量與偏差常量;γ為正則化參數(shù);ek∈R為誤差變量.利用拉格朗日法求解上述優(yōu)化問題,首先構(gòu)造優(yōu)化問題對(duì)應(yīng)的拉格朗日函數(shù),如下所示:

    (2)

    其中,α=[α1,…,αN]T為拉格朗日乘子向量;然后分別求解拉格朗日函數(shù)(α,ω,b,e)對(duì)αk、ω、b、ek的偏導(dǎo)數(shù),從而將優(yōu)化問題(1)轉(zhuǎn)化為求解一組線性方程,如下:

    (3)

    其中,1N為N階列向量,其中的元素都為1;為IN階單位矩陣;Ω為N×N階矩陣,其元素為Ωk=K(x,xk)=?(x)T?(xk),K(x,xk)為滿足Mercer條件的核函數(shù);Y=[y1,…,yN]T.由式(3)求出a和b后,可得到LSSVM回歸模型為

    b.

    (4)

    選取泛化能力強(qiáng)的RBF核函數(shù)作為L(zhǎng)SSVM核函數(shù),則式(4)轉(zhuǎn)化為

    ‖x-xk‖2/σ2)]+b,

    (5)

    其中,σ為核函數(shù)參數(shù).

    2 自適應(yīng)動(dòng)態(tài)粒子群優(yōu)化算法

    慣性權(quán)重δ和加速系數(shù)c1、c2是影響PSO算法性能的2個(gè)重要控制參數(shù),它們的取值直接影響到算法的搜索能力和收斂速度,其中,δ描述了粒子上一步迭代的速度對(duì)當(dāng)前速度影響的程度.δ取值越大,粒子的速度越大,有利于粒子發(fā)現(xiàn)新的解域,算法的全局搜索能力越強(qiáng);δ取值越小,粒子的速度越小,有利于粒子在當(dāng)前空間搜索更優(yōu)解,算法的局部搜索能力越強(qiáng).目前采用較多的慣性權(quán)重是線性遞減策略,其值隨著迭代次數(shù)的增加而線性減小.線性遞減策略雖可一定程度上平衡全局搜索能力和局部搜索能力,但粒子的實(shí)際搜索過程是一個(gè)非線性過程,該策略難以正確反映粒子的真實(shí)搜索過程,而且線性遞減策略無法考慮適應(yīng)度函數(shù)提供的信息,將導(dǎo)致粒子搜索方向的啟發(fā)性不強(qiáng),使得算法的搜索速度較慢而易陷入局部最優(yōu)值[14].c1、c2反映了粒子自身認(rèn)知信息與社會(huì)認(rèn)知信息的交流程度.迭代前期階段,希望c1取值較大、c2取值較小,便于粒子進(jìn)行全局尋優(yōu)、避免陷入局部極值;迭代后期階段,應(yīng)有c1取值較小、c2取值較大,使粒子迅速、準(zhǔn)確地收斂于全局最優(yōu)解.傳統(tǒng)PSO算法中,c1、c2取相同的固定值,粒子的自身認(rèn)知能力與社會(huì)認(rèn)知能力相同,算法在迭代初期雖可保持較快的收斂速度,但到迭代后期,由于粒子逐步統(tǒng)一化,容易陷入局部最優(yōu).

    慣性權(quán)重、加速系數(shù)兩參數(shù)相輔相成,共同影響著PSO算法的尋優(yōu)能力.若對(duì)它們單獨(dú)進(jìn)行調(diào)整,將削弱粒子搜索過程的統(tǒng)一性,難以適應(yīng)復(fù)雜非線性問題的優(yōu)化.ADPSO算法采用自適應(yīng)動(dòng)態(tài)調(diào)整策略,使慣性權(quán)重與加速系數(shù)同時(shí)隨適應(yīng)度值變化,以增強(qiáng)粒子的全局搜索能力與局部搜索能力、獲得更優(yōu)的搜索結(jié)果.

    假設(shè)在D維目標(biāo)搜索空間中,有一個(gè)由m個(gè)粒子組成的群體,每個(gè)粒子代表所優(yōu)化問題的一個(gè)潛在解,第i個(gè)粒子的位置記為Si=(si1,…,siD),其速度記為Vi=(vi1,…,viD).在整個(gè)種群中,每個(gè)粒子通過目標(biāo)函數(shù)的適應(yīng)度值搜索個(gè)體最優(yōu)位置Pibest=(si1best,…,siDbest)和整個(gè)粒子種群的最優(yōu)位置,也就是全局最優(yōu)位置Gbest=(s1best,…,sDbest).每個(gè)粒子的速度與位置按以下公式更新:

    (6)

    (7)

    其中,δmax、δmin分別為慣性權(quán)重的最大值與最小值,F(xiàn)為粒子的當(dāng)前適應(yīng)度值,F(xiàn)avg、Fmin分別為粒子群的平均適應(yīng)度值與最小適應(yīng)度值.由式(7)可知,當(dāng)F高于Favg時(shí),應(yīng)取較大的δ值,使粒子更快地飛向更優(yōu)的目標(biāo)搜索空間;當(dāng)F低于Favg時(shí),應(yīng)取較小的δ值,使粒子在當(dāng)前目標(biāo)搜索空間中尋找更優(yōu)位置.

    c1和c2為加速系數(shù),這里設(shè)為隨著適應(yīng)度值進(jìn)行動(dòng)態(tài)更新,更新公式為

    (8)

    c2(t)=4-c1(t),

    (9)

    其中,F(xiàn)avg(t)為粒子第t步迭代的平均適應(yīng)度值;Fg(t)為種群第t步迭代的最優(yōu)適應(yīng)度值;c1(t)、c2(t)的取值在[0,4]之間.根據(jù)式(8)和式(9),迭代前期階段,F(xiàn)avg(t)與Fg(t)的差異較大,可獲得較大的c1(t)值、較小的c2(t)值,便于粒子進(jìn)行全局搜索、免于陷入局部最優(yōu)值;迭代后期階段,F(xiàn)avg(t)與Fg(t)的差異較小,可獲得較小的c1(t)值、較大的c2(t)值,有利于增加粒子的多樣性,提高算法的收斂速度與搜索精度.

    3 LSSVM-ADPSO模型

    由式(3)和式(5)可知,基于RBF核函數(shù)的LSSVM有2個(gè)參數(shù)需要確定:正則化參數(shù)γ和核參數(shù)σ,γ用于協(xié)調(diào)LSSVM的誤差與復(fù)雜程度,σ反映了數(shù)據(jù)樣本的范圍特性.參數(shù)組合(γ,σ)的合理選取對(duì)提高LSSVM的學(xué)習(xí)能力與泛化能力至關(guān)重要.傳統(tǒng)參數(shù)確定方法,如經(jīng)驗(yàn)法、試湊法、網(wǎng)格分割法、交叉驗(yàn)證法等需要進(jìn)行大量計(jì)算,且具有人為選擇的盲目性,將使得LSSVM收斂速度減慢,難以得到最優(yōu)解.本文采用ADPSO算法選擇LSSVM最優(yōu)參數(shù).

    利用ADPSO算法選擇LSSVM最優(yōu)參數(shù)(記作LSSVM-ADPSO模型),本質(zhì)上是將LSSVM的構(gòu)建過程與ADPSO算法的尋優(yōu)過程有機(jī)融合,以每個(gè)粒子代表LSSVM的一組參數(shù)組合(γ,σ),在(γ,σ)構(gòu)成的二維目標(biāo)搜索空間中,粒子群通過目標(biāo)函數(shù)確定的適應(yīng)度值來搜尋全局最優(yōu)解.LSSVM-ADPSO模型構(gòu)建及預(yù)測(cè)的具體步驟如下:

    步驟1 數(shù)據(jù)標(biāo)準(zhǔn)化處理.對(duì)原始樣本數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理,計(jì)算公式如下:

    (10)

    其中,θmax、θmin分別為原始樣本數(shù)據(jù)的最大、最小值;θ、θ′分別為標(biāo)準(zhǔn)化前、后的數(shù)據(jù)樣本.將θ′劃分為訓(xùn)練樣本集L與驗(yàn)證樣本集Z2部分.

    步驟2 粒子群初始化.設(shè)置粒子群的各種參數(shù),如種群規(guī)模數(shù)m,慣性權(quán)重的最大值與最小值δmax、δmin,最大迭代步數(shù)tmax等.隨機(jī)給出粒子的初始速度與初始位置,粒子位置對(duì)應(yīng)LSSVM的參數(shù)組合.

    步驟3 適應(yīng)度函數(shù)定義.將L劃分為訓(xùn)練樣本子集U和檢驗(yàn)樣本子集V,定義適應(yīng)度函數(shù)為L(zhǎng)SSVM的訓(xùn)練誤差與檢驗(yàn)誤差之和,如下:

    (11)

    步驟4 粒子群最優(yōu)位置更新.根據(jù)式(11)計(jì)算粒子i的適應(yīng)度值,比較粒子i的當(dāng)前適應(yīng)度值與自身最優(yōu)適應(yīng)度值,若更優(yōu),更新粒子的個(gè)體最優(yōu)位置為當(dāng)前位置;比較每個(gè)粒子的個(gè)體最優(yōu)位置適應(yīng)度值與種群最優(yōu)位置適應(yīng)度值,若更優(yōu),更新種群的最優(yōu)位置為該粒子的個(gè)體最優(yōu)位置.慣性權(quán)重δ按式(7)更新,加速系數(shù)c1(t)和c2(t)分別按式(8)與式(9)更新.

    步驟5 終止條件判斷.若滿足期望誤差或預(yù)設(shè)的迭代步數(shù),則停止搜索,粒子群的全局最優(yōu)位置即為L(zhǎng)SSVM的最優(yōu)參數(shù)組合(γ*,σ*);若不滿足要求,轉(zhuǎn)回步驟3,進(jìn)行新一輪搜索.

    步驟6 模型構(gòu)建與預(yù)測(cè).根據(jù)最優(yōu)參數(shù)組合(γ*,σ*)建立LSSVM并預(yù)測(cè)高速公路交通量,再將預(yù)測(cè)值轉(zhuǎn)化為原始高速公路交通量預(yù)測(cè)值,轉(zhuǎn)化公式如下:

    (12)

    4 實(shí)例計(jì)算

    4.1 指標(biāo)選取

    以某市高速公路交通量為例,驗(yàn)證LSSVM-ADPSO模型的有效性.根據(jù)文獻(xiàn)[15],將高速公路車流量作為高速公路交通量的量化指標(biāo),選取如下11個(gè)社會(huì)經(jīng)濟(jì)指標(biāo)作為高速公路交通量的經(jīng)濟(jì)影響因素:地區(qū)生產(chǎn)總值(GDP)、地區(qū)固定資產(chǎn)投資總額、地區(qū)進(jìn)出口總額、地區(qū)工業(yè)生產(chǎn)總值、地區(qū)第一產(chǎn)業(yè)總產(chǎn)值、地區(qū)第二產(chǎn)業(yè)總產(chǎn)值、地區(qū)第三產(chǎn)業(yè)總產(chǎn)值、地區(qū)總?cè)丝跀?shù)、地區(qū)社會(huì)消費(fèi)品總額、地區(qū)貨物運(yùn)輸量、地區(qū)鐵路貨物運(yùn)輸量.選取1993—2014年的以上指標(biāo)數(shù)據(jù)進(jìn)行實(shí)例計(jì)算,相應(yīng)數(shù)據(jù)來源于《國(guó)家統(tǒng)計(jì)局網(wǎng)站》.

    4.2 模型訓(xùn)練與預(yù)測(cè)

    將標(biāo)準(zhǔn)化到[0,1]區(qū)間的樣本數(shù)據(jù)劃分為訓(xùn)練樣本集L(1993—2008年數(shù)據(jù))和驗(yàn)證樣本集Z(2009—2014年數(shù)據(jù)).ADPSO算法優(yōu)化LSSVM時(shí),其自身控制參數(shù)設(shè)定如下:群體規(guī)模m設(shè)為10;最大迭代步數(shù)tmax設(shè)為20;慣性權(quán)重的最大值δmax、最小值δmin分別設(shè)為0.9和0.1.由于ADPSO算法為隨機(jī)優(yōu)化算法,每次優(yōu)化的結(jié)果可能在一定范圍內(nèi)產(chǎn)生偏差.為減少ADPSO算法優(yōu)化偏差對(duì)LSSVM參數(shù)選擇的影響,通過ADPSO算法連續(xù)優(yōu)化LSSVM參數(shù)多次,選擇最小訓(xùn)練誤差對(duì)應(yīng)的參數(shù)為最優(yōu)參數(shù)組合(γ*,σ*).根據(jù)(γ*,σ*)構(gòu)建LSSVM并預(yù)測(cè)2009—2014年的高速公路交通量,再將預(yù)測(cè)值轉(zhuǎn)化為原始高速公路交通量的預(yù)測(cè)值.

    基于相同數(shù)據(jù)樣本集,分別利用AIWPSO算法優(yōu)化參數(shù)的LSSVM(記作LSSVM-AIWPSO模型)、DACPSO算法優(yōu)化參數(shù)的LSSVM(記作LSSVM-DACPSO模型)預(yù)測(cè)高速公路交通量.AIWPSO算法優(yōu)化LSSVM時(shí),其自身控制參數(shù)設(shè)定如下:群體規(guī)模m設(shè)為10;最大迭代步數(shù)tmax設(shè)為20;慣性權(quán)重的最大值δmax、最小值δmin分別設(shè)為0.9和0.1;加速系數(shù)c1、c2均設(shè)為2.DACPSO算法優(yōu)化LSSVM時(shí),其自身控制參數(shù)設(shè)定如下:群體規(guī)模m設(shè)為10;最大迭代步數(shù)tmax設(shè)為20;慣性權(quán)重的最大值δmax、最小值δmin分別設(shè)為0.9和0.1.將LSSVM-AIWPSO模型和LSSVM-DACPSO模型的預(yù)測(cè)結(jié)果與LSSVM-ADPSO模型進(jìn)行比較.

    4.3 結(jié)果分析

    選用4種損失函數(shù)評(píng)價(jià)各模型的預(yù)測(cè)性能:均方根誤差(root mean squared error,RMSE)、平均絕對(duì)誤差(mean absolute error,MAE)、平均百分比誤差(mean percentage error,MPE),以及西爾統(tǒng)計(jì)量(THEIL),它們分別定義如下:

    (13)

    (14)

    (15)

    (16)

    由表1可知,LSSVM-ADPSO模型的預(yù)測(cè)性能優(yōu)于LSSVM-AIWPSO模型和LSSVM-DACPSO模型,其最大、最小預(yù)測(cè)相對(duì)誤差,以及4個(gè)損失函數(shù)值均小于LSSVM-AIWPSO模型和LSSVM-DACPSO模型的對(duì)應(yīng)值.這主要是由于ADPSO算法依據(jù)適應(yīng)度值同時(shí)更新慣性權(quán)重和加速系數(shù)的值,顯著增強(qiáng)了算法的尋優(yōu)能力,能正確搜索到LSSVM的全局最優(yōu)參數(shù),進(jìn)而提高了LSSVM的預(yù)測(cè)性能.此外,LSSVM-ADPSO模型2009—2010年的預(yù)測(cè)相對(duì)誤差為-3.27%和3.19%,2011—2012年的預(yù)測(cè)相對(duì)誤差有所減小,分別為1.08%和0.24%,2013—2014年的預(yù)測(cè)相對(duì)誤差又明顯增大,分別4.15%和-4.57%,但仍明顯高于2009—2010年的預(yù)測(cè)相對(duì)誤差,表明LSSVM-ADPSO模型適合于高速公路交通量的短期預(yù)測(cè).

    通過比較LSSVM-AIWPSO模型和LSSVM-DACPSO模型可知,LSSVM-AIWPSO模型的最大預(yù)測(cè)相對(duì)誤差稍大于LSSVM-DACPSO模型的對(duì)應(yīng)值、最小預(yù)測(cè)相對(duì)誤差小于LSSVM-AIWPSO模型的對(duì)應(yīng)值;LSSVM-AIWPSO模型的MAE、MPE值小于LSSVM-DACPSO模型的對(duì)應(yīng)值,而RMSE、THEIL值又大于LSSVM-DACPSO模型的對(duì)應(yīng)值.因此,暫時(shí)無法確定兩模型預(yù)測(cè)性能的優(yōu)劣,需進(jìn)一步研究.

    表1 不同模型預(yù)測(cè)結(jié)果比較

    10TV代表 10 thousand vehicles.

    圖1 不同模型預(yù)測(cè)值比較Fig.1 Comparison of forecasting values of different models

    由圖1可看出,3種模型均較準(zhǔn)確地預(yù)測(cè)出高速公路交通量的逐年增長(zhǎng)趨勢(shì),在整個(gè)預(yù)測(cè)樣本期內(nèi),LSSVM-ADPSO模型更好地揭示出高速公路交通量的變化特征,特別是在變化幅度比較大的2012年和2014年,其預(yù)測(cè)值比其他兩模型更接近于實(shí)際值;而LSSVM-AIWPSO模型和LSSVM-DACPSO模型2014年的預(yù)測(cè)值明顯偏離實(shí)際值.

    5 結(jié)論

    提出高速公路交通量的LSSVM-ADPSO預(yù)測(cè)模型,LSSVM用于預(yù)測(cè)高速公路交通量,ADPSO算法用作LSSVM最優(yōu)參數(shù)的選擇方法.對(duì)某市高速公路交通實(shí)例計(jì)算表明,LSSVM-ADPSO模型的預(yù)測(cè)性能優(yōu)于LSSVM-AIWPSO模型和LSSVM-DACPSO模型,對(duì)變化幅度較大的高速公路交通量具有更優(yōu)的預(yù)測(cè)效果,適用于高速公路交通量的短期預(yù)測(cè).LSSVM-AIWPSO模型與LSSVM-DACPSO模型預(yù)測(cè)性能的優(yōu)劣有待進(jìn)一步研究.

    [1] 徐建閩,尹宏賓.廣珠高速公路交通量預(yù)測(cè)[J].系統(tǒng)工程理論與實(shí)踐,1998,18(4):122-125.DOI:10.3321/j.issn:1000-6788.1998.04.021. XU J M,YIN H B.Guangzhou-Zhuhai expressway traffic volume forecast[J].Systems Engineering-Theory & Practice,1998,18(4):122-125.DOI:10.3321/j.issn:1000-6788.1998.04.021.

    [2] 朱順應(yīng),王紅,李關(guān)壽.路段上短時(shí)間區(qū)段內(nèi)交通量預(yù)測(cè)ARIMA模型[J].重慶交通大學(xué)學(xué)報(bào)(自然科學(xué)版),2003,22(1):76-77,95.DOI:10.3969/j.issn.1674-0696.2003.01.019. ZHU S Y,WANG H,LI G S.The ARIMA model used in forecasting of traffic volume in short interval on the link[J].Journal of Chongqing Jiaotong University,2003,22(1):76-77,95.DOI:10.3969/j.issn.1674-0696.2003.01.019.

    [3] 高連生,易誕,毛娜,等.基于改進(jìn)的GM(1,1)的長(zhǎng)期交通量預(yù)測(cè)模型[J].鐵道科學(xué)與工程學(xué)報(bào),2015,12(1):203-207.DOI:10.3969/j.issn.1672-7029.2015.01.033. GAO L S,YI D,MAO N,LI L.A forecast model for long -term traffic volu me based on improved GM(1,1)model[J].Journal of Railway Science and Engineering,2015,12(1):203-207.DOI:10.3969/j.issn.1672-7029.2015.01.033.

    [4] 王炎,王華.基于BP-NN和遺傳算法的高速公路交通量預(yù)測(cè)[J].計(jì)算機(jī)工程與應(yīng)用,2006(4):226-228.DOI:10.3321/j.issn:1002-8331.2006.04.070. WANG Y,WANG H.Prediction of highway traffic based on BP-neural network and genetic algorithms[J].Computer Engineering and Applications,2006(4):226-228.DOI:10.3321/j.issn:1002-8331.2006.04.070.

    [5] VAPNIK V N.An overview of statistical learning theory[J].IEEE Transactions on Neural Networks,1999,l0(5):988-999.DOI:10.1109/72.788640.

    [6] HU W B,YAN L P,LIU K Z,et al.A Short-term traffic flow forecasting method based on the hybrid PSO-SVR[J].Neural Processing Letters,2016,43(1):155-172.DOI:10.1007/s11063-015-9409-6.

    [7] AHN J,KO E,KIM E Y.Highway traffic flow prediction using support vector regression and Bayesian classifier[Z].The IEEE International Conference on Big Data and Smart Computing (BigComp),HongKong,2016.DOI:10.1109/BIGCOMP.2016.7425919.

    [8] SUYKENS J A K,VANDEVALLE J.Least squares support vector machine classifiers[J].Neural Processing Letters,1999,9(3):293-300.DOI:10.1023/A:1018628609742.

    [9] 戢小輝.基于灰色關(guān)聯(lián)的LS-SVM道路交通事故預(yù)測(cè)[J].計(jì)算機(jī)應(yīng)用研究,2016,33(3):806-809.DOI:10.3969/j.issn.1001-3695.2016.03.037. JI X H.Forecast model of road traffic accidents based on LS-SVMwith grey correlation analysis[J].Application Research of Computers,2016,33(3):806-809.DOI:10.3969/j.issn.1001-3695.2016.03.037.

    [10] 鄧佳佳,劉爽.基于免疫模糊聚類的LSSVM在短期負(fù)荷預(yù)測(cè)中的應(yīng)用[J].河北大學(xué)學(xué)報(bào)(自然科學(xué)版),2012,32(3):234-239.DOI:10.3969/j.issn.1000-1565.2012.03.003. DENG J J,LIU S.Application of LSSVM based on immune fuzzy clustering algorithm to short-term load forecasting[J].Journal of Hebei University(Natural Science Edition),2012,32(3):234-239.DOI:10.3969/j.issn.1000-1565.2012.03.003.

    [11] POLI R,KENNEDY J,BLACKWELL T.Particle swarm optimization:An overview[J].Swarm Intelligence,2007,1(1):33-57.DOI:10.1007/s11721-007-0002-0.

    [12] SHI Y,EBERHART R C.Empirical study of particle swarm optimization[Z].The IEEE International Congress on Evolutionary Computation,Piscataway,NJ,USA,1999.DOI:10.1109/CEC.1999.785511.

    [13] 紀(jì)震,廖惠連,吳青華.粒子群算法及應(yīng)用[M].北京:科學(xué)出版社,2009.

    [14] 姜長(zhǎng)弘,張永恒,王盛慧.基于自適應(yīng)粒子群算法的篦冷機(jī)電液位置飼服系統(tǒng)PID參數(shù)優(yōu)化[J].液壓與氣動(dòng),2016(7):44-49.DOI:10.11832/j.issn.1000-4858.2016.07.009. JIANG C H,ZHANG Y H,WANG S H.PID parameter optimization of the gate cooling electro-hydraulic position servo system based on adaptive particle swarm optimization algorithm[J].Chinese Hydraulics & Pneumatics,2016(7):44-49.DOI:10.11832/j.issn.1000-4858.2016.07.009.

    [15] 趙澤輝.高速公路交通量預(yù)測(cè)及山區(qū)高速公路后評(píng)價(jià)研究[D].大連:大連理工大學(xué),2015. ZHAO Z H.Highway traffic forecasting and post-evaluation research of mountain highway[D].Dalian:Dalian University of Technology,2015.

    (責(zé)任編輯:孟素蘭)

    A forecasting method of highway traffic flow using LSSVM optimized by ADPSO algorithm

    SI Wenjing1,FENG Xibo2,GENG Liyan3,4,ZHANG Zhanfu5
    (1.Construction Engineering Department,North China Institute of Aerospace Engineering,Langfang 065000,China;2.Hebei Province Expressway Langfang Beisanxian County Management Department,Langfang 065000,China;3.School of Economics and Management,Shijiazhuang Tiedao University,Shijiazhuang 050043,China;4.Business School,Manchester Metropolitan University,Manchester,M15 6BH,UK;5.Sifang College,Shijiazhuang Tiedao University,Shijiazhuang 051132,China)

    There is a complex nonlinear relationship between highway traffic flow and its influencing factors.Combing least squares support vector machines (LSSVM) with adaptive dynamic particle swarm optimization (ADPSO) algorithm,this paper proposed a new highway traffic flow forecasting method based on LSSVM optimized by ADPSO algorithm.Highway traffic flow was forecasted by LSSVM with the advantages of easy modeling and high precision.And the optimal parameters of LSSVM were selected based on the good optimization ability of ADPSO algorithm.An example analysis on the highway traffic flow in a city was performed to test the effectiveness of LSSVM-ADPSO model.The results indicate that the proposed method has better highway traffic flow forecasting performance and is suitable for short-term highway traffic flow forecasting.

    highway;traffic flow forecasting;adaptive dynamic particle swarm optimization algorithm;least squares support vector machines

    2016-12-25

    國(guó)家自然科學(xué)基金青年項(xiàng)目(61503261);河北省交通運(yùn)輸廳科技計(jì)劃項(xiàng)目(Y-2010024);北華航天工業(yè)學(xué)院科研基金項(xiàng)目(KY-2015-09); 河北省軟科學(xué)研究計(jì)劃項(xiàng)目(15456106D);河北省高等學(xué)校青年拔尖人才計(jì)劃項(xiàng)目(BJ2014097);河北省社會(huì)科學(xué)發(fā)展重點(diǎn)課題(2015020206);國(guó)家留學(xué)基金委(CSC)公派留學(xué)地方合作項(xiàng)目(201608130165);河北省高校人文社會(huì)科學(xué)重點(diǎn)研究基地石家莊鐵道大學(xué)工程建設(shè)管理研究中心資助項(xiàng)目;河北省軟科學(xué)工程建設(shè)管理研究基地資助項(xiàng)目;河北省重點(diǎn)學(xué)科管理科學(xué)與工程資助項(xiàng)目

    司文靜(1980—),女,河北遷安人,北華航天工業(yè)學(xué)院講師,主要從事道路與橋梁工程方面的研究. E-mail:siwenjing1021@sina.com

    封喜波(1979—),男,河北平山人,河北省高速公路廊坊北三縣管理處高級(jí)工程師,主要從事交通工程方面的研究. E-mail:6511129@qq.com

    10.3969/j.issn.1000-1565.2017.03.013

    U491.14

    A

    1000-1565(2017)03-0302-07

    猜你喜歡
    交通量適應(yīng)度粒子
    改進(jìn)的自適應(yīng)復(fù)制、交叉和突變遺傳算法
    基于ETC門架數(shù)據(jù)的高速公路交通量轉(zhuǎn)換探究
    基于粒子群優(yōu)化的橋式起重機(jī)模糊PID控制
    基于動(dòng)態(tài)差法的交通量監(jiān)測(cè)技術(shù)應(yīng)用
    基于粒子群優(yōu)化極點(diǎn)配置的空燃比輸出反饋控制
    高速公路補(bǔ)償交通量模型研究
    基于四階段法的公路交通量預(yù)測(cè)研究
    北方交通(2016年12期)2017-01-15 13:52:51
    基于空調(diào)導(dǎo)風(fēng)板成型工藝的Kriging模型適應(yīng)度研究
    基于Matlab的α粒子的散射實(shí)驗(yàn)?zāi)M
    物理與工程(2014年4期)2014-02-27 11:23:08
    基于兩粒子糾纏態(tài)隱形傳送四粒子GHZ態(tài)
    国产高清三级在线| 日韩欧美精品免费久久| videos熟女内射| 免费av毛片视频| 狠狠精品人妻久久久久久综合| 一级a做视频免费观看| 国产永久视频网站| 亚州av有码| 欧美日韩视频高清一区二区三区二| 王馨瑶露胸无遮挡在线观看| 亚洲欧美日韩卡通动漫| 国产探花极品一区二区| 日韩 亚洲 欧美在线| 免费黄色在线免费观看| 简卡轻食公司| 乱系列少妇在线播放| 99久久精品国产国产毛片| 中国国产av一级| 亚洲欧美成人精品一区二区| 一区二区三区免费毛片| 一级a做视频免费观看| 少妇丰满av| 国产精品无大码| 80岁老熟妇乱子伦牲交| 狂野欧美白嫩少妇大欣赏| 欧美性感艳星| 又粗又硬又长又爽又黄的视频| 视频中文字幕在线观看| 卡戴珊不雅视频在线播放| 日日撸夜夜添| 成人亚洲精品一区在线观看 | 日本wwww免费看| 日韩欧美精品v在线| 精品久久久久久久久av| 2021少妇久久久久久久久久久| 91精品伊人久久大香线蕉| 1000部很黄的大片| 国产 一区 欧美 日韩| 日本一二三区视频观看| 成人特级av手机在线观看| 日韩成人av中文字幕在线观看| 国产成人福利小说| 日韩人妻高清精品专区| 美女被艹到高潮喷水动态| 午夜激情福利司机影院| 一级黄片播放器| 天堂网av新在线| 久久久精品免费免费高清| 极品教师在线视频| 国产精品精品国产色婷婷| 2021天堂中文幕一二区在线观| 国产乱人视频| 国产午夜福利久久久久久| 2022亚洲国产成人精品| 三级男女做爰猛烈吃奶摸视频| 亚洲精品国产av蜜桃| 搡老乐熟女国产| 久久久久久久久大av| 黄片wwwwww| 少妇人妻 视频| 我的老师免费观看完整版| 欧美一级a爱片免费观看看| 免费黄色在线免费观看| 高清欧美精品videossex| a级毛色黄片| 国产精品久久久久久久久免| 午夜老司机福利剧场| www.色视频.com| 全区人妻精品视频| 日韩国内少妇激情av| 国产av码专区亚洲av| 在线观看免费高清a一片| 国产乱来视频区| 女的被弄到高潮叫床怎么办| 人妻少妇偷人精品九色| 有码 亚洲区| 中文字幕av成人在线电影| 麻豆久久精品国产亚洲av| 成年版毛片免费区| 在线天堂最新版资源| 国产成人免费观看mmmm| 色视频在线一区二区三区| 日韩成人av中文字幕在线观看| 亚洲欧美精品专区久久| 欧美激情久久久久久爽电影| 国产精品99久久久久久久久| 国产综合懂色| 国产精品不卡视频一区二区| 午夜福利在线在线| 久久99热这里只有精品18| 免费不卡的大黄色大毛片视频在线观看| 啦啦啦在线观看免费高清www| 边亲边吃奶的免费视频| 国产精品.久久久| 久久久久久久久久久丰满| 亚洲人与动物交配视频| 免费黄色在线免费观看| 美女视频免费永久观看网站| 少妇 在线观看| 国产成人免费无遮挡视频| 国产精品人妻久久久影院| 色吧在线观看| 亚洲精品日韩av片在线观看| eeuss影院久久| 国产熟女欧美一区二区| 欧美变态另类bdsm刘玥| 国产黄色视频一区二区在线观看| 久久久久九九精品影院| 免费观看无遮挡的男女| 国产免费一区二区三区四区乱码| 观看美女的网站| 黄片无遮挡物在线观看| 夜夜看夜夜爽夜夜摸| 亚洲熟女精品中文字幕| 精品国产三级普通话版| 日韩三级伦理在线观看| 午夜激情久久久久久久| 国产国拍精品亚洲av在线观看| 免费播放大片免费观看视频在线观看| 三级国产精品欧美在线观看| 亚洲av电影在线观看一区二区三区 | 免费高清在线观看视频在线观看| 日本熟妇午夜| 乱码一卡2卡4卡精品| 丝袜脚勾引网站| 亚洲精品日本国产第一区| 在线免费观看不下载黄p国产| 少妇的逼水好多| 我的老师免费观看完整版| 美女主播在线视频| 色哟哟·www| 欧美zozozo另类| 久久97久久精品| 精品久久久久久久久av| 2018国产大陆天天弄谢| 日韩大片免费观看网站| 久久久久久国产a免费观看| 国产一区亚洲一区在线观看| 久久久久久久久大av| 亚洲一级一片aⅴ在线观看| 美女xxoo啪啪120秒动态图| 欧美 日韩 精品 国产| 女人久久www免费人成看片| 草草在线视频免费看| 大香蕉97超碰在线| 久久久久久久久久成人| 丰满人妻一区二区三区视频av| 国产片特级美女逼逼视频| 日韩在线高清观看一区二区三区| 亚洲av.av天堂| 亚洲国产精品成人久久小说| 伊人久久精品亚洲午夜| 少妇裸体淫交视频免费看高清| 成年免费大片在线观看| 亚洲第一区二区三区不卡| 啦啦啦在线观看免费高清www| 日韩一区二区三区影片| 国产黄片美女视频| 亚洲欧洲日产国产| 亚洲美女搞黄在线观看| 久久久久久九九精品二区国产| 观看美女的网站| 新久久久久国产一级毛片| 两个人的视频大全免费| 久久久久久久精品精品| 久久久a久久爽久久v久久| 久久久亚洲精品成人影院| 日本爱情动作片www.在线观看| 亚洲美女搞黄在线观看| 欧美精品一区二区大全| 久久久久久久精品精品| 观看免费一级毛片| 午夜福利在线在线| 国产成人91sexporn| 国产亚洲午夜精品一区二区久久 | 自拍偷自拍亚洲精品老妇| 中文天堂在线官网| 91aial.com中文字幕在线观看| 搡老乐熟女国产| 国产精品99久久99久久久不卡 | 欧美高清成人免费视频www| 夫妻性生交免费视频一级片| 亚洲国产欧美人成| 精品亚洲乱码少妇综合久久| 白带黄色成豆腐渣| 五月伊人婷婷丁香| 久久久成人免费电影| 日韩av在线免费看完整版不卡| 国产高潮美女av| 国产免费一级a男人的天堂| 日韩人妻高清精品专区| 亚洲一级一片aⅴ在线观看| 国语对白做爰xxxⅹ性视频网站| 成年av动漫网址| 欧美区成人在线视频| 91久久精品国产一区二区三区| 欧美3d第一页| 日韩成人伦理影院| 在线a可以看的网站| 女人久久www免费人成看片| 在线a可以看的网站| 十八禁网站网址无遮挡 | 各种免费的搞黄视频| 大片免费播放器 马上看| 欧美极品一区二区三区四区| 新久久久久国产一级毛片| 狠狠精品人妻久久久久久综合| 国产精品99久久99久久久不卡 | 精品一区在线观看国产| 亚洲精品456在线播放app| 亚洲第一区二区三区不卡| 免费电影在线观看免费观看| 欧美成人一区二区免费高清观看| 男人狂女人下面高潮的视频| 亚洲av中文字字幕乱码综合| 亚洲精品自拍成人| 丰满人妻一区二区三区视频av| xxx大片免费视频| 一本色道久久久久久精品综合| 亚洲精品色激情综合| 久久精品综合一区二区三区| 国产免费一区二区三区四区乱码| 成人二区视频| 久久鲁丝午夜福利片| 亚洲无线观看免费| 国产成人freesex在线| 99久久九九国产精品国产免费| 国产精品久久久久久久电影| 男女啪啪激烈高潮av片| 国产男人的电影天堂91| 欧美丝袜亚洲另类| 1000部很黄的大片| 18禁裸乳无遮挡免费网站照片| 午夜免费观看性视频| 精品国产一区二区三区久久久樱花 | 国产真实伦视频高清在线观看| 欧美精品人与动牲交sv欧美| 看非洲黑人一级黄片| 精品午夜福利在线看| 欧美zozozo另类| 久久久久久久大尺度免费视频| 直男gayav资源| av在线播放精品| 中文乱码字字幕精品一区二区三区| 午夜精品一区二区三区免费看| 五月伊人婷婷丁香| 成年免费大片在线观看| 日本黄大片高清| 亚洲色图综合在线观看| 国产精品不卡视频一区二区| 2018国产大陆天天弄谢| 日韩一区二区三区影片| 美女被艹到高潮喷水动态| 特大巨黑吊av在线直播| 久久99热6这里只有精品| 在线免费观看不下载黄p国产| 亚洲丝袜综合中文字幕| 人人妻人人看人人澡| 大话2 男鬼变身卡| 熟妇人妻不卡中文字幕| 亚洲av欧美aⅴ国产| videossex国产| 亚洲在线观看片| 日本-黄色视频高清免费观看| 国产精品av视频在线免费观看| 少妇被粗大猛烈的视频| 哪个播放器可以免费观看大片| 我的女老师完整版在线观看| 国产视频首页在线观看| 日本一二三区视频观看| 新久久久久国产一级毛片| 全区人妻精品视频| 久久久国产一区二区| 亚洲精品日本国产第一区| 极品教师在线视频| 国产在线男女| 搞女人的毛片| 欧美激情国产日韩精品一区| 神马国产精品三级电影在线观看| 国产精品三级大全| 久久99热6这里只有精品| 成人高潮视频无遮挡免费网站| 边亲边吃奶的免费视频| 精品国产三级普通话版| 欧美日韩视频高清一区二区三区二| 国产乱来视频区| 麻豆成人av视频| 高清午夜精品一区二区三区| 永久网站在线| 欧美日韩在线观看h| 视频区图区小说| 最新中文字幕久久久久| 一级二级三级毛片免费看| 久久精品久久久久久久性| 欧美xxⅹ黑人| 国产成人免费无遮挡视频| 亚洲性久久影院| 亚洲国产欧美人成| 欧美一级a爱片免费观看看| 精品一区二区免费观看| 99久久人妻综合| 国产黄色免费在线视频| 一级二级三级毛片免费看| a级毛色黄片| 国产大屁股一区二区在线视频| 午夜免费男女啪啪视频观看| 99re6热这里在线精品视频| av在线观看视频网站免费| 亚洲国产精品成人综合色| 水蜜桃什么品种好| 久久久精品免费免费高清| 欧美日韩视频精品一区| 婷婷色av中文字幕| 国产午夜精品久久久久久一区二区三区| 国产成人aa在线观看| av天堂中文字幕网| 精品久久久久久电影网| 五月玫瑰六月丁香| 啦啦啦在线观看免费高清www| a级毛片免费高清观看在线播放| 亚洲精品久久午夜乱码| 国产黄色视频一区二区在线观看| av在线老鸭窝| 亚洲精品自拍成人| 成人特级av手机在线观看| 午夜激情久久久久久久| 亚洲精品色激情综合| 国产精品熟女久久久久浪| 欧美变态另类bdsm刘玥| 蜜臀久久99精品久久宅男| av一本久久久久| 18禁动态无遮挡网站| 深夜a级毛片| 99热6这里只有精品| 国产免费一级a男人的天堂| 免费观看的影片在线观看| 人体艺术视频欧美日本| 美女cb高潮喷水在线观看| 99视频精品全部免费 在线| 国产成人a∨麻豆精品| 黄色一级大片看看| 久久韩国三级中文字幕| 波野结衣二区三区在线| 插阴视频在线观看视频| 久久久久久久久久人人人人人人| 精品久久久久久久末码| 国产av码专区亚洲av| av国产久精品久网站免费入址| 久久久久久久久久久免费av| 大片免费播放器 马上看| 看十八女毛片水多多多| 美女内射精品一级片tv| 最近中文字幕2019免费版| 亚洲高清免费不卡视频| 99视频精品全部免费 在线| 国模一区二区三区四区视频| 夜夜看夜夜爽夜夜摸| 亚洲va在线va天堂va国产| 夫妻性生交免费视频一级片| 日韩免费高清中文字幕av| 精品国产露脸久久av麻豆| 亚洲高清免费不卡视频| 熟女电影av网| 色婷婷久久久亚洲欧美| 日韩人妻高清精品专区| 三级国产精品片| 久久99精品国语久久久| 啦啦啦啦在线视频资源| 黄片wwwwww| 成人亚洲精品av一区二区| 亚洲精品成人久久久久久| 欧美成人精品欧美一级黄| 精品人妻视频免费看| 香蕉精品网在线| 少妇被粗大猛烈的视频| 伊人久久国产一区二区| 国产精品人妻久久久影院| 精品久久久噜噜| 大片电影免费在线观看免费| 看黄色毛片网站| 国产高清国产精品国产三级 | www.av在线官网国产| 亚洲aⅴ乱码一区二区在线播放| 天天一区二区日本电影三级| 六月丁香七月| 亚洲第一区二区三区不卡| 男人和女人高潮做爰伦理| 国产免费一级a男人的天堂| 简卡轻食公司| 亚洲真实伦在线观看| 日本黄色片子视频| eeuss影院久久| 99热这里只有精品一区| 成人国产av品久久久| 18禁裸乳无遮挡动漫免费视频 | 欧美性猛交╳xxx乱大交人| 日日摸夜夜添夜夜添av毛片| 久久99热这里只频精品6学生| 国产成人freesex在线| 一二三四中文在线观看免费高清| av福利片在线观看| 男人和女人高潮做爰伦理| 男的添女的下面高潮视频| 久久久久国产精品人妻一区二区| 国产av码专区亚洲av| 免费观看的影片在线观看| 午夜福利网站1000一区二区三区| 99热全是精品| eeuss影院久久| 亚洲图色成人| 大码成人一级视频| 少妇的逼水好多| 国产精品女同一区二区软件| 狂野欧美激情性bbbbbb| 日日啪夜夜撸| 极品教师在线视频| 久久久久久久久大av| 国产成人精品一,二区| 国产黄a三级三级三级人| 国产又色又爽无遮挡免| 日本爱情动作片www.在线观看| 白带黄色成豆腐渣| 毛片女人毛片| 亚洲激情五月婷婷啪啪| 香蕉精品网在线| 国产精品秋霞免费鲁丝片| tube8黄色片| 午夜激情久久久久久久| 尾随美女入室| 久久韩国三级中文字幕| 亚洲欧美一区二区三区黑人 | 午夜免费鲁丝| 成人免费观看视频高清| 国产精品熟女久久久久浪| 国内少妇人妻偷人精品xxx网站| 一区二区三区乱码不卡18| 小蜜桃在线观看免费完整版高清| av在线天堂中文字幕| 看免费成人av毛片| 婷婷色综合www| 成人毛片a级毛片在线播放| 国产精品一区www在线观看| 精品国产一区二区三区久久久樱花 | 亚洲精品久久久久久婷婷小说| 免费观看的影片在线观看| 中文天堂在线官网| 又爽又黄无遮挡网站| 草草在线视频免费看| 国产成人a区在线观看| 欧美另类一区| 久久久久国产网址| 国产亚洲av嫩草精品影院| 国内少妇人妻偷人精品xxx网站| 少妇人妻久久综合中文| 乱码一卡2卡4卡精品| 一级片'在线观看视频| 欧美高清成人免费视频www| 日本爱情动作片www.在线观看| 国产精品一区二区在线观看99| 婷婷色综合www| 久久久亚洲精品成人影院| 欧美3d第一页| 国产成人91sexporn| 联通29元200g的流量卡| 欧美日韩视频高清一区二区三区二| 看十八女毛片水多多多| 国产老妇伦熟女老妇高清| 亚洲真实伦在线观看| 亚洲精品国产av蜜桃| 各种免费的搞黄视频| 少妇的逼水好多| av在线播放精品| 一级毛片黄色毛片免费观看视频| 久久精品夜色国产| 免费播放大片免费观看视频在线观看| 国产成人精品一,二区| 狂野欧美激情性xxxx在线观看| 十八禁网站网址无遮挡 | 久久久久久久精品精品| 全区人妻精品视频| 韩国av在线不卡| 日韩精品有码人妻一区| 亚洲在线观看片| 欧美性猛交╳xxx乱大交人| 网址你懂的国产日韩在线| 老司机影院成人| 亚洲精品,欧美精品| 精品99又大又爽又粗少妇毛片| 亚洲av免费在线观看| 日韩在线高清观看一区二区三区| 大片免费播放器 马上看| 纵有疾风起免费观看全集完整版| 亚洲精品456在线播放app| 我的女老师完整版在线观看| 免费av观看视频| 亚洲欧美日韩卡通动漫| 欧美激情久久久久久爽电影| 久久精品久久精品一区二区三区| 在线精品无人区一区二区三 | 中国三级夫妇交换| av线在线观看网站| 少妇人妻 视频| 精品99又大又爽又粗少妇毛片| 日本午夜av视频| 超碰97精品在线观看| 国产真实伦视频高清在线观看| 又爽又黄无遮挡网站| 街头女战士在线观看网站| 国产69精品久久久久777片| 日本与韩国留学比较| 国产69精品久久久久777片| 日本与韩国留学比较| 日本三级黄在线观看| 一级a做视频免费观看| 亚洲国产日韩一区二区| 日韩亚洲欧美综合| 黄片无遮挡物在线观看| tube8黄色片| 欧美一区二区亚洲| 欧美日韩综合久久久久久| 免费黄色在线免费观看| 菩萨蛮人人尽说江南好唐韦庄| 国产伦精品一区二区三区四那| 最近中文字幕高清免费大全6| 日本-黄色视频高清免费观看| 久久99热6这里只有精品| 日日撸夜夜添| 国产伦理片在线播放av一区| 亚洲精品国产av成人精品| 日日啪夜夜撸| 中文字幕人妻熟人妻熟丝袜美| 狂野欧美激情性xxxx在线观看| 波野结衣二区三区在线| 亚洲精品成人av观看孕妇| 国产精品秋霞免费鲁丝片| 欧美变态另类bdsm刘玥| videossex国产| 午夜福利在线在线| 乱码一卡2卡4卡精品| 国产精品精品国产色婷婷| 人体艺术视频欧美日本| 日本黄色片子视频| 国产一区亚洲一区在线观看| 国产成人freesex在线| 中文乱码字字幕精品一区二区三区| 精华霜和精华液先用哪个| 国产色婷婷99| 亚洲aⅴ乱码一区二区在线播放| 中文天堂在线官网| 九九爱精品视频在线观看| av一本久久久久| 色网站视频免费| 韩国高清视频一区二区三区| 国产欧美另类精品又又久久亚洲欧美| 伊人久久精品亚洲午夜| av女优亚洲男人天堂| 久久久久精品久久久久真实原创| 国产在线男女| 能在线免费看毛片的网站| 男女边摸边吃奶| 亚洲精品一区蜜桃| 欧美性感艳星| 身体一侧抽搐| 午夜精品国产一区二区电影 | tube8黄色片| 天美传媒精品一区二区| 亚洲精品aⅴ在线观看| 久久久久久久大尺度免费视频| 日本猛色少妇xxxxx猛交久久| 国产白丝娇喘喷水9色精品| 免费黄色在线免费观看| 成人国产av品久久久| 听说在线观看完整版免费高清| 国产av码专区亚洲av| 女人被狂操c到高潮| 国产男人的电影天堂91| 伦理电影大哥的女人| 久久久久国产精品人妻一区二区| 七月丁香在线播放| 内射极品少妇av片p| 五月天丁香电影| 国内揄拍国产精品人妻在线| 男人添女人高潮全过程视频| 国产精品久久久久久av不卡| 亚洲一区二区三区欧美精品 | 噜噜噜噜噜久久久久久91| 在线天堂最新版资源| 国产精品久久久久久久久免| 精品久久国产蜜桃| 大陆偷拍与自拍| 国产精品久久久久久精品古装| 日本猛色少妇xxxxx猛交久久| 老司机影院毛片| 超碰97精品在线观看| 精品久久久久久久久av| 久久久久久久久久人人人人人人| 欧美3d第一页| 国产精品三级大全| 最近中文字幕2019免费版| 一级av片app| 欧美少妇被猛烈插入视频| 夜夜看夜夜爽夜夜摸| 国产欧美日韩一区二区三区在线 | 女人久久www免费人成看片| 久久人人爽人人片av| 精品午夜福利在线看| 成人亚洲精品一区在线观看 | av黄色大香蕉| 国产黄片视频在线免费观看| 亚洲色图av天堂| 一级黄片播放器| 亚洲自偷自拍三级| 国产黄片美女视频| 干丝袜人妻中文字幕| 97在线视频观看|