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

    基于隔代映射算子的差分進(jìn)化算法

    2016-06-27 04:22:25符純明陳光宋
    中國機(jī)械工程 2016年11期

    符純明 姜 潮 陳光宋 吉 磊

    1.湖南大學(xué)汽車車身先進(jìn)設(shè)計(jì)制造國家重點(diǎn)實(shí)驗(yàn)室,長沙,4100822.南京理工大學(xué),南京,210094

    基于隔代映射算子的差分進(jìn)化算法

    符純明1姜潮1陳光宋2吉磊2

    1.湖南大學(xué)汽車車身先進(jìn)設(shè)計(jì)制造國家重點(diǎn)實(shí)驗(yàn)室,長沙,4100822.南京理工大學(xué),南京,210094

    摘要:提出一種基于隔代映射算子的差分進(jìn)化算法以求解優(yōu)化問題,該方法在保證解的精度的同時(shí)具有較快的收斂速度。在經(jīng)典的差分進(jìn)化算法基礎(chǔ)上,采用反向?qū)W習(xí)策略產(chǎn)生初始種群,并采用兩種差分變異策略產(chǎn)生變異個(gè)體,以增加種群的多樣性;利用隔代映射算子產(chǎn)生三個(gè)新個(gè)體替換當(dāng)前進(jìn)化種群中最差的三個(gè)個(gè)體,以實(shí)現(xiàn)精英策略提升算法的收斂性;為了保持種群的多樣性和避免獲得局部解,利用探測(cè)算子策略產(chǎn)生新個(gè)體加入進(jìn)化種群。采用11個(gè)單峰、多峰測(cè)試函數(shù)和兩個(gè)工程實(shí)例驗(yàn)證了該方法的有效性。

    關(guān)鍵詞:差分進(jìn)化算法; 隔代映射算子; 反向?qū)W習(xí); 探測(cè)算子

    0引言

    在實(shí)際工程問題中,經(jīng)常遇到一類較為復(fù)雜的優(yōu)化問題,其目標(biāo)函數(shù)不連續(xù)或不可導(dǎo),導(dǎo)致傳統(tǒng)基于梯度的優(yōu)化方法難以有效求解。而基于種群的進(jìn)化算法不需要目標(biāo)函數(shù)的梯度信息,只需采用適應(yīng)度函數(shù)即可對(duì)復(fù)雜目標(biāo)函數(shù)進(jìn)行優(yōu)化求解,為求解該類問題提供了一種有效的途徑。差分進(jìn)化(differential evolution,DE)算法是進(jìn)化算法(evolutionary algorithms,EAs)中一種簡單有效的隨機(jī)搜索算法,最早由Storn等[1]于1995年提出,并用于求解切比雪夫多項(xiàng)式擬合問題。因DE操作簡單,控制參數(shù)少,具有較強(qiáng)的全局優(yōu)化性和魯棒性,故被廣泛應(yīng)用于機(jī)械優(yōu)化設(shè)計(jì)、機(jī)器人控制、車間節(jié)能調(diào)度等領(lǐng)域[2-4]。

    DE算法的性能主要取決于試驗(yàn)向量產(chǎn)生策略(交叉和變異算子)和控制參數(shù)(種群大小、縮放因子、交叉概率)。在改進(jìn)試驗(yàn)向量產(chǎn)生策略和控制參數(shù)方面,國內(nèi)外學(xué)者進(jìn)行了一系列研究和探索。如Fan等[5]提出了一種三角變異算子用于加速DE算法收斂性能,并同時(shí)引入一參數(shù)控制三角變異算子的使用頻率以保持種群的多樣性。Zhang等[6]利用當(dāng)前群體中最好解和較好解的信息產(chǎn)生變異個(gè)體加入進(jìn)化種群,并將已淘汰的個(gè)體儲(chǔ)存于一外部種群中參與隨后的進(jìn)化,其縮放因子和交叉概率分別通過標(biāo)準(zhǔn)正態(tài)分布和柯西分布產(chǎn)生。但該方法需要較大的儲(chǔ)存空間保存淘汰的個(gè)體。Das等[7]提出了一種新型的局部搜索模型,采用目標(biāo)向量鄰域中的最優(yōu)解產(chǎn)生變異個(gè)體,該方法動(dòng)態(tài)地利用臨近個(gè)體產(chǎn)生更優(yōu)良個(gè)體。Wang等[8]針對(duì)當(dāng)前DE交叉算子普遍存在的缺陷,提出了一種正交交叉策略,有效地提升了DE算法搜索性能。Ronkkonen等[9]采用大量測(cè)試函數(shù)分析了初始種群大小和縮放因子對(duì)傳統(tǒng)的DE算法性能的影響。

    通常,DE算法求解優(yōu)化問題時(shí),對(duì)控制參數(shù)較為敏感,找到一組合理的控制參數(shù)通常需要反復(fù)調(diào)試,因此,國內(nèi)外學(xué)者發(fā)展了適應(yīng)性和自適應(yīng)性策略調(diào)節(jié)其控制參數(shù)[10-11]。適應(yīng)性策略是指搜索過程中不斷利用反饋信息來調(diào)節(jié)參數(shù),而自適應(yīng)性策略是指將控制參數(shù)編碼于個(gè)體參與進(jìn)化過程。Zhu等[10]提出了一種動(dòng)態(tài)調(diào)整種群大小的策略,可有效地篩除進(jìn)化種群中冗余個(gè)體從而減少計(jì)算量。Brest等[11]提出了一種新型自適應(yīng)算法,將縮放因子和交叉概率進(jìn)行編碼并參與進(jìn)化,有利于增加種群的多樣性。Wang等[12]建立了試驗(yàn)向量產(chǎn)生策略知識(shí)庫和控制參數(shù)設(shè)置知識(shí)庫,并提出組合差分進(jìn)化(composite DE,CoDE)算法來求解優(yōu)化問題,有效地利用多組策略和參數(shù)產(chǎn)生不同的變異個(gè)體,保證了該算法求解不同問題時(shí)具有優(yōu)良的性能。車林仙等[13]提出了一種自適應(yīng)逃逸差分進(jìn)化算法,并用于求解并聯(lián)機(jī)器人主驅(qū)動(dòng)機(jī)構(gòu)位置的優(yōu)化問題。楊曉明等[14]利用改進(jìn)的差分進(jìn)化算法優(yōu)化設(shè)計(jì)了盤式制動(dòng)器。王前等[15]提出了一種自適應(yīng)差分進(jìn)化方法,并將其應(yīng)用于輪胎模型的參數(shù)辨識(shí)中。

    一系列改進(jìn)的DE算法被應(yīng)用,獲得了較好的效果,但求解某些優(yōu)化問題時(shí)這些算法存在收斂速度慢且易于陷入局部最優(yōu)的缺陷。為此,本文將連續(xù)兩代中的最優(yōu)個(gè)體用于構(gòu)造搜索方向以便快速有效地獲得更優(yōu)個(gè)體,提出一種基于隔代映射算子的差分進(jìn)化(intergeneration projection DE,IPDE)算法,相比現(xiàn)有方法,該方法在獲得全局最優(yōu)解的同時(shí)具有較快的收斂速度。

    1差分進(jìn)化算法

    通常的優(yōu)化問題可定義為

    (1)

    式中,f(x)為目標(biāo)函數(shù);xk、xkL和xkU分別為設(shè)計(jì)變量、設(shè)計(jì)變量的下界和上界;n為設(shè)計(jì)變量的個(gè)數(shù)。

    DE算法[1]的基本思想是:對(duì)當(dāng)前種群中的個(gè)體進(jìn)行變異和交叉操作產(chǎn)生新個(gè)體,然后采用貪婪準(zhǔn)則在兩個(gè)種群中選出較優(yōu)個(gè)體組成下一代進(jìn)化種群,不斷迭代直至滿足收斂條件,其流程如圖1所示。

    圖1 DE算法流程圖

    首先,DE算法在可行域空間隨機(jī)產(chǎn)生N個(gè)個(gè)體組成初始種群[1]:

    (2)

    i=1,2,…,N;j=1,2,…,n

    (3)

    (4)

    DE/rand/1策略收斂慢但具有較強(qiáng)的全局搜索能力,利于保持種群的多樣性。DE/best/1策略則具有較快的收斂速度,但容易陷入局部最優(yōu)。

    再次,采用二項(xiàng)式交叉算子產(chǎn)生試驗(yàn)向量ui,其中

    (5)

    式中,Cr為交叉概率;prj為第j維向量在[0,1]中隨機(jī)產(chǎn)生的一實(shí)數(shù)。

    最后,在目標(biāo)向量和試驗(yàn)向量中選擇較優(yōu)個(gè)體保留至第t+1代:

    (6)

    上述操作不斷迭代直至滿足收斂條件。

    2IPDE算法及構(gòu)造

    傳統(tǒng)DE算法求解某些較為復(fù)雜優(yōu)化問題時(shí),其后期收斂效率低或易于陷入局部最優(yōu)。為克服以上缺點(diǎn),本文提出的IPDE算法采用反向?qū)W習(xí)策略[16]產(chǎn)生初始種群,并利用進(jìn)化種群中連續(xù)兩代的最優(yōu)個(gè)體來構(gòu)造搜索方向,以便快速有效地獲得更優(yōu)個(gè)體,即隔代映射策略。同時(shí)為了保持種群的多樣性,每隔T代由當(dāng)前種群中的最優(yōu)個(gè)體和最差個(gè)體產(chǎn)生兩個(gè)新個(gè)體,擇優(yōu)選擇一個(gè)個(gè)體進(jìn)入下一代進(jìn)化種群。為此,IPDE算法求解優(yōu)化問題時(shí),在保證解的精度的同時(shí)并具有快速收斂的特點(diǎn)。

    2.1隔代映射策略

    基于種群的進(jìn)化算法中,其核心之一是如何快速有效地產(chǎn)生優(yōu)良個(gè)體,從而加速算法收斂性能。本文引入一種隔代映射策略[17],以快速收斂到局部最優(yōu)解,從而加速算法收斂性能。如圖2所示,隔代映射策略采用連續(xù)兩代種群中的最優(yōu)個(gè)體來構(gòu)造搜索方向,從而在兩代最優(yōu)個(gè)體的搜索方向上產(chǎn)生三個(gè)新個(gè)體。因新個(gè)體在連續(xù)兩代的最優(yōu)個(gè)體的方向上產(chǎn)生,有助于快速生成接近局部最優(yōu)解的個(gè)體。

    圖2 隔代映射策略示意圖(2變量)

    (7)

    2.2算法構(gòu)造

    本文提出的IPDE算法同時(shí)結(jié)合反向?qū)W習(xí)[16]和隔代映射策略的優(yōu)點(diǎn)求解優(yōu)化問題,在保證精度的同時(shí)具有快速收斂的特點(diǎn),其迭代過程如下:

    (1)在搜索空間隨機(jī)產(chǎn)生N個(gè)個(gè)體,并利用反向?qū)W習(xí)策略產(chǎn)生N個(gè)個(gè)體便于均勻分布于搜索空間,最后在2N個(gè)個(gè)體中根據(jù)目標(biāo)函數(shù)值選出較優(yōu)的N個(gè)個(gè)體構(gòu)成初始種群。

    (2)隨機(jī)執(zhí)行DE/best/1或DE/rand/1差分變異策略產(chǎn)生變異個(gè)體。首先,隨機(jī)產(chǎn)生一個(gè)[0, 1]之間的數(shù),當(dāng)該隨機(jī)數(shù)小于0.5時(shí),執(zhí)行DE/rand/1策略,否則執(zhí)行DE/best/1策略;其次,采用二項(xiàng)式交叉策略生成N個(gè)試驗(yàn)向量個(gè)體;最后,在試驗(yàn)向量和目標(biāo)向量中選出N個(gè)個(gè)體進(jìn)入下一代。

    (3)當(dāng)代數(shù)t>1時(shí)執(zhí)行隔代映射算子策略[17]產(chǎn)生三個(gè)新個(gè)體,替換進(jìn)化種群中最差的三個(gè)個(gè)體。

    (4)判斷當(dāng)前進(jìn)化代數(shù)t是否為參數(shù)T的整數(shù)倍,如果是則執(zhí)行探測(cè)算子策略產(chǎn)生兩個(gè)體,擇優(yōu)選擇一個(gè)個(gè)體加入下一代種群,以增加進(jìn)化種群的多樣性,否則不執(zhí)行。

    (5)判斷目標(biāo)函數(shù)調(diào)用次數(shù)是否大于最大目標(biāo)函數(shù)調(diào)用次數(shù),是則輸出當(dāng)前最優(yōu)解,否則轉(zhuǎn)到步驟(2)。

    其算法流程如圖3所示。

    圖3 IPDE算法流程

    在步驟(1)中,反向?qū)W習(xí)策略使初始種群能更好地分布于搜索空間,便于提高算法的全局優(yōu)化性能。步驟(2)采用隨機(jī)策略較好地避免了單獨(dú)采用DE/rand/1策略收斂速度慢的缺點(diǎn),同時(shí)也避免了單獨(dú)采用DE/best/1策略收斂速度較快可能導(dǎo)致局部收斂的缺陷。在步驟(3)中,采用隔代映射算子較好地利用了當(dāng)前代和上一代種群中的最優(yōu)個(gè)體信息,有利于快速有效地指導(dǎo)進(jìn)化種群產(chǎn)生較優(yōu)個(gè)體。

    3測(cè)試函數(shù)及工程應(yīng)用

    3.1測(cè)試函數(shù)及參數(shù)設(shè)置

    測(cè)試函數(shù)采用2005年IEEE進(jìn)化計(jì)算大會(huì)公布的測(cè)試函數(shù)中的部分測(cè)試函數(shù)[19]:單峰函數(shù)F1~F5和多峰函數(shù)F6、F8~F12。在11個(gè)測(cè)試函數(shù)中,設(shè)計(jì)變量n為30。IPDE算法的參數(shù)設(shè)置為:N取70,F(xiàn)取0.8,Cr取0.9,執(zhí)行探測(cè)算子的代數(shù)T取20,隔代映射算子中參數(shù)γ=s=λ=0.6。為了便于對(duì)比各算法的魯棒性,IPDE算法對(duì)每個(gè)測(cè)試函數(shù)獨(dú)立運(yùn)算25次,統(tǒng)計(jì)其函數(shù)誤差值(f(x′)-f(x*))的平均值和標(biāo)準(zhǔn)差,其中x*為測(cè)試函數(shù)已知的最優(yōu)解。當(dāng)算法最大函數(shù)評(píng)價(jià)次數(shù)達(dá)到300 000次后,算法迭代結(jié)束并輸出最優(yōu)解x′。

    3.2算法性能比較

    3.2.1與兩種改進(jìn)的DE算法比較

    基于差分進(jìn)化模式的SaDE(DEwithstrategyadaptation)算法[20]和EPSDE(ensembleofmutationstrategiesandcontrolparameterswithDE)算法[21]求解優(yōu)化問題時(shí),具有較優(yōu)異的綜合性能。SaDE算法[20]利用先前搜索過程中積累的經(jīng)驗(yàn),以自適應(yīng)的方式同時(shí)選擇試驗(yàn)向量產(chǎn)生策略和控制參數(shù)。EPSDE算法[21]則采用三種試驗(yàn)向量產(chǎn)生策略和一組參數(shù)的隨機(jī)組合產(chǎn)生變異個(gè)體。為此,將IPDE算法與SaDE算法和EPSDE算法進(jìn)行比較,其中,SaDE算法和EPSDE算法求解測(cè)試函數(shù)的仿真實(shí)驗(yàn)統(tǒng)計(jì)數(shù)據(jù)分別通過參考文獻(xiàn)[20-21]獲得,具體統(tǒng)計(jì)結(jié)果如表1所示。表1中,“-”,“+”和“≈”分別表示相應(yīng)算法統(tǒng)計(jì)性能差于、優(yōu)于和相似于IPDE。m和s分別表示25次獨(dú)立運(yùn)行的函數(shù)誤差值的平均值和標(biāo)準(zhǔn)差。

    表1 IPDE與SaDE、EPSDE、CLPSO和GL-25在11個(gè)測(cè)試函數(shù)上的統(tǒng)計(jì)結(jié)果

    (a)F1收斂曲線

    (b)F2收斂曲線

    (c)F3收斂曲線

    (d)F4收斂曲線圖4 F1、F2、F3和F4收斂曲線

    從表1中得出,IPDE算法相比于SaDE算法和EPSDE算法在整體上具有較優(yōu)異的全局優(yōu)化性能。具體地,IPDE算法分別在5個(gè)和6個(gè)測(cè)試函數(shù)上性能優(yōu)于SaDE算法和EPSDE算法,在3個(gè)測(cè)試函數(shù)上性能相似于SaDE算法和EPSDE算法,分別在3個(gè)和2個(gè)函數(shù)上劣于SaDE算法和EPSDE算法。IPDE算法、SaDE算法和EPSDE算法求解11個(gè)測(cè)試函數(shù)收斂曲線如圖4~圖6所示。由圖4~圖6可以看出,IPDE算法在整體上一致收斂到全局最優(yōu)解,其收斂性能在函數(shù)F2~F8、F10和F11上顯著優(yōu)于SaDE算法,在F2~F4、F8、F11和F12上顯著優(yōu)于EPSDE算法。3.2.2與兩種非DE模式的進(jìn)化算法比較

    為進(jìn)一步驗(yàn)證該算法的有效性,采用非差分進(jìn)化模式的CLPSO(comprehensive learning particle swarm optimizer)算法[22]和GL-25(global and local real-coded genetic algorithms)算法[23]進(jìn)行對(duì)比分析。CLPSO算法[22]以粒子群算法為基礎(chǔ),其核心為一個(gè)粒子采用其他所有粒子的最好信息更新其速度。GL-25算法[23]是一類基于混合實(shí)數(shù)編碼的遺傳算法。兩種算法求解測(cè)試函數(shù)的仿真結(jié)果統(tǒng)計(jì)于表1。從表1中得出,IPDE算法在整體上顯著優(yōu)于CLPSO算法和GL-25算法。就11個(gè)測(cè)試函數(shù)的統(tǒng)計(jì)特性而言,IPDE算法分別在8個(gè)和10個(gè)測(cè)試函數(shù)上優(yōu)于CLPSO算法和GL-25算法。與IPDE算法相比,CLPSO算法僅在2個(gè)測(cè)試函數(shù)(F1和F9)上占優(yōu),而GL-25算法在11個(gè)測(cè)試函數(shù)上均不能占優(yōu)。

    (a)F5收斂曲線

    (b)F6收斂曲線

    (c)F8收斂曲線

    (d)F9收斂曲線圖5 F5、F6、F8和F9收斂曲線

    由上述5種算法求解11個(gè)測(cè)試函數(shù)的仿真結(jié)果和收斂曲線得出:相比于對(duì)比算法,IPDE算法在獲得最優(yōu)解的同時(shí)具有較好的全局優(yōu)化和快速收斂性能。

    按照13號(hào)線停站方案,運(yùn)行圖鋪畫結(jié)果顯示(快車-慢車-慢車-慢車-快車的發(fā)車間隔為3 min-3 min-2 min-2 min),遠(yuǎn)期早高峰慢車需在白芒站、羅租站、同觀路站、東周路站和長春北站待避。

    3.3工程實(shí)例

    3.3.1火炮身管結(jié)構(gòu)優(yōu)化設(shè)計(jì)

    火炮的射擊精度在很大程度上取決于射彈的初始條件,而火炮發(fā)射時(shí)身管的振動(dòng)對(duì)射彈初始條件有較大的影響,因此,優(yōu)化設(shè)計(jì)身管的結(jié)構(gòu)來調(diào)節(jié)其固有頻率具有重要的意義。本文采用文獻(xiàn)[24]所簡化的身管模型,如圖7所示,其模型主要由6根Timoshenko梁組成,其中m1和m2為簡化的炮口制退器和炮尾的集中質(zhì)量,圖中兩支撐點(diǎn)為搖架的前后支撐點(diǎn),L為身管的總長,xi(i=1,2,…,5)表示各個(gè)等截面Timoshenko梁的長度。

    (a)F10收斂曲線

    (b)F11收斂曲線

    (c)F12收斂曲線圖6 F10、F11和F12收斂曲線

    圖7 簡化的火炮身管模型

    身管模型的一階固有頻率的高低可表征其剛度的大小,所以將最大化身管的一階固有頻率ω作為優(yōu)化目標(biāo)函數(shù),其設(shè)計(jì)變量為身管的各軸向長度xk(k=1,2,…,5),其數(shù)學(xué)模型表示為

    (8)

    某汽車車架結(jié)構(gòu)[25]如圖8所示,由2根縱梁和8根橫梁組成,需優(yōu)化其橫梁的布置使車架在y方向上具有最大剛度。bi(i=1,2,…,8)表示8根橫梁,xi(i=1,2,3,4)表示各橫梁之間的跨距。車架為汽車的基座,大多數(shù)關(guān)鍵部件如發(fā)動(dòng)機(jī)、懸架、油箱和駕駛室等都固定于車架,而這些部件通過連接件對(duì)車架產(chǎn)生載荷,通過簡化處理,獲得車架的靜力學(xué)模型如圖9所示。圖中Q1、Q2、Q3和Q4分別表示駕駛室、發(fā)動(dòng)機(jī)、油箱和載重的等價(jià)質(zhì)量作用于車架的載荷,分別取值為19 600N,39 200N,4410N,343 000N;三角形表示不同方向的固定約束。車架材料的密度為ρ=7.86×103kg/m3,彈性模量E=2.0×108Pa,泊松比ν=0.3。

    圖8 某車架簡化模型(mm)

    圖9 車架有限元模型及加載

    橫梁b1、b6、b7、b8固定,其他橫梁間的跨距需優(yōu)化。因?yàn)檐嚰茏冃卧趛向上的最大位移可表征其剛度的大小,故將其作為目標(biāo)函數(shù),建立如下的優(yōu)化問題:

    (9)

    式中,dmax為車架變形后在y向上的最大位移。

    建立車架的有限元模型,采用SHELL63單元?jiǎng)澐志W(wǎng)格。采用IPDE算法優(yōu)化,初始種群N為30,收斂準(zhǔn)則為連續(xù)兩代獲得最大位移誤差小于10-3mm,其他參數(shù)設(shè)置與求解測(cè)試函數(shù)的參數(shù)設(shè)置一致。該算法優(yōu)化求解獲得的最小最大位移為1.682mm,其對(duì)應(yīng)的最優(yōu)設(shè)計(jì)向量為(410mm, 1580mm, 1030mm, 791mm),降低了原模型的最大位移,使車架結(jié)構(gòu)設(shè)計(jì)更為合理。工程設(shè)計(jì)師可根據(jù)優(yōu)化結(jié)果結(jié)合工程經(jīng)驗(yàn)指導(dǎo)實(shí)際的工程設(shè)計(jì)。

    4結(jié)語

    提出了一種基于隔代映射算子的差分進(jìn)化算法來求解優(yōu)化問題,在保證解的精度的同時(shí)具有較快的收斂速度。11個(gè)單峰和多峰測(cè)試函數(shù)驗(yàn)證了該算法具有優(yōu)秀的全局優(yōu)化與快速收斂性。其統(tǒng)計(jì)性能指標(biāo)表明,IPDE算法相比于SaDE算法和EPSDE算法分別在5個(gè)和6個(gè)測(cè)試函數(shù)中優(yōu)于SaDE算法和EPSDE算法,并分別在8個(gè)和11個(gè)測(cè)試函數(shù)上顯著優(yōu)于非差分模式的CLPSO算法和GL-25算法。最后,算法被應(yīng)用于兩個(gè)實(shí)際工程問題的求解,獲得了比優(yōu)化前更好的設(shè)計(jì)方案。

    參考文獻(xiàn):

    [1]StornR,PriceK.DifferentialEvolution-aSimpleandEfficientHeuristicforGlobalOptimizationoverContinuousSpaces[J].JournalofGlobalOptimization, 1997, 11(4): 341-359.

    [2]Mezura-MontesE,CoelloCAC,VelZquez-ReyesJ,etal.MultipleTrialVectorsinDifferentialEvolutionforEngineeringDesign[J].EngineeringOptimization, 2007, 39(5): 567-589.

    [3]陳勇, 吳云翔, 王亞良, 等. 訂單不確定下雙資源約束多裝配線魯棒調(diào)度[J]. 中國機(jī)械工程, 2014, 25(12): 1567-1573.

    ChenYong,WuYunxiang,WangYaliang,etal.Multi-assemblyLineRobustSchedulingofDoubleResourceConstrainsunderUncertainOrders[J].JournalofMechanicalEngineering, 2014, 25(12): 1567-1573.

    [4]趙燕偉, 張立萍, 張景玲, 等. 加工裝配式流水車間節(jié)能調(diào)度建模與優(yōu)化[J]. 中國機(jī)械工程, 2014, 25(16): 2196-2203.

    ZhaoYanwei,ZhangLiping,ZhangJingling,etal.ModelingandOptimizationofProcess-assembly-typeFlow-shopSchedulingProblemwithEnergySaving[J].JournalofMechanicalEngineering, 2014, 25(16): 2196-2203.

    [5]FanHY,LampinenJ.ATrigonometricMutationOperationtoDifferentialEvolution[J].JournalofGlobalOptimization, 2003, 27(1): 105-129.

    [6]ZhangJ,SandersonAC.JADE:AdaptiveDifferentialEvolutionwithOptionalExternalArchive[J].IEEETransactionsonEvolutionaryComputation, 2009, 13(5): 945-958.

    [7]DasS,AbrahamA,ChakrabortyUK,etal.DifferentialEvolutionUsingaNeighborhood-basedMutationOperator[J].IEEETransactionsonEvolutionaryComputation, 2009, 13(3): 526-553.

    [8]WangY,CaiZ,ZhangQ.EnhancingtheSearchAbilityofDifferentialEvolutionthroughOrthogonalCrossover[J].InformationSciences, 2012, 185(1): 153-177.

    [9]RonkkonenJ,KukkonenS,PriceKV.Real-parameterOptimizationwithDifferentialEvolution[C]//ProceedingsoftheIEEECongressonEvolutionaryComputation(CEC'2005),Piscataway,NJ:IEEEPress,2005:506-513.

    [10]ZhuW,TangY,FangJA,etal.AdaptivePopulationTuningSchemeforDifferentialEvolution[J].InformationSciences, 2013, 223: 164-191.

    [11]BrestJ,GreinerS,BoskovicB,etal.Self-adaptingControlParametersinDifferentialEvolution:aComparativeStudyonNumericalBenchmarkProblems[J].IEEETransactionsonEvolutionaryComputation, 2006, 10(6): 646-657.

    [12]WangY,CaiZ,ZhangQ.DifferentialEvolutionwithCompositeTrialVectorGenerationStrategiesandControlParameters[J].IEEETransactionsonEvolutionaryComputation, 2011, 15(1): 55-66.

    [13]車林仙, 何兵, 程志紅. 6-CRS并聯(lián)機(jī)器人機(jī)構(gòu)及其位置分析[J]. 中國機(jī)械工程, 2010,21(14): 1669-1675.

    CheLinxian,HeBing,ChengZhihong.A6-CRSParallelManipulatorandItsPositionalAnalysis[J].JournalofMechanicalEngineering, 2010,21(14): 1669-1675.

    [14]楊曉明, 邱清盈, 馮培恩, 等. 盤式制動(dòng)器的全性能優(yōu)化設(shè)計(jì)[J]. 中國機(jī)械工程, 2005, 16(7): 630-633.

    YangXiaoming,QiuQingying,FengPei’en,etal.OptimalDesignforOverallPerformanceofDiskBrake[J].JournalofMechanicalEngineering,2005, 16(7): 630-633.

    [15]王前, 楊志堅(jiān), 丁康. 基于新自適應(yīng)差分進(jìn)化算法的MagicFormula輪胎模型參數(shù)辨識(shí)方法[J]. 機(jī)械工程學(xué)報(bào), 2014, 50(6): 120-128.

    WangQian,YangZhijian,DingKang.MethodinIdentifyingtheParametersofMagicFormulaTireModelBasedonNewSelf-adaptiveDifferentialEvolution[J].JournalofMechanicalEngineering, 2014, 50(6): 120-128.

    [16]RahnamayanS,TizhooshHR,SalamaMM.Opposition-basedDifferentialEvolution[J].IEEETransactionsonEvolutionaryComputation, 2008, 12(1): 64-79.

    [17]Xu Y, Li G, Wu Z. A Novel Hybrid Genetic Algorithm Using Local Optimizer Based on Heuristic Pattern Move[J]. Applied Artificial Intelligence, 2001, 15(7): 601-631.

    [18]劉桂萍. 基于微型遺傳算法的多目標(biāo)優(yōu)化方法及應(yīng)用研究[D].長沙: 湖南大學(xué), 2007.

    [19]Suganthan P N, Hansen N, Liang J J, et al. Problem Definitions and Evaluation Criteria for the CEC 2005 Special Session on Real-parameter Optimization[R]. Singapore:Nanyang Technological University,2005.

    [20]Qin A K, Huang V L, Suganthan P N. Differential Evolution Algorithm with Strategy Adaptation for Global Numerical Optimization[J]. IEEE Transactions on Evolutionary Computation, 2009, 13(2): 398-417.

    [21]Mallipeddi R, Suganthan P N, Pan Q K, et al.Differential Evolution Algorithm with Ensemble of Parameters and Mutation Strategies[J]. Applied Soft Computing, 2011, 11(2): 1679-1696.

    [22]Liang J J, Qin A K, Suganthan P N, et al. Comprehensive Learning Particle Swarm Optimizer for Global Optimization of Multimodal Functions[J]. IEEE Transactions on Evolutionary Computation, 2006, 10(3): 281-295.

    [23]Garcia-Martnez C, Lozano M, Herrera F, et al. Global and Local Real-coded Genetic Algorithms Based on Parent-centric Crossover Operators[J]. European Journal of Operational Research, 2008, 185(3): 1088-1113.

    [24]陳光宋, 錢林方, 徐亞棟, 等. 身管橫向固有振動(dòng)的半解析解法[J]. 兵工學(xué)報(bào), 2012, 33(10): 1168-1172.

    Chen Guangsong, Qian Linfang, Xu Yadong, et al.Semi-analysis Solution of Nature Frequency of Transverse Vibration of a Barrel[J]. Acta Armamentarii, 2012, 33(10): 1168-1172.

    [25]Jiang C, Han X, Lu G Y, et al. Correlation Analysis of Non-probabilistic Convex Model and Corresponding Structural Reliability Technique[J]. Computer Methods in Applied Mechanics and Engineering, 2011, 200(33): 2528-2546.

    (編輯王艷麗)

    Differential Evolution Algorithm with Intergeneration Projection Operator

    Fu Chunming1Jiang Chao1Chen Guangsong2Ji Lei2

    1.State Key Laboratory of Advanced Design and Manufacturing for Vehicle Body,Hunan University,Changsha,410082 2.Nanjing University of Science and Technology,Nanjing,210094

    Abstract:A DE based on intergeneration projection operator with good optimum and fast convergence performance was proposed to solve optimization problems. The proposed method based on the classical differential evolution mainly included the following characteristics. Firstly, for improving the diversity of population, opposition learning was employed to generate initial population and two different strategies were randomly selected to generate new mutant individuals. Secondly, an intergeneration projection operator was designed to generate three offsprings to substitute for the three worst individuals into the next generation. Thirdly, the exploratory operator was introduced to generate the new individuals into the next generation for keeping the diversity of evolutionary population and avoiding to obtain local solution. Finally, the performances of IPDE algorithm were verified by the eleven single-and multi-modal benchmark tests and two practical engineering problems.

    Key words:differential evolution(DE) algorithm; intergeneration projection(IP) operator; opposition learning; explorative operator

    收稿日期:2015-07-27

    基金項(xiàng)目:國家自然科學(xué)基金資助項(xiàng)目(11172096);教育部全國百篇優(yōu)秀博士論文資助項(xiàng)目(201235);湖南省杰出青年基金資助項(xiàng)目(14JJ1016)

    中圖分類號(hào):TP18

    DOI:10.3969/j.issn.1004-132X.2016.11.019

    作者簡介:符純明,男,1987年生。湖南大學(xué)機(jī)械與運(yùn)載工程學(xué)院博士研究生。主要研究方向?yàn)橹悄芩惴岸嗄繕?biāo)優(yōu)化。姜潮(通信作者),男,1978年生。湖南大學(xué)機(jī)械與運(yùn)載工程學(xué)院教授、博士研究生導(dǎo)師。陳光宋,男,1987年生。南京理工大學(xué)機(jī)械工程學(xué)院博士研究生。吉磊,男,1990年生。南京理工大學(xué)機(jī)械工程學(xué)院博士研究生。

    他把我摸到了高潮在线观看| 国产精品自产拍在线观看55亚洲 | 五月开心婷婷网| 日本欧美视频一区| av天堂在线播放| 天堂俺去俺来也www色官网| 亚洲精品成人av观看孕妇| 中文字幕精品免费在线观看视频| 丰满迷人的少妇在线观看| 国产人伦9x9x在线观看| 免费久久久久久久精品成人欧美视频| 波多野结衣一区麻豆| 亚洲av电影在线进入| 国产亚洲欧美在线一区二区| 精品久久久久久,| 日本黄色视频三级网站网址 | 在线永久观看黄色视频| 波多野结衣av一区二区av| 亚洲少妇的诱惑av| 天堂俺去俺来也www色官网| 丝袜在线中文字幕| 纯流量卡能插随身wifi吗| www.自偷自拍.com| 精品人妻1区二区| 波多野结衣av一区二区av| 国产高清激情床上av| 无遮挡黄片免费观看| 99久久人妻综合| 午夜视频精品福利| 国产一区二区三区视频了| 又黄又爽又免费观看的视频| 欧美国产精品一级二级三级| 欧美日韩福利视频一区二区| 欧美日韩福利视频一区二区| 国产极品粉嫩免费观看在线| 亚洲熟女精品中文字幕| 最新的欧美精品一区二区| a在线观看视频网站| 亚洲人成电影观看| 日韩大码丰满熟妇| 亚洲精品自拍成人| 免费女性裸体啪啪无遮挡网站| 久久99一区二区三区| 久久99一区二区三区| 满18在线观看网站| 久久中文字幕一级| 叶爱在线成人免费视频播放| 岛国毛片在线播放| 久久香蕉国产精品| 99精品欧美一区二区三区四区| 成人免费观看视频高清| 大陆偷拍与自拍| 最近最新中文字幕大全免费视频| 69av精品久久久久久| 久久精品91无色码中文字幕| 这个男人来自地球电影免费观看| 国产黄色免费在线视频| 嫁个100分男人电影在线观看| 黑丝袜美女国产一区| 久久精品人人爽人人爽视色| 国产免费现黄频在线看| 极品教师在线免费播放| 91av网站免费观看| 满18在线观看网站| 日韩大码丰满熟妇| 波多野结衣一区麻豆| 精品无人区乱码1区二区| 国产亚洲欧美98| 国产精品久久电影中文字幕 | 国产免费男女视频| 久久精品亚洲精品国产色婷小说| 美女国产高潮福利片在线看| 少妇 在线观看| 国产精品一区二区精品视频观看| 国产精品98久久久久久宅男小说| 国产精品亚洲一级av第二区| 久久香蕉精品热| 亚洲精品成人av观看孕妇| 一进一出抽搐gif免费好疼 | 老司机靠b影院| 老熟妇仑乱视频hdxx| 十八禁高潮呻吟视频| videos熟女内射| 久久久精品免费免费高清| 男女床上黄色一级片免费看| 久久中文看片网| 啦啦啦 在线观看视频| 亚洲av电影在线进入| 狠狠狠狠99中文字幕| 国精品久久久久久国模美| 国精品久久久久久国模美| 欧美乱色亚洲激情| 欧美日韩精品网址| 动漫黄色视频在线观看| 天天躁夜夜躁狠狠躁躁| 午夜福利视频在线观看免费| 日韩欧美一区视频在线观看| 亚洲精品一二三| 国产亚洲欧美在线一区二区| 亚洲av日韩在线播放| 国产精华一区二区三区| 操出白浆在线播放| 久久久久视频综合| 国产激情久久老熟女| 亚洲 欧美一区二区三区| 午夜福利影视在线免费观看| 水蜜桃什么品种好| 婷婷成人精品国产| 波多野结衣av一区二区av| 欧美激情 高清一区二区三区| 天堂√8在线中文| 日本a在线网址| 欧美精品人与动牲交sv欧美| 黄片大片在线免费观看| 久久香蕉激情| 亚洲免费av在线视频| 欧美大码av| 黄色毛片三级朝国网站| 欧美 亚洲 国产 日韩一| 不卡av一区二区三区| 999精品在线视频| 亚洲国产精品sss在线观看 | 亚洲精品成人av观看孕妇| 午夜精品国产一区二区电影| 亚洲第一av免费看| av欧美777| 深夜精品福利| 欧美日韩亚洲综合一区二区三区_| 在线观看午夜福利视频| 人人妻人人澡人人看| av网站在线播放免费| 精品第一国产精品| 久久久久国产精品人妻aⅴ院 | 最新的欧美精品一区二区| 很黄的视频免费| 国产人伦9x9x在线观看| 热re99久久国产66热| 两性午夜刺激爽爽歪歪视频在线观看 | 97人妻天天添夜夜摸| 成人av一区二区三区在线看| 久久国产精品男人的天堂亚洲| 国产av一区二区精品久久| 女性生殖器流出的白浆| 老鸭窝网址在线观看| 两性夫妻黄色片| 久久精品国产99精品国产亚洲性色 | 精品久久蜜臀av无| 精品久久久久久久久久免费视频 | 极品人妻少妇av视频| 国内久久婷婷六月综合欲色啪| 国产人伦9x9x在线观看| 黑人操中国人逼视频| 曰老女人黄片| 飞空精品影院首页| 欧美人与性动交α欧美软件| 午夜福利影视在线免费观看| 精品电影一区二区在线| avwww免费| 欧美黄色片欧美黄色片| 亚洲,欧美精品.| 日韩欧美一区二区三区在线观看 | 久久久久国内视频| 久99久视频精品免费| 亚洲成人免费av在线播放| 18禁观看日本| 黄色女人牲交| 一区二区三区激情视频| 女人高潮潮喷娇喘18禁视频| 久久亚洲精品不卡| a级片在线免费高清观看视频| 亚洲精品美女久久av网站| 久久人妻熟女aⅴ| 亚洲一区二区三区欧美精品| 日本撒尿小便嘘嘘汇集6| 高清在线国产一区| 日日爽夜夜爽网站| 亚洲国产精品合色在线| 91字幕亚洲| 国产黄色免费在线视频| 十分钟在线观看高清视频www| 色播在线永久视频| 久久久久久人人人人人| 免费看a级黄色片| 热re99久久精品国产66热6| 国产免费av片在线观看野外av| 狠狠狠狠99中文字幕| 久久久国产一区二区| 在线十欧美十亚洲十日本专区| 亚洲国产精品sss在线观看 | 天天影视国产精品| 人成视频在线观看免费观看| 欧美黄色淫秽网站| 日韩免费高清中文字幕av| 三上悠亚av全集在线观看| 色94色欧美一区二区| 一进一出抽搐动态| 一级毛片精品| 久久久久久免费高清国产稀缺| 免费少妇av软件| 另类亚洲欧美激情| 久久精品亚洲熟妇少妇任你| 久久久久久免费高清国产稀缺| 在线视频色国产色| 亚洲熟妇中文字幕五十中出 | 国产亚洲一区二区精品| 精品亚洲成a人片在线观看| 亚洲精品一二三| 亚洲精品一二三| 色94色欧美一区二区| 欧美激情高清一区二区三区| 免费观看人在逋| 两性夫妻黄色片| 中文字幕色久视频| 99精品在免费线老司机午夜| 亚洲专区国产一区二区| 亚洲欧美色中文字幕在线| 国产一区二区激情短视频| 一级作爱视频免费观看| av天堂久久9| 97人妻天天添夜夜摸| 日本五十路高清| 欧美激情 高清一区二区三区| cao死你这个sao货| 99热网站在线观看| 少妇猛男粗大的猛烈进出视频| 黑人巨大精品欧美一区二区mp4| 亚洲成人免费av在线播放| 美女高潮到喷水免费观看| 久久精品91无色码中文字幕| 一级片'在线观看视频| 亚洲国产欧美网| 欧美日韩黄片免| av在线播放免费不卡| 成人亚洲精品一区在线观看| 欧美精品一区二区免费开放| 欧美日韩亚洲国产一区二区在线观看 | 超色免费av| 免费在线观看亚洲国产| av国产精品久久久久影院| 妹子高潮喷水视频| 1024香蕉在线观看| 国产精品久久久久久人妻精品电影| 久久久精品国产亚洲av高清涩受| 久久人妻av系列| 岛国在线观看网站| 91成人精品电影| 大香蕉久久网| 久久九九热精品免费| 精品久久久久久久久久免费视频 | 亚洲色图综合在线观看| 亚洲熟女毛片儿| 热99久久久久精品小说推荐| 一夜夜www| 高清av免费在线| 久久久久国产精品人妻aⅴ院 | 一区在线观看完整版| 夫妻午夜视频| 美女国产高潮福利片在线看| 日本五十路高清| 极品教师在线免费播放| 久99久视频精品免费| 91在线观看av| 村上凉子中文字幕在线| 成人亚洲精品一区在线观看| 久久这里只有精品19| 亚洲欧美一区二区三区久久| av天堂在线播放| 亚洲人成77777在线视频| 日韩有码中文字幕| 亚洲国产欧美网| 亚洲精品在线观看二区| 欧美日韩亚洲综合一区二区三区_| 国产成人免费观看mmmm| 咕卡用的链子| 新久久久久国产一级毛片| 99精国产麻豆久久婷婷| 一区在线观看完整版| 亚洲色图综合在线观看| 亚洲 国产 在线| 天堂√8在线中文| 国产极品粉嫩免费观看在线| 久久人人爽av亚洲精品天堂| 国产99白浆流出| 日韩 欧美 亚洲 中文字幕| 啪啪无遮挡十八禁网站| 国产亚洲欧美精品永久| 免费在线观看视频国产中文字幕亚洲| 纯流量卡能插随身wifi吗| 俄罗斯特黄特色一大片| 成熟少妇高潮喷水视频| 日日夜夜操网爽| 在线av久久热| 美女福利国产在线| 一级a爱片免费观看的视频| 99国产精品一区二区蜜桃av | 成人影院久久| 亚洲欧美日韩另类电影网站| 伊人久久大香线蕉亚洲五| 国产av精品麻豆| 精品国产一区二区久久| 亚洲欧美一区二区三区久久| 天天躁狠狠躁夜夜躁狠狠躁| 久久热在线av| 精品高清国产在线一区| av中文乱码字幕在线| 亚洲av熟女| 欧美中文综合在线视频| 水蜜桃什么品种好| 高清视频免费观看一区二区| 欧美日韩视频精品一区| 啦啦啦在线免费观看视频4| 不卡一级毛片| 国产三级黄色录像| 亚洲免费av在线视频| 老司机深夜福利视频在线观看| 午夜福利欧美成人| 老鸭窝网址在线观看| 91老司机精品| 老熟女久久久| 亚洲欧美一区二区三区黑人| 正在播放国产对白刺激| 亚洲熟女毛片儿| 国产视频一区二区在线看| 成年动漫av网址| 久9热在线精品视频| 午夜福利影视在线免费观看| 日韩欧美三级三区| 亚洲精品乱久久久久久| 国产成人欧美| 91成年电影在线观看| 丰满的人妻完整版| 久久香蕉激情| 曰老女人黄片| 国产成人精品在线电影| 多毛熟女@视频| 伊人久久大香线蕉亚洲五| 欧美激情极品国产一区二区三区| 男人的好看免费观看在线视频 | 动漫黄色视频在线观看| 精品久久久久久久毛片微露脸| 色尼玛亚洲综合影院| av中文乱码字幕在线| 黄频高清免费视频| 一级片'在线观看视频| 变态另类成人亚洲欧美熟女 | 久久国产精品影院| 久久久精品免费免费高清| 51午夜福利影视在线观看| 国产成人av教育| 脱女人内裤的视频| 建设人人有责人人尽责人人享有的| 国产又色又爽无遮挡免费看| 99精品久久久久人妻精品| 午夜福利影视在线免费观看| 1024香蕉在线观看| 后天国语完整版免费观看| 丝袜美足系列| 日韩欧美一区视频在线观看| 国产成人影院久久av| 亚洲av成人一区二区三| 久久久久久久国产电影| 国产日韩欧美亚洲二区| 两个人免费观看高清视频| 天堂动漫精品| 亚洲av电影在线进入| 亚洲一码二码三码区别大吗| 欧美日韩成人在线一区二区| 在线观看免费午夜福利视频| e午夜精品久久久久久久| 法律面前人人平等表现在哪些方面| 亚洲国产中文字幕在线视频| 色播在线永久视频| 成年版毛片免费区| 999精品在线视频| 女人高潮潮喷娇喘18禁视频| 日韩精品免费视频一区二区三区| 91麻豆精品激情在线观看国产 | 天天躁狠狠躁夜夜躁狠狠躁| 在线永久观看黄色视频| 99re6热这里在线精品视频| 99精国产麻豆久久婷婷| 色精品久久人妻99蜜桃| 女警被强在线播放| 可以免费在线观看a视频的电影网站| 欧美中文综合在线视频| 一级毛片女人18水好多| 欧美乱色亚洲激情| 日韩中文字幕欧美一区二区| 男人舔女人的私密视频| 无遮挡黄片免费观看| 他把我摸到了高潮在线观看| 99在线人妻在线中文字幕 | 久久人妻av系列| 午夜91福利影院| av免费在线观看网站| 伊人久久大香线蕉亚洲五| 一区福利在线观看| 亚洲精品美女久久av网站| 狠狠狠狠99中文字幕| 如日韩欧美国产精品一区二区三区| 欧美丝袜亚洲另类 | 涩涩av久久男人的天堂| 精品国产美女av久久久久小说| 久久精品人人爽人人爽视色| 亚洲精华国产精华精| 久99久视频精品免费| 欧美日韩亚洲国产一区二区在线观看 | 国产午夜精品久久久久久| 亚洲一码二码三码区别大吗| 中文字幕人妻熟女乱码| 国产亚洲欧美精品永久| 亚洲精品av麻豆狂野| 色老头精品视频在线观看| 久久久久久人人人人人| 欧美黑人精品巨大| 久久国产精品大桥未久av| 青草久久国产| 美女扒开内裤让男人捅视频| 一区二区三区国产精品乱码| 水蜜桃什么品种好| av天堂在线播放| 天堂动漫精品| 夜夜爽天天搞| e午夜精品久久久久久久| 国产成人av激情在线播放| 亚洲国产精品合色在线| 伊人久久大香线蕉亚洲五| 亚洲成人国产一区在线观看| 叶爱在线成人免费视频播放| 制服诱惑二区| 国产精品永久免费网站| 天天躁日日躁夜夜躁夜夜| 国产亚洲一区二区精品| 啦啦啦免费观看视频1| 19禁男女啪啪无遮挡网站| 久久久久国内视频| 婷婷丁香在线五月| 黑人操中国人逼视频| 精品国产超薄肉色丝袜足j| 免费少妇av软件| 午夜福利在线免费观看网站| 两性夫妻黄色片| 黑人猛操日本美女一级片| 无限看片的www在线观看| 悠悠久久av| 国产精品国产av在线观看| а√天堂www在线а√下载 | 国产精品98久久久久久宅男小说| 精品亚洲成国产av| 日韩熟女老妇一区二区性免费视频| 亚洲国产欧美一区二区综合| 欧美亚洲 丝袜 人妻 在线| 日本vs欧美在线观看视频| 日日夜夜操网爽| 久久久久久免费高清国产稀缺| 一级片免费观看大全| 制服人妻中文乱码| 久久久久久人人人人人| 久久人妻熟女aⅴ| 搡老乐熟女国产| 一边摸一边抽搐一进一出视频| 十八禁网站免费在线| 在线观看免费日韩欧美大片| 91大片在线观看| 国产亚洲精品久久久久久毛片 | 久久亚洲真实| 欧美日韩亚洲国产一区二区在线观看 | e午夜精品久久久久久久| 国产成人精品久久二区二区免费| 国产精品综合久久久久久久免费 | 精品乱码久久久久久99久播| 久久久久国内视频| 国产高清videossex| 精品少妇一区二区三区视频日本电影| 亚洲精品粉嫩美女一区| 看免费av毛片| 宅男免费午夜| 成年女人毛片免费观看观看9 | av视频免费观看在线观看| 视频区图区小说| 久久久久国产精品人妻aⅴ院 | 成人影院久久| 午夜免费成人在线视频| 国产黄色免费在线视频| 午夜免费鲁丝| 欧美 日韩 精品 国产| 成人免费观看视频高清| 精品无人区乱码1区二区| 丝袜美足系列| 久久久国产欧美日韩av| 亚洲国产精品一区二区三区在线| 黄色视频,在线免费观看| 看片在线看免费视频| 五月开心婷婷网| 久久精品国产综合久久久| 免费观看人在逋| 成人三级做爰电影| 欧美日韩精品网址| 另类亚洲欧美激情| 久久久久国内视频| 怎么达到女性高潮| 大型黄色视频在线免费观看| 极品教师在线免费播放| 精品久久蜜臀av无| 9热在线视频观看99| 国产91精品成人一区二区三区| 午夜精品在线福利| 成人精品一区二区免费| 午夜福利视频在线观看免费| 男女之事视频高清在线观看| 亚洲综合色网址| 亚洲欧美日韩另类电影网站| 中文字幕人妻丝袜一区二区| bbb黄色大片| 搡老熟女国产l中国老女人| 国产一区在线观看成人免费| 午夜福利在线观看吧| xxx96com| 又黄又粗又硬又大视频| 国产在线一区二区三区精| 天堂中文最新版在线下载| 免费av中文字幕在线| 久久精品亚洲熟妇少妇任你| 一进一出抽搐gif免费好疼 | 久久精品亚洲精品国产色婷小说| 搡老岳熟女国产| 国产精品久久视频播放| 精品一区二区三区av网在线观看| 亚洲专区中文字幕在线| 久久中文字幕人妻熟女| svipshipincom国产片| 国产亚洲一区二区精品| 90打野战视频偷拍视频| 色综合欧美亚洲国产小说| 日本黄色视频三级网站网址 | 高清av免费在线| 亚洲中文字幕日韩| 久久久国产精品麻豆| 久久精品国产99精品国产亚洲性色 | 99精品久久久久人妻精品| 女人爽到高潮嗷嗷叫在线视频| avwww免费| 99精品久久久久人妻精品| 国产精品久久久久成人av| 80岁老熟妇乱子伦牲交| 精品熟女少妇八av免费久了| 一本大道久久a久久精品| 91麻豆精品激情在线观看国产 | 国内毛片毛片毛片毛片毛片| 亚洲精品国产色婷婷电影| 一进一出抽搐gif免费好疼 | 日韩欧美一区二区三区在线观看 | 免费在线观看亚洲国产| 九色亚洲精品在线播放| 亚洲精品国产精品久久久不卡| 99久久综合精品五月天人人| 中文亚洲av片在线观看爽 | 国产精品永久免费网站| 欧美+亚洲+日韩+国产| 亚洲五月天丁香| 免费在线观看日本一区| 久久狼人影院| 精品熟女少妇八av免费久了| 成年人午夜在线观看视频| 啦啦啦 在线观看视频| 精品久久蜜臀av无| 飞空精品影院首页| 人成视频在线观看免费观看| 久久 成人 亚洲| 日韩精品免费视频一区二区三区| 精品熟女少妇八av免费久了| 国产乱人伦免费视频| 久久精品人人爽人人爽视色| 国产精品永久免费网站| 午夜免费成人在线视频| 亚洲五月婷婷丁香| 在线十欧美十亚洲十日本专区| 日本精品一区二区三区蜜桃| 精品午夜福利视频在线观看一区| 国产99白浆流出| 精品福利永久在线观看| 亚洲av成人一区二区三| 欧美黄色片欧美黄色片| 国产成人欧美在线观看 | 狠狠婷婷综合久久久久久88av| 精品第一国产精品| 国产精品久久久av美女十八| 日韩有码中文字幕| 久久香蕉精品热| 成在线人永久免费视频| 看免费av毛片| 91麻豆精品激情在线观看国产 | av网站免费在线观看视频| 久久国产精品人妻蜜桃| 国产精品香港三级国产av潘金莲| 久久香蕉激情| 成人特级黄色片久久久久久久| 国产成人精品久久二区二区91| 午夜日韩欧美国产| 老汉色av国产亚洲站长工具| 国产精品免费大片| 色在线成人网| 人人妻人人澡人人看| 热re99久久精品国产66热6| 母亲3免费完整高清在线观看| 一级毛片高清免费大全| 黄色女人牲交| 美女扒开内裤让男人捅视频| 亚洲成人免费av在线播放| 视频区图区小说| 女人被躁到高潮嗷嗷叫费观| 超色免费av| 在线十欧美十亚洲十日本专区| 国产男女超爽视频在线观看| 国产高清国产精品国产三级|