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

    一種基于截?cái)鄼C(jī)制的穩(wěn)態(tài)優(yōu)化算法求解多目標(biāo)優(yōu)化問題

    2018-11-30 01:51:48荊東星張清安
    關(guān)鍵詞:支配穩(wěn)態(tài)種群

    荊東星 張清安

    (湘西民族職業(yè)技術(shù)學(xué)院計(jì)算機(jī)系 湖南 湘西 416000)

    0 引 言

    許多實(shí)際優(yōu)化問題,通常需要同時(shí)求解2~3個(gè)高度復(fù)雜、非線性并且?guī)в袥_突性的目標(biāo)[1]。這類問題通常被定義為多目標(biāo)優(yōu)化問題(MOPs)。而多目標(biāo)優(yōu)化算法(MOEAs)非常適合于求解此類問題。因此,在過去的幾十年里,MOEAs得到了飛速發(fā)展,并在工程領(lǐng)域得到廣泛應(yīng)用[2-4]。

    對(duì)于MOPs來說,MOEAs的一個(gè)關(guān)鍵優(yōu)勢(shì)是基于種群特性,它允許種群中的個(gè)體在單個(gè)執(zhí)行過程中同時(shí)收斂到Pareto面的不同區(qū)域。一般來說,可以將MOEAs的優(yōu)化目標(biāo)簡(jiǎn)單歸納為以下兩個(gè)方面[5]:

    1) 收斂:將種群到MOP的Pareto面的距離最小化。

    2) 分布:種群廣泛且均勻的分散在MOP的Pareto面上。

    為了實(shí)現(xiàn)以上目標(biāo),近年來,研究者提出了大量有效的MOEAs[6]。而根據(jù)其每次進(jìn)化過程中產(chǎn)生的子代個(gè)數(shù)可以將現(xiàn)有的MOEAs大致分為兩類:一類是基于種群進(jìn)化的MOEA[7];另一類是基于個(gè)體進(jìn)化的穩(wěn)態(tài)MOEA[8]。第一類中,具有代表性的算法有Zitzler等提出的SPEA[9]及其改進(jìn)版本SPEA2[10],Srinivas等[11]提出的NSGA,以及Deb等[7]在此基礎(chǔ)上提出的NSGA-II等。而第二類的代表性算法則有Deb等[8]提出的ε-MOEA?;诜N群的MOEA一直備受關(guān)注,而基于個(gè)體的穩(wěn)態(tài)MOEA自從2006年被提出后,在多目標(biāo)優(yōu)化領(lǐng)域鮮有文獻(xiàn)出現(xiàn)。

    因此,本文首先分析ε-MOEA所存在的問題,然后提出一種基于截?cái)鄼C(jī)制的穩(wěn)態(tài)MOEA。該算法首先采用Pareto支配控制進(jìn)化種群和歸檔種群的收斂性,然后采用一種新的截?cái)鄼C(jī)制維持歸檔集種群的分布性能。另外,本文的算法除默認(rèn)參數(shù)外(如:種群大小,變異率,交叉率等),無需任何參數(shù)。通過實(shí)驗(yàn)比較分析可以發(fā)現(xiàn),本文算法在求解2~3維MOPs時(shí),能夠獲得良好性能(收斂性和分布性)的解集。

    1 基本概念

    1.1 問題描述及相關(guān)性質(zhì)

    為了便于描述,本文所采用的優(yōu)化問題為最小化問題。一般地,最大化問題通常也可以通過公式轉(zhuǎn)換為最小化問題,問題的表達(dá)形式可以定義如下:

    minimizeF(x)=f1(x),…,fm(x))T

    (1)

    s.t.gi(x)≥0j=1,2,…,J

    hk(x)=0k=1,2,…,K

    x∈Ω

    在多目標(biāo)優(yōu)化中,以下概念得到了很好的定義和廣泛應(yīng)用。

    1) Pareto支配:對(duì)于兩個(gè)不同的解,x1,x2∈X,如果滿足式(2),則認(rèn)為x1支配x2,并可以用x1x2表示。

    (2)

    2) Pareto最優(yōu)解集:對(duì)于單個(gè)解x*∈X,如果不存在其他解x′∈X滿足x′x*,那么它將被視為Pareto最優(yōu)解。所有Pareto最優(yōu)解構(gòu)成了Pareto最優(yōu)解集。

    1.2 穩(wěn)態(tài)MOEA模型

    2003年由Deb等人首次提出穩(wěn)態(tài)MOEA[8],為多目標(biāo)優(yōu)化研究提供了一條新的思路。穩(wěn)態(tài)MOEA中,兩個(gè)種群(進(jìn)化種群EA和歸檔集Archive)同時(shí)進(jìn)行優(yōu)化。首先通過隨機(jī)初始化生成初始EA種群,將該種群的非支配解集存入到Archive內(nèi),然后分別從EA和Archive中隨機(jī)挑選一個(gè)個(gè)體(如圖1中的p和e)進(jìn)行交叉變異操作生成一個(gè)新的子代個(gè)體c。最后根據(jù)EA保留機(jī)制和Archive保留機(jī)制分別判斷子代個(gè)體能否被EA種群或Archive種群接納。其示意圖如圖1所示。

    圖1 穩(wěn)態(tài)算法模型

    一般而言,EA保留機(jī)制采用的是Pareto支配策略,通過Pareto支配策略將收斂性好的個(gè)體保留下來。而Archive保留機(jī)制不僅需要考慮種群的收斂性,還要保證種群能夠均勻且廣泛地分布在整個(gè)Pareto面上。例如,ε-MOEA的Archive保留機(jī)制采用的是ε-支配策略來保證種群的收斂性和分布性。

    2 ε-支配的缺陷

    ε-支配方法是由Laumanns等[12]提出,主要思想是根據(jù)值將目標(biāo)空間劃分為若干個(gè)網(wǎng)格,每個(gè)網(wǎng)格內(nèi)只允許一個(gè)個(gè)體存在。因此,當(dāng)兩個(gè)解p、q滿足下式條件時(shí),認(rèn)為pε-支配。

    (1-εi)·fi(p)≤fi(q) ?i∈{1,2,…,m}

    (3)

    式中:p和q為目標(biāo)空間的兩個(gè)解;m為目標(biāo)的維度。

    圖2給出了ε-支配的示意圖,個(gè)體p的ε-支配范圍為ABCDA所圍繞的區(qū)域,而原始Pareto支配范圍為PECFP構(gòu)成的區(qū)域??梢园l(fā)現(xiàn),ε-支配的區(qū)域要明顯大于Pareto支配的區(qū)域。在ABCDA區(qū)域的個(gè)體都將被p支配掉。另外,當(dāng)同一個(gè)網(wǎng)格內(nèi)存在多個(gè)解時(shí),如圖2中的1和2,ε-支配需要從中選擇一個(gè)性能最好的解保留下來。首先判斷這些解的Pareto支配關(guān)系,只考慮其中的非Pareto支配解。然后從中挑選距離所在網(wǎng)格的圓點(diǎn)最近的個(gè)體。在圖2中,解1離網(wǎng)格圓點(diǎn)更近,因此它將被保留。

    圖2 ε-支配

    進(jìn)一步發(fā)現(xiàn),種群的大小以及性能都依賴于值的設(shè)置,當(dāng)ε值設(shè)置過大,種群的數(shù)量將會(huì)過小,MOPs的某些區(qū)域可能很難搜索到;而當(dāng)值設(shè)置過小,種群的數(shù)量則會(huì)過于龐大而限制算法的優(yōu)化速度。如何為特定問題設(shè)置合適的值沒有一個(gè)很好的解決方案。另外,ε-支配很難保證種群的廣泛性,因其特有的性質(zhì),位于Pareto面邊界的解很容易被支配掉[13]。如圖3所示,目標(biāo)空間被劃分為400個(gè)均勻網(wǎng)格,網(wǎng)格的寬度等于ε。根據(jù)ε的性質(zhì),在目標(biāo)空間中最大能夠容納20個(gè)非ε支配解,然而ε-支配策略實(shí)際只能允許12個(gè)非ε支配解在目標(biāo)空間中。因?yàn)檫@些解將會(huì)把與所在網(wǎng)格平行或垂直的其他網(wǎng)格中的所有解支配掉。從圖3中解的分布可以發(fā)現(xiàn),靠近邊界的解的分布要比中間位置的稀疏。

    圖3 非支配解集

    為了使得Archive中的解的數(shù)量可控,且具有良好的性能。本文對(duì)ε-MOEA的Archive保留機(jī)制進(jìn)行了修改,并提出了一種新的穩(wěn)態(tài)優(yōu)化算法。

    3 一種新的穩(wěn)態(tài)算法

    在本文提出的算法(PTEA)中,采用了Pareto支配以及截?cái)鄼C(jī)制來共同維持Archive種群的性能。Pareto支配維持Archive種群的收斂性,截?cái)鄼C(jī)制維持Archive種群的分布性。

    在截?cái)鄼C(jī)制中,首先采用歐氏距離計(jì)算各個(gè)個(gè)體之間的距離,然后選擇最短距離的兩個(gè)個(gè)體,再判斷它們第二短的距離。如果第二短的距離也相等再判斷它們第三短的距離,依次類推,直到找到距離不一樣為止,然后刪除距離短的個(gè)體。其與SPEA2的修剪類似,唯一區(qū)別是SPEA2只考慮了兩次比較。而當(dāng)維度達(dá)到3維時(shí),在目標(biāo)空間中可能存在多個(gè)第二短距離相等的情況。

    PTEA的算法步驟如下:

    Step1初始化參數(shù)。設(shè)置種群大小N,最大進(jìn)化代數(shù)MAXGEN,交叉率Pc,變異率Pm。

    Step2初始化種群。采用隨機(jī)初始化方法初始化大小為N的EA種群。

    Step3構(gòu)建初始Archive。將EA種群進(jìn)行Pareto非支配排序,將非支配解集復(fù)制到Archive中。

    Step4生成單個(gè)子個(gè)體。從EA種群和Archive種群中隨機(jī)挑選一個(gè)個(gè)體進(jìn)行交叉變異產(chǎn)生新個(gè)體。

    Step5更新EA種群。如果子代個(gè)體被EA種群中的個(gè)體支配,不考慮更新。如果EA種群中不存在個(gè)體支配新的子代,判斷EA種群是否存在個(gè)體被子代支配。如果存在則隨機(jī)替換一個(gè)被支配的個(gè)體,如果不存在,則隨機(jī)替換EA種群中的個(gè)體。

    Step6更新Archive種群。首先判斷子代與Archive種群中個(gè)體的支配關(guān)系,如果子代被支配,不更新Archive種群。如果Archive種群中存在個(gè)體被子代支配,從中隨機(jī)挑選一個(gè)被子代替換。當(dāng)Archive種群中的解與子代互為非支配解時(shí),需考慮兩種情況:(1) Archive種群數(shù)量小于N-1,此時(shí)直接將子代保留到Archive中。(2) Archive種群數(shù)量大于N-1,此時(shí)引入截?cái)鄼C(jī)制替換Archive中的解。這樣可以保證種群大小為一個(gè)固定值。

    Step7判斷終止條件。本文的終止條件是進(jìn)化次數(shù)t=MAXGEN。如果未滿足條件則返回到Step 4,否則終止條件。

    一直以來,林語堂始終堅(jiān)持充分肯定女性的生命存在和人生價(jià)值,他呼吁:“我愿意看見新時(shí)代的女子,——她要無愧的標(biāo)立、表現(xiàn)、發(fā)揮女性的不同,建造新女性于別個(gè)的女性之上?!盵2]150可見,他的女性觀是進(jìn)步的,是符合時(shí)代潮流的。作者將木蘭作為貫穿小說的中心人物,既賦予她超凡脫俗的氣質(zhì)和才華,又使她扮演著賢妻良母的傳統(tǒng)角色,把經(jīng)過其文化道德篩選的全部女性美揉合于他的女主人公木蘭一身,使木蘭達(dá)到形體美、心靈美、人性美的高度統(tǒng)一,繼而升華到理想美的境界。[3]

    Step8輸出最終解集,算法結(jié)束。

    4 性能評(píng)價(jià)指標(biāo)及測(cè)試問題

    4.1 評(píng)價(jià)指標(biāo)

    為了更好地分析算法的性能,本文采用的評(píng)價(jià)指標(biāo)為反向世代距離(IGD)[14]。IGD作為綜合評(píng)價(jià)指標(biāo),通過計(jì)算Pareto面到種群的距離的平均值來判斷種群的分布性和收斂性。其計(jì)算公式表示為:

    (4)

    式中:P為最終種群;|P*|為Pareto面上參考點(diǎn)集P*的數(shù)量;d(x*,P)為參考點(diǎn)x*到離種群中各個(gè)解的最小距離。由式(4)可知,當(dāng)IGD值越小時(shí),種群的性能越好。

    4.2 測(cè)試問題

    本文選擇了8個(gè)具有代表性的測(cè)試問題[15-16]來驗(yàn)證算法的性能,分別為:

    1) 規(guī)整簡(jiǎn)單型Pareto面問題:DTLZ1、DTLZ2以及凸面DTLZ2(CDTLZ2);

    2) 倒轉(zhuǎn)Pareto面問題:IDTLZ1、IDTLZ2;

    3) 非連續(xù)問題:DTLZ7;

    5) 目標(biāo)范圍不一致問題:SDTLZ2。

    本文還對(duì)IGD數(shù)據(jù)進(jìn)行了顯著性比較[17],其中+、-、=分別表示比PETA顯著好、顯著差以及沒有顯著區(qū)別。

    5 實(shí)驗(yàn)結(jié)果

    本文挑選了3種具有代表性的多目標(biāo)優(yōu)化算法(ε-MOEA[8]、NSGA-II[7]以及PESA-II[18])進(jìn)行對(duì)比實(shí)驗(yàn),得到的3維優(yōu)化問題的最終解集。其中為穩(wěn)態(tài)算法,在進(jìn)化種群中采用Pareto支配維持進(jìn)化種群的收斂性,在歸檔集中采用ε-支配維持歸檔集種群的收斂和分布性。ε-MOEA最終目的是要獲得一組良好性能的歸檔集種群。NSGA-II采用非支配排序策略將種群劃分為若干個(gè)非支配層,再根據(jù)層序?qū)⒏复鷤€(gè)體逐層保留到下一代,當(dāng)保留下來的個(gè)體大于種群大小時(shí),則在當(dāng)前非支配層(臨界層)中挑選部分個(gè)體保留下來。NSGA-II在臨界層中采用聚集距離維持算法的分布性,從該層中挑選分布性好的個(gè)體保留下來。SPEA-II是在PESA[19]基礎(chǔ)上的改進(jìn),它將目標(biāo)空間劃分為若干個(gè)超網(wǎng)格,PESA-II的過程包括選擇超網(wǎng)格以及在該網(wǎng)格中隨機(jī)選擇一個(gè)個(gè)體保留。

    圖4給出了PTEA與其他3種算法在3目標(biāo)DTLZ1測(cè)試問題上的最終解集。從直觀上說,PTEA和ε-MOEA的分布性表現(xiàn)得最好,然而ε-MOEA很難搜索到Pareto面的邊界。 通過圖4可以看出ε-MOEA的廣泛性沒有PTEA的好。而其他兩個(gè)算法的最終解集的分布明顯不如前面的兩個(gè)算法。

    (a) PTEA (b) ε-MOEA

    (c) NSGA-II (d) PESA-II圖4 各算法在3目標(biāo)DTLZ1測(cè)試問題上的最終解集

    圖5給出的是各算法在3目標(biāo)DTLZ2測(cè)試問題上的最終解集,不管是均勻性,還是廣泛性,PTEA的表現(xiàn)都要優(yōu)于其他3種算法。PTEA的解集均勻地覆蓋了整個(gè)Pareto面。而對(duì)于其他3種算法來說,ε-MOEA和PESA-II的性能差不多,都要好于NSGA-II。NSGA-II的解集雖然能夠廣泛地分布在Pareto面上,但是它的均勻性不是很理想。

    (a) PTEA (b) ε-MOEA

    (c) NSGA-II (d) PESA-II圖5 各算法在3目標(biāo)DTLZ2測(cè)試問題上的最終解集

    在降維問題上,如圖6所示,PTEA已經(jīng)覆蓋到了整個(gè)Pareto面,而在Pareto面的兩端解的分布明顯稀疏與中間部分。而其他兩個(gè)算法同樣有小部分區(qū)域沒有解的分布。比如圖6(c)的Pareto面的上半部分區(qū)域以及圖6(d)的Pareto面的上半部分和中部區(qū)域相對(duì)來說都比較稀疏。

    (a) PTEA (b) ε-MOEA

    (c) NSGA-II (d) PESA-II圖6 各算法在3目標(biāo)DTLZ5測(cè)試問題上的最終解集

    圖7給出了4種算法在非連續(xù)問題(DTLZ7)的優(yōu)化結(jié)果,該問題的Pareto面由4個(gè)區(qū)域組成??梢园l(fā)現(xiàn),ε-MOEA的分布性要稍好于PTEA,但是廣泛性沒有本文算法好。在DTLZ7問題上,本文算法的分布性與NSGA-II和PESA-II沒有太大的區(qū)別。

    (a) PTEA (b) ε-MOEA

    (c) NSGA-II (d) PESA-II圖7 各算法在3目標(biāo)DTLZ7測(cè)試問題上的最終解集

    對(duì)于CDTLZ1測(cè)試問題,它是DTLZ1測(cè)試問題的變體,其Pareto面為凸面。如圖8所示,表現(xiàn)最好的是PTEA,其次是PESA-II。雖然ε-MOEA的解集看似均勻分布,但它沒有搜索到整個(gè)Pareto面,而是大部分解集都聚集到了Pareto面的中部位置。而與ε-MOEA的解集分布相反的是NSGA-II的大部分解集聚集在了Pareto面的邊界。

    (a) PTEA (b) ε-MOEA

    (c) NSGA-II (d) PESA-II圖8 各算法在3目標(biāo)CDTLZ1測(cè)試問題上的最終解集

    圖9和圖10給出了各算法在倒轉(zhuǎn)Pareto面問題上的解集分布,在這兩個(gè)問題上,PTEA的性能同樣也是最好的。

    (a) PTEA (b) ε-MOEA

    (c) NSGA-II (d) PESA-II圖9 各算法在3目標(biāo)IDTLZ1測(cè)試問題上的最終解集

    (a) PTEA (b) ε-MOEA

    (c) NSGA-II (d) PESA-II圖10 各算法在3目標(biāo)IDTLZ2測(cè)試問題上的最終解集

    最后是SDTLZ1問題,該問題的各個(gè)目標(biāo)的取值范圍都不相同,本文將各個(gè)維度i的取值范圍設(shè)置為[0,2i]。從圖11中解集的分布可以看出,PTEA和ε-MOEA的分布最均勻,而4種算法的分布廣泛性相當(dāng)。

    (a) PTEA (b) ε-MOEA

    (c) NSGA-II (d) PESA-II圖11 各算法在3目標(biāo)SDTLZ2測(cè)試問題上的最終解集

    因此,從分布均勻以及分布廣度兩個(gè)角度觀察,本文算法都要優(yōu)于其他3種算法。同時(shí),本文算法避免了額外的參數(shù)設(shè)置,極大提高了算法的實(shí)用性。

    本文除了從直觀上體現(xiàn)算法性能外,還給出了定量分析。如表1給出了4種算法在不同測(cè)試函數(shù)上綜合性能的數(shù)據(jù)結(jié)果,表中的數(shù)據(jù)為算法運(yùn)行30次統(tǒng)計(jì)得到的IGD平均值。 表1中的A、B、C、D、E、F、G、H分別表示DTLZ1、DTLZ2、DTLZ5、DTLZ7、CDTLZ2、IDTLZ1、IDTLZ2、SDTLZ2測(cè)試問題。

    表1 算法在各個(gè)測(cè)試問題上的IGD值

    續(xù)表1

    從表1中可以看出,總共16個(gè)測(cè)試問題,PTEA有14個(gè)表現(xiàn)得最好,僅3維的DTLZ7和3維的SDTLZ2分別稍差于NSGA-II和ε-MOEA算法。而就顯著性而言,本文算法有14個(gè)測(cè)試問題顯著優(yōu)于ε-MOEA,15個(gè)測(cè)試問題顯著優(yōu)于NSGA-II和PESA-II。

    6 結(jié) 語

    本文針對(duì)穩(wěn)態(tài)多目標(biāo)優(yōu)化算法,提出了一種基于截?cái)鄼C(jī)制的穩(wěn)態(tài)優(yōu)化算法(PTEA)。本文還通過一系列對(duì)比實(shí)驗(yàn)驗(yàn)證了該算法的性能。從算法的最終解集的分布以及IGD統(tǒng)計(jì)數(shù)據(jù)可以發(fā)現(xiàn),PTEA能有效地平衡種群的分布性和收斂性。

    雖然PTEA能夠很好地優(yōu)化2~3維的優(yōu)化問題,但是它很難在高維優(yōu)化問題上獲得良好的解集。其主要原因是PTEA依賴Pareto支配機(jī)制引導(dǎo)種群收斂,而Pareto支配隨著維度的增加其作用將急劇減弱。因此,引入一種新的收斂機(jī)制將是后續(xù)的研究方向。

    猜你喜歡
    支配穩(wěn)態(tài)種群
    邢氏水蕨成功繁衍并建立種群 等
    山西省發(fā)現(xiàn)刺五加種群分布
    可變速抽水蓄能機(jī)組穩(wěn)態(tài)運(yùn)行特性研究
    碳化硅復(fù)合包殼穩(wěn)態(tài)應(yīng)力與失效概率分析
    電廠熱力系統(tǒng)穩(wěn)態(tài)仿真軟件開發(fā)
    煤氣與熱力(2021年4期)2021-06-09 06:16:54
    被貧窮生活支配的恐懼
    意林(2021年9期)2021-05-28 20:26:14
    元中期歷史劇對(duì)社會(huì)穩(wěn)態(tài)的皈依與維護(hù)
    中華戲曲(2020年1期)2020-02-12 02:28:18
    跟蹤導(dǎo)練(四)4
    基于決策空間變換最近鄰方法的Pareto支配性預(yù)測(cè)
    隨心支配的清邁美食探店記
    Coco薇(2016年8期)2016-10-09 00:02:56
    欧美精品av麻豆av| 99热国产这里只有精品6| 欧美xxⅹ黑人| 精品福利永久在线观看| 国产av精品麻豆| 久久精品国产鲁丝片午夜精品| 精品一区二区三卡| 制服丝袜香蕉在线| 有码 亚洲区| 一区二区三区乱码不卡18| 亚洲精品国产av蜜桃| 中国美白少妇内射xxxbb| 伊人亚洲综合成人网| 免费不卡的大黄色大毛片视频在线观看| 日日爽夜夜爽网站| 熟妇人妻不卡中文字幕| 看免费成人av毛片| 亚洲天堂av无毛| 久久热在线av| 超碰97精品在线观看| 久久毛片免费看一区二区三区| 亚洲av日韩在线播放| 七月丁香在线播放| 久久精品国产综合久久久 | av免费观看日本| 黑人欧美特级aaaaaa片| 久热久热在线精品观看| www.av在线官网国产| 国产精品偷伦视频观看了| 亚洲欧洲日产国产| 三级国产精品片| 精品一品国产午夜福利视频| 亚洲精品久久久久久婷婷小说| 欧美少妇被猛烈插入视频| 丝袜脚勾引网站| 国产免费又黄又爽又色| 国产福利在线免费观看视频| 亚洲精品国产色婷婷电影| 亚洲国产精品成人久久小说| 午夜免费男女啪啪视频观看| 男女边摸边吃奶| 国产男女内射视频| 丝袜脚勾引网站| 美女中出高潮动态图| 日本黄大片高清| 亚洲精品色激情综合| 国产成人午夜福利电影在线观看| 中文字幕人妻丝袜制服| 久久久久视频综合| 午夜福利影视在线免费观看| 国产精品久久久av美女十八| 一级片免费观看大全| 亚洲内射少妇av| 亚洲精品日本国产第一区| 国产高清三级在线| 母亲3免费完整高清在线观看 | 视频中文字幕在线观看| 精品99又大又爽又粗少妇毛片| 不卡视频在线观看欧美| 99香蕉大伊视频| 高清不卡的av网站| 久久精品国产a三级三级三级| 国产在线视频一区二区| 色婷婷av一区二区三区视频| 18禁在线无遮挡免费观看视频| 午夜精品国产一区二区电影| 国产在线免费精品| 日韩大片免费观看网站| a级毛片黄视频| 亚洲国产毛片av蜜桃av| 视频在线观看一区二区三区| 人人澡人人妻人| 美女视频免费永久观看网站| 免费高清在线观看日韩| av又黄又爽大尺度在线免费看| 国产精品久久久久成人av| 国产精品久久久av美女十八| 久久久亚洲精品成人影院| 男女国产视频网站| 女人被躁到高潮嗷嗷叫费观| 国产av精品麻豆| 18禁裸乳无遮挡动漫免费视频| 亚洲高清免费不卡视频| 亚洲丝袜综合中文字幕| 蜜臀久久99精品久久宅男| 中文字幕最新亚洲高清| 欧美成人精品欧美一级黄| 亚洲国产日韩一区二区| 18禁裸乳无遮挡动漫免费视频| 欧美国产精品va在线观看不卡| 欧美97在线视频| 亚洲成人av在线免费| 成人国语在线视频| 秋霞在线观看毛片| 看免费成人av毛片| 自拍欧美九色日韩亚洲蝌蚪91| 精品一区二区三区视频在线| 99国产精品免费福利视频| 亚洲天堂av无毛| 国产1区2区3区精品| 99久久人妻综合| 香蕉丝袜av| 国产男人的电影天堂91| 久久久精品区二区三区| 国产免费现黄频在线看| 有码 亚洲区| 欧美日韩综合久久久久久| 亚洲国产av新网站| 大片免费播放器 马上看| 日本91视频免费播放| 在线看a的网站| 有码 亚洲区| 欧美成人午夜精品| 久久鲁丝午夜福利片| 少妇人妻 视频| 99九九在线精品视频| 97精品久久久久久久久久精品| 欧美日韩av久久| 天堂中文最新版在线下载| 90打野战视频偷拍视频| 捣出白浆h1v1| 久久久久久久久久久免费av| 满18在线观看网站| 日韩熟女老妇一区二区性免费视频| videosex国产| av在线老鸭窝| 午夜视频国产福利| 亚洲成色77777| 亚洲一区二区三区欧美精品| 两性夫妻黄色片 | 国精品久久久久久国模美| 成人无遮挡网站| 一二三四在线观看免费中文在 | 久久免费观看电影| 亚洲内射少妇av| 午夜福利视频精品| 免费人成在线观看视频色| 一级片免费观看大全| 久久午夜综合久久蜜桃| 乱人伦中国视频| 国产成人精品久久久久久| 2018国产大陆天天弄谢| 国产精品久久久久久久久免| 国产精品久久久久久久电影| 国产片内射在线| 侵犯人妻中文字幕一二三四区| 国产成人午夜福利电影在线观看| 日本av免费视频播放| 热re99久久国产66热| 国产色婷婷99| 波多野结衣一区麻豆| 日日爽夜夜爽网站| 久久久精品免费免费高清| 国产精品无大码| 视频中文字幕在线观看| 伊人久久国产一区二区| 丝袜人妻中文字幕| 亚洲综合精品二区| 欧美+日韩+精品| 久久精品久久久久久噜噜老黄| 我要看黄色一级片免费的| 午夜免费观看性视频| 国产熟女午夜一区二区三区| 久久国内精品自在自线图片| 国产成人精品一,二区| 亚洲av日韩在线播放| 丰满迷人的少妇在线观看| 熟妇人妻不卡中文字幕| 丝袜人妻中文字幕| 国产亚洲精品久久久com| 久久久欧美国产精品| 国产极品天堂在线| 亚洲精品国产av成人精品| 久久毛片免费看一区二区三区| 欧美日韩成人在线一区二区| 亚洲av福利一区| 一级片免费观看大全| 成年动漫av网址| 久久精品国产a三级三级三级| 亚洲人成网站在线观看播放| 中文字幕人妻丝袜制服| 99香蕉大伊视频| 看非洲黑人一级黄片| 午夜福利网站1000一区二区三区| 午夜福利网站1000一区二区三区| 久久这里只有精品19| 久久人人爽av亚洲精品天堂| 亚洲成人手机| 国产日韩欧美视频二区| 国产日韩欧美在线精品| 不卡视频在线观看欧美| 高清黄色对白视频在线免费看| 亚洲av成人精品一二三区| 久久久精品94久久精品| 国产老妇伦熟女老妇高清| 久久久久久久大尺度免费视频| 久久亚洲国产成人精品v| 伊人亚洲综合成人网| 国产一区亚洲一区在线观看| 日韩免费高清中文字幕av| 黑人巨大精品欧美一区二区蜜桃 | 大码成人一级视频| 女性被躁到高潮视频| 中文字幕亚洲精品专区| 久久久久久久久久久久大奶| 最后的刺客免费高清国语| 少妇的逼水好多| 一区二区av电影网| 亚洲精品中文字幕在线视频| 美女大奶头黄色视频| 免费观看无遮挡的男女| 免费在线观看完整版高清| 国产成人91sexporn| 一边摸一边做爽爽视频免费| 男女免费视频国产| 99久久综合免费| kizo精华| 少妇的逼好多水| 18禁观看日本| 午夜av观看不卡| 亚洲精品日本国产第一区| 久久鲁丝午夜福利片| 国产亚洲欧美精品永久| 高清在线视频一区二区三区| 中文字幕人妻丝袜制服| 亚洲美女黄色视频免费看| 亚洲丝袜综合中文字幕| 亚洲成av片中文字幕在线观看 | 国产一区二区激情短视频 | 午夜福利网站1000一区二区三区| 午夜久久久在线观看| 丰满少妇做爰视频| 亚洲精品第二区| 男女无遮挡免费网站观看| 欧美日韩一区二区视频在线观看视频在线| 亚洲综合色惰| 欧美精品av麻豆av| 我要看黄色一级片免费的| 中文字幕制服av| 国产成人一区二区在线| 久久久久久久久久久久大奶| 久久久国产一区二区| 9热在线视频观看99| 午夜影院在线不卡| 插逼视频在线观看| 各种免费的搞黄视频| 9191精品国产免费久久| 欧美 日韩 精品 国产| 热99久久久久精品小说推荐| 日韩av在线免费看完整版不卡| 欧美国产精品一级二级三级| a级毛片黄视频| 久热这里只有精品99| 国产精品国产三级专区第一集| 色视频在线一区二区三区| 日本欧美视频一区| 最近最新中文字幕免费大全7| 欧美xxxx性猛交bbbb| 波野结衣二区三区在线| 黄色毛片三级朝国网站| 亚洲精品美女久久久久99蜜臀 | 亚洲成av片中文字幕在线观看 | 日日摸夜夜添夜夜爱| 午夜福利乱码中文字幕| 啦啦啦在线观看免费高清www| 久久久久久人妻| 丝袜人妻中文字幕| 水蜜桃什么品种好| 亚洲精品aⅴ在线观看| 一二三四中文在线观看免费高清| 久久av网站| 又粗又硬又长又爽又黄的视频| 男人舔女人的私密视频| 国产熟女午夜一区二区三区| 亚洲成人一二三区av| 国产永久视频网站| 国产精品久久久久久精品电影小说| 人人妻人人澡人人爽人人夜夜| 中国国产av一级| 国产色婷婷99| 日韩中文字幕视频在线看片| 一个人免费看片子| 日本wwww免费看| 在线观看www视频免费| 成人毛片60女人毛片免费| 99九九在线精品视频| 不卡视频在线观看欧美| 综合色丁香网| 亚洲久久久国产精品| 日本猛色少妇xxxxx猛交久久| 男女午夜视频在线观看 | 午夜影院在线不卡| 各种免费的搞黄视频| 高清毛片免费看| 国产精品成人在线| 97精品久久久久久久久久精品| 天堂中文最新版在线下载| 国产精品久久久久久av不卡| 在线精品无人区一区二区三| av一本久久久久| 日韩视频在线欧美| 国产精品人妻久久久影院| 熟妇人妻不卡中文字幕| 国产国语露脸激情在线看| 超碰97精品在线观看| 亚洲精品日韩在线中文字幕| 最近的中文字幕免费完整| 肉色欧美久久久久久久蜜桃| 一本大道久久a久久精品| 波野结衣二区三区在线| 精品卡一卡二卡四卡免费| 精品少妇黑人巨大在线播放| 亚洲精品日本国产第一区| 在线天堂中文资源库| 国产片内射在线| 免费人成在线观看视频色| 亚洲成国产人片在线观看| 日韩制服丝袜自拍偷拍| 亚洲欧美成人综合另类久久久| 欧美亚洲日本最大视频资源| 国产精品无大码| 少妇熟女欧美另类| 日产精品乱码卡一卡2卡三| a级毛片黄视频| 午夜福利在线观看免费完整高清在| av国产久精品久网站免费入址| 成年人免费黄色播放视频| 国产片内射在线| 久久这里有精品视频免费| 又粗又硬又长又爽又黄的视频| 丝袜喷水一区| av免费在线看不卡| 久久久国产欧美日韩av| 日韩三级伦理在线观看| 久久国产精品男人的天堂亚洲 | 蜜桃在线观看..| 国产精品免费大片| 80岁老熟妇乱子伦牲交| 狂野欧美激情性xxxx在线观看| 亚洲国产看品久久| 人人妻人人澡人人看| 免费日韩欧美在线观看| 人人妻人人爽人人添夜夜欢视频| 男人操女人黄网站| 黄色毛片三级朝国网站| 巨乳人妻的诱惑在线观看| 欧美人与性动交α欧美精品济南到 | 国产无遮挡羞羞视频在线观看| 日本-黄色视频高清免费观看| 精品一区在线观看国产| 涩涩av久久男人的天堂| 蜜臀久久99精品久久宅男| 高清不卡的av网站| 91成人精品电影| 亚洲精品日韩在线中文字幕| 黄色毛片三级朝国网站| 国产又爽黄色视频| 一区二区三区精品91| 国产免费一区二区三区四区乱码| 成人综合一区亚洲| 2022亚洲国产成人精品| 亚洲国产最新在线播放| 国产精品久久久av美女十八| 超色免费av| 99九九在线精品视频| 赤兔流量卡办理| 欧美人与性动交α欧美精品济南到 | 亚洲精品aⅴ在线观看| 青青草视频在线视频观看| 在线观看人妻少妇| 十分钟在线观看高清视频www| 精品人妻一区二区三区麻豆| 爱豆传媒免费全集在线观看| 成年女人在线观看亚洲视频| 多毛熟女@视频| 考比视频在线观看| 国产精品女同一区二区软件| 国产成人一区二区在线| 亚洲,一卡二卡三卡| 国产欧美亚洲国产| 午夜91福利影院| 99精国产麻豆久久婷婷| 久久国内精品自在自线图片| 欧美精品人与动牲交sv欧美| 一级黄片播放器| 亚洲激情五月婷婷啪啪| 黄网站色视频无遮挡免费观看| 超色免费av| 另类精品久久| 国产极品粉嫩免费观看在线| 国产乱人偷精品视频| av免费在线看不卡| 精品视频人人做人人爽| 久久久精品94久久精品| 亚洲,欧美,日韩| 色视频在线一区二区三区| 欧美日本中文国产一区发布| 99久久综合免费| 精品国产一区二区三区久久久樱花| 亚洲综合精品二区| 交换朋友夫妻互换小说| 一边摸一边做爽爽视频免费| 哪个播放器可以免费观看大片| 国产精品蜜桃在线观看| 最新中文字幕久久久久| 国产乱人偷精品视频| 国产精品国产三级国产av玫瑰| tube8黄色片| 国产成人精品婷婷| 人妻 亚洲 视频| 91久久精品国产一区二区三区| 久久亚洲国产成人精品v| 综合色丁香网| 性高湖久久久久久久久免费观看| 一本大道久久a久久精品| 男女下面插进去视频免费观看 | 精品国产乱码久久久久久小说| 内地一区二区视频在线| 18禁裸乳无遮挡动漫免费视频| 亚洲三级黄色毛片| 成人影院久久| 亚洲人成网站在线观看播放| 爱豆传媒免费全集在线观看| 日韩欧美精品免费久久| 亚洲精品久久午夜乱码| 极品少妇高潮喷水抽搐| 国产高清不卡午夜福利| 两个人看的免费小视频| 在线观看人妻少妇| 2022亚洲国产成人精品| 毛片一级片免费看久久久久| av视频免费观看在线观看| 少妇人妻久久综合中文| 妹子高潮喷水视频| 国产精品 国内视频| 中文字幕精品免费在线观看视频 | 午夜激情久久久久久久| 久久精品国产鲁丝片午夜精品| 69精品国产乱码久久久| 日韩一区二区三区影片| 亚洲av欧美aⅴ国产| 国产精品久久久久成人av| 人人妻人人澡人人爽人人夜夜| 久久久久久久精品精品| 亚洲精华国产精华液的使用体验| 一级片'在线观看视频| 久久久国产欧美日韩av| 亚洲欧美成人综合另类久久久| 美女视频免费永久观看网站| 日本与韩国留学比较| 最近的中文字幕免费完整| 天堂中文最新版在线下载| 美女内射精品一级片tv| 最近最新中文字幕大全免费视频 | 欧美日韩国产mv在线观看视频| 久久久久久久久久人人人人人人| 亚洲国产色片| 香蕉国产在线看| 欧美成人午夜免费资源| 啦啦啦在线观看免费高清www| 国产又色又爽无遮挡免| 久久国产亚洲av麻豆专区| 你懂的网址亚洲精品在线观看| 亚洲国产色片| 国产 一区精品| 亚洲精品一区蜜桃| 一级a做视频免费观看| 欧美激情国产日韩精品一区| 国语对白做爰xxxⅹ性视频网站| 中文天堂在线官网| 黄色 视频免费看| 久久精品国产自在天天线| videos熟女内射| 国国产精品蜜臀av免费| 日韩,欧美,国产一区二区三区| 最后的刺客免费高清国语| 亚洲av欧美aⅴ国产| 亚洲av综合色区一区| 亚洲国产av影院在线观看| 亚洲国产看品久久| 国产熟女午夜一区二区三区| 久久精品国产亚洲av天美| 国产成人精品婷婷| 久久久久久久久久成人| av国产久精品久网站免费入址| 交换朋友夫妻互换小说| 亚洲国产欧美日韩在线播放| 日韩三级伦理在线观看| 伦精品一区二区三区| 男男h啪啪无遮挡| 纵有疾风起免费观看全集完整版| 亚洲精品456在线播放app| 亚洲av.av天堂| 赤兔流量卡办理| 久久 成人 亚洲| 高清黄色对白视频在线免费看| 妹子高潮喷水视频| 亚洲精品乱久久久久久| 久久av网站| 成人无遮挡网站| 亚洲精品乱码久久久久久按摩| 在线亚洲精品国产二区图片欧美| 欧美精品国产亚洲| 国产一区二区三区综合在线观看 | 免费高清在线观看视频在线观看| av线在线观看网站| 97在线视频观看| 大香蕉久久网| 一本久久精品| 国产黄频视频在线观看| 国产亚洲av片在线观看秒播厂| www.熟女人妻精品国产 | 中文乱码字字幕精品一区二区三区| av播播在线观看一区| 国产亚洲精品第一综合不卡 | 精品一品国产午夜福利视频| 国国产精品蜜臀av免费| 国产精品女同一区二区软件| 三上悠亚av全集在线观看| 国产 一区精品| 妹子高潮喷水视频| 哪个播放器可以免费观看大片| 少妇人妻精品综合一区二区| 纵有疾风起免费观看全集完整版| 中文乱码字字幕精品一区二区三区| 成年美女黄网站色视频大全免费| 美国免费a级毛片| 成人亚洲欧美一区二区av| 好男人视频免费观看在线| 精品国产露脸久久av麻豆| 亚洲欧美清纯卡通| 99视频精品全部免费 在线| 丝袜美足系列| 久久久国产精品麻豆| 91精品伊人久久大香线蕉| 精品人妻熟女毛片av久久网站| 久久久久久久精品精品| 亚洲精品美女久久av网站| 免费观看无遮挡的男女| 人人妻人人澡人人爽人人夜夜| av国产久精品久网站免费入址| 国产熟女欧美一区二区| 99久国产av精品国产电影| 国产免费一级a男人的天堂| 97人妻天天添夜夜摸| 午夜日本视频在线| 国产亚洲一区二区精品| 免费观看a级毛片全部| 中文字幕免费在线视频6| 少妇人妻久久综合中文| 一个人免费看片子| 色哟哟·www| 久久亚洲国产成人精品v| 日本av手机在线免费观看| 性高湖久久久久久久久免费观看| 成人无遮挡网站| 亚洲第一区二区三区不卡| 蜜桃国产av成人99| 自拍欧美九色日韩亚洲蝌蚪91| 不卡视频在线观看欧美| 亚洲欧美中文字幕日韩二区| 精品国产一区二区久久| 欧美老熟妇乱子伦牲交| 国产精品久久久久成人av| av免费观看日本| 秋霞伦理黄片| 久久精品国产综合久久久 | 51国产日韩欧美| 搡女人真爽免费视频火全软件| 精品人妻在线不人妻| 18禁在线无遮挡免费观看视频| 亚洲精品久久久久久婷婷小说| 91成人精品电影| 日本爱情动作片www.在线观看| 永久免费av网站大全| 久久久久久久国产电影| 亚洲欧洲精品一区二区精品久久久 | 2021少妇久久久久久久久久久| 黄片无遮挡物在线观看| 男女无遮挡免费网站观看| av卡一久久| 一区二区三区精品91| 国产又爽黄色视频| 美女福利国产在线| 只有这里有精品99| 成人综合一区亚洲| av卡一久久| 日本与韩国留学比较| 日韩视频在线欧美| 亚洲成人手机| 欧美精品一区二区免费开放| 久久精品熟女亚洲av麻豆精品| 大片免费播放器 马上看| 91国产中文字幕| 男女无遮挡免费网站观看| 亚洲婷婷狠狠爱综合网| 一区二区三区精品91| 精品国产一区二区三区四区第35| 国产成人a∨麻豆精品| 大香蕉久久成人网| 国产淫语在线视频| 亚洲成人手机| 黄片无遮挡物在线观看| 亚洲精品久久成人aⅴ小说| 日韩欧美精品免费久久| 亚洲天堂av无毛| 国产午夜精品一二区理论片| 欧美人与性动交α欧美软件 | 欧美97在线视频| 美女内射精品一级片tv| 国产一区二区三区av在线| 久热这里只有精品99| 春色校园在线视频观看| 中文字幕人妻熟女乱码| 亚洲精品一二三|