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

    基于M-H采樣的快速反向微分進化算法

    2014-06-07 05:53:22涂維維葛洪偉楊金龍袁運浩
    計算機工程 2014年11期
    關鍵詞:優(yōu)化

    涂維維,葛洪偉,楊金龍,袁運浩

    (江南大學物聯(lián)網(wǎng)工程學院,江蘇無錫214211)

    基于M-H采樣的快速反向微分進化算法

    涂維維,葛洪偉,楊金龍,袁運浩

    (江南大學物聯(lián)網(wǎng)工程學院,江蘇無錫214211)

    反向微分進化(ODE)算法基于反向優(yōu)化對種群進行初始化更新以保持種群多樣性。但該算法中反向個體容易偏離全局最優(yōu)個體,不能很快達到全局最優(yōu),在函數(shù)優(yōu)化過程中收斂速度慢且容易陷入局部最優(yōu)。為此,提出一種基于M-H采樣的快速反向微分進化算法。M-H采樣用于ODE算法的變異操作,滿足馬爾可夫鏈可逆條件。馬爾可夫鏈的一步轉移概率根據(jù)個體等級分配的選擇概率進行計算,既能選擇最優(yōu)個體,又能尋找優(yōu)化方向并保持種群多樣性。仿真結果表明,M-H采樣得到的個體具有馬爾可夫鏈平穩(wěn)分布特性,該算法在單峰函數(shù)和多峰函數(shù)優(yōu)化中都能快速收斂,全局和局部搜索性能達到平衡,具有較高的搜索精度及較好的魯棒性。

    微分進化算法;反向微分進化算法;轉移概率;平穩(wěn)分布;馬爾可夫鏈蒙特卡洛;反向學習

    1 概述

    進化算法是一種模擬生物進化過程求解優(yōu)化問題的啟發(fā)式自適應人工智能技術。1995年Storn和Price等人提出微分進化(Differential Evolution, DE)[1-2]算法,其主要特點是算法簡單、收斂速度快、可調參數(shù)少、魯棒性好,相對于其他優(yōu)化算法有較強的全局收斂能力和穩(wěn)定性。變異操作是微分進化的核心操作,對微分進化的影響很大,文獻[3]是關于微分進化算法的參數(shù)和變異策略的綜述和應用。近年來DE成功應用于不同領域,如:工程優(yōu)化設計,數(shù)字濾波器設計,圖像處理,數(shù)據(jù)挖掘,多傳感器融合等。但是微分進化算法在高維多峰函數(shù)優(yōu)化中容易陷入局部最優(yōu),收斂速度慢、優(yōu)化性能不穩(wěn)定。針對DE算法的缺陷,許多學者提出了很多改進方法,如文獻[4]采用自適應控制參數(shù)來提高微分進化算法的收斂性能,但是收斂速度依然較慢;文獻[5]在文獻[4]的基礎上增加了種群規(guī)模減少機制和3個變異策略,改善收斂速度和DE算法的魯棒性,但是算法性能不穩(wěn)定,易于陷入局部最優(yōu);文獻[6]引入基于排序的變異,將種群個體進行適應度排序后,進行迭代更新,維持局部搜索和全局搜索的平衡;文獻[7]將種群分成多個子群,動態(tài)調整每個子種群的個體數(shù)目,增加子種群間個體信息交換,并且采用局部搜索和自適應方法,然而該算法參數(shù)較多,選擇困難,且常由于參數(shù)選擇不當導致性能不穩(wěn)定;文獻[8]提出一種分階段的二次變異,提高算法的全局搜索能力,但增加了時間復雜度;文獻[9]采用鄰域變異方法,在某個設定的領域內取得變異的個體,這種方法易于陷入局部最優(yōu);文獻[10]通過求反向種群來保持種群多樣性,但是反向個體容易偏離全局最優(yōu)個體,更新速率較慢,不能很快地達到全局最優(yōu)。

    本文采用馬爾可夫鏈蒙特卡洛(Markov Chain Monte Carlo,MCMC)思想,提出基于 Metropolis-Hastings(M-H)采樣的算法用于反向微分進化(Opposition Differential Evolution,ODE)[10]算法(MHODE)的變異操作。M-H算法具有馬爾可夫鏈的平穩(wěn)性,所采樣的個體具有平穩(wěn)分布的特性,根據(jù)該采樣進行變異操作,能使M-HODE算法的收斂速度加快,收斂趨于平穩(wěn)。

    2 微分進化和反向微分進化算法

    2.1 微分進化算法

    圖1 DE算法流程

    微分進化算法流程:

    (1)初始化種群

    在決策空間X內,用式(1)隨機產(chǎn)生初始向量:

    (2)變異操作

    微分進化算法的差分向量與縮放因子相乘后跟基向量進行向量合成,一般采用DE/rand/1變異策略,變異操作的公式為:

    (3)交叉操作

    對變異產(chǎn)生的新個體和當前個體進行交叉操作,以增加種群個體的多樣性。經(jīng)過交叉操作得到實驗向量:

    其交叉操作公式如下:

    其中,j=1,2,…,D;CR∈[0,1]為交叉概率;jrand∈[1,D]為均勻選取的隨機整數(shù)。

    (4)選擇操作

    DE算法通過選擇操作,對實驗個體和當前個體進行適應度評價,根據(jù)式(4)決定候選個體是否取代當前個體。

    2.2 反向微分進化算法

    2.2.1 反向學習理論

    反向微分進化算法是基于反向學習(Oppositionbased Learning,OBL)理論的微分進化算法。OBL理論思想如下:

    定義1(反向個體) 在多維空間R,p=(x1,x2,…,xD)是D維空間的一個個體,x1,x2,…,xD∈R,xi∈[a,b]?i∈{1,2,…,D},式(5)求反向個體=。

    2.2.2 ODE算法步驟

    ODE算法步驟具體如下:

    步驟1 基于反向學習的種群初始化。

    (1)隨機產(chǎn)生均勻分布的種群Po,大小為NP。

    (2)用定義1中的式(5)來計算Po中每一個個體的反向個體opoi,j=aj+bj-poi,j,使用反向優(yōu)化方法從集合{po,Opo}中選NP個適應度最好的個體作為初始種群。

    步驟2 根據(jù)迭代條件,對種群個體進行微分進化算法的變異、交叉、選擇操作。

    步驟3 隨機反向代跳操作,若隨機的rand(0, 1)<Jr(Jr是跳轉概率),MINPj和MAXP

    j分別是當前種群P的區(qū)間最小和最大值,用式(6)求出反向個體,然后從集合{P,OP}選擇Np個最優(yōu)個體作為當前代種群P。

    opi,j=MINpj+MAXp

    j-pi,j(6)

    步驟4 判斷是否達到迭代終止條件,否則轉向步驟2。

    3 基于M-H采樣的反向微分進化算法

    3.1 M-H采樣方法

    3.1.1 馬爾可夫鏈基本思想

    M-H算法是出現(xiàn)較早一般化的馬爾可夫鏈蒙特卡洛[11-12]采樣方法,下面先介紹馬爾可夫鏈(Markov Chain)。

    定義3(馬爾可夫鏈) 假定隨機序列{X0,X1,…}滿足馬爾可夫性,即在任意時刻t≥0時,序列中某時刻t+1的狀態(tài)Xt+1由條件分布p(Xt+1|Xt)產(chǎn)生,它只依賴時刻t的狀態(tài),與之前的狀態(tài){X0,X1,…,Xt-1}無關,這樣的一個隨機序列被稱為馬爾可夫鏈,馬爾可夫鏈的轉移核表示如下:

    p(x,x′)=p(x→x′)=p(xt+1=x′|xt=x) (7)

    定義4(馬爾可夫鏈可逆) 若π(dx)P(x,dy)= π(dy)p(y,dx),x,y∈X則狀態(tài)空間X上的馬爾可夫鏈關于π(·)可逆。

    定義5(平穩(wěn)分布) 設πj(t)=π{X(t)=j},j∈X,如果關于π(x)可逆的馬爾可夫鏈{X(t),t≥0}滿足: πj=lim

    t→∞πj(t),j∈X,則稱π(x)為該馬爾可夫鏈的平穩(wěn)分布,平穩(wěn)分布馬爾科夫鏈的狀態(tài)轉移與初始值狀態(tài)無關,只與時間間隔有關。

    3.1.2 M-H采樣思想

    MCMC[13-14]基本原理是建立一個平穩(wěn)分布,采樣得到的馬爾可夫鏈樣本是一個π(x)樣本,其基本步驟可概括為:

    (1)構造馬爾可夫鏈,在空間X的樣本上選擇一個合適的馬爾可夫鏈,轉移核設為p(*,*),使其收斂到平穩(wěn)分布π(x)。

    (2)樣本提取過程:由空間X的某一點出發(fā)X0,用步驟(1)構造的馬爾可夫鏈進行抽樣模擬,產(chǎn)生序列X1,X2,…,Xn。

    (3)蒙特卡洛積分:對任意整數(shù)m和任意足夠大的整數(shù)n,任一個目標函數(shù)的f(x)期望估計為: f(x(t)) (8)

    M-H方法最早由Metropolis等人提出,之后由Hastings加以推廣,形成Metropolis-Hastings方法。

    原理如下:假設馬爾可夫鏈第t個時刻的狀態(tài)為xt,π(x)是求解問題的目標分布,W(x,xt)是對稱的轉移函數(shù)。Metropolis-Hasting算法通過以下2步得到t+1時刻的狀態(tài):

    (1)從轉移分布W(x,x′)中得到x′,這里W(x, x′)是對稱的概率轉移函數(shù),即W(x,x′)=W(x′,x)。

    (2)取隨機數(shù)U,使U服從(0,1)均勻隨機分布,令r=min(1,π(x′)W(x′,x) π(x)W(x,x′) E[f(x)]= 1 n-m

    n

    t=m+1

    ),如果U≤r,則令xt+1=x′;否則令xt+1=x。

    3.2 M-HODE算法描述

    微分進化算法中最核心的操作是變異操作,MHODE算法提出一種,將Metropolis-Hastings采樣方法用于ODE的變異操作,M-HODE在初始化種群后求反向種群,合并初始種群和反向種群并計算每個種群的適應度值,根據(jù)適應度值大小進行升序排列,選取前NP個個體作為初始種群。在變異操作中用M-H采樣方法選擇基向量和差分向量,采樣個體組成的馬爾可夫鏈滿足馬爾可夫平穩(wěn)分布。在進行變異,交叉,選擇操作后,隨機的進行反向種群更新,保持了種群的多樣性,利于 M-HODE算法的全部優(yōu)化。

    3.3 M-HODE算法步驟

    M-HODE算法步驟具體如下:

    步驟1 隨機產(chǎn)生均勻分布的種群P0,種群大小為NP,用式(5)計算得到種群個體的反向個體opoi,j=aj+bj-poi,j。從集合 P0,OP0),則可得r=min(1, π(x′) π(x) { }中選NP個個體作為初始種群。

    步驟2 將種群里所有個體的按適應度值進行升序排列,計算排列后每一個個體的等級 Ri,公式為:

    Ri=Np-i,i=1,2,…,Np (9)

    步驟3 對每個向量個體進行等級分配后,計算每個個體的選擇概率,這里用到已提出的平方模型式(10)[15],向量xi的選擇概率計算公式為:

    pi= RiNp ,i=1,2,…,Np (10)

    步驟4 使用Metropolis-Hastings算法進行抽樣參加變異操作的個體。以DE/rand/1策略作為實例,選取xr0,xr1,xr2,xr3為馬爾可夫鏈。

    M-H采樣方法具體如下:

    (1)當rand(0,1)>pr0且r0≠i條件滿足時,隨機選擇r0∈{1,Np}。/*xr0作為馬爾可夫鏈的初始

    ■■

    2

    ■■向量*/

    (3)如果U≤r,則x1=xr1,否則x1=xr0。/*選取x1為基向量*/

    (5)如果U≤r,則x2=xr2,否則x2=xr1并且x1=xr2。/*x2作為差分向量*/

    (6)r3≠r1、r3≠r0和r3≠i條件滿足時隨機選取r3∈{1,NP}。/*隨機選取x3*/

    步驟5 當函數(shù)評價次數(shù)小于最大評價次數(shù)以及迭代次數(shù)小于最大迭代次數(shù)2個條件滿足時,進行差分進化算法的變異、交叉、選擇操作得到種群P。

    步驟6 隨機反向代跳操作:如果跳轉概率Jr大于(0,1)間的一個隨機數(shù),那么根據(jù)式(6)求出種群P的反向種群OP,然后從集合{P,OP}中選擇NP個個體作為當前代種群P。

    步驟7 重復迭代到算法停止迭代的條件為止。

    4 仿真實驗與結果分析

    為測試M-HODE算法的有效性和正確性,這里用11個常用的標準測試函數(shù)進行仿真評價,將MHODE,DE,JDE及ODE算法進行比較。這11個測試函數(shù)中既有高維單峰也有高維多峰函數(shù)。其中,f4,f5,f73個函數(shù)是特殊維數(shù),分別是10維、2維、2維,其他函數(shù)均是30維,另外,f6是一個噪聲函數(shù)。測試函數(shù)的部分參數(shù)設置如表1所示。f5的最優(yōu)值fmin=-1,其他函數(shù)的最優(yōu)值fmin都為0。

    本文實驗仿真環(huán)境為Matlab2012,運行于Intel (R)Pentium(R)4,CPU 2.93 GHz,1 GB內存的聯(lián)想臺式電腦。仿真時設置最大迭代次數(shù)MAXIter= 1 000,函數(shù)適應度最大評價次數(shù)Max_NFFES= 10 000×D,縮放因子和交叉概率初始值分別設置為F=0.5和CR=0.8,反向操作的跳轉概率Jr按照ODE算法的Jr取一樣值。為驗證該算法的精確性,將每個算法獨立運行50次,取50次實驗結果的平均值和標準方差作為評價算法的主要性能指標,為驗證算法的快速收斂,在實驗條件的基礎上加上目標值VTR=1.0E-10的條件,計算50次所用平均時間。

    表1 測試函數(shù)及相關參數(shù)

    表2是DE,JDE,ODE及M-HODE算法的實驗結果,加粗數(shù)據(jù)表示最優(yōu)值。統(tǒng)計表2中的測試函數(shù)數(shù)據(jù),f0~f4這5個高維單峰函數(shù)中有4個函數(shù)M-HODE均值優(yōu)于DE,JDE,ODE算法;2維函數(shù)f5的實驗結果都是一樣的,針對噪聲函數(shù)f6,M-HODE方差比其他3個算法方差更小,說明M-HODE具有更好的抗噪性;多峰函數(shù)f7~f10中有3個函數(shù) MHODE均值優(yōu)于DE,JDE,ODE算法,M-HODE方差比其他3個算法較小,該優(yōu)越性歸功于M-HODE算法中M-H采樣的平穩(wěn)性??梢?針對一些高維多峰函數(shù)優(yōu)化,M-HODE具有明顯優(yōu)勢。

    圖2~圖4給出了f1,f2和f103個具有典型高維測試函數(shù)的迭代曲線,迭代次數(shù)為1 000次,按一定比例用線條標出,從圖2~圖4可以看出,M-HODE算法、DE算法、JDE算法、ODE算法曲線趨勢相似,但是M-HODE算法有更快的收斂速度和更優(yōu)的精度。圖4對于f10函數(shù),M-HODE算法表現(xiàn)出更加明顯的收斂速度和收斂精度。這歸功于M-HODE采樣方法的穩(wěn)定性和反向進化的種群多樣性,促進了算法的快速收斂。針對其余的測試函數(shù),實驗也獲得了相似的進化曲線。

    表2 DE,JDE,ODE,M-HODE算法均值和方差

    圖2 f1函數(shù)迭代曲線

    圖3 f2函數(shù)迭代曲線

    圖4f10函數(shù)迭代曲線

    圖5是DE,JDE,ODE和M-HODE算法運行時間比較,仿真設置為VTR=1.0E-10,最大迭代次數(shù)MaxIter=1 000,函數(shù)適應度最大評價次數(shù)Max_NFFES=10 000D。從圖5可以看出,在運行時間上M-HODE優(yōu)于DE,JDE,ODE這3個算法,針對有些函數(shù)M-HODE算法所用迭代時間甚至可以快到上百倍。M-H采樣很大程度上不僅可以提高差微分進化函數(shù)的優(yōu)化精度,而且可以加快函數(shù)優(yōu)化的收斂速度,使優(yōu)化很快趨于平穩(wěn),并減少時間復雜度。

    圖5 算法運行時間比較

    5 結束語

    本文提出一種基于M-H采樣的快速反向微分進化算法,將M-H采樣策略用于反向微分進化算法中,M-H采樣個體組成平穩(wěn)分布的馬爾可夫鏈,能改進ODE算法的變異操作。通過對11個典型的Benchmark函數(shù)進行測試,結果表明,M-HODE算法不僅在高維單峰函數(shù)和高維多峰函數(shù)優(yōu)化中具有較高的精度,而且M-H采樣能加快M-HODE算法的收斂速度,避免陷入局部最優(yōu),從而實現(xiàn)M-H快速反向微分進化算法。

    [1] Storn R,Price K.Differential Evolution——A Simple and Efficient Adaptive Scheme for Global Optimization over Continuous Spaces[R].Berkley,USA:ICSI,Technical Report:TR-95-012,1995.

    [2] Storn R,Price K.Differential Evolution——A Simple and Efficient Heuristic for Global Optimization over Continuous Spaces[J].Journal of Global Optimization,1997,11(4): 341-359.

    [3] Mallipeddi R,Suganthan P N,Pan Q K,etal.Differential Evolution Algorithm With Ensembleof Parameters and Mutate on Strategies[J].Applied Soft Computing,2011,11(2):1679-1696.

    [4] Brest J,Greiner S,Boskovic B,et al.Self-adapting Control Parameters in Differential Evolution: A Comparative Study on Numerical Benchmark Problems[J].IEEE Transactions on Evolutionary Computation, 2006,10(6):646-657.

    [5] Brest J,Maucec M S.Self-adaptive Differential Evolution Algorithm Using Population Size Reduction and Three Stragies[J].Soft Computing,2011,15(11):2157-2174.

    [6] Gong Wenyin,Cai Zhihua.Differential Evolution with Ranking-based Mutation Operators[J].IEEE Transactions on Evolutionary Computation,2013,43(6):1-16.

    [7] 張雪霞,陳維榮,戴朝華.帶局部搜索的動態(tài)多群體自適應差分進化算法及函數(shù)優(yōu)化[J].電子學報,2010, 38(8):1825-1830.

    [8] 王筱珍,李 鵬,俞國燕.分階段二次變異的多目標混沌差分進化算法[J].控制與決策,2011,26(3):457-463.

    [9] Qu B Y,Suganthan P N,Liang J J.Differential Evolution with Neighborhood Mutation for Multimodal Optimization[J].IEEE Transactions on Evolutionary Computation,2012,16(5):601-614.

    [10] Rahnamayan S,Tizhoosh H R,Salama M M A.Oppositionbased Differential Evolution[J].IEEE Transactions on Evolutionary Computation,2008,12(1):64-79.

    [11] Karandikar R L.On the Markov Chain Monte Carlo (MCMC)Method[J].Sadhana,2006,31(2):81-104.

    [12] Gallagher K,Charvin K,Nielsen S,et al.Markov Chain Monte Carlo(MCMC)Sampling Methods to Determine Optimal Models,Model Resolution and Model Choice for Earth Science Problems[J].Marine and Petroleum Geology,2009,26(4):525-535.

    [13] Chen Peide.Hasting-metropolis Algorithms and Reference Measures[J].Elsevier Statistics&Probability Letters,1998, 38(4):323-328.

    [14] Eidsvik J,Tjelmeland H.On Directional Metropolis——Hastings Algorithms[J].Statistics and Computing, 2006,16(1):93-106.

    [15] Kalelo P,Ali M M.Differential Evolution Algorithm Using Hybrid Mutation[J].Computational Optimization and Application,2007,37(2):231-246.

    編輯 陸燕菲

    Fast Opposition Differential Evolution Algorithm Based on M-H Sampling

    TU Weiwei,GE Hongwei,YANG Jinlong,YUAN Yunhao
    (School of IOT Engineering,Jiangnan University,Wuxi 214211,China)

    In Differential Evolution(DE)algorithm,the population initialization is updated by using opposition optimization rule,so as to maintain the population diversity.However,the reverse individuals are easy to deviate from the global optimal solution,has slow convergence speed and easy to fall into local optimum in function optimization.This paper proposes a fast Opposition Differential Evolution(ODE)algorithm based on M-H(Metropolis-Hastings)sampling method.M-H sampling is used in the mutation operation of ODE.M-H sampling satisfies Markov Chain reversible conditions.One step transition probability of Markov Chain is calculated according to the selecting probability of individual ranking-assignment,not only chooses the best individual,but also searches for the optimal direction and keeps the population diversity.Simulation results show these individuals from M-H sampling have Markov stationary distribution property.The algorithm can accelerate the speed of convergence in unimodal functions and multimodal functions,balance the performance of global searching and local searching,and has higher precision and better robustness.

    Differential Evolution(DE)algorithm;Opposition Differential Evolution(ODE)algorithm;transition probability;stationary distribution;Markov Chain Monte Carlo(MCMC);Opposition-based Learning(OBL)

    1000-3428(2014)11-0155-05

    A

    TP301.6

    10.3969/j.issn.1000-3428.2014.11.031

    國家自然科學基金資助項目(61305017);江蘇省自然科學基金資助項目(BK20130154);江蘇高校優(yōu)勢學科建設工程基金資助項目。

    涂維維(1987-),女,碩士研究生,主研方向:微分進化,人工智能,模式識別;葛洪偉,教授、博士;楊金龍、袁運浩,副教授、博士。

    2013-12-13

    2014-01-13E-mail:tuweiweier@163.com

    中文引用格式:涂維維,葛洪偉,楊金龍,等.基于M-H采樣的快速反向微分進化算法[J].計算機工程,2014,40(11): 155-159.

    英文引用格式:Tu Weiwei,Ge Hongwei,Yang Jinlong,et al.Fast Opposition Differential Evolution Algorithm Based on M-H Sampling[J].Computer Engineering,2014,40(11):155-159.

    猜你喜歡
    優(yōu)化
    超限高層建筑結構設計與優(yōu)化思考
    PEMFC流道的多目標優(yōu)化
    能源工程(2022年1期)2022-03-29 01:06:28
    民用建筑防煙排煙設計優(yōu)化探討
    關于優(yōu)化消防安全告知承諾的一些思考
    一道優(yōu)化題的幾何解法
    由“形”啟“數(shù)”優(yōu)化運算——以2021年解析幾何高考題為例
    圍繞“地、業(yè)、人”優(yōu)化產(chǎn)業(yè)扶貧
    事業(yè)單位中固定資產(chǎn)會計處理的優(yōu)化
    消費導刊(2018年8期)2018-05-25 13:20:08
    4K HDR性能大幅度優(yōu)化 JVC DLA-X8 18 BC
    幾種常見的負載均衡算法的優(yōu)化
    電子制作(2017年20期)2017-04-26 06:57:45
    久久久久久久精品吃奶| 久久久久国产一级毛片高清牌| 国产真人三级小视频在线观看| 99久久久亚洲精品蜜臀av| 一级黄色大片毛片| 国产精品电影一区二区三区| 色综合欧美亚洲国产小说| 亚洲五月婷婷丁香| 国产精品久久久久久人妻精品电影| 女性生殖器流出的白浆| 十八禁网站免费在线| 18禁国产床啪视频网站| 男女下面进入的视频免费午夜 | 欧美精品啪啪一区二区三区| 日日爽夜夜爽网站| 中文在线观看免费www的网站 | 精品久久久久久久末码| 亚洲精品美女久久av网站| 久久人人精品亚洲av| 在线十欧美十亚洲十日本专区| 久久久久久久久中文| 国内久久婷婷六月综合欲色啪| 欧美色欧美亚洲另类二区| 日本撒尿小便嘘嘘汇集6| 国产日本99.免费观看| 精品欧美国产一区二区三| 1024手机看黄色片| 中文字幕人妻熟女乱码| 日韩高清综合在线| 精品欧美国产一区二区三| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品美女特级片免费视频播放器 | 成人午夜高清在线视频 | 日韩一卡2卡3卡4卡2021年| 最近在线观看免费完整版| 国产精品 国内视频| 国产亚洲av嫩草精品影院| 午夜福利高清视频| 黄色片一级片一级黄色片| cao死你这个sao货| 亚洲天堂国产精品一区在线| 国产熟女午夜一区二区三区| 99久久无色码亚洲精品果冻| 国产成人精品无人区| www.www免费av| 丝袜人妻中文字幕| 免费人成视频x8x8入口观看| 一级片免费观看大全| 又大又爽又粗| 久9热在线精品视频| 亚洲欧美激情综合另类| 午夜福利成人在线免费观看| 午夜成年电影在线免费观看| 中文字幕最新亚洲高清| 久久青草综合色| 真人做人爱边吃奶动态| 国产精品亚洲美女久久久| 无限看片的www在线观看| 国产亚洲av高清不卡| ponron亚洲| 亚洲精品一区av在线观看| 午夜福利欧美成人| 草草在线视频免费看| 国产精品二区激情视频| 黄色视频,在线免费观看| 亚洲中文av在线| 亚洲九九香蕉| 国产精品综合久久久久久久免费| 国产免费男女视频| 亚洲国产精品合色在线| 妹子高潮喷水视频| 大型黄色视频在线免费观看| 亚洲成人精品中文字幕电影| 国产一区在线观看成人免费| 宅男免费午夜| 18禁国产床啪视频网站| 国产私拍福利视频在线观看| 亚洲激情在线av| 国产91精品成人一区二区三区| 宅男免费午夜| 色综合站精品国产| 男女之事视频高清在线观看| 国产精品久久久久久人妻精品电影| 18美女黄网站色大片免费观看| 日韩有码中文字幕| 少妇的丰满在线观看| 国产久久久一区二区三区| 亚洲人成网站高清观看| 色播在线永久视频| 曰老女人黄片| 午夜两性在线视频| 伊人久久大香线蕉亚洲五| 波多野结衣av一区二区av| 久久狼人影院| 天堂动漫精品| 欧美日韩亚洲综合一区二区三区_| 午夜亚洲福利在线播放| 国产精品久久久久久人妻精品电影| 午夜福利视频1000在线观看| 18禁黄网站禁片午夜丰满| 老司机午夜福利在线观看视频| 午夜久久久在线观看| 色尼玛亚洲综合影院| 九色国产91popny在线| 国产1区2区3区精品| 久久久久免费精品人妻一区二区 | 国产亚洲av高清不卡| 欧美三级亚洲精品| 2021天堂中文幕一二区在线观 | 一进一出抽搐动态| 日韩一卡2卡3卡4卡2021年| 91麻豆av在线| 亚洲欧洲精品一区二区精品久久久| 午夜免费鲁丝| av电影中文网址| 亚洲熟女毛片儿| 精品人妻1区二区| 亚洲片人在线观看| 亚洲精品美女久久久久99蜜臀| 久久久水蜜桃国产精品网| 国产精品av久久久久免费| 色播亚洲综合网| 久久久久九九精品影院| 非洲黑人性xxxx精品又粗又长| 法律面前人人平等表现在哪些方面| 我的亚洲天堂| 欧美 亚洲 国产 日韩一| 变态另类成人亚洲欧美熟女| 成人国产一区最新在线观看| 美女午夜性视频免费| 天堂动漫精品| 搡老妇女老女人老熟妇| 午夜福利视频1000在线观看| 97超级碰碰碰精品色视频在线观看| 在线av久久热| 99久久久亚洲精品蜜臀av| 欧美又色又爽又黄视频| 欧美黑人精品巨大| 久久久水蜜桃国产精品网| 国产精品久久久久久人妻精品电影| 欧美日韩一级在线毛片| 午夜福利成人在线免费观看| 午夜免费鲁丝| 中文字幕精品免费在线观看视频| 国产一级毛片七仙女欲春2 | 国产极品粉嫩免费观看在线| 欧美中文日本在线观看视频| 男女那种视频在线观看| 亚洲国产欧美网| 国产精品99久久99久久久不卡| 大型黄色视频在线免费观看| 搡老岳熟女国产| 国产成年人精品一区二区| 一夜夜www| 香蕉国产在线看| 成人亚洲精品av一区二区| 99国产精品99久久久久| 国产乱人伦免费视频| 亚洲精品一区av在线观看| 波多野结衣高清作品| 欧美激情高清一区二区三区| 男人操女人黄网站| 法律面前人人平等表现在哪些方面| 亚洲av中文字字幕乱码综合 | 中文字幕另类日韩欧美亚洲嫩草| 99在线视频只有这里精品首页| 成人亚洲精品av一区二区| 中文资源天堂在线| 欧美日韩黄片免| 在线观看一区二区三区| 国产精品自产拍在线观看55亚洲| 禁无遮挡网站| 美女扒开内裤让男人捅视频| 一级毛片高清免费大全| 又紧又爽又黄一区二区| 午夜福利在线在线| 美女高潮喷水抽搐中文字幕| 国产成人精品久久二区二区91| 久久久精品国产亚洲av高清涩受| 我的亚洲天堂| 亚洲精品一卡2卡三卡4卡5卡| 国产野战对白在线观看| 国产精品日韩av在线免费观看| 禁无遮挡网站| 色av中文字幕| 看片在线看免费视频| 亚洲三区欧美一区| 一级a爱片免费观看的视频| 欧美激情极品国产一区二区三区| 精品久久久久久久末码| 精品久久蜜臀av无| 国产成人精品久久二区二区免费| 精品一区二区三区视频在线观看免费| 国产精品香港三级国产av潘金莲| 91成年电影在线观看| 久9热在线精品视频| 亚洲国产精品久久男人天堂| 黑人操中国人逼视频| 欧美亚洲日本最大视频资源| 美女大奶头视频| www.999成人在线观看| 99国产精品一区二区蜜桃av| 在线av久久热| av电影中文网址| 成人国产一区最新在线观看| 欧美黑人欧美精品刺激| a级毛片在线看网站| 亚洲第一欧美日韩一区二区三区| 香蕉久久夜色| 久久精品人妻少妇| 国产精品一区二区三区四区久久 | 中出人妻视频一区二区| 波多野结衣巨乳人妻| 亚洲精品国产区一区二| 他把我摸到了高潮在线观看| 亚洲av第一区精品v没综合| 欧美中文日本在线观看视频| 成人精品一区二区免费| 亚洲专区国产一区二区| 久久人人精品亚洲av| 国产91精品成人一区二区三区| 成人午夜高清在线视频 | 少妇粗大呻吟视频| 久久中文看片网| 狠狠狠狠99中文字幕| xxxwww97欧美| 真人做人爱边吃奶动态| 久久精品91蜜桃| 男女做爰动态图高潮gif福利片| 两个人视频免费观看高清| e午夜精品久久久久久久| 久久青草综合色| 亚洲电影在线观看av| 国产av一区二区精品久久| 亚洲精品中文字幕一二三四区| 免费在线观看成人毛片| 最新美女视频免费是黄的| 中文字幕精品免费在线观看视频| 亚洲精品一区av在线观看| 黄色视频,在线免费观看| 中文字幕精品亚洲无线码一区 | 国产成人欧美| 99精品欧美一区二区三区四区| 91国产中文字幕| 日本一本二区三区精品| 成年版毛片免费区| 久久亚洲精品不卡| 国产精品久久视频播放| 国产精品久久电影中文字幕| 久久精品亚洲精品国产色婷小说| 男人舔奶头视频| 在线观看一区二区三区| 村上凉子中文字幕在线| 成在线人永久免费视频| 两人在一起打扑克的视频| 亚洲av片天天在线观看| 麻豆久久精品国产亚洲av| 美女高潮喷水抽搐中文字幕| a在线观看视频网站| 国内久久婷婷六月综合欲色啪| 色婷婷久久久亚洲欧美| 欧美在线黄色| 亚洲aⅴ乱码一区二区在线播放 | 女人被狂操c到高潮| 久久午夜综合久久蜜桃| 曰老女人黄片| 婷婷六月久久综合丁香| 中文字幕另类日韩欧美亚洲嫩草| 99久久精品国产亚洲精品| 黄色毛片三级朝国网站| 身体一侧抽搐| 国产一区二区在线av高清观看| 后天国语完整版免费观看| 国内精品久久久久精免费| 老熟妇仑乱视频hdxx| 岛国视频午夜一区免费看| 午夜福利在线观看吧| 日本在线视频免费播放| 亚洲一区中文字幕在线| 午夜福利欧美成人| 日日干狠狠操夜夜爽| 首页视频小说图片口味搜索| 欧美 亚洲 国产 日韩一| 亚洲美女黄片视频| 久久中文字幕人妻熟女| e午夜精品久久久久久久| 女人爽到高潮嗷嗷叫在线视频| 国产亚洲精品第一综合不卡| 国产成人精品无人区| 亚洲国产精品久久男人天堂| 亚洲一区二区三区色噜噜| 长腿黑丝高跟| 97超级碰碰碰精品色视频在线观看| 欧美性长视频在线观看| 亚洲精品久久成人aⅴ小说| 1024视频免费在线观看| svipshipincom国产片| 国产精品亚洲一级av第二区| 99在线人妻在线中文字幕| x7x7x7水蜜桃| 欧美性长视频在线观看| 黄色女人牲交| 欧美人与性动交α欧美精品济南到| 啦啦啦观看免费观看视频高清| 精品人妻1区二区| 母亲3免费完整高清在线观看| 黄片播放在线免费| 成人永久免费在线观看视频| 午夜福利在线在线| 亚洲黑人精品在线| 我的亚洲天堂| 国产精品久久久久久亚洲av鲁大| 免费看美女性在线毛片视频| 法律面前人人平等表现在哪些方面| 日本成人三级电影网站| 国产精品二区激情视频| 日韩一卡2卡3卡4卡2021年| 熟妇人妻久久中文字幕3abv| 男女做爰动态图高潮gif福利片| 在线观看免费日韩欧美大片| 波多野结衣巨乳人妻| 久久国产精品男人的天堂亚洲| 国产精品 国内视频| 久久青草综合色| 免费看a级黄色片| 亚洲黑人精品在线| 国内精品久久久久精免费| 免费看a级黄色片| 99久久无色码亚洲精品果冻| 日本三级黄在线观看| 一区二区三区激情视频| 国产一区二区激情短视频| 久久草成人影院| 亚洲五月色婷婷综合| 国产亚洲欧美在线一区二区| 性欧美人与动物交配| 变态另类成人亚洲欧美熟女| 美女高潮到喷水免费观看| www国产在线视频色| 国产亚洲av高清不卡| 国产精品日韩av在线免费观看| 国产单亲对白刺激| 午夜两性在线视频| 91成年电影在线观看| 国产精品久久电影中文字幕| 好男人电影高清在线观看| 免费观看人在逋| 久久国产精品人妻蜜桃| 国产精品久久久久久亚洲av鲁大| 色播亚洲综合网| 91大片在线观看| 一进一出抽搐gif免费好疼| 免费搜索国产男女视频| 亚洲在线自拍视频| 亚洲国产精品999在线| 亚洲成av人片免费观看| 亚洲国产欧美网| 国产精品久久久人人做人人爽| 成年人黄色毛片网站| 国产亚洲精品一区二区www| 国产在线观看jvid| 国产欧美日韩一区二区精品| 亚洲国产欧美网| 丁香欧美五月| 国产熟女午夜一区二区三区| 搞女人的毛片| 免费在线观看完整版高清| 性欧美人与动物交配| 一区二区三区激情视频| 制服丝袜大香蕉在线| 黄片播放在线免费| 久久人人精品亚洲av| 午夜精品在线福利| 国产伦一二天堂av在线观看| 香蕉av资源在线| 国产亚洲欧美98| 亚洲一区二区三区色噜噜| 国产亚洲精品久久久久5区| 国产aⅴ精品一区二区三区波| 国产精品影院久久| 男女视频在线观看网站免费 | 脱女人内裤的视频| 国产真人三级小视频在线观看| 99精品久久久久人妻精品| 女性被躁到高潮视频| 99在线人妻在线中文字幕| 亚洲avbb在线观看| 亚洲成人久久性| 日本a在线网址| 黄片播放在线免费| 亚洲美女黄片视频| 日韩大码丰满熟妇| www.熟女人妻精品国产| 久久久久国产精品人妻aⅴ院| 免费在线观看完整版高清| 欧美中文日本在线观看视频| 精品国产亚洲在线| 亚洲av熟女| 可以在线观看毛片的网站| 一a级毛片在线观看| 人人妻人人澡人人看| 老司机深夜福利视频在线观看| 他把我摸到了高潮在线观看| 99热6这里只有精品| 亚洲 欧美一区二区三区| 人人妻人人澡欧美一区二区| av超薄肉色丝袜交足视频| 日本在线视频免费播放| 日韩欧美 国产精品| 中亚洲国语对白在线视频| 国产成+人综合+亚洲专区| 免费高清视频大片| 亚洲第一欧美日韩一区二区三区| 国产99白浆流出| 免费搜索国产男女视频| 亚洲国产欧美网| 亚洲熟妇中文字幕五十中出| 国产久久久一区二区三区| 91成人精品电影| 男人操女人黄网站| 国产精品国产高清国产av| 成人免费观看视频高清| 在线播放国产精品三级| 2021天堂中文幕一二区在线观 | 一边摸一边做爽爽视频免费| 亚洲七黄色美女视频| 一区二区三区国产精品乱码| 91麻豆精品激情在线观看国产| 午夜福利成人在线免费观看| 亚洲专区中文字幕在线| 99riav亚洲国产免费| 又黄又粗又硬又大视频| 国产精品九九99| 国产在线精品亚洲第一网站| 男女做爰动态图高潮gif福利片| 美女 人体艺术 gogo| 亚洲中文字幕日韩| 啪啪无遮挡十八禁网站| 好看av亚洲va欧美ⅴa在| 国产亚洲欧美在线一区二区| 最近最新中文字幕大全电影3 | 久久久久国内视频| 国产成人av教育| 女同久久另类99精品国产91| 欧美黑人巨大hd| 日本熟妇午夜| a级毛片a级免费在线| 中文字幕av电影在线播放| 国产三级黄色录像| 狂野欧美激情性xxxx| 亚洲国产精品sss在线观看| 午夜福利在线观看吧| 777久久人妻少妇嫩草av网站| 国产高清有码在线观看视频 | 国产单亲对白刺激| 亚洲精品国产区一区二| 嫁个100分男人电影在线观看| 国产免费男女视频| 欧美日韩中文字幕国产精品一区二区三区| 日韩欧美国产在线观看| www日本黄色视频网| 在线观看免费午夜福利视频| 在线观看66精品国产| 日日干狠狠操夜夜爽| 精品国产国语对白av| 天堂√8在线中文| 亚洲人成网站高清观看| 久久精品夜夜夜夜夜久久蜜豆 | 大型av网站在线播放| 国产亚洲精品av在线| 日韩三级视频一区二区三区| 日韩欧美国产在线观看| 亚洲国产看品久久| 在线看三级毛片| ponron亚洲| 1024视频免费在线观看| 亚洲一区二区三区不卡视频| 90打野战视频偷拍视频| 日日夜夜操网爽| 88av欧美| 亚洲男人的天堂狠狠| 一级a爱片免费观看的视频| 久久这里只有精品19| 999久久久国产精品视频| 国产色视频综合| 欧美成人午夜精品| 亚洲aⅴ乱码一区二区在线播放 | 成年版毛片免费区| 一夜夜www| 欧美乱码精品一区二区三区| 色综合站精品国产| 可以在线观看毛片的网站| 波多野结衣高清作品| 久久久久精品国产欧美久久久| avwww免费| 人妻久久中文字幕网| 亚洲成人精品中文字幕电影| x7x7x7水蜜桃| 欧美+亚洲+日韩+国产| 在线观看免费午夜福利视频| 亚洲成人久久性| 丝袜人妻中文字幕| 亚洲欧美精品综合一区二区三区| 国产精品久久久久久人妻精品电影| 亚洲va日本ⅴa欧美va伊人久久| 欧美性猛交黑人性爽| 日韩大尺度精品在线看网址| 国产欧美日韩一区二区三| www.精华液| 一区福利在线观看| 色综合欧美亚洲国产小说| 久久国产乱子伦精品免费另类| 久久精品国产99精品国产亚洲性色| 亚洲真实伦在线观看| 在线看三级毛片| 国产精品影院久久| 在线看三级毛片| 国产男靠女视频免费网站| e午夜精品久久久久久久| 亚洲真实伦在线观看| 非洲黑人性xxxx精品又粗又长| 亚洲全国av大片| 女警被强在线播放| 大型黄色视频在线免费观看| 国产真实乱freesex| ponron亚洲| 无人区码免费观看不卡| 亚洲全国av大片| 精品一区二区三区四区五区乱码| 午夜福利高清视频| 女警被强在线播放| 97碰自拍视频| 久久久精品国产亚洲av高清涩受| 伊人久久大香线蕉亚洲五| 国产成人一区二区三区免费视频网站| 精品电影一区二区在线| 可以在线观看的亚洲视频| 国内少妇人妻偷人精品xxx网站 | 久久人人精品亚洲av| 色尼玛亚洲综合影院| 久热这里只有精品99| 国产av一区在线观看免费| 欧美午夜高清在线| 午夜久久久在线观看| 少妇被粗大的猛进出69影院| 欧美大码av| 男女床上黄色一级片免费看| 我的亚洲天堂| 亚洲精品粉嫩美女一区| 国产精品野战在线观看| 国产av不卡久久| 中文字幕高清在线视频| 国产精品九九99| 99久久99久久久精品蜜桃| 黄频高清免费视频| 婷婷亚洲欧美| 视频在线观看一区二区三区| 亚洲五月色婷婷综合| 这个男人来自地球电影免费观看| 国产亚洲av高清不卡| 成人手机av| 国产爱豆传媒在线观看 | 成年女人毛片免费观看观看9| 色在线成人网| 亚洲一区高清亚洲精品| 日韩 欧美 亚洲 中文字幕| √禁漫天堂资源中文www| 国产免费av片在线观看野外av| 久久久久精品国产欧美久久久| 亚洲,欧美精品.| 精品国内亚洲2022精品成人| 欧美性猛交黑人性爽| 国产熟女xx| cao死你这个sao货| 精品不卡国产一区二区三区| 嫁个100分男人电影在线观看| 老司机靠b影院| 美女 人体艺术 gogo| 波多野结衣av一区二区av| 成熟少妇高潮喷水视频| 他把我摸到了高潮在线观看| 国产在线精品亚洲第一网站| 亚洲精品美女久久久久99蜜臀| 最新美女视频免费是黄的| АⅤ资源中文在线天堂| 亚洲av成人一区二区三| 日本a在线网址| 777久久人妻少妇嫩草av网站| 视频在线观看一区二区三区| 国产精品亚洲一级av第二区| 久久香蕉激情| 午夜福利在线在线| 人妻久久中文字幕网| 男女下面进入的视频免费午夜 | 在线观看66精品国产| 久久九九热精品免费| 久久欧美精品欧美久久欧美| www国产在线视频色| 在线永久观看黄色视频| 黄色a级毛片大全视频| 国产午夜精品久久久久久| 国产精品亚洲美女久久久| 免费人成视频x8x8入口观看| 动漫黄色视频在线观看| 18禁黄网站禁片午夜丰满| 久久人人精品亚洲av| 亚洲国产精品久久男人天堂| 亚洲熟妇中文字幕五十中出| 精品国产国语对白av| av视频在线观看入口| 一二三四在线观看免费中文在| 禁无遮挡网站| 久久精品人妻少妇| 色综合欧美亚洲国产小说| 日韩欧美三级三区|