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

    年最大洪峰流量群居蜘蛛優(yōu)化投影尋蹤預(yù)測(cè)模型

    2017-06-05 22:02:54王文川劉惠敏雷冠軍劉寬邱林
    南水北調(diào)與水利科技 2017年2期
    關(guān)鍵詞:預(yù)測(cè)

    王文川 劉惠敏 雷冠軍 劉寬 邱林

    摘要:年最大洪峰流量預(yù)測(cè),受較多的復(fù)雜因素的影響,不確定性較強(qiáng),用常規(guī)統(tǒng)計(jì)方法做出準(zhǔn)確預(yù)報(bào)具有較大困難。從水文序列本身出發(fā),提出將投影回歸模型應(yīng)用于年最大洪峰流量預(yù)測(cè),為了更好獲得投影尋蹤模型參數(shù)和預(yù)測(cè)精度,提出了運(yùn)用延遲相關(guān)系數(shù)法確定回歸預(yù)測(cè)因子、群居蜘蛛算法優(yōu)化投影尋蹤模型最佳投影方向參數(shù)n、利用最小二乘法確定多項(xiàng)式的權(quán)系數(shù)c、嶺函數(shù)個(gè)數(shù)M的群居蜘蛛優(yōu)化投影尋蹤年最大洪峰流量預(yù)測(cè)模型,結(jié)合長(zhǎng)江宜昌站(1882年-2004年)的年最大洪峰流量資料進(jìn)行實(shí)例預(yù)測(cè),訓(xùn)練階段平均絕對(duì)相對(duì)誤差為8.61%,預(yù)測(cè)階段平均絕對(duì)相對(duì)誤差為10.5 1%,該模型預(yù)測(cè)效果較好,模型結(jié)果穩(wěn)定,可有效應(yīng)用于年最大洪峰流量預(yù)測(cè)。

    關(guān)鍵詞:SSO算法;參數(shù)投影尋蹤;混合智能;年最大洪峰流量;預(yù)測(cè)

    中圖分類(lèi)號(hào):TV 124 文獻(xiàn)標(biāo)識(shí)碼:A 文章編號(hào):1672-1683(2017)02-0009-06

    鑒于年最大洪峰流量在防洪渡汛工作中具有重要作用,所以其預(yù)測(cè)是中長(zhǎng)期水文預(yù)報(bào)的重點(diǎn)。但由于其涉及較大的時(shí)空尺度范圍,并受到大氣運(yùn)。王渺林在趨勢(shì)分析與周期分析的基礎(chǔ)上對(duì)年最大洪峰流量的演變規(guī)律進(jìn)行研究。金鑫家運(yùn)用小波分析對(duì)水文序列進(jìn)行由粗到細(xì)的分析,分析了年最大洪峰流量的周期和預(yù)報(bào):王文圣等人在小波變換序列的基礎(chǔ)上提出了最近鄰抽樣回歸組合預(yù)報(bào),結(jié)果表明提出的方法好于諧波預(yù)報(bào)法:宋荷花等利用模糊模式識(shí)別模型對(duì)湘江湘潭站年最大洪峰流量進(jìn)行預(yù)報(bào):李紅霞等探討了基于貝葉斯正則化方法的神經(jīng)網(wǎng)絡(luò)年最大洪峰流量預(yù)測(cè)模型:孫樹(shù)青等將統(tǒng)計(jì)相關(guān)性分析與模糊方法相結(jié)合,建立了基于BP神經(jīng)網(wǎng)絡(luò)的最大洪峰流量進(jìn)行預(yù)報(bào)模型,這些工作作為對(duì)最大洪峰流量預(yù)測(cè)建模的有益嘗試,促進(jìn)了其發(fā)展,但人工神經(jīng)網(wǎng)絡(luò)在解決預(yù)報(bào)問(wèn)題方面仍存在網(wǎng)絡(luò)結(jié)構(gòu)選擇困難、全局搜索能力差、易陷入局部最優(yōu)值點(diǎn)且模擬結(jié)果不穩(wěn)定等問(wèn)題,在一定程度上限制了其在實(shí)際生產(chǎn)中的推廣應(yīng)用。趙璀將灰色系統(tǒng)理論運(yùn)用到中長(zhǎng)期水文預(yù)報(bào)中,取得了較好的結(jié)果。張漢雄認(rèn)為年最大洪峰流量序列是一個(gè)純隨機(jī)現(xiàn)象,因而選擇采用馬爾科夫鏈給出預(yù)測(cè)值的大致范圍,但該方法的結(jié)果受序列長(zhǎng)度的影響較明顯。胡致強(qiáng)在自適應(yīng)控制理論發(fā)展的基礎(chǔ)上,運(yùn)用Kalman濾波對(duì)年最大洪峰流量進(jìn)行預(yù)報(bào),如果模型和水文系統(tǒng)不匹配,則會(huì)使預(yù)報(bào)誤差較大。

    投影尋蹤將應(yīng)用數(shù)學(xué)、現(xiàn)代統(tǒng)計(jì)與計(jì)算機(jī)技術(shù)相結(jié)合進(jìn)行數(shù)據(jù)處理,在解決樣本數(shù)量少和維數(shù)較大等問(wèn)題時(shí),具有明顯的優(yōu)勢(shì),在評(píng)估、聚類(lèi)分析、洪水預(yù)報(bào)等領(lǐng)域得到了很好的應(yīng)用。

    投影尋蹤建模過(guò)程是利用線性投影的多個(gè)多元回歸的線性組合,優(yōu)選投影方向a、多項(xiàng)式系數(shù)c和嶺函數(shù)個(gè)數(shù)M的最優(yōu)組合。根據(jù)niedman和Stu-etzle提出的投影尋蹤技術(shù)多重平滑實(shí)現(xiàn)法,核心是采用多層分組迭代交替優(yōu)化方法,選取m組參數(shù)以及嶺函數(shù)的最優(yōu)項(xiàng)數(shù)目。把參數(shù)分為若干組,選擇其中一組作為參數(shù)尋優(yōu)的對(duì)象,其他組都給定一初值。求得結(jié)果后,把尋優(yōu)所得的參數(shù)極值點(diǎn)作初值,另選一組參數(shù)在這一初值下尋優(yōu),多次反復(fù)尋優(yōu)直到最后的一組參數(shù)使得目標(biāo)函數(shù)達(dá)到最小。

    傳統(tǒng)的投影尋蹤方法尋優(yōu)速度較慢,且尋優(yōu)結(jié)果與初始值的設(shè)定有關(guān)。為了更好地尋找模型的參數(shù),文獻(xiàn)分別采用不同的尋優(yōu)算法與投影尋蹤算法相結(jié)合,確定最優(yōu)的參數(shù)組合,以此來(lái)進(jìn)行預(yù)測(cè)。本文提出了通過(guò)延遲相關(guān)系數(shù)法確定回歸預(yù)測(cè)因子、群居蜘蛛算法f Social Spider Optimization,SSO)優(yōu)化投影尋蹤(Projeetion pursuit regres-sion,PPR)模型最佳投影方向參數(shù)a,最小二乘法確定多項(xiàng)式權(quán)系數(shù)c,合格率控制參數(shù)M個(gè)數(shù)相結(jié)合的混合智能投影尋蹤年最大洪峰流量預(yù)測(cè)模型,結(jié)合長(zhǎng)江宜昌站(1882年-2004年)的年最大洪峰流量資料進(jìn)行了實(shí)例分析。

    1基于參數(shù)投影尋蹤回歸的年最大洪峰流量建模原理

    投影尋蹤回歸法是Friedman和Stuetzle思想是尋找能反映高維數(shù)據(jù)的結(jié)構(gòu)特征的若干個(gè)投影方向,將高維數(shù)據(jù)投影到一維空間,通過(guò)優(yōu)選出的嶺函數(shù)的代數(shù)和去逼近回歸函數(shù),以此構(gòu)造模型。其回歸預(yù)測(cè)建模的原理主要如下。

    設(shè)x為p維的預(yù)測(cè)因子,a為p維投影向量,f(x)為預(yù)測(cè)的年最大洪峰流量,為了解決多元線性回歸不能有效反映實(shí)際非線性的問(wèn)題,投影尋蹤回歸模型用一系列的嶺函數(shù)g(aTx)的代數(shù)和去逼近回歸函數(shù)f(x),因此,預(yù)測(cè)量與預(yù)測(cè)因子的投影回歸關(guān)系能被表達(dá)為

    (1)式中:M為嶺函數(shù)的個(gè)數(shù);g為嶺函數(shù)。當(dāng)M,g等于1時(shí),式(l)就變?yōu)槌R?guī)的多元函數(shù)。

    在Friedman和Stuetzle最初的研究中提出的采用龐大的簡(jiǎn)單函數(shù)集合,通過(guò)分層組迭代交替優(yōu)化技術(shù)去光滑逼近,實(shí)現(xiàn)難度較大,在很大程度上限制了該模型的推廣與應(yīng)用。Hwang等人研究表明Hermite多項(xiàng)式具有較強(qiáng)的光滑、趨勢(shì)、動(dòng)態(tài)、非線性擬合能力,且多項(xiàng)式次數(shù)的增加會(huì)增強(qiáng)起擬合復(fù)雜程度的能力,因此,可采用可變階遞推的正交Hermite多項(xiàng)式擬合一維嶺函數(shù)。則投影回歸關(guān)系式就變?yōu)?/p>

    (2)式中:n表示輸入樣本的個(gè)數(shù);zi表示第i個(gè)輸入樣本在投影方向。上的投影,其表示達(dá)式為

    (3)

    j為多項(xiàng)式的階數(shù);c為多項(xiàng)式的系數(shù);h為正交Hermite多項(xiàng)式,h的數(shù)學(xué)表達(dá)式如下:

    (4)式中;j!代表多項(xiàng)式階數(shù)i的階乘;φ(z)=

    (5)

    由此回歸預(yù)測(cè)建模就可轉(zhuǎn)化為求解式(6)的最小化問(wèn)題。

    (6)

    (7)

    (8)

    根據(jù)水文情報(bào)預(yù)報(bào)規(guī)范(GB/T 22482-2008),合格預(yù)報(bào)次數(shù)與預(yù)報(bào)總次數(shù)之比的百分?jǐn)?shù)為合格率QR,其計(jì)算公式為:

    (9)式中QR為合格率(%);nh為合格預(yù)報(bào)的次數(shù):n為預(yù)報(bào)的總次數(shù)。

    所以,投影尋蹤回歸預(yù)測(cè)的關(guān)鍵就轉(zhuǎn)化為優(yōu)選投影回歸模型的投影方向參數(shù)a、多項(xiàng)式權(quán)系數(shù)c和嶺函數(shù)的個(gè)數(shù)M,參數(shù)優(yōu)選的結(jié)果直接影響模型的泛化能力、模擬及預(yù)測(cè)精度。

    2群居蜘蛛優(yōu)化算法

    群居蜘蛛優(yōu)化算法(SSO)是由Cuevas等人基于對(duì)群居蜘蛛的協(xié)作行為,提出的一種新的隨機(jī)全局優(yōu)化技術(shù)。群居型蜘蛛是群居性昆蟲(chóng)的典型例子,是群體成員能夠完成一系列復(fù)雜的協(xié)同行為的蜘蛛。在群居蜘蛛群落中,每一個(gè)個(gè)體,由于性別的差異,完成一系列任務(wù)如:捕食、交配、織網(wǎng)和相互交流的任務(wù)。蜘蛛的網(wǎng)是群落重要的一部分,它不僅是全體成員生存環(huán)境的一部分,而且是它們之間的交流渠道。因此重要的信息(例如捕食或者是交配)由網(wǎng)通過(guò)震動(dòng)的形式傳播。這些信息,被當(dāng)作局部的知識(shí),每個(gè)成員利用這些知識(shí)來(lái)指導(dǎo)它們各自的協(xié)作行為,同時(shí)影響了群落的社會(huì)規(guī)范。

    群居蜘蛛生物群落有兩個(gè)基本的組成部分:群落成員和公共的網(wǎng)。群落成員被分為兩類(lèi):雌蜘蛛、雄蜘蛛。群居蜘蛛是高度的母系氏種群。一些研究表明雄蜘蛛的數(shù)量差不多達(dá)到全部群落成員數(shù)目的30%。群落成員之間有直接和間接的交流。直接的交流意味著身體接觸或者是體液的交換諸如交配。對(duì)于間接的交流,公共的網(wǎng)作為一個(gè)交流媒介用來(lái)傳達(dá)重要的信息,這些信息每一個(gè)群落成員都能夠感受得到。這些被編碼成震動(dòng)的信息是成員之間集體協(xié)作的一個(gè)關(guān)鍵的方面。群落成員將震動(dòng)解碼成幾條信息諸如被捕獲的動(dòng)物的大小,鄰居成員的特征等等。這種震動(dòng)的強(qiáng)度依賴(lài)于激發(fā)他們的蜘蛛的重量和距離。

    它們相互作用的方式是依賴(lài)于蜘蛛的性別的。雌蜘蛛展示了群居式的主要傾向,由于震動(dòng)依賴(lài)于激發(fā)他們的蜘蛛的重量和距離,很強(qiáng)有力的震動(dòng)或者是大蜘蛛激發(fā)的,或者是距離較近的蜘蛛發(fā)出的。與雌蜘蛛不同的是,雄蜘蛛的行為是生殖導(dǎo)向的。雄性蜘蛛將他們自己視為雄性組織的領(lǐng)袖,他們控制著群落的資源。因此,雄性個(gè)體被分為兩類(lèi):優(yōu)勢(shì)者和劣勢(shì)者。優(yōu)勢(shì)的雄性蜘蛛和劣勢(shì)的雄性蜘蛛比起來(lái)有更好的適應(yīng)度。在一個(gè)典型的行為中,優(yōu)勢(shì)蜘蛛被公共網(wǎng)上的最近的雌性蜘蛛所吸引。相反,劣勢(shì)雄性蜘蛛圍繞在優(yōu)勢(shì)雄性群體的周?chē)脙?yōu)勢(shì)雄性種群浪費(fèi)的資源。

    交配是一個(gè)重要的生理行為,該行為不僅保證了種群的延續(xù),而且能夠使信息在種群之間交流。群居蜘蛛的交配發(fā)生在優(yōu)勢(shì)雄性蜘蛛和雌性蜘蛛之間。在這種情況下,當(dāng)一個(gè)優(yōu)勢(shì)蜘蛛在一只或是更多的雌性蜘蛛的特定的范圍內(nèi),它們和所有的蜘蛛交配為了產(chǎn)生后代。

    作為優(yōu)化技術(shù)其求解優(yōu)化問(wèn)題最小值的步驟可參考文獻(xiàn),在此不再贅述。

    3年最大洪峰流量混合智能預(yù)測(cè)建模步驟

    根據(jù)前面介紹投影尋蹤回歸建模的原理和群居蜘蛛優(yōu)化算法求解最小值優(yōu)化問(wèn)題的步驟,年最大洪峰流量混合智能預(yù)測(cè)建模的具體步驟如下。

    (1)數(shù)據(jù)歸一化。為消除數(shù)據(jù)在量綱和標(biāo)準(zhǔn)差數(shù)值水平上的差異,利用下式對(duì)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理:

    (10)式中:qi表示年最大洪峰流量序列;xi表示歸一化后的流量序列;n表示序列的容量。

    (2)確定預(yù)測(cè)因子。利用相關(guān)技術(shù)確定徑流序列的預(yù)測(cè)因子。時(shí)序xi延遲k步的自相關(guān)系數(shù)Rk可通過(guò)下式獲得:

    (11)

    (12)

    (13)則推斷時(shí)序xi延遲k步相依性顯著,將xi-k作為xi的預(yù)測(cè)因子。μα12從正態(tài)分布表中查得,在本文研究中取80%的置信水平。

    (3)生成初始投影方向。設(shè)定蜘蛛的種群規(guī)模Ⅳ,按約束條件式(7)生成N組PP模型投影方向a的初始值。

    (4)計(jì)算多項(xiàng)式權(quán)系數(shù)c。根據(jù)生成的投影方向,利用式(3)計(jì)算投影值z(mì),利用式(4)計(jì)算r階Hermite多項(xiàng)式hr(z),然后最小二乘法計(jì)算多項(xiàng)式權(quán)系數(shù)c。

    (5)計(jì)算相應(yīng)的適應(yīng)度值,評(píng)價(jià)蜘蛛個(gè)體。根據(jù)第4步得到權(quán)系數(shù)c,利用式(2)計(jì)算回歸值,利用式(6)計(jì)算適應(yīng)度值,進(jìn)入群居蜘蛛優(yōu)化算法的步驟2。經(jīng)過(guò)群居蜘蛛優(yōu)化算法一系列的機(jī)制,更新個(gè)體,進(jìn)行迭代計(jì)算,獲得最優(yōu)的投影方向參數(shù)。和多項(xiàng)式權(quán)系數(shù)c,第一個(gè)嶺函數(shù)優(yōu)化結(jié)束。

    (6)模型優(yōu)化終止和結(jié)果輸出。根據(jù)優(yōu)選得到的投影方向參數(shù)a和多項(xiàng)式權(quán)系數(shù)c,計(jì)算擬合殘差和合格率,如果滿(mǎn)足合格率終止準(zhǔn)則,輸出預(yù)測(cè)值和相應(yīng)的參數(shù),用εi代替yi,轉(zhuǎn)入(3)進(jìn)行下一個(gè)嶺函數(shù)的優(yōu)化。為了避免程序進(jìn)入死循環(huán),模型優(yōu)化終止準(zhǔn)則采用合格率和最大嶺函數(shù)個(gè)數(shù)相結(jié)合的方式。

    4實(shí)例驗(yàn)證

    選取長(zhǎng)江宜昌站1882年-2004年共123年的年最大洪峰流量序列{q1|i=1,2.…,123},應(yīng)用前113年(1882年-1994年)訓(xùn)練投影尋蹤混合智能預(yù)測(cè)模型確定模型參數(shù),用后10年(1995年-2004年)的年最大洪峰流量進(jìn)行檢驗(yàn)。通過(guò)計(jì)算計(jì)算該序列前30階自相關(guān)系數(shù)Rk和與之相應(yīng)的上、下限R1k,R2k值,其中置信水平取80%,計(jì)算結(jié)果表明R1,R2,R3,R20,R21,R22,R25,R28,R30的相依性在置信水平80%的條件下是顯著的。所以,對(duì)qi預(yù)測(cè)的因子取為qi-1,qi-2,qi-3,qi-20,qi-21,qi-22,qq-25,qqi-28,qqi-30,則前113年(1882年-1994年)可得到83組訓(xùn)練數(shù)據(jù)。訓(xùn)練過(guò)程,a的取值范圍為[-l,1],Hermite多項(xiàng)式階數(shù)為r=6,群居蜘蛛優(yōu)化算法的種群規(guī)模N=50,最大迭代次數(shù)Gmax為500;嶺函數(shù)的最大個(gè)數(shù)Mmax=3。經(jīng)群居蜘蛛算法和最小二乘方法混合優(yōu)化,得到長(zhǎng)江宜昌站年最大洪峰流量的預(yù)測(cè)模型為

    (14)

    根據(jù)得到的混合智能投影尋蹤年最大洪峰流量預(yù)測(cè)模型,計(jì)算1995年-2004年的年最大洪峰流量預(yù)測(cè)值。圖1給出了本文提出模型的(1912年-1994年)訓(xùn)練模擬值和(1995年-2004年)預(yù)測(cè)值與宜昌站年最大洪峰流量觀測(cè)值的對(duì)比圖。訓(xùn)練階段和預(yù)測(cè)階段相關(guān)誤差統(tǒng)計(jì)分析結(jié)果見(jiàn)表1。根據(jù)水文情報(bào)預(yù)報(bào)規(guī)范(GB/T 22482-2008),合格率超過(guò)85%的為預(yù)報(bào)等級(jí)為甲等。根據(jù)表1中統(tǒng)計(jì)分析結(jié)果可以看出,在訓(xùn)練階段的合格率為94.0%,在檢驗(yàn)階段的合格率為90%,預(yù)報(bào)等級(jí)為甲等。同時(shí)也可以看出,在訓(xùn)練階段相對(duì)誤差絕對(duì)值在[0,15]區(qū)間內(nèi)的百分比達(dá)84 30%,在[0,10]區(qū)間內(nèi)的百分比達(dá)67.5%:在檢驗(yàn)階段相對(duì)誤差絕對(duì)值在[0,15]區(qū)間內(nèi)的百分比達(dá)70.0%,在[0,10]區(qū)間內(nèi)的百分比達(dá)60%:訓(xùn)練階段的平均絕對(duì)相對(duì)誤差僅有8.61;檢驗(yàn)階段的最大相對(duì)誤差絕對(duì)值為22.44%,最小相對(duì)誤差絕對(duì)值為1.76%,平均絕對(duì)相對(duì)誤差僅有10.51。所以提出混合智能投影尋蹤年最大洪峰流量預(yù)測(cè)模型宜昌站的年最大洪峰流量的模擬預(yù)測(cè)獲得了非常好的效果。

    為了測(cè)試預(yù)報(bào)結(jié)果的穩(wěn)定性,在相同的計(jì)算機(jī)系統(tǒng)和環(huán)境下,不改變種群規(guī)模和迭代次數(shù)條件下,運(yùn)行10次,優(yōu)選出的模型參數(shù)沒(méi)有變化,說(shuō)明提出的混合智能投影尋蹤預(yù)測(cè)模型穩(wěn)定性強(qiáng)。

    5結(jié)論

    對(duì)于年最大洪峰流量預(yù)測(cè),傳統(tǒng)的統(tǒng)計(jì)分析方法很難取得滿(mǎn)意的效果,這從對(duì)其預(yù)測(cè)的文獻(xiàn)檢索也可以印證,其研究成果明顯少于其他中長(zhǎng)期水文預(yù)報(bào)(如年、月徑流量)的預(yù)報(bào)。為此,本文嘗試性地將參數(shù)投影尋蹤模型引入年最大洪峰流量預(yù)報(bào),提出了運(yùn)用延遲相關(guān)系數(shù)法確定回歸預(yù)測(cè)因子、群居蜘蛛算法優(yōu)化投影尋蹤模型參數(shù)的混合智能投影尋蹤年最大洪峰流量預(yù)測(cè)模型。以長(zhǎng)江宜昌站(1882年-2004年)的年最大洪峰流量資料為依據(jù)進(jìn)行了實(shí)例分析,結(jié)果表明,建立模型能獲得很好的精度,檢驗(yàn)階段預(yù)報(bào)合格率達(dá)到90%,明顯高于文獻(xiàn)中用諧波預(yù)報(bào)法所得的檢驗(yàn)階段合格率68.2%以及基于小波變換序列的最近鄰抽樣回歸組合預(yù)報(bào)所得的檢驗(yàn)階段合格率81.8%。同時(shí)提出模型智能化程度高,預(yù)測(cè)因子,模型參數(shù)可以完全有計(jì)算機(jī)完成,且結(jié)果穩(wěn)定性好,可有效應(yīng)用于年最大洪峰流量預(yù)測(cè),值得在實(shí)際中推廣應(yīng)用。

    猜你喜歡
    預(yù)測(cè)
    無(wú)可預(yù)測(cè)
    黃河之聲(2022年10期)2022-09-27 13:59:46
    基于PCC-CNN-GRU的短期風(fēng)電功率預(yù)測(cè)
    選修2—2期中考試預(yù)測(cè)卷(A卷)答案與提示
    選修2-2期中考試預(yù)測(cè)卷(A卷)
    選修2-2期中考試預(yù)測(cè)卷(B卷)
    選修2—2期中考試預(yù)測(cè)卷(B卷)
    選修2—2期中考試預(yù)測(cè)卷(A卷)
    選修2—2期中考試預(yù)測(cè)卷(B卷)答案與提示
    “預(yù)測(cè)”得準(zhǔn)
    不可預(yù)測(cè)
    男人舔女人的私密视频| av在线蜜桃| 国产视频一区二区在线看| 身体一侧抽搐| www.自偷自拍.com| 国产爱豆传媒在线观看| 女人高潮潮喷娇喘18禁视频| 久久天堂一区二区三区四区| 岛国视频午夜一区免费看| 国产精品久久久人人做人人爽| 久久精品夜夜夜夜夜久久蜜豆| 国产蜜桃级精品一区二区三区| h日本视频在线播放| 网址你懂的国产日韩在线| 日韩欧美精品v在线| 一个人免费在线观看的高清视频| 国产精品美女特级片免费视频播放器 | 一本精品99久久精品77| 变态另类成人亚洲欧美熟女| 欧美午夜高清在线| 十八禁人妻一区二区| 亚洲精品乱码久久久v下载方式 | 久久亚洲精品不卡| cao死你这个sao货| 欧美国产日韩亚洲一区| 欧美日韩精品网址| 久久精品国产清高在天天线| 在线观看午夜福利视频| www.精华液| 欧美日韩瑟瑟在线播放| 精品国产三级普通话版| 亚洲aⅴ乱码一区二区在线播放| 性欧美人与动物交配| netflix在线观看网站| 国内少妇人妻偷人精品xxx网站 | 成人国产综合亚洲| 久久久成人免费电影| 欧美在线一区亚洲| 美女扒开内裤让男人捅视频| 国产精品一区二区精品视频观看| 宅男免费午夜| 好男人在线观看高清免费视频| 国产又黄又爽又无遮挡在线| 女同久久另类99精品国产91| 欧美黄色片欧美黄色片| 男人和女人高潮做爰伦理| 国产 一区 欧美 日韩| 免费av不卡在线播放| 成人三级黄色视频| 亚洲av熟女| 久久精品国产99精品国产亚洲性色| 黑人操中国人逼视频| 欧美+亚洲+日韩+国产| 国产精品日韩av在线免费观看| 别揉我奶头~嗯~啊~动态视频| 精品久久久久久久久久免费视频| 国产精品综合久久久久久久免费| 露出奶头的视频| 女警被强在线播放| 国产精品亚洲美女久久久| 国产av在哪里看| 成在线人永久免费视频| 熟女少妇亚洲综合色aaa.| 欧美精品啪啪一区二区三区| 色综合亚洲欧美另类图片| 欧美黑人巨大hd| 国产激情久久老熟女| 激情在线观看视频在线高清| 午夜免费观看网址| 亚洲成人久久爱视频| 国产欧美日韩精品一区二区| 嫩草影视91久久| 露出奶头的视频| 夜夜夜夜夜久久久久| 一本综合久久免费| 国产主播在线观看一区二区| 国产亚洲精品久久久com| 久久国产精品人妻蜜桃| 一进一出抽搐gif免费好疼| 99在线人妻在线中文字幕| 久久久国产成人免费| 99热精品在线国产| 嫩草影视91久久| 国产精品自产拍在线观看55亚洲| 国产69精品久久久久777片 | 麻豆成人午夜福利视频| 欧美乱码精品一区二区三区| 97人妻精品一区二区三区麻豆| 搡老妇女老女人老熟妇| 午夜成年电影在线免费观看| 身体一侧抽搐| 久久中文字幕人妻熟女| 法律面前人人平等表现在哪些方面| 国产精品1区2区在线观看.| 香蕉av资源在线| 亚洲五月婷婷丁香| 中文字幕高清在线视频| 亚洲国产精品sss在线观看| 18禁美女被吸乳视频| 欧美日韩亚洲国产一区二区在线观看| 精品一区二区三区视频在线观看免费| 可以在线观看毛片的网站| 12—13女人毛片做爰片一| 久久久久久久午夜电影| 婷婷丁香在线五月| 久久精品91蜜桃| 亚洲精品美女久久久久99蜜臀| 真人一进一出gif抽搐免费| 午夜精品一区二区三区免费看| 国产一区二区三区视频了| 变态另类丝袜制服| 免费在线观看成人毛片| 两性午夜刺激爽爽歪歪视频在线观看| 一进一出好大好爽视频| 无人区码免费观看不卡| 国产伦在线观看视频一区| 18禁黄网站禁片免费观看直播| 久久久久久久精品吃奶| 啦啦啦韩国在线观看视频| www.自偷自拍.com| 中亚洲国语对白在线视频| 99热这里只有是精品50| 午夜视频精品福利| 午夜亚洲福利在线播放| 精品一区二区三区av网在线观看| 女人高潮潮喷娇喘18禁视频| 麻豆av在线久日| 亚洲在线观看片| 亚洲专区中文字幕在线| 成人午夜高清在线视频| 伦理电影免费视频| 国内少妇人妻偷人精品xxx网站 | 国产精品98久久久久久宅男小说| 女生性感内裤真人,穿戴方法视频| 亚洲欧美激情综合另类| 三级男女做爰猛烈吃奶摸视频| 在线免费观看的www视频| 精品熟女少妇八av免费久了| 国产伦一二天堂av在线观看| 日韩免费av在线播放| 亚洲美女视频黄频| 中文字幕人妻丝袜一区二区| 免费大片18禁| 99在线视频只有这里精品首页| 熟女少妇亚洲综合色aaa.| 99精品在免费线老司机午夜| 成人国产综合亚洲| 中文亚洲av片在线观看爽| 欧美+亚洲+日韩+国产| 国产私拍福利视频在线观看| 精品免费久久久久久久清纯| 国产成人aa在线观看| 久久久久久久久免费视频了| 久久久久国产精品人妻aⅴ院| 国产精品久久久久久亚洲av鲁大| 中文字幕人妻丝袜一区二区| 最近最新免费中文字幕在线| 少妇的丰满在线观看| 成人精品一区二区免费| 99久久久亚洲精品蜜臀av| 日韩欧美三级三区| 性色avwww在线观看| 欧洲精品卡2卡3卡4卡5卡区| 九九在线视频观看精品| 黄色视频,在线免费观看| 色av中文字幕| 日本免费一区二区三区高清不卡| 免费人成视频x8x8入口观看| 99国产精品99久久久久| 美女免费视频网站| 美女cb高潮喷水在线观看 | 黄色视频,在线免费观看| 国产精品久久电影中文字幕| 18禁观看日本| 欧美xxxx黑人xx丫x性爽| 99久久精品一区二区三区| 精品久久久久久,| 中亚洲国语对白在线视频| 99热精品在线国产| 国产伦精品一区二区三区视频9 | 婷婷丁香在线五月| 国内少妇人妻偷人精品xxx网站 | 嫩草影视91久久| 亚洲av熟女| 亚洲欧美日韩高清在线视频| 午夜a级毛片| 好男人在线观看高清免费视频| 国产午夜精品论理片| 一进一出好大好爽视频| 女人被狂操c到高潮| 日本成人三级电影网站| 成人特级av手机在线观看| 欧美另类亚洲清纯唯美| 国产乱人视频| 国产精品久久久久久人妻精品电影| 99热这里只有是精品50| 久久精品夜夜夜夜夜久久蜜豆| 欧美色视频一区免费| netflix在线观看网站| 日韩欧美一区二区三区在线观看| 99热6这里只有精品| 一二三四在线观看免费中文在| 19禁男女啪啪无遮挡网站| 亚洲av电影不卡..在线观看| av天堂在线播放| 国产精品久久视频播放| 特级一级黄色大片| 国产成人福利小说| 高清在线国产一区| 亚洲国产精品sss在线观看| 美女cb高潮喷水在线观看 | 99riav亚洲国产免费| 久久久久久久久中文| 国产精品av视频在线免费观看| 别揉我奶头~嗯~啊~动态视频| 国产精品乱码一区二三区的特点| 久久精品综合一区二区三区| 最新美女视频免费是黄的| 精品人妻1区二区| 99riav亚洲国产免费| 香蕉久久夜色| 久久久久久久午夜电影| 亚洲在线观看片| 美女免费视频网站| 操出白浆在线播放| 精品国产亚洲在线| 一夜夜www| 99热精品在线国产| 精品乱码久久久久久99久播| 欧洲精品卡2卡3卡4卡5卡区| 国产野战对白在线观看| tocl精华| 日本成人三级电影网站| 久久久久久久精品吃奶| 搡老熟女国产l中国老女人| 亚洲av美国av| 欧美国产日韩亚洲一区| 人人妻人人澡欧美一区二区| 国产真实乱freesex| 岛国在线观看网站| 偷拍熟女少妇极品色| xxx96com| 午夜福利欧美成人| 成人鲁丝片一二三区免费| 全区人妻精品视频| 日本免费a在线| 怎么达到女性高潮| 午夜福利免费观看在线| 国产高清三级在线| 久久久成人免费电影| 日本黄色片子视频| 国产乱人视频| 国产淫片久久久久久久久 | 成人18禁在线播放| 巨乳人妻的诱惑在线观看| 国产淫片久久久久久久久 | 伦理电影免费视频| 99re在线观看精品视频| xxxwww97欧美| 国产精品久久久久久亚洲av鲁大| 精品熟女少妇八av免费久了| 亚洲一区二区三区色噜噜| 国产成人精品无人区| 中文字幕av在线有码专区| 国产蜜桃级精品一区二区三区| 日本黄色片子视频| 精品人妻1区二区| 麻豆国产av国片精品| 黄色日韩在线| 欧美zozozo另类| 波多野结衣巨乳人妻| 亚洲激情在线av| 老司机午夜福利在线观看视频| 午夜福利高清视频| 精品欧美国产一区二区三| 国产精品国产高清国产av| 淫秽高清视频在线观看| 亚洲人成网站高清观看| 18禁国产床啪视频网站| 亚洲av电影不卡..在线观看| 18禁黄网站禁片免费观看直播| 99久久国产精品久久久| 亚洲中文字幕日韩| 嫩草影院入口| 熟妇人妻久久中文字幕3abv| or卡值多少钱| 色综合婷婷激情| 美女高潮喷水抽搐中文字幕| 日韩欧美国产在线观看| www日本黄色视频网| 国产成人一区二区三区免费视频网站| 亚洲 欧美 日韩 在线 免费| 久久久久九九精品影院| 免费无遮挡裸体视频| 国产精品av久久久久免费| 精品国内亚洲2022精品成人| 老汉色av国产亚洲站长工具| 国产综合懂色| 国产精品久久久久久精品电影| 女人被狂操c到高潮| 又粗又爽又猛毛片免费看| 午夜精品在线福利| 欧美日韩黄片免| 色精品久久人妻99蜜桃| 网址你懂的国产日韩在线| 观看美女的网站| 9191精品国产免费久久| 国产黄a三级三级三级人| xxx96com| 国产成年人精品一区二区| 午夜成年电影在线免费观看| 热99在线观看视频| 国产精品 欧美亚洲| 天天躁狠狠躁夜夜躁狠狠躁| 久久精品国产99精品国产亚洲性色| 欧美乱妇无乱码| 9191精品国产免费久久| 嫩草影院精品99| 日日摸夜夜添夜夜添小说| 欧美黑人欧美精品刺激| 日本一本二区三区精品| 黄片小视频在线播放| 日本三级黄在线观看| 美女免费视频网站| 精品日产1卡2卡| 亚洲成人精品中文字幕电影| 女人被狂操c到高潮| 午夜福利在线在线| 国产日本99.免费观看| 久久久久久人人人人人| 成年人黄色毛片网站| 三级男女做爰猛烈吃奶摸视频| 69av精品久久久久久| 国产伦在线观看视频一区| 亚洲av电影在线进入| 日日夜夜操网爽| 十八禁网站免费在线| 天堂影院成人在线观看| 97超视频在线观看视频| 国产 一区 欧美 日韩| 日韩欧美一区二区三区在线观看| 级片在线观看| 岛国视频午夜一区免费看| 天堂影院成人在线观看| 97超视频在线观看视频| 亚洲欧美日韩高清在线视频| 九色成人免费人妻av| 亚洲av成人精品一区久久| 精品国产美女av久久久久小说| www.精华液| 最新中文字幕久久久久 | 国产精品影院久久| 99久久精品国产亚洲精品| 久久欧美精品欧美久久欧美| xxxwww97欧美| 国产成人系列免费观看| 亚洲av电影不卡..在线观看| 国产伦精品一区二区三区视频9 | 麻豆一二三区av精品| 长腿黑丝高跟| 一级毛片精品| 九九热线精品视视频播放| 波多野结衣巨乳人妻| 国产亚洲欧美在线一区二区| 99国产精品99久久久久| 99国产综合亚洲精品| 亚洲av片天天在线观看| 日韩欧美免费精品| 日韩精品中文字幕看吧| 午夜福利在线观看吧| 老熟妇乱子伦视频在线观看| 嫩草影院精品99| 一本精品99久久精品77| 精品乱码久久久久久99久播| 亚洲第一电影网av| 欧美av亚洲av综合av国产av| 中文字幕精品亚洲无线码一区| 日本免费一区二区三区高清不卡| 国内精品久久久久久久电影| 国产野战对白在线观看| 特大巨黑吊av在线直播| avwww免费| www.精华液| 黄色女人牲交| 最近最新中文字幕大全电影3| 美女黄网站色视频| 免费观看精品视频网站| 视频区欧美日本亚洲| 日日夜夜操网爽| 久久久久久国产a免费观看| 成人精品一区二区免费| 日韩欧美免费精品| 曰老女人黄片| 免费看十八禁软件| 亚洲电影在线观看av| 成年人黄色毛片网站| 黄色 视频免费看| 日本一二三区视频观看| 免费在线观看影片大全网站| 无人区码免费观看不卡| 天天躁狠狠躁夜夜躁狠狠躁| 法律面前人人平等表现在哪些方面| 久久99热这里只有精品18| 欧美日韩亚洲国产一区二区在线观看| 免费av毛片视频| 日韩欧美国产在线观看| 国产av不卡久久| 99国产精品99久久久久| 我要搜黄色片| 日韩国内少妇激情av| 免费人成视频x8x8入口观看| 亚洲国产高清在线一区二区三| 欧美乱妇无乱码| 欧美一区二区精品小视频在线| 午夜激情欧美在线| 91av网一区二区| 99热这里只有精品一区 | 99久国产av精品| 国产一区二区在线观看日韩 | 一级毛片精品| 免费大片18禁| 国产精品一区二区精品视频观看| 亚洲成av人片免费观看| 精品电影一区二区在线| 久久久国产成人精品二区| 99re在线观看精品视频| 国产又黄又爽又无遮挡在线| 国产不卡一卡二| 国产激情久久老熟女| 国产精品一及| 丁香六月欧美| 偷拍熟女少妇极品色| 亚洲国产精品合色在线| 黄片大片在线免费观看| 久久久久久九九精品二区国产| 女警被强在线播放| 成人三级黄色视频| av黄色大香蕉| 国产精品九九99| 欧美日韩国产亚洲二区| 欧美黑人巨大hd| 欧美三级亚洲精品| 国产麻豆成人av免费视频| 亚洲乱码一区二区免费版| 老司机午夜福利在线观看视频| 草草在线视频免费看| 制服人妻中文乱码| 亚洲精品国产精品久久久不卡| 9191精品国产免费久久| 啪啪无遮挡十八禁网站| 欧美三级亚洲精品| 很黄的视频免费| 国产真人三级小视频在线观看| 亚洲欧洲精品一区二区精品久久久| 伊人久久大香线蕉亚洲五| 热99re8久久精品国产| 国产野战对白在线观看| 久久久久国产一级毛片高清牌| 久久亚洲精品不卡| 啦啦啦观看免费观看视频高清| 90打野战视频偷拍视频| 成人国产综合亚洲| 青草久久国产| 波多野结衣高清无吗| 亚洲第一欧美日韩一区二区三区| 免费电影在线观看免费观看| 中文字幕av在线有码专区| 亚洲av第一区精品v没综合| 日本成人三级电影网站| 一进一出抽搐动态| 日韩 欧美 亚洲 中文字幕| 国产aⅴ精品一区二区三区波| 亚洲精品粉嫩美女一区| 99国产精品一区二区三区| 欧美黑人巨大hd| av在线蜜桃| 两人在一起打扑克的视频| 午夜a级毛片| 91老司机精品| 亚洲性夜色夜夜综合| 亚洲色图av天堂| 亚洲色图 男人天堂 中文字幕| 嫩草影视91久久| 亚洲av电影不卡..在线观看| www.999成人在线观看| 亚洲av片天天在线观看| 母亲3免费完整高清在线观看| 国内精品美女久久久久久| 国产视频内射| 国产欧美日韩一区二区三| 99热精品在线国产| 男女之事视频高清在线观看| 国产精品久久久久久久电影 | 国产亚洲精品综合一区在线观看| 免费搜索国产男女视频| 久久久久久人人人人人| 欧美成人免费av一区二区三区| 久久久久久久久久黄片| 91在线观看av| 一边摸一边抽搐一进一小说| 又粗又爽又猛毛片免费看| 成人亚洲精品av一区二区| 国产精品,欧美在线| 看免费av毛片| 精品国产乱码久久久久久男人| 国产毛片a区久久久久| 禁无遮挡网站| 亚洲真实伦在线观看| 欧洲精品卡2卡3卡4卡5卡区| 女生性感内裤真人,穿戴方法视频| 日本熟妇午夜| 给我免费播放毛片高清在线观看| 色在线成人网| xxx96com| 观看美女的网站| 夜夜看夜夜爽夜夜摸| x7x7x7水蜜桃| 九九久久精品国产亚洲av麻豆 | 99久久精品热视频| 精品久久久久久久人妻蜜臀av| 99国产精品一区二区蜜桃av| 12—13女人毛片做爰片一| 亚洲第一欧美日韩一区二区三区| 国产高清视频在线观看网站| 国产精品亚洲美女久久久| 嫁个100分男人电影在线观看| 99国产极品粉嫩在线观看| 亚洲第一电影网av| 中文字幕熟女人妻在线| 国产精品1区2区在线观看.| 国产美女午夜福利| 久久这里只有精品中国| 国产av麻豆久久久久久久| cao死你这个sao货| 日韩免费av在线播放| 亚洲avbb在线观看| 成年女人永久免费观看视频| 网址你懂的国产日韩在线| 一进一出抽搐动态| 99久久精品热视频| 9191精品国产免费久久| 日韩高清综合在线| 亚洲美女黄片视频| 国产成人av激情在线播放| 搡老妇女老女人老熟妇| 两个人的视频大全免费| 美女扒开内裤让男人捅视频| 欧美一区二区国产精品久久精品| 一二三四在线观看免费中文在| 国产成人av教育| 国产高清激情床上av| 小蜜桃在线观看免费完整版高清| 色精品久久人妻99蜜桃| av天堂在线播放| 国产 一区 欧美 日韩| 在线观看美女被高潮喷水网站 | 色尼玛亚洲综合影院| 亚洲无线在线观看| 久久国产精品影院| 日本黄色片子视频| 十八禁网站免费在线| 色老头精品视频在线观看| 日韩有码中文字幕| 久久久成人免费电影| 91九色精品人成在线观看| 日本一二三区视频观看| 搡老岳熟女国产| 亚洲 欧美一区二区三区| 免费人成视频x8x8入口观看| 国产三级中文精品| 国产aⅴ精品一区二区三区波| 91av网站免费观看| 两个人的视频大全免费| 欧美中文综合在线视频| 噜噜噜噜噜久久久久久91| 精品不卡国产一区二区三区| 欧美色欧美亚洲另类二区| 深夜精品福利| 在线十欧美十亚洲十日本专区| 欧美日韩一级在线毛片| 日韩欧美一区二区三区在线观看| 搡老熟女国产l中国老女人| 成年人黄色毛片网站| 色噜噜av男人的天堂激情| 看黄色毛片网站| 午夜福利在线观看吧| 一区二区三区高清视频在线| 国产激情久久老熟女| 最近视频中文字幕2019在线8| 精品一区二区三区视频在线观看免费| 成人三级做爰电影| a级毛片a级免费在线| 久久精品91蜜桃| 99久国产av精品| 欧美性猛交黑人性爽| 一区二区三区激情视频| 老司机福利观看| 欧美性猛交黑人性爽| 一区二区三区激情视频| 人妻久久中文字幕网| 久久久久久久午夜电影| 全区人妻精品视频| 色av中文字幕| 中文字幕熟女人妻在线| 999久久久精品免费观看国产| 亚洲avbb在线观看| 成在线人永久免费视频| 麻豆久久精品国产亚洲av| 国产激情久久老熟女| 久久中文字幕人妻熟女| 国产高潮美女av| 色尼玛亚洲综合影院| 午夜亚洲福利在线播放| 伊人久久大香线蕉亚洲五| 国产在线精品亚洲第一网站|