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

    改進(jìn)的基于差分進(jìn)化的群集蜘蛛優(yōu)化算法*

    2021-06-26 01:57:38魯海燕胡士娟沈莞薔
    傳感器與微系統(tǒng) 2021年6期
    關(guān)鍵詞:雌性雄性蜘蛛

    向 蕾,魯海燕,2,胡士娟,沈莞薔,2

    (1.江南大學(xué) 理學(xué)院,江蘇 無(wú)錫 214122;2.無(wú)錫市生物計(jì)算工程技術(shù)研究中心,江蘇 無(wú)錫 214122)

    0 引 言

    近些年,新型群智能體算法在函數(shù)優(yōu)化、工程應(yīng)用等方面得到了快速發(fā)展,基于生物自然特性的群智能算法不斷被提出,例如蟻群[1](artificial bee colony,ABC)算法、粒子群優(yōu)化[2](particle swarm optimization,PSO)算法、蝙蝠算法[3](bat algorithm,BA)、螢火蟲(chóng)優(yōu)化[4](glowworm swarm optimization,GSO)算法等。2013年Cuevas E提出了一種模擬自然界蜘蛛互相協(xié)作行為的新型群體智能優(yōu)化算法—群集蜘蛛優(yōu)化(social spider optimization,SSO)算法[5],其研究表明,與PSO和ABC算法相比,SSO算法具有更快的收斂速度和更高的求解精度。目前,國(guó)內(nèi)外許多學(xué)者已把SSO算法應(yīng)用于約束優(yōu)化[6]、經(jīng)濟(jì)負(fù)荷調(diào)度[7]、路徑規(guī)劃[8]問(wèn)題中。雖然SSO算法在一些實(shí)際問(wèn)題中已經(jīng)得到了成功的應(yīng)用,但與一些先進(jìn)算法相比,仍存在求解精度低、收斂速度慢、易陷入局部極值等缺陷。因此,一些學(xué)者提出了多種改進(jìn)的算法。王艷嬌等人[9]提出一種基于動(dòng)態(tài)學(xué)習(xí)策略的群集蜘蛛優(yōu)化算法,該算法利用云模型對(duì)蜘蛛位置更新公式中的隨機(jī)部分進(jìn)行改進(jìn),同時(shí)加入動(dòng)態(tài)學(xué)習(xí)策略,提高了算法的收斂速度;趙汝鑫等人[10]采用反向?qū)W習(xí)和精英選擇機(jī)制進(jìn)行種群初始化,使算法的全局搜索能力得到了提高;朱林波等人[11]將群集蜘蛛算法和差分進(jìn)化算法相結(jié)合,提出了基于差分進(jìn)化變異策略的群集蜘蛛算法,該算法對(duì)交配后的新一代蜘蛛進(jìn)行變異,有效地提高了種群多樣性。

    本文提出了帶有差分變異算子和粒子群慣性權(quán)重及學(xué)習(xí)因子的改進(jìn)SSO算法。利用差分進(jìn)化(differential evolution,DE)算法的變異操作,更新部分雌蜘蛛的位置,并引入一組非線性的慣性權(quán)重及反余弦非對(duì)稱的學(xué)習(xí)因子,使得算法前期具有較快的收斂速度,后期具有較強(qiáng)的開(kāi)發(fā)能力。

    1 SSO算法

    在SSO算法中,整個(gè)搜索空間被視作蜘蛛運(yùn)動(dòng)的蜘蛛網(wǎng),而每只蜘蛛的位置則被作為目標(biāo)優(yōu)化函數(shù)的候選解,且每只蜘蛛的權(quán)重大小由蜘蛛的適應(yīng)值大小決定。蜘蛛個(gè)體之間通過(guò)分工合作、信息交流和繁衍后代等行為實(shí)現(xiàn)對(duì)問(wèn)題最優(yōu)解進(jìn)行搜索。

    1.1 種群初始化

    設(shè)種群數(shù)目設(shè)為N,雌性蜘蛛數(shù)目為Nf=[(0.9-0.25·rand)·N],雄性蜘蛛數(shù)目為Nm=N-Nf,雌性蜘蛛種群為F={f1,f2,…,fNf},雄性蜘蛛種群為M={m1,m2,…,mNm}。

    雌性和雄性蜘蛛初始化種群公式如下

    i∈{1,2,…,Nf},h∈(1,2,…,Nm},

    j∈{1,2,…,D}

    (1)

    式中 [·]為向下取整函數(shù),rand為[0,1]之間隨機(jī)數(shù),fij和mij分別為第i個(gè)雌性和雄性個(gè)體的第j維變量,fjmax和fjmin為雌性蜘蛛個(gè)體第j維變量的上限和下限,mjmax和mjmin為雄性蜘蛛個(gè)體第j維變量的上限和下限。

    1.2 雌性蜘蛛位置更新

    雌性蜘蛛根據(jù)自身振動(dòng)的大小吸引或排斥其他蜘蛛,其位置更新模式如下

    (2)

    式中Sc為距離個(gè)體蜘蛛fi最近且體重比它重的蜘蛛fc的體重,Sb為蜘蛛群體中體重最重的蜘蛛fb的體重,Vibci為個(gè)體蜘蛛fi和個(gè)體蜘蛛fc之間的振動(dòng)信息,PF為概率因子,α,β,δ,rand均為[0,1]之間的隨機(jī)數(shù)。

    雌性蜘蛛振動(dòng)信息公式如下

    (3)

    (4)

    式中J(Si)為第i個(gè)個(gè)體的目標(biāo)適應(yīng)值,dij為個(gè)體fi與fj的歐氏距離。

    1.3 雄性蜘蛛位置更新

    以中間個(gè)體(中間個(gè)體指的是按權(quán)重大小對(duì)蜘蛛排序處于中間位置的個(gè)體)的權(quán)重ωNf+n為準(zhǔn)線,SSO算法把雄性蜘蛛分成兩類,其中個(gè)體權(quán)重大于ωNf+n的雄性蜘蛛為支配雄蜘蛛,否則為非支配雄蜘蛛。支配雄蜘蛛吸引與之靠近的雌蜘蛛,而非支配雄蜘蛛向著雄蜘蛛群的中間個(gè)體靠近。雄蜘蛛位置更新模式如下

    (5)

    1.4 交配行為

    在雄性蜘蛛中,支配雄蜘蛛會(huì)與其交配半徑r內(nèi)的雌蜘蛛進(jìn)行交配,交配半徑定義如下

    (6)

    交配過(guò)程中,交配蜘蛛越重對(duì)新的子代蜘蛛個(gè)體影響越大,每只蜘蛛的影響可能性為Psk,根據(jù)輪盤賭的方法確定,公式如下

    (7)

    雌性和雄性蜘蛛交配完成后,若新的子代優(yōu)于上一代中體重最輕的個(gè)體,則由新的蜘蛛代替上一代中體重最輕的個(gè)體;否則,保持原蜘蛛種群不變。

    2 改進(jìn)的群集蜘蛛算法

    2.1 差分變異策略

    DE算法是一種基于種群個(gè)體之間合作與競(jìng)爭(zhēng)行為的群體進(jìn)化算法,主要包括變異、交叉和選擇三種操作。首先,DE算法對(duì)父代種群進(jìn)行變異和交叉操作,產(chǎn)生子代種群,然后利用選擇操作,選取適應(yīng)值較優(yōu)的個(gè)體作為新一代種群。變異操作中的變異個(gè)體由下式產(chǎn)生,vi=xr1+F·(xr2-xr3),i∈{1,2,…,Np},r1≠r2≠r3≠i。其中,xr2-xr3為偏差變量;F為變異算子,控制偏差變量的放大作用。由上式可知,DE算法中的變異個(gè)體由基向量以及父代種群中的兩個(gè)個(gè)體的偏差變量組成。因?yàn)樵诟复N群中隨機(jī)選取兩個(gè)個(gè)體的組合繁多,所以算法在迭代過(guò)程中能維持良好的多樣性。變異算子的大小決定算法向偏差向量學(xué)習(xí)的能力大小,因而合理設(shè)置變異算子能夠使得算法具有良好的勘探能力。

    本文借鑒DE算法中的變異操作,對(duì)原群集蜘蛛算法中的部分雌蜘蛛進(jìn)行變異操作,具體如下:

    設(shè)置閾值m,隨機(jī)產(chǎn)生一個(gè)數(shù)Z∈(0,1),當(dāng)Z≤m時(shí),雌蜘蛛變異公式如下

    (8)

    式中fp,fq由雌性蜘蛛種群F={f1,f2,…,fNf}中隨機(jī)產(chǎn)生,且p≠q;T為變異算子,這里令T=g/G,g為當(dāng)前迭代次數(shù),G為總迭代次數(shù)??梢钥闯?,變異算子T隨著迭代次數(shù)增加動(dòng)態(tài)變大,在迭代初期,T值比較小,偏差變量對(duì)變異影響小,蜘蛛的變異較大程度受群體最優(yōu)的蜘蛛影響,蜘蛛會(huì)較快地往群體最優(yōu)個(gè)體靠近;在算法運(yùn)行后期,蜘蛛聚集在最優(yōu)的蜘蛛周圍,T值比較大,此時(shí)偏差變量對(duì)變異的影響較大,使得種群在最優(yōu)蜘蛛周圍進(jìn)行局部的精細(xì)搜索,提高算法搜索精度。

    2.2 非線性慣性權(quán)重策略與反余弦非對(duì)稱學(xué)習(xí)因子策略

    本文借鑒文獻(xiàn)[12]中粒子群算法的反余弦非對(duì)稱的學(xué)習(xí)因子的搜索機(jī)制,為蜘蛛位置更新設(shè)計(jì)了一組新的慣性權(quán)重k及學(xué)習(xí)因子c,其公式如下

    (9)

    本文設(shè)置kmax=1.1,kmin=0.4使得慣性權(quán)重k在0.2~0.9之間非線性增大,這時(shí)種群向自身學(xué)習(xí)力度最好;c1,c2的選取與文獻(xiàn)[12]一致,c1max=2.75,c1min=1.25,c2max=2.25,c2min=0.5;i為當(dāng)前迭代次數(shù),T為最大迭代次數(shù)。

    本文所提出的慣性權(quán)重k及學(xué)習(xí)因子c隨迭代次數(shù)增加的變化曲線如圖1所示。

    圖1 慣性權(quán)重k及加速因子c1和c2的迭代曲線

    當(dāng)設(shè)置kmax=1.1,kmin=0.4時(shí),如圖1(a)所示,慣性權(quán)重k在0.2~0.9之間隨著迭代次數(shù)的增加非線性增大,這里將非線性增大的慣性權(quán)重加入基本群集蜘蛛算法的位置更新中。由圖1(b)可以看出,c1和c2呈現(xiàn)非對(duì)稱的增大與減小,使得種群在算法初期更多的獲得全局最優(yōu)個(gè)體的信息,從而以較快的速度進(jìn)入局部搜索,在算法后期,種群更多的獲得局部最優(yōu)個(gè)體的信息并保持一定的搜索速度,從而進(jìn)行局部的精細(xì)搜索。根據(jù)上述策略,以下給出雌雄蜘蛛的位置更新公式。

    雌蜘蛛和雄性蜘蛛位置更新公式如下

    (10)

    (11)

    式中k,c1,c2為式(9)給出的非線性慣性權(quán)重及反余弦非對(duì)稱學(xué)習(xí)因子。

    2.3 本文算法步驟

    1) 初始化參數(shù):初始化種群數(shù)目N;初始解的維度D;權(quán)重系數(shù)的最大值kmax和最小值kmin;學(xué)習(xí)因子的最大值cmax和最小值cmin。2) 隨機(jī)產(chǎn)生一組個(gè)數(shù)為N的初始解,分配雌性蜘蛛和雄性蜘蛛種群。3) 計(jì)算每一只蜘蛛相應(yīng)的適應(yīng)值。4) 隨機(jī)產(chǎn)生一個(gè)Z(Z∈(0,1)),若Z≤0.5,按式(8)更新雌性蜘蛛個(gè)體位置;否則,按式(10)更新雌性蜘蛛個(gè)體位置,并保留新老雌性蜘蛛個(gè)體中較為優(yōu)秀的個(gè)體作為下一代的蜘蛛個(gè)體。5) 按式(11)更新雄性蜘蛛個(gè)體的位置,并保留新老雄性蜘蛛個(gè)體中較為優(yōu)秀的50%作為進(jìn)入下一代的蜘蛛個(gè)體。6) 雌性蜘蛛進(jìn)行交配,選取適應(yīng)值優(yōu)的新個(gè)體保留。7) 判斷是否滿足結(jié)束條件,若滿足,直接輸出結(jié)果;否則返回步驟(3)。

    3 實(shí)驗(yàn)結(jié)果與分析

    3.1 基本測(cè)試函數(shù)

    本文選取了8個(gè)測(cè)試函數(shù)對(duì)本文算法的性能進(jìn)行對(duì)比測(cè)試。測(cè)試函數(shù)f1~f8依次為Sphere,Schwefel’s P2.22,Schwefel’s P1.2,Bentcigar ,Alpine,Griewank,Rastrigin,Ackley。其中,f1~f4為單峰函數(shù)和f5~f8為多峰函數(shù),理論上的最優(yōu)值均為0。

    3.2 仿真實(shí)驗(yàn)結(jié)果分析

    3.2.1 與相關(guān)改進(jìn)算法的比較

    將本文算法與一種基于差分進(jìn)化變異策略的改進(jìn)群集蜘蛛優(yōu)化(SSO-DE)算法[13]進(jìn)行對(duì)比。文獻(xiàn)[13]中所提出的SSO-DE算法和本文所提出的算法都利用了差分變異策略對(duì)雌性蜘蛛進(jìn)行了變異操作,不同點(diǎn)在于SSO-DE是利用兩個(gè)偏差變量誘導(dǎo)差分變異,本文算法則是利用一個(gè)偏差變量誘導(dǎo)差分變異。為了方便比較,本文將兩個(gè)算法的參數(shù)設(shè)置與文獻(xiàn)[13]保持一致,由于文獻(xiàn)[13]僅提供5個(gè)測(cè)試函數(shù)的數(shù)值,下面將兩種算法在文獻(xiàn)[13]中所涉及的5個(gè)經(jīng)典函數(shù)上進(jìn)行尋優(yōu),比較最優(yōu)解的平均值(mean)和標(biāo)準(zhǔn)方差(Std)來(lái)體現(xiàn)算法的優(yōu)越性。具體實(shí)驗(yàn)結(jié)果分別由表1給出。

    由表1所示SSO-DE和本文算法的實(shí)驗(yàn)結(jié)果可以看出,本文算法在單峰函數(shù)Sphere,Schwefel’s P2.2和Quartic上的求解精度和穩(wěn)定性高于SSO-DE算法。兩個(gè)算法在多峰函數(shù)Rastrigin上都尋找到了最優(yōu)解,在多峰函數(shù)Ackley上的求解精度幾乎相同。由上可知,本文算法的求解性能優(yōu)于SSO-DE算法。

    表1 SSO-DE和本文算法實(shí)驗(yàn)結(jié)果

    3.2.2 與其他幾種群智能算法的比較

    將DESSOcw算法與PSO、BA、SSO、隨機(jī)慣性權(quán)重的粒子群算法RandPSO(random weighting particle swarm optimization algorithm)算法進(jìn)行對(duì)比實(shí)驗(yàn),PSO算法參數(shù)設(shè)置參照文獻(xiàn)[2],BA算法參數(shù)設(shè)置參照文獻(xiàn)[3],RandPSO算法的參數(shù)參照文獻(xiàn)[14],SSO算法參數(shù)設(shè)置參照文獻(xiàn)[5],本文算法參數(shù)設(shè)置參照本文2.1,2.2節(jié)。實(shí)驗(yàn)分成3組,維數(shù)分別設(shè)置為10維、30維和50維,所有算法的種群群體規(guī)模為NP=50,最大迭代次數(shù)為1 000。這里記錄5種算法在8個(gè)經(jīng)典測(cè)試函數(shù)上分別獨(dú)立運(yùn)行30次的最優(yōu)解的平均值(Mean)和標(biāo)準(zhǔn)方差(Std)(見(jiàn)表2),同時(shí)列出以上5種算法求解問(wèn)題為30維時(shí)的收斂圖(見(jiàn)圖2)。

    表2 5種算法的測(cè)試結(jié)果

    圖2 5種算法在不同測(cè)試函數(shù)上的收斂曲線

    由表2的數(shù)據(jù)對(duì)比可以得出,與PSO,BA,SSO,Rand-PSO算法相比,無(wú)論在10維、30維還是50維的函數(shù)上,本文算法的求解精度和穩(wěn)定性都有了很大提高。針對(duì)4個(gè)單峰函數(shù),本文算法的求解精度十分接近理論最優(yōu)值,隨著維度的提高,本文算法的求解精度并沒(méi)有明顯下降。對(duì)于4個(gè)多峰函數(shù),本文算法在f6和f7兩個(gè)函數(shù)上可以求解到理論最優(yōu)值0,說(shuō)明本文提出的算法在求解函數(shù)f6和f7時(shí)顯示出良好的求解性能。同時(shí),在三種維度下,本文算法在8個(gè)測(cè)試函數(shù)上的標(biāo)準(zhǔn)方差最小,尤其在f1,f4,f6~f8函數(shù)上的標(biāo)準(zhǔn)方差都為0,說(shuō)明算法具有更強(qiáng)的穩(wěn)定性??傮w來(lái)說(shuō),本文提出的本文算法的求解精度和穩(wěn)定性均優(yōu)于其他4種算法,顯示了極高的求解性能。

    算法的收斂曲線可以直觀地對(duì)比各算法的收斂性能。由圖2(a)~(d)可知,在求解這4個(gè)單峰函數(shù)時(shí),PSO算法、RandPSO算法、BA算法、SSO算法的收斂曲線平緩,說(shuō)明這4種算法在收斂過(guò)程中多樣性差,易陷入局部最優(yōu)解,而本文算法的收斂曲線較為陡峭,尋找到的最優(yōu)解均優(yōu)于其他4個(gè)算法。由圖2(f)和圖2(g)可知,在求解Griewank和Rastrigin這兩個(gè)多峰函數(shù)時(shí),本文算法收斂速度和求解精度均優(yōu)于其他4種算法,不僅收斂到最優(yōu)解0,而且收斂速度極快。綜上所述,本文算法的收斂速度均優(yōu)于其他4種算法。

    為了便于觀察解集的分布情況,給出5個(gè)算法在4個(gè)經(jīng)典函數(shù)上獨(dú)立運(yùn)行30次的最優(yōu)解集箱須圖(見(jiàn)圖3)。從圖3(a)~(h)中可以看出,對(duì)于4個(gè)單峰函數(shù)和4個(gè)多峰函數(shù),本文算法得到的數(shù)據(jù)整體都比其他4種算法更接近橫坐標(biāo)軸,說(shuō)明本文算法得到的解均優(yōu)于其他算法。另外,本文算法所得的數(shù)據(jù)范圍均小于其他4種算法,且未出現(xiàn)異常解,說(shuō)明本文算法不僅收斂精度高,而且具有極高的穩(wěn)定性。

    圖3 5種算法在不同測(cè)試函數(shù)上最優(yōu)解的方差分析

    4 結(jié) 論

    本文算法利用差分進(jìn)化算法的變異策略,對(duì)部分雌蜘蛛進(jìn)行變異,從而增加了種群的多樣性;同時(shí),在個(gè)體更新過(guò)程中加入一組非線性的慣性權(quán)重及學(xué)習(xí)因子,以更好地平衡算法局部搜索和全局搜索能力。仿真實(shí)驗(yàn)表明,與相關(guān)算法相比,本文給出的改進(jìn)算法具有更高的求解精度和更好的穩(wěn)定性,有一定實(shí)際應(yīng)用的潛力。

    猜你喜歡
    雌性雄性蜘蛛
    來(lái)的是誰(shuí)?
    麥穗魚(yú)(雄性)
    垂釣(2023年11期)2024-01-21 16:07:04
    大鰭鱊(雄性)
    垂釣(2023年9期)2023-12-10 19:39:30
    連續(xù)超促排卵致腎精不足伴生育力低下雌性小鼠模型制備和比較研究
    小蜘蛛凍僵了,它在哪兒呢?
    河南一種雌性蚜蠅首次記述
    蜘蛛
    萌物
    飛碟探索(2016年5期)2016-05-10 23:44:30
    慢性焦慮刺激對(duì)成年雌性大鼠性激素水平的影響
    大蜘蛛
    国产单亲对白刺激| 大又大粗又爽又黄少妇毛片口| 嫩草影院入口| 色哟哟·www| 成人特级av手机在线观看| 18禁裸乳无遮挡免费网站照片| 亚洲国产精品成人久久小说| 久久久精品94久久精品| 国产免费视频播放在线视频 | 亚洲av免费高清在线观看| 亚洲综合精品二区| 又爽又黄a免费视频| 日韩欧美国产在线观看| 三级经典国产精品| 看非洲黑人一级黄片| 欧美一区二区亚洲| 18禁在线无遮挡免费观看视频| 免费在线观看成人毛片| 亚洲av中文字字幕乱码综合| 天堂中文最新版在线下载 | 午夜久久久久精精品| 国产一区有黄有色的免费视频 | 亚洲成色77777| 村上凉子中文字幕在线| 国产国拍精品亚洲av在线观看| 精品久久久久久久末码| 国产又黄又爽又无遮挡在线| 天堂中文最新版在线下载 | 波野结衣二区三区在线| 亚洲av男天堂| 久久久久九九精品影院| 丝袜喷水一区| 亚洲欧美精品自产自拍| 日韩欧美三级三区| 男人的好看免费观看在线视频| 蜜桃亚洲精品一区二区三区| 成人高潮视频无遮挡免费网站| 久久热精品热| 七月丁香在线播放| 午夜免费激情av| 国产在线一区二区三区精 | 亚州av有码| 亚洲自拍偷在线| 国产美女午夜福利| 亚洲欧美中文字幕日韩二区| 精品无人区乱码1区二区| 午夜激情欧美在线| 国产乱来视频区| 亚洲精品456在线播放app| 亚洲精品一区蜜桃| 小蜜桃在线观看免费完整版高清| 特大巨黑吊av在线直播| 一二三四中文在线观看免费高清| av卡一久久| 午夜激情福利司机影院| a级毛色黄片| 亚洲久久久久久中文字幕| 亚洲一区高清亚洲精品| 国产老妇伦熟女老妇高清| 六月丁香七月| 欧美成人a在线观看| 少妇裸体淫交视频免费看高清| 狂野欧美激情性xxxx在线观看| 丰满乱子伦码专区| 久久久久久久久久黄片| 国产成人免费观看mmmm| 日日撸夜夜添| 亚洲精品自拍成人| 日本免费在线观看一区| 亚洲第一区二区三区不卡| ponron亚洲| 国产精品一区二区三区四区免费观看| 国产精品麻豆人妻色哟哟久久 | 亚洲精品国产av成人精品| 国产精品久久视频播放| 久久久精品94久久精品| 久久精品人妻少妇| 成人毛片a级毛片在线播放| 午夜福利高清视频| 亚洲中文字幕一区二区三区有码在线看| 黄色配什么色好看| 天堂中文最新版在线下载 | 超碰av人人做人人爽久久| 噜噜噜噜噜久久久久久91| 日韩强制内射视频| 亚洲精品自拍成人| 国产成人精品久久久久久| 一级爰片在线观看| 国产成人a区在线观看| 99热6这里只有精品| 国产淫片久久久久久久久| 亚洲欧美一区二区三区国产| av卡一久久| 观看美女的网站| 欧美又色又爽又黄视频| 久久精品久久久久久噜噜老黄 | 亚洲人成网站在线播| 狠狠狠狠99中文字幕| 床上黄色一级片| 身体一侧抽搐| 精品久久久久久久人妻蜜臀av| 不卡视频在线观看欧美| 欧美高清成人免费视频www| 永久网站在线| 午夜a级毛片| 亚洲内射少妇av| 两个人视频免费观看高清| 中文字幕精品亚洲无线码一区| 日韩亚洲欧美综合| 国产精品久久久久久精品电影小说 | 最近中文字幕2019免费版| 亚洲av二区三区四区| 晚上一个人看的免费电影| 高清av免费在线| eeuss影院久久| 亚洲一级一片aⅴ在线观看| 天天躁日日操中文字幕| 免费在线观看成人毛片| 国产一区有黄有色的免费视频 | 亚洲图色成人| 亚洲人成网站在线观看播放| 久久精品91蜜桃| 久久国产乱子免费精品| 国产精品一区二区三区四区久久| 国国产精品蜜臀av免费| 插逼视频在线观看| 色尼玛亚洲综合影院| 国产一级毛片在线| 亚洲欧美精品综合久久99| 日本黄色片子视频| 桃色一区二区三区在线观看| 欧美高清性xxxxhd video| 亚洲18禁久久av| 成人av在线播放网站| 国产精品日韩av在线免费观看| 亚洲乱码一区二区免费版| 汤姆久久久久久久影院中文字幕 | 亚洲欧美日韩高清专用| 青春草国产在线视频| 国产一级毛片在线| 如何舔出高潮| 亚洲三级黄色毛片| 黄色配什么色好看| 我的老师免费观看完整版| 精品久久久久久久久亚洲| 五月玫瑰六月丁香| 久久久久久大精品| 日韩欧美精品免费久久| 午夜激情欧美在线| 亚洲国产精品国产精品| 久久亚洲国产成人精品v| 亚洲欧美精品自产自拍| 亚洲自拍偷在线| av视频在线观看入口| 欧美高清性xxxxhd video| av.在线天堂| 淫秽高清视频在线观看| 午夜福利成人在线免费观看| 大话2 男鬼变身卡| 视频中文字幕在线观看| 日韩欧美国产在线观看| 自拍偷自拍亚洲精品老妇| av在线播放精品| 国产乱人偷精品视频| 两个人的视频大全免费| 女人十人毛片免费观看3o分钟| 日韩一区二区三区影片| 成年版毛片免费区| 少妇的逼好多水| 亚洲人与动物交配视频| 少妇高潮的动态图| 久久久久久国产a免费观看| 成年免费大片在线观看| 天堂中文最新版在线下载 | 1000部很黄的大片| 一级毛片久久久久久久久女| 国产又黄又爽又无遮挡在线| 精品人妻视频免费看| 欧美成人a在线观看| 免费人成在线观看视频色| 亚洲精品久久久久久婷婷小说 | 久久这里有精品视频免费| 黄片无遮挡物在线观看| 亚洲五月天丁香| 色尼玛亚洲综合影院| 高清毛片免费看| 黑人高潮一二区| 欧美一区二区国产精品久久精品| 精品久久国产蜜桃| 日韩av在线大香蕉| 午夜老司机福利剧场| 亚洲欧美中文字幕日韩二区| 国产69精品久久久久777片| 午夜久久久久精精品| 国产探花极品一区二区| 麻豆久久精品国产亚洲av| eeuss影院久久| 男女国产视频网站| 亚洲国产色片| 高清在线视频一区二区三区 | 亚洲最大成人中文| 中文字幕av在线有码专区| 春色校园在线视频观看| 国产黄色视频一区二区在线观看 | 日日摸夜夜添夜夜添av毛片| 国产精品国产三级专区第一集| 亚洲精品色激情综合| 亚洲av电影在线观看一区二区三区 | 天堂中文最新版在线下载 | 女人十人毛片免费观看3o分钟| 最近视频中文字幕2019在线8| 中文乱码字字幕精品一区二区三区 | 我要搜黄色片| 亚洲最大成人手机在线| 91精品国产九色| 国产 一区 欧美 日韩| 国产精品99久久久久久久久| 国产乱人视频| 麻豆国产97在线/欧美| 日日摸夜夜添夜夜爱| 51国产日韩欧美| 成人美女网站在线观看视频| 国产午夜福利久久久久久| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲av中文字字幕乱码综合| 亚洲综合色惰| 国产精品美女特级片免费视频播放器| 成人高潮视频无遮挡免费网站| 亚洲av免费在线观看| 国产精品美女特级片免费视频播放器| 三级男女做爰猛烈吃奶摸视频| АⅤ资源中文在线天堂| 精品久久久久久久久久久久久| 久久热精品热| 亚洲四区av| 黄色配什么色好看| 国产黄a三级三级三级人| 老司机福利观看| 国产精品综合久久久久久久免费| 久久精品久久精品一区二区三区| 99热这里只有是精品50| 婷婷六月久久综合丁香| 日韩一区二区三区影片| 在线观看av片永久免费下载| 国产亚洲精品久久久com| 一个人看视频在线观看www免费| 美女国产视频在线观看| 欧美成人免费av一区二区三区| 亚洲国产高清在线一区二区三| 听说在线观看完整版免费高清| 尤物成人国产欧美一区二区三区| 亚洲欧美日韩无卡精品| 欧美xxxx性猛交bbbb| 中文亚洲av片在线观看爽| 免费看美女性在线毛片视频| 精品久久久久久久人妻蜜臀av| 成人国产麻豆网| 国产精华一区二区三区| 亚洲经典国产精华液单| 一区二区三区四区激情视频| 国产成人一区二区在线| 狂野欧美激情性xxxx在线观看| 国产av一区在线观看免费| 亚洲乱码一区二区免费版| 日本欧美国产在线视频| 国产熟女欧美一区二区| 亚洲精品影视一区二区三区av| 中文字幕久久专区| 搡女人真爽免费视频火全软件| 欧美xxxx性猛交bbbb| 你懂的网址亚洲精品在线观看 | 少妇高潮的动态图| 最近中文字幕高清免费大全6| 一级爰片在线观看| 日韩欧美三级三区| 老师上课跳d突然被开到最大视频| 亚洲精品日韩在线中文字幕| 国产精品一区www在线观看| 小说图片视频综合网站| 午夜福利在线观看吧| 亚洲欧美日韩东京热| 久久久精品94久久精品| 天堂av国产一区二区熟女人妻| 婷婷六月久久综合丁香| 久久久久久久国产电影| 看非洲黑人一级黄片| 日韩一区二区三区影片| 久久欧美精品欧美久久欧美| 在现免费观看毛片| 亚洲中文字幕一区二区三区有码在线看| 狠狠狠狠99中文字幕| 国产精品不卡视频一区二区| 99在线人妻在线中文字幕| 尤物成人国产欧美一区二区三区| 色噜噜av男人的天堂激情| 亚洲伊人久久精品综合 | 欧美成人免费av一区二区三区| 日韩国内少妇激情av| 亚洲成人中文字幕在线播放| 欧美一区二区亚洲| 如何舔出高潮| 人体艺术视频欧美日本| 亚州av有码| 国产精品国产三级国产专区5o | 一级二级三级毛片免费看| 欧美成人免费av一区二区三区| 亚洲欧美日韩东京热| 免费av毛片视频| 乱系列少妇在线播放| 国产午夜精品论理片| 日韩精品有码人妻一区| 最近手机中文字幕大全| 午夜视频国产福利| 国产午夜精品久久久久久一区二区三区| 久久久久久久久中文| 精品熟女少妇av免费看| 国产成人精品婷婷| 亚洲欧美中文字幕日韩二区| 久久这里有精品视频免费| 中文欧美无线码| 午夜福利在线观看吧| 中文字幕久久专区| 国产午夜精品久久久久久一区二区三区| 亚洲av一区综合| 亚洲久久久久久中文字幕| av.在线天堂| 青青草视频在线视频观看| 精品人妻一区二区三区麻豆| 男的添女的下面高潮视频| 久久99热这里只频精品6学生 | 精品国内亚洲2022精品成人| 国产亚洲5aaaaa淫片| 一个人观看的视频www高清免费观看| 99热全是精品| 一边摸一边抽搐一进一小说| 午夜精品国产一区二区电影 | 人人妻人人澡欧美一区二区| 色哟哟·www| 精品午夜福利在线看| 国产淫片久久久久久久久| 中文资源天堂在线| 精品人妻偷拍中文字幕| 97在线视频观看| 欧美极品一区二区三区四区| 三级国产精品欧美在线观看| 国产精品女同一区二区软件| 亚洲天堂国产精品一区在线| 又粗又爽又猛毛片免费看| 99国产精品一区二区蜜桃av| 成年版毛片免费区| 一区二区三区免费毛片| 美女脱内裤让男人舔精品视频| 69人妻影院| 久久久久久久久久成人| 一区二区三区四区激情视频| 国产白丝娇喘喷水9色精品| 69av精品久久久久久| 国产免费男女视频| 亚洲欧美日韩卡通动漫| 国产精品1区2区在线观看.| 国产精品久久久久久av不卡| 最近中文字幕2019免费版| 成人无遮挡网站| 国产精品一区二区在线观看99 | 国语对白做爰xxxⅹ性视频网站| 欧美激情久久久久久爽电影| av卡一久久| 亚洲av二区三区四区| 欧美成人a在线观看| 网址你懂的国产日韩在线| 嫩草影院新地址| 日本熟妇午夜| 久久韩国三级中文字幕| 插逼视频在线观看| 大香蕉97超碰在线| 成人毛片a级毛片在线播放| 欧美色视频一区免费| 1000部很黄的大片| 免费黄色在线免费观看| 国产一区二区三区av在线| 免费看av在线观看网站| 欧美bdsm另类| 国产乱人偷精品视频| 99久久精品热视频| 99热这里只有是精品50| 国产亚洲午夜精品一区二区久久 | 天堂影院成人在线观看| av又黄又爽大尺度在线免费看 | 精品国产三级普通话版| 亚洲精品一区蜜桃| 五月玫瑰六月丁香| 成人午夜精彩视频在线观看| 免费看美女性在线毛片视频| 九草在线视频观看| 日韩欧美精品免费久久| 国产免费一级a男人的天堂| 国产真实乱freesex| 中文字幕av成人在线电影| av女优亚洲男人天堂| 午夜福利高清视频| 久久精品久久久久久久性| 久久欧美精品欧美久久欧美| 久久草成人影院| 少妇丰满av| 亚洲国产精品久久男人天堂| 国产精品国产三级国产av玫瑰| 国产大屁股一区二区在线视频| 热99在线观看视频| 亚洲天堂国产精品一区在线| 久久精品影院6| 欧美变态另类bdsm刘玥| 亚洲av一区综合| 黄色配什么色好看| 插阴视频在线观看视频| 男的添女的下面高潮视频| 秋霞在线观看毛片| 精华霜和精华液先用哪个| 欧美成人午夜免费资源| 爱豆传媒免费全集在线观看| 赤兔流量卡办理| 欧美性猛交╳xxx乱大交人| 熟妇人妻久久中文字幕3abv| 91久久精品国产一区二区成人| 熟妇人妻久久中文字幕3abv| 男女视频在线观看网站免费| 国产极品精品免费视频能看的| 在现免费观看毛片| 2022亚洲国产成人精品| 国产片特级美女逼逼视频| 黄色日韩在线| 亚洲综合色惰| 97人妻精品一区二区三区麻豆| 观看美女的网站| 一级毛片aaaaaa免费看小| 久久精品影院6| 亚洲最大成人手机在线| 在线观看66精品国产| 亚洲国产精品成人综合色| 免费观看精品视频网站| 日韩大片免费观看网站 | 精品久久久久久成人av| 久久热精品热| 国产淫片久久久久久久久| 永久网站在线| 欧美zozozo另类| 亚洲欧美日韩东京热| 久久午夜福利片| 永久网站在线| 国产视频内射| 免费播放大片免费观看视频在线观看 | 赤兔流量卡办理| 观看美女的网站| 国产精品1区2区在线观看.| 国产亚洲av嫩草精品影院| 亚洲第一区二区三区不卡| 熟女人妻精品中文字幕| 日本午夜av视频| 亚洲国产精品久久男人天堂| 日本色播在线视频| 天天躁日日操中文字幕| 亚洲欧美清纯卡通| 日本熟妇午夜| 欧美一区二区国产精品久久精品| 美女高潮的动态| 精品人妻偷拍中文字幕| 日本欧美国产在线视频| 免费av不卡在线播放| 嫩草影院新地址| 国产伦理片在线播放av一区| 美女国产视频在线观看| 成年av动漫网址| 色噜噜av男人的天堂激情| 欧美3d第一页| 日本五十路高清| 少妇人妻精品综合一区二区| 欧美区成人在线视频| av卡一久久| 网址你懂的国产日韩在线| 亚洲欧美日韩卡通动漫| 国产免费男女视频| 国国产精品蜜臀av免费| 亚洲三级黄色毛片| 国产成人福利小说| 看免费成人av毛片| 老司机影院成人| 汤姆久久久久久久影院中文字幕 | 好男人在线观看高清免费视频| 日韩欧美 国产精品| 亚洲一区高清亚洲精品| 九九久久精品国产亚洲av麻豆| 日韩欧美国产在线观看| 精品久久国产蜜桃| www.av在线官网国产| 日韩高清综合在线| 久久欧美精品欧美久久欧美| 别揉我奶头 嗯啊视频| 床上黄色一级片| 国产成人精品一,二区| 我要看日韩黄色一级片| 一级黄色大片毛片| 99在线视频只有这里精品首页| 午夜福利在线在线| 免费黄网站久久成人精品| 久久精品影院6| 两性午夜刺激爽爽歪歪视频在线观看| 成人午夜高清在线视频| 人妻夜夜爽99麻豆av| 亚洲国产色片| 久久精品国产自在天天线| 国产色婷婷99| 国产黄片视频在线免费观看| 亚洲欧美精品自产自拍| 亚洲精品456在线播放app| 精品一区二区免费观看| 亚洲av中文字字幕乱码综合| 毛片一级片免费看久久久久| ponron亚洲| 中文字幕制服av| 亚洲av福利一区| 欧美日韩国产亚洲二区| 国产在视频线精品| 黄片无遮挡物在线观看| 99热这里只有是精品50| 一个人免费在线观看电影| 日本免费在线观看一区| 免费看av在线观看网站| 国产一区亚洲一区在线观看| 麻豆成人av视频| 国产精华一区二区三区| 欧美日韩在线观看h| 国产精品麻豆人妻色哟哟久久 | 久久精品久久久久久久性| 色吧在线观看| 久久久久久久久久成人| 精品欧美国产一区二区三| 亚洲精品日韩av片在线观看| 成人综合一区亚洲| 亚洲欧美日韩卡通动漫| ponron亚洲| 亚洲美女搞黄在线观看| 综合色av麻豆| 人人妻人人澡人人爽人人夜夜 | 久久久久久久国产电影| 亚洲国产欧美人成| 在线观看66精品国产| 少妇熟女欧美另类| 七月丁香在线播放| 国产伦一二天堂av在线观看| 大香蕉97超碰在线| www日本黄色视频网| 青春草视频在线免费观看| 久久6这里有精品| 午夜福利高清视频| 亚洲精品乱码久久久v下载方式| 长腿黑丝高跟| 国内精品一区二区在线观看| 水蜜桃什么品种好| 最新中文字幕久久久久| 午夜免费男女啪啪视频观看| 插阴视频在线观看视频| 少妇人妻一区二区三区视频| 99热精品在线国产| 六月丁香七月| 国产三级在线视频| 成人综合一区亚洲| 男女视频在线观看网站免费| 日日啪夜夜撸| 亚洲人成网站高清观看| av视频在线观看入口| 久久久亚洲精品成人影院| 尤物成人国产欧美一区二区三区| 九九热线精品视视频播放| av天堂中文字幕网| 日本黄色片子视频| 又粗又硬又长又爽又黄的视频| 国产午夜精品论理片| 日韩亚洲欧美综合| 久久精品国产自在天天线| 小蜜桃在线观看免费完整版高清| 久久6这里有精品| 亚洲精品自拍成人| 国产色爽女视频免费观看| 狠狠狠狠99中文字幕| 日韩一本色道免费dvd| 久久精品国产鲁丝片午夜精品| 综合色av麻豆| 久久精品91蜜桃| 国产精品蜜桃在线观看| 精品久久国产蜜桃| 欧美日韩综合久久久久久| 又爽又黄无遮挡网站| 精品久久国产蜜桃| 九九在线视频观看精品| 如何舔出高潮| 淫秽高清视频在线观看| 亚洲欧美精品自产自拍| 欧美高清性xxxxhd video| 欧美性猛交╳xxx乱大交人| 三级国产精品欧美在线观看| 亚洲图色成人| 视频中文字幕在线观看| 在线观看av片永久免费下载| 国产又黄又爽又无遮挡在线| 久久6这里有精品| 亚洲精品久久久久久婷婷小说 | 久久精品夜夜夜夜夜久久蜜豆| 99热精品在线国产| 丰满少妇做爰视频| 免费观看a级毛片全部| 天堂网av新在线| 亚洲精品影视一区二区三区av| 波多野结衣高清无吗| 亚洲在线观看片| 国产精品久久视频播放| 成人三级黄色视频|