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

    一種基于狀態(tài)空間模型的進(jìn)化算法

    2014-08-08 23:47:59李茂軍賈玲

    李茂軍+賈玲

    收稿日期:2013-09-18

    作者簡(jiǎn)介:李茂軍(1964—),男,湖南寧鄉(xiāng)人,教授,博士,研究方向:智能控制與智能計(jì)算。

    文章編號(hào):1003-6199(2014)02-0085-04

    摘 要:傳統(tǒng)進(jìn)化算法主要通過(guò)選擇、重組和變異這三種遺傳操作實(shí)現(xiàn)種群的進(jìn)化。在進(jìn)化過(guò)程中通常需要設(shè)定群體規(guī)模、交叉概率和變異概率等參數(shù),而且它們的值會(huì)直接影響計(jì)算結(jié)果及精度。為了簡(jiǎn)化操作過(guò)程,設(shè)計(jì)一種基于離散系統(tǒng)狀態(tài)空間模型的進(jìn)化算法,這種算法采用實(shí)數(shù)編碼方式,構(gòu)造一個(gè)狀態(tài)進(jìn)化矩陣來(lái)實(shí)現(xiàn)重組和變異的功能,提高算法的可操作性和可靠性。并將該算法應(yīng)用于求解無(wú)約束全局優(yōu)化問(wèn)題,對(duì)幾種典型的測(cè)試函數(shù)進(jìn)行仿真,結(jié)果表明:這種新的進(jìn)化算法具有搜索能力強(qiáng)、收斂速度快、計(jì)算精度高、操作簡(jiǎn)單等優(yōu)點(diǎn),對(duì)相關(guān)研究有參考作用。

    關(guān)鍵詞:進(jìn)化算法;狀態(tài)空間模型;實(shí)數(shù)編碼;狀態(tài)進(jìn)化矩陣

    中圖分類(lèi)號(hào):TP301.6文獻(xiàn)標(biāo)識(shí)碼:A

    

    An Evolutionary Algorithm Based on StatespaceModel

    

    LI Maojun,JIA Ling

    (College of Electrical and Information Engineering, Changsha University of Science & Technology, Changsha,Hunan 410114,China)

    Abstract:The traditional evolutionary algorithm primarily through three genetic operators: selection, recombination and mutation operations, to achieve the evolution of the population. In the process of evolution, it usually needs to set the crossover probability and mutation probability, which will directly affect the results and precision. In order to simplify the procedure, we design a new evolutionary algorithm, which based on discrete state-space model system and using real-encoding method. The algorithm constructs a state evolution matrix to achieve the function of recombination and mutation, and improve the operability and reliability of the algorithm. We do some simulation based on several typical test functions, the results shows that: this new evolutionary algorithm has many advantages, such as strong search capability, rapid convergence, high precision, simple operation, etc. It has useful reference for relevant studies.

    Key words:evolutionary algorithm; state-space model ; real-encoding; state evolution matrix

    1 引 言

    進(jìn)化算法(EA)是一類(lèi)模擬生物進(jìn)化機(jī)制的智能優(yōu)化方法,如遺傳算法(GA)[1]、蟻群算法(ACO)[2]、模擬退火算法(SA)[3]等。同傳統(tǒng)的梯度法、牛頓法、窮舉法等優(yōu)化算法相比,進(jìn)化計(jì)算具有自組織、自適應(yīng)、自學(xué)習(xí)、不受問(wèn)題性質(zhì)限制的優(yōu)點(diǎn),因此進(jìn)化算法常用來(lái)解決復(fù)雜的工程優(yōu)化問(wèn)題[4]。隨著科學(xué)的發(fā)展和應(yīng)用需求的增加,傳統(tǒng)進(jìn)化算法已不能滿足工程應(yīng)用需要。幾十年來(lái),許多學(xué)者嘗試了很多方法來(lái)更好地解決優(yōu)化問(wèn)題,如對(duì)傳統(tǒng)進(jìn)化算法進(jìn)行改進(jìn)、引入新的理論、結(jié)合兩種或兩種以上進(jìn)化算法等來(lái)處理優(yōu)化問(wèn)題,取得了一定的效果[5-7]。

    文獻(xiàn)[8]提出一種改進(jìn)的遺傳算法,為了避免連續(xù)函數(shù)優(yōu)化過(guò)程中的早熟收斂和搜索遲鈍,在簡(jiǎn)單遺傳算法基礎(chǔ)上提出了劃分尋優(yōu)區(qū)間、基于排序和最佳保留的輪盤(pán)賭選擇算子, 并采用擇優(yōu)交叉算子和二元變異算子,提高了算法的運(yùn)行效率和收斂速度,并可避免陷入局部最優(yōu);文獻(xiàn)[9] 針對(duì)粒子群算法(PSO)算法存在進(jìn)化后期收斂速度慢、易陷入局部最優(yōu)點(diǎn)的缺點(diǎn),提出了一種多向?qū)W習(xí)型的粒子群優(yōu)化算法,該算法中粒子通過(guò)同時(shí)追隨自己找到的最優(yōu)解、隨機(jī)的其他粒子同維度的最優(yōu)解和整個(gè)群的最優(yōu)解來(lái)完成速度更新,通過(guò)判別區(qū)域邊界來(lái)完成位置優(yōu)化更新,通過(guò)對(duì)全局最優(yōu)位置進(jìn)行小范圍擾動(dòng),以增強(qiáng)算法跳出局部最優(yōu)的能力。明顯改善了全局搜索能力,并且能夠有效避免早熟收斂問(wèn)題。文獻(xiàn)[10]提出了一種結(jié)合免疫克隆算子的量子遺傳算法(QGA),采用免疫克隆操作及交叉策略提高抗體成熟力及親和性,增強(qiáng)抗體群分布的多樣性及穩(wěn)定性,有效克服了量子遺傳算法容易陷于局部最優(yōu)及計(jì)算緩慢的不足。

    傳統(tǒng)進(jìn)化算法存在的問(wèn)題有:一、編程過(guò)程比較復(fù)雜,算法開(kāi)始要先對(duì)所求問(wèn)題進(jìn)行編碼,最后對(duì)找到最優(yōu)解還要進(jìn)行解碼;二、遺傳算子和初始種群的選擇對(duì)解的品質(zhì)影響很大,大部分需要依靠經(jīng)驗(yàn)來(lái)選擇;三、搜索速度比較慢,要得到精確度高的解需要花很長(zhǎng)時(shí)間;四、容易出現(xiàn)早熟收斂現(xiàn)象,陷入局部最優(yōu)解而無(wú)法跳出。

    針對(duì)傳統(tǒng)進(jìn)化算法的存在的問(wèn)題,本文提出一種基于狀態(tài)空間模型的進(jìn)化算法(SEA),這種算法采用實(shí)數(shù)編碼方式,構(gòu)造一個(gè)狀態(tài)進(jìn)化矩陣來(lái)實(shí)現(xiàn)種群進(jìn)化,并通過(guò)選種池中的選擇操作實(shí)現(xiàn)優(yōu)勝劣汰的自然選擇機(jī)制。通過(guò)對(duì)幾種典型函數(shù)的測(cè)試結(jié)果表明,該算法具有很強(qiáng)的搜索能力和很高的搜索精度,能快速地找到問(wèn)題的全局最優(yōu)解。

    2 基于狀態(tài)空間模型的進(jìn)化算法

    算法基于離散系統(tǒng)狀態(tài)空間模型,引入進(jìn)化計(jì)算的基本思想,構(gòu)造一種基于狀態(tài)空間模型X'(k+1)=GX(k)(其中X(k)為第k個(gè)采樣時(shí)刻的狀態(tài)向量,G為狀態(tài)進(jìn)化矩陣)的進(jìn)化算法。在這種算法中,進(jìn)化算法的群體表示為狀態(tài)向量X(k),X(k)表示第k代群體,狀態(tài)向量X(k)包含N個(gè)分量,每個(gè)分量均表示1個(gè)個(gè)體(這里的個(gè)體是按傳統(tǒng)進(jìn)化算法中的實(shí)數(shù)編碼方法而得到的),每個(gè)個(gè)體包含M個(gè)變量。在這里,狀態(tài)向量X(k)實(shí)際上是一個(gè)N×M矩陣,該矩陣的每一行表示一個(gè)個(gè)體,每一個(gè)元素是變量的實(shí)數(shù)值。群體進(jìn)化通過(guò)狀態(tài)進(jìn)化矩陣G實(shí)現(xiàn),G是一個(gè)N×N的矩陣??苫谶M(jìn)化算法中群體進(jìn)化的基本方法來(lái)構(gòu)造狀態(tài)進(jìn)化矩陣G,也可以通過(guò)其它途徑來(lái)構(gòu)造狀態(tài)進(jìn)化矩陣G。本文主要基于進(jìn)化計(jì)算中群體進(jìn)化的基本思想來(lái)構(gòu)造狀態(tài)進(jìn)化矩陣G。其計(jì)算模型如圖1所示。

    計(jì)算技術(shù)與自動(dòng)化2014年6月

    第33卷第2期李茂軍等:一種基于狀態(tài)空間模型的進(jìn)化算法

    圖1 基于狀態(tài)空間模型的進(jìn)化算法計(jì)算模型

    2.1 初始種群

    設(shè)種群規(guī)模為N,待優(yōu)化問(wèn)題的變量有M個(gè),則第k代群體X(k)為

    X(k)=x11x12…x1Mx21x22…x2M…xN1xN2…xNM (1) 

    其中Xi=(xi1,xi2,…,xiM),i=1,2,…,N為優(yōu)化問(wèn)題的一組可行解,xij∈[αj,βj],i=1,2,…,N,j=1,2,…,M。

    在算法初始化階段,用隨機(jī)函數(shù)產(chǎn)生N組、每組M個(gè)分別在αj,βj(j=1,2,…,M)上服從均勻分布的實(shí)數(shù)構(gòu)成初始種群X0。

    2.2 狀態(tài)進(jìn)化矩陣

    本算法采用一個(gè)狀態(tài)進(jìn)化矩陣G來(lái)模擬進(jìn)化算法中的重組與變異操作過(guò)程,這使得算法操作起來(lái)非常簡(jiǎn)單,因此,狀態(tài)進(jìn)化矩陣的構(gòu)造是算法的重點(diǎn)。本文按照式(2)構(gòu)造狀態(tài)進(jìn)化矩陣G

    G=g11g12…g1Ng21g22…g2N…gN1gN2…gNN (2)

    其中0≤gij≤1,i,j=1,2,…,N,且∑nj=1gij≤1 ,矩陣G中元素gij(i,j=1,2,…,N)的值是隨機(jī)確定的,本文中G為滿足上述條件的隨機(jī)常數(shù)矩陣,具體構(gòu)造流程如圖2所示。

    2.3 選種池和選擇操作

    同傳統(tǒng)進(jìn)化算法類(lèi)似,選擇操作是一個(gè)擇優(yōu)的過(guò)程,體現(xiàn)了優(yōu)勝劣汰的進(jìn)化思想,本文按最小化優(yōu)化問(wèn)題尋優(yōu)。如圖1所示,從初始群體X(0)開(kāi)始,按照X'(k+1)=GX(k)迭代,可依次得到一系列群體X'(1),X'(2),X'(3),……。群體X'(k+1)和X(k)(k= 0,1,2,…)同時(shí)進(jìn)入選種池,按照適應(yīng)度函數(shù)f()計(jì)算選種池中每個(gè)群體的適應(yīng)度值Y(k),按照從小到大排列,選擇前N個(gè)個(gè)體組成新一代群體X(k+1),如此反復(fù),直到算法滿足停機(jī)條件。

    圖2 狀態(tài)進(jìn)化矩陣構(gòu)造流程圖

    2.4 算法描述

    Step 1. 初始化種群 X(0),并計(jì)算初始種群的適應(yīng)度值;

    Step 2. 根據(jù)圖2流程得到狀態(tài)進(jìn)化矩陣G;

    Step 3. 進(jìn)行進(jìn)化操作,X'(k+1)=G?X(k),并計(jì)算種群X'(k+1)適應(yīng)度值;

    Step 4. 將種群X(k),X'(k+1)同時(shí)放入選種池,按適應(yīng)度值從小到大排列,選擇前N個(gè)個(gè)體作為新一代的群體X(k+1);

    Step 5. 是否滿足終止條件:若是,則結(jié)束;否則,轉(zhuǎn)到Step 3。

    3 仿真實(shí)驗(yàn)與性能分析

    實(shí)驗(yàn)仿真平臺(tái)為Matlab 7.1。本實(shí)驗(yàn)采用了3個(gè)經(jīng)典測(cè)試函數(shù),這些測(cè)試函數(shù)為高維多模函數(shù),具有大量的局部最優(yōu)點(diǎn),是優(yōu)化領(lǐng)域中公認(rèn)的較難優(yōu)化的函數(shù)。

    1)Shubert函數(shù)

    f1=∑5i=1icos [(i+1)x+i]?

    ∑5i=1icos [(i+1)y+i]-10≤x,y≤10 (3)

    有多個(gè)極小值點(diǎn),但只有一個(gè)全局最小值-186.73。

    2)Schwefel's函數(shù)

    f2=-xsin (x)-ysin (y)-500≤x,y≤500(4)

    是一個(gè)具有典型欺騙問(wèn)題的函數(shù),有1個(gè)全局極小值點(diǎn)取值近似為-837.9658,在(420.96 87,420. 9687)處,距離另一個(gè)局部最優(yōu)點(diǎn)很遠(yuǎn),因此如果陷入局部最優(yōu)就很難跳出。

    (3)GoldsteinPrice函數(shù)

    f3=[1+(x+y+1)2(19-14x+3x2-14y+6xy+3y2)]?[30+(2x-3y)2(18-32x+12x2+48y-36xy+27y2)]-2≤x,y≤2(5)

    這是一個(gè)多模函數(shù),有多個(gè)極小值點(diǎn),但只有一個(gè)全局最小值3,極小值點(diǎn)為(0,-1)。

    在測(cè)試中,將種群大小設(shè)置為500,每次運(yùn)行代數(shù)為60代,每個(gè)測(cè)試函數(shù)運(yùn)行測(cè)試50次。表1顯示了函數(shù)f1~f3運(yùn)行50次的實(shí)驗(yàn)結(jié)果,包括平均最優(yōu)值、標(biāo)準(zhǔn)差和平均收斂代數(shù)。全局最優(yōu)值為f1~f3理論最優(yōu)值,平均最優(yōu)值fmav按照公式(6)計(jì)算得到(其中fmi為第i次運(yùn)用本算法求解f1~f3優(yōu)化問(wèn)題得到的最優(yōu)解對(duì)應(yīng)的適應(yīng)度函數(shù)值),標(biāo)準(zhǔn)差為平均最優(yōu)值與全局最優(yōu)值之差,平均收斂代數(shù)為每次實(shí)驗(yàn)找到最優(yōu)解所需要的最少代數(shù)的平均值。

    fmav=∑50i=1fmi/50,i=1,2,…,50(6)

    表1 算法對(duì)3個(gè)測(cè)試函數(shù)50次實(shí)驗(yàn)的結(jié)果

    函數(shù)

    全局最優(yōu)值

    平均最優(yōu)值

    標(biāo)準(zhǔn)差

    平均收斂代數(shù)

    f1

    -186.7316

    -186.7316

    0

    9

    f2

    -837.9658

    -837.9658

    0

    12

    f3

    3.0000

    3.0000

    0

    8

    圖3~5分別顯示了運(yùn)用本算法求解函數(shù)f1~f3優(yōu)化問(wèn)題所得到的最優(yōu)解對(duì)應(yīng)的適應(yīng)度函數(shù)值隨進(jìn)化代數(shù)變化曲線。

    圖3 算法求解f1優(yōu)化問(wèn)題適應(yīng)度函數(shù)值變化曲線 

    圖4 算法求解f2優(yōu)化問(wèn)題適應(yīng)度函數(shù)值變化曲線

    圖5 算法求解f3優(yōu)化問(wèn)題適應(yīng)度函數(shù)值變化曲線

    從表1和圖3~5可以看出,本文提出的算法在3個(gè)函數(shù)優(yōu)化中取得了良好的效果,找到最優(yōu)解的成功率高,收斂速度快。通過(guò)50次實(shí)驗(yàn)結(jié)果可以看出,基于狀態(tài)空間模型的進(jìn)化算法穩(wěn)定性好,成功率高。

    4 結(jié) 論

    本文設(shè)計(jì)了一種基于離散系統(tǒng)狀態(tài)空間模型的進(jìn)化算法,以基本進(jìn)化算法的思想為基礎(chǔ),借鑒矩陣?yán)碚摵蛯?shí)數(shù)編碼方法,通過(guò)構(gòu)造一個(gè)狀態(tài)進(jìn)化矩陣,來(lái)實(shí)現(xiàn)種群的進(jìn)化,提高了算法的可操作性和可靠性。通過(guò)對(duì)3個(gè)經(jīng)典優(yōu)化測(cè)試函數(shù)進(jìn)行優(yōu)化實(shí)驗(yàn),結(jié)果表明:本文提出的算法增強(qiáng)了防止陷入局部最優(yōu)的能力,能夠較大地提高收斂速度和精度,而且穩(wěn)定性能很好。

    雖然算法的仿真結(jié)果是令人滿意的,但是由于實(shí)驗(yàn)條件的局限性和實(shí)驗(yàn)次數(shù)的有限性,算法收斂性和遍歷性還需要通過(guò)進(jìn)一步的數(shù)學(xué)證明。如何構(gòu)造更好的狀態(tài)進(jìn)化矩陣,是需要進(jìn)一步研究的問(wèn)題。

    參考文獻(xiàn)

    [1] 劉鯖潔,陳桂明,劉小方. 基于矩陣編碼的遺傳算法研究[J].計(jì)算機(jī)工程:2011, 37(13):160-162.

    [2] IRINA CIORNEI, ELIAS KYRIAKIDES. Hybrid Ant Colony-Genetic Algorithm (GAAPI) for Global Continuous Optimization [J]. Systems, Man, and Cybernetics, part B, 2012, 42(1), 34-245.

    [3] 董麗麗,龔光紅,李妮,等.基于云模型的自適應(yīng)并行模擬退火遺傳算法[J].北京航空航天大學(xué)學(xué)報(bào):2011,37 (9):1132-1136.

    [4] 劉正龍,楊艷梅.基于遺傳算法的數(shù)值優(yōu)化約束問(wèn)題的研究[J].計(jì)算機(jī)系統(tǒng)應(yīng)用:2013,22(5):139-142.

    [5] 王巍,趙文紅,王宇平.一種有效的解無(wú)約束全局優(yōu)化的進(jìn)化算法[J].控制理論與應(yīng)用:2010,27 (5):570-574.

    [6] 陳曉峰,楊廣明,黃明.一種實(shí)數(shù)編碼的量子差分進(jìn)化算法[J].小型微型計(jì)算機(jī)系統(tǒng):2013,34(5):1141-1146.

    [7] YUSUKE T,YOSHIFUMIO,SHINJIW, et al. BinaryBased Topology Optimization of Magnetostatic Shielding by a Hybrid Evolutionary Algorithm Combining Genetic Algorithm and ExtendedCompact Genetic Algorithm [J]. Magnetics:2013,49 (5) :2093-2096.

    [8] 王越,許全文,黃麗豐.基于改進(jìn)遺傳算法的連續(xù)函數(shù)優(yōu)化[J].重慶理工大學(xué)學(xué)報(bào):2011,25(2):62-67.

    [9] 闞超豪.多向?qū)W習(xí)自適應(yīng)的粒子群算法[J].計(jì)算機(jī)工程與應(yīng)用:2013,49(6):23-28.

    [10]徐雪松,王四春.基于免疫量子遺傳算法的多峰函數(shù)尋優(yōu)[J].計(jì)算機(jī)應(yīng)用:2012,32(6):1674-1677.

    計(jì)算技術(shù)與自動(dòng)化2014年6月

    第33卷第2期李茂軍等:一種基于狀態(tài)空間模型的進(jìn)化算法

    圖1 基于狀態(tài)空間模型的進(jìn)化算法計(jì)算模型

    2.1 初始種群

    設(shè)種群規(guī)模為N,待優(yōu)化問(wèn)題的變量有M個(gè),則第k代群體X(k)為

    X(k)=x11x12…x1Mx21x22…x2M…xN1xN2…xNM (1) 

    其中Xi=(xi1,xi2,…,xiM),i=1,2,…,N為優(yōu)化問(wèn)題的一組可行解,xij∈[αj,βj],i=1,2,…,N,j=1,2,…,M。

    在算法初始化階段,用隨機(jī)函數(shù)產(chǎn)生N組、每組M個(gè)分別在αj,βj(j=1,2,…,M)上服從均勻分布的實(shí)數(shù)構(gòu)成初始種群X0。

    2.2 狀態(tài)進(jìn)化矩陣

    本算法采用一個(gè)狀態(tài)進(jìn)化矩陣G來(lái)模擬進(jìn)化算法中的重組與變異操作過(guò)程,這使得算法操作起來(lái)非常簡(jiǎn)單,因此,狀態(tài)進(jìn)化矩陣的構(gòu)造是算法的重點(diǎn)。本文按照式(2)構(gòu)造狀態(tài)進(jìn)化矩陣G

    G=g11g12…g1Ng21g22…g2N…gN1gN2…gNN (2)

    其中0≤gij≤1,i,j=1,2,…,N,且∑nj=1gij≤1 ,矩陣G中元素gij(i,j=1,2,…,N)的值是隨機(jī)確定的,本文中G為滿足上述條件的隨機(jī)常數(shù)矩陣,具體構(gòu)造流程如圖2所示。

    2.3 選種池和選擇操作

    同傳統(tǒng)進(jìn)化算法類(lèi)似,選擇操作是一個(gè)擇優(yōu)的過(guò)程,體現(xiàn)了優(yōu)勝劣汰的進(jìn)化思想,本文按最小化優(yōu)化問(wèn)題尋優(yōu)。如圖1所示,從初始群體X(0)開(kāi)始,按照X'(k+1)=GX(k)迭代,可依次得到一系列群體X'(1),X'(2),X'(3),……。群體X'(k+1)和X(k)(k= 0,1,2,…)同時(shí)進(jìn)入選種池,按照適應(yīng)度函數(shù)f()計(jì)算選種池中每個(gè)群體的適應(yīng)度值Y(k),按照從小到大排列,選擇前N個(gè)個(gè)體組成新一代群體X(k+1),如此反復(fù),直到算法滿足停機(jī)條件。

    圖2 狀態(tài)進(jìn)化矩陣構(gòu)造流程圖

    2.4 算法描述

    Step 1. 初始化種群 X(0),并計(jì)算初始種群的適應(yīng)度值;

    Step 2. 根據(jù)圖2流程得到狀態(tài)進(jìn)化矩陣G;

    Step 3. 進(jìn)行進(jìn)化操作,X'(k+1)=G?X(k),并計(jì)算種群X'(k+1)適應(yīng)度值;

    Step 4. 將種群X(k),X'(k+1)同時(shí)放入選種池,按適應(yīng)度值從小到大排列,選擇前N個(gè)個(gè)體作為新一代的群體X(k+1);

    Step 5. 是否滿足終止條件:若是,則結(jié)束;否則,轉(zhuǎn)到Step 3。

    3 仿真實(shí)驗(yàn)與性能分析

    實(shí)驗(yàn)仿真平臺(tái)為Matlab 7.1。本實(shí)驗(yàn)采用了3個(gè)經(jīng)典測(cè)試函數(shù),這些測(cè)試函數(shù)為高維多模函數(shù),具有大量的局部最優(yōu)點(diǎn),是優(yōu)化領(lǐng)域中公認(rèn)的較難優(yōu)化的函數(shù)。

    1)Shubert函數(shù)

    f1=∑5i=1icos [(i+1)x+i]?

    ∑5i=1icos [(i+1)y+i]-10≤x,y≤10 (3)

    有多個(gè)極小值點(diǎn),但只有一個(gè)全局最小值-186.73。

    2)Schwefel's函數(shù)

    f2=-xsin (x)-ysin (y)-500≤x,y≤500(4)

    是一個(gè)具有典型欺騙問(wèn)題的函數(shù),有1個(gè)全局極小值點(diǎn)取值近似為-837.9658,在(420.96 87,420. 9687)處,距離另一個(gè)局部最優(yōu)點(diǎn)很遠(yuǎn),因此如果陷入局部最優(yōu)就很難跳出。

    (3)GoldsteinPrice函數(shù)

    f3=[1+(x+y+1)2(19-14x+3x2-14y+6xy+3y2)]?[30+(2x-3y)2(18-32x+12x2+48y-36xy+27y2)]-2≤x,y≤2(5)

    這是一個(gè)多模函數(shù),有多個(gè)極小值點(diǎn),但只有一個(gè)全局最小值3,極小值點(diǎn)為(0,-1)。

    在測(cè)試中,將種群大小設(shè)置為500,每次運(yùn)行代數(shù)為60代,每個(gè)測(cè)試函數(shù)運(yùn)行測(cè)試50次。表1顯示了函數(shù)f1~f3運(yùn)行50次的實(shí)驗(yàn)結(jié)果,包括平均最優(yōu)值、標(biāo)準(zhǔn)差和平均收斂代數(shù)。全局最優(yōu)值為f1~f3理論最優(yōu)值,平均最優(yōu)值fmav按照公式(6)計(jì)算得到(其中fmi為第i次運(yùn)用本算法求解f1~f3優(yōu)化問(wèn)題得到的最優(yōu)解對(duì)應(yīng)的適應(yīng)度函數(shù)值),標(biāo)準(zhǔn)差為平均最優(yōu)值與全局最優(yōu)值之差,平均收斂代數(shù)為每次實(shí)驗(yàn)找到最優(yōu)解所需要的最少代數(shù)的平均值。

    fmav=∑50i=1fmi/50,i=1,2,…,50(6)

    表1 算法對(duì)3個(gè)測(cè)試函數(shù)50次實(shí)驗(yàn)的結(jié)果

    函數(shù)

    全局最優(yōu)值

    平均最優(yōu)值

    標(biāo)準(zhǔn)差

    平均收斂代數(shù)

    f1

    -186.7316

    -186.7316

    0

    9

    f2

    -837.9658

    -837.9658

    0

    12

    f3

    3.0000

    3.0000

    0

    8

    圖3~5分別顯示了運(yùn)用本算法求解函數(shù)f1~f3優(yōu)化問(wèn)題所得到的最優(yōu)解對(duì)應(yīng)的適應(yīng)度函數(shù)值隨進(jìn)化代數(shù)變化曲線。

    圖3 算法求解f1優(yōu)化問(wèn)題適應(yīng)度函數(shù)值變化曲線 

    圖4 算法求解f2優(yōu)化問(wèn)題適應(yīng)度函數(shù)值變化曲線

    圖5 算法求解f3優(yōu)化問(wèn)題適應(yīng)度函數(shù)值變化曲線

    從表1和圖3~5可以看出,本文提出的算法在3個(gè)函數(shù)優(yōu)化中取得了良好的效果,找到最優(yōu)解的成功率高,收斂速度快。通過(guò)50次實(shí)驗(yàn)結(jié)果可以看出,基于狀態(tài)空間模型的進(jìn)化算法穩(wěn)定性好,成功率高。

    4 結(jié) 論

    本文設(shè)計(jì)了一種基于離散系統(tǒng)狀態(tài)空間模型的進(jìn)化算法,以基本進(jìn)化算法的思想為基礎(chǔ),借鑒矩陣?yán)碚摵蛯?shí)數(shù)編碼方法,通過(guò)構(gòu)造一個(gè)狀態(tài)進(jìn)化矩陣,來(lái)實(shí)現(xiàn)種群的進(jìn)化,提高了算法的可操作性和可靠性。通過(guò)對(duì)3個(gè)經(jīng)典優(yōu)化測(cè)試函數(shù)進(jìn)行優(yōu)化實(shí)驗(yàn),結(jié)果表明:本文提出的算法增強(qiáng)了防止陷入局部最優(yōu)的能力,能夠較大地提高收斂速度和精度,而且穩(wěn)定性能很好。

    雖然算法的仿真結(jié)果是令人滿意的,但是由于實(shí)驗(yàn)條件的局限性和實(shí)驗(yàn)次數(shù)的有限性,算法收斂性和遍歷性還需要通過(guò)進(jìn)一步的數(shù)學(xué)證明。如何構(gòu)造更好的狀態(tài)進(jìn)化矩陣,是需要進(jìn)一步研究的問(wèn)題。

    參考文獻(xiàn)

    [1] 劉鯖潔,陳桂明,劉小方. 基于矩陣編碼的遺傳算法研究[J].計(jì)算機(jī)工程:2011, 37(13):160-162.

    [2] IRINA CIORNEI, ELIAS KYRIAKIDES. Hybrid Ant Colony-Genetic Algorithm (GAAPI) for Global Continuous Optimization [J]. Systems, Man, and Cybernetics, part B, 2012, 42(1), 34-245.

    [3] 董麗麗,龔光紅,李妮,等.基于云模型的自適應(yīng)并行模擬退火遺傳算法[J].北京航空航天大學(xué)學(xué)報(bào):2011,37 (9):1132-1136.

    [4] 劉正龍,楊艷梅.基于遺傳算法的數(shù)值優(yōu)化約束問(wèn)題的研究[J].計(jì)算機(jī)系統(tǒng)應(yīng)用:2013,22(5):139-142.

    [5] 王巍,趙文紅,王宇平.一種有效的解無(wú)約束全局優(yōu)化的進(jìn)化算法[J].控制理論與應(yīng)用:2010,27 (5):570-574.

    [6] 陳曉峰,楊廣明,黃明.一種實(shí)數(shù)編碼的量子差分進(jìn)化算法[J].小型微型計(jì)算機(jī)系統(tǒng):2013,34(5):1141-1146.

    [7] YUSUKE T,YOSHIFUMIO,SHINJIW, et al. BinaryBased Topology Optimization of Magnetostatic Shielding by a Hybrid Evolutionary Algorithm Combining Genetic Algorithm and ExtendedCompact Genetic Algorithm [J]. Magnetics:2013,49 (5) :2093-2096.

    [8] 王越,許全文,黃麗豐.基于改進(jìn)遺傳算法的連續(xù)函數(shù)優(yōu)化[J].重慶理工大學(xué)學(xué)報(bào):2011,25(2):62-67.

    [9] 闞超豪.多向?qū)W習(xí)自適應(yīng)的粒子群算法[J].計(jì)算機(jī)工程與應(yīng)用:2013,49(6):23-28.

    [10]徐雪松,王四春.基于免疫量子遺傳算法的多峰函數(shù)尋優(yōu)[J].計(jì)算機(jī)應(yīng)用:2012,32(6):1674-1677.

    計(jì)算技術(shù)與自動(dòng)化2014年6月

    第33卷第2期李茂軍等:一種基于狀態(tài)空間模型的進(jìn)化算法

    圖1 基于狀態(tài)空間模型的進(jìn)化算法計(jì)算模型

    2.1 初始種群

    設(shè)種群規(guī)模為N,待優(yōu)化問(wèn)題的變量有M個(gè),則第k代群體X(k)為

    X(k)=x11x12…x1Mx21x22…x2M…xN1xN2…xNM (1) 

    其中Xi=(xi1,xi2,…,xiM),i=1,2,…,N為優(yōu)化問(wèn)題的一組可行解,xij∈[αj,βj],i=1,2,…,N,j=1,2,…,M。

    在算法初始化階段,用隨機(jī)函數(shù)產(chǎn)生N組、每組M個(gè)分別在αj,βj(j=1,2,…,M)上服從均勻分布的實(shí)數(shù)構(gòu)成初始種群X0。

    2.2 狀態(tài)進(jìn)化矩陣

    本算法采用一個(gè)狀態(tài)進(jìn)化矩陣G來(lái)模擬進(jìn)化算法中的重組與變異操作過(guò)程,這使得算法操作起來(lái)非常簡(jiǎn)單,因此,狀態(tài)進(jìn)化矩陣的構(gòu)造是算法的重點(diǎn)。本文按照式(2)構(gòu)造狀態(tài)進(jìn)化矩陣G

    G=g11g12…g1Ng21g22…g2N…gN1gN2…gNN (2)

    其中0≤gij≤1,i,j=1,2,…,N,且∑nj=1gij≤1 ,矩陣G中元素gij(i,j=1,2,…,N)的值是隨機(jī)確定的,本文中G為滿足上述條件的隨機(jī)常數(shù)矩陣,具體構(gòu)造流程如圖2所示。

    2.3 選種池和選擇操作

    同傳統(tǒng)進(jìn)化算法類(lèi)似,選擇操作是一個(gè)擇優(yōu)的過(guò)程,體現(xiàn)了優(yōu)勝劣汰的進(jìn)化思想,本文按最小化優(yōu)化問(wèn)題尋優(yōu)。如圖1所示,從初始群體X(0)開(kāi)始,按照X'(k+1)=GX(k)迭代,可依次得到一系列群體X'(1),X'(2),X'(3),……。群體X'(k+1)和X(k)(k= 0,1,2,…)同時(shí)進(jìn)入選種池,按照適應(yīng)度函數(shù)f()計(jì)算選種池中每個(gè)群體的適應(yīng)度值Y(k),按照從小到大排列,選擇前N個(gè)個(gè)體組成新一代群體X(k+1),如此反復(fù),直到算法滿足停機(jī)條件。

    圖2 狀態(tài)進(jìn)化矩陣構(gòu)造流程圖

    2.4 算法描述

    Step 1. 初始化種群 X(0),并計(jì)算初始種群的適應(yīng)度值;

    Step 2. 根據(jù)圖2流程得到狀態(tài)進(jìn)化矩陣G;

    Step 3. 進(jìn)行進(jìn)化操作,X'(k+1)=G?X(k),并計(jì)算種群X'(k+1)適應(yīng)度值;

    Step 4. 將種群X(k),X'(k+1)同時(shí)放入選種池,按適應(yīng)度值從小到大排列,選擇前N個(gè)個(gè)體作為新一代的群體X(k+1);

    Step 5. 是否滿足終止條件:若是,則結(jié)束;否則,轉(zhuǎn)到Step 3。

    3 仿真實(shí)驗(yàn)與性能分析

    實(shí)驗(yàn)仿真平臺(tái)為Matlab 7.1。本實(shí)驗(yàn)采用了3個(gè)經(jīng)典測(cè)試函數(shù),這些測(cè)試函數(shù)為高維多模函數(shù),具有大量的局部最優(yōu)點(diǎn),是優(yōu)化領(lǐng)域中公認(rèn)的較難優(yōu)化的函數(shù)。

    1)Shubert函數(shù)

    f1=∑5i=1icos [(i+1)x+i]?

    ∑5i=1icos [(i+1)y+i]-10≤x,y≤10 (3)

    有多個(gè)極小值點(diǎn),但只有一個(gè)全局最小值-186.73。

    2)Schwefel's函數(shù)

    f2=-xsin (x)-ysin (y)-500≤x,y≤500(4)

    是一個(gè)具有典型欺騙問(wèn)題的函數(shù),有1個(gè)全局極小值點(diǎn)取值近似為-837.9658,在(420.96 87,420. 9687)處,距離另一個(gè)局部最優(yōu)點(diǎn)很遠(yuǎn),因此如果陷入局部最優(yōu)就很難跳出。

    (3)GoldsteinPrice函數(shù)

    f3=[1+(x+y+1)2(19-14x+3x2-14y+6xy+3y2)]?[30+(2x-3y)2(18-32x+12x2+48y-36xy+27y2)]-2≤x,y≤2(5)

    這是一個(gè)多模函數(shù),有多個(gè)極小值點(diǎn),但只有一個(gè)全局最小值3,極小值點(diǎn)為(0,-1)。

    在測(cè)試中,將種群大小設(shè)置為500,每次運(yùn)行代數(shù)為60代,每個(gè)測(cè)試函數(shù)運(yùn)行測(cè)試50次。表1顯示了函數(shù)f1~f3運(yùn)行50次的實(shí)驗(yàn)結(jié)果,包括平均最優(yōu)值、標(biāo)準(zhǔn)差和平均收斂代數(shù)。全局最優(yōu)值為f1~f3理論最優(yōu)值,平均最優(yōu)值fmav按照公式(6)計(jì)算得到(其中fmi為第i次運(yùn)用本算法求解f1~f3優(yōu)化問(wèn)題得到的最優(yōu)解對(duì)應(yīng)的適應(yīng)度函數(shù)值),標(biāo)準(zhǔn)差為平均最優(yōu)值與全局最優(yōu)值之差,平均收斂代數(shù)為每次實(shí)驗(yàn)找到最優(yōu)解所需要的最少代數(shù)的平均值。

    fmav=∑50i=1fmi/50,i=1,2,…,50(6)

    表1 算法對(duì)3個(gè)測(cè)試函數(shù)50次實(shí)驗(yàn)的結(jié)果

    函數(shù)

    全局最優(yōu)值

    平均最優(yōu)值

    標(biāo)準(zhǔn)差

    平均收斂代數(shù)

    f1

    -186.7316

    -186.7316

    0

    9

    f2

    -837.9658

    -837.9658

    0

    12

    f3

    3.0000

    3.0000

    0

    8

    圖3~5分別顯示了運(yùn)用本算法求解函數(shù)f1~f3優(yōu)化問(wèn)題所得到的最優(yōu)解對(duì)應(yīng)的適應(yīng)度函數(shù)值隨進(jìn)化代數(shù)變化曲線。

    圖3 算法求解f1優(yōu)化問(wèn)題適應(yīng)度函數(shù)值變化曲線 

    圖4 算法求解f2優(yōu)化問(wèn)題適應(yīng)度函數(shù)值變化曲線

    圖5 算法求解f3優(yōu)化問(wèn)題適應(yīng)度函數(shù)值變化曲線

    從表1和圖3~5可以看出,本文提出的算法在3個(gè)函數(shù)優(yōu)化中取得了良好的效果,找到最優(yōu)解的成功率高,收斂速度快。通過(guò)50次實(shí)驗(yàn)結(jié)果可以看出,基于狀態(tài)空間模型的進(jìn)化算法穩(wěn)定性好,成功率高。

    4 結(jié) 論

    本文設(shè)計(jì)了一種基于離散系統(tǒng)狀態(tài)空間模型的進(jìn)化算法,以基本進(jìn)化算法的思想為基礎(chǔ),借鑒矩陣?yán)碚摵蛯?shí)數(shù)編碼方法,通過(guò)構(gòu)造一個(gè)狀態(tài)進(jìn)化矩陣,來(lái)實(shí)現(xiàn)種群的進(jìn)化,提高了算法的可操作性和可靠性。通過(guò)對(duì)3個(gè)經(jīng)典優(yōu)化測(cè)試函數(shù)進(jìn)行優(yōu)化實(shí)驗(yàn),結(jié)果表明:本文提出的算法增強(qiáng)了防止陷入局部最優(yōu)的能力,能夠較大地提高收斂速度和精度,而且穩(wěn)定性能很好。

    雖然算法的仿真結(jié)果是令人滿意的,但是由于實(shí)驗(yàn)條件的局限性和實(shí)驗(yàn)次數(shù)的有限性,算法收斂性和遍歷性還需要通過(guò)進(jìn)一步的數(shù)學(xué)證明。如何構(gòu)造更好的狀態(tài)進(jìn)化矩陣,是需要進(jìn)一步研究的問(wèn)題。

    參考文獻(xiàn)

    [1] 劉鯖潔,陳桂明,劉小方. 基于矩陣編碼的遺傳算法研究[J].計(jì)算機(jī)工程:2011, 37(13):160-162.

    [2] IRINA CIORNEI, ELIAS KYRIAKIDES. Hybrid Ant Colony-Genetic Algorithm (GAAPI) for Global Continuous Optimization [J]. Systems, Man, and Cybernetics, part B, 2012, 42(1), 34-245.

    [3] 董麗麗,龔光紅,李妮,等.基于云模型的自適應(yīng)并行模擬退火遺傳算法[J].北京航空航天大學(xué)學(xué)報(bào):2011,37 (9):1132-1136.

    [4] 劉正龍,楊艷梅.基于遺傳算法的數(shù)值優(yōu)化約束問(wèn)題的研究[J].計(jì)算機(jī)系統(tǒng)應(yīng)用:2013,22(5):139-142.

    [5] 王巍,趙文紅,王宇平.一種有效的解無(wú)約束全局優(yōu)化的進(jìn)化算法[J].控制理論與應(yīng)用:2010,27 (5):570-574.

    [6] 陳曉峰,楊廣明,黃明.一種實(shí)數(shù)編碼的量子差分進(jìn)化算法[J].小型微型計(jì)算機(jī)系統(tǒng):2013,34(5):1141-1146.

    [7] YUSUKE T,YOSHIFUMIO,SHINJIW, et al. BinaryBased Topology Optimization of Magnetostatic Shielding by a Hybrid Evolutionary Algorithm Combining Genetic Algorithm and ExtendedCompact Genetic Algorithm [J]. Magnetics:2013,49 (5) :2093-2096.

    [8] 王越,許全文,黃麗豐.基于改進(jìn)遺傳算法的連續(xù)函數(shù)優(yōu)化[J].重慶理工大學(xué)學(xué)報(bào):2011,25(2):62-67.

    [9] 闞超豪.多向?qū)W習(xí)自適應(yīng)的粒子群算法[J].計(jì)算機(jī)工程與應(yīng)用:2013,49(6):23-28.

    [10]徐雪松,王四春.基于免疫量子遺傳算法的多峰函數(shù)尋優(yōu)[J].計(jì)算機(jī)應(yīng)用:2012,32(6):1674-1677.

    精品欧美一区二区三区在线| 日韩一卡2卡3卡4卡2021年| 操出白浆在线播放| 交换朋友夫妻互换小说| 在线观看舔阴道视频| 日韩免费高清中文字幕av| 日本av免费视频播放| 国产精品九九99| 人人妻人人澡人人爽人人夜夜| av有码第一页| 精品国产乱码久久久久久男人| 国产免费现黄频在线看| 国产免费福利视频在线观看| 欧美精品人与动牲交sv欧美| 亚洲激情五月婷婷啪啪| 免费观看a级毛片全部| 精品国产乱码久久久久久小说| 日本一区二区免费在线视频| 亚洲伊人久久精品综合| 国产成人免费无遮挡视频| 狂野欧美激情性bbbbbb| 考比视频在线观看| 免费高清在线观看视频在线观看| 伦理电影免费视频| 一边摸一边做爽爽视频免费| 久久久欧美国产精品| 精品福利永久在线观看| 国产精品 欧美亚洲| 欧美精品啪啪一区二区三区 | 日韩制服丝袜自拍偷拍| 午夜福利在线免费观看网站| 99热网站在线观看| 高清欧美精品videossex| 自线自在国产av| 19禁男女啪啪无遮挡网站| 欧美精品一区二区免费开放| 日日爽夜夜爽网站| 纯流量卡能插随身wifi吗| 人人妻人人爽人人添夜夜欢视频| 嫁个100分男人电影在线观看| 欧美另类亚洲清纯唯美| av免费在线观看网站| 亚洲国产看品久久| 亚洲av国产av综合av卡| 欧美大码av| 精品少妇一区二区三区视频日本电影| 国产真人三级小视频在线观看| 久久女婷五月综合色啪小说| 国产成人欧美在线观看 | 伦理电影免费视频| 麻豆乱淫一区二区| 在线永久观看黄色视频| 欧美成狂野欧美在线观看| 91大片在线观看| 精品少妇一区二区三区视频日本电影| 欧美老熟妇乱子伦牲交| 亚洲黑人精品在线| 中文字幕人妻丝袜制服| 18禁黄网站禁片午夜丰满| 91麻豆av在线| 久久久久精品国产欧美久久久 | 午夜两性在线视频| 国产xxxxx性猛交| 亚洲欧美精品自产自拍| 丰满人妻熟妇乱又伦精品不卡| 欧美 亚洲 国产 日韩一| 无遮挡黄片免费观看| 99国产综合亚洲精品| 国产一区二区在线观看av| 王馨瑶露胸无遮挡在线观看| 男女午夜视频在线观看| 永久免费av网站大全| 一本色道久久久久久精品综合| 丝袜美足系列| 欧美精品高潮呻吟av久久| 久久久欧美国产精品| 亚洲欧美精品自产自拍| 国产一区二区在线观看av| 丰满迷人的少妇在线观看| 久久精品亚洲熟妇少妇任你| 新久久久久国产一级毛片| avwww免费| 国产精品.久久久| 51午夜福利影视在线观看| 两性夫妻黄色片| 女人精品久久久久毛片| 国产精品 欧美亚洲| 亚洲av片天天在线观看| 亚洲五月色婷婷综合| 日韩制服丝袜自拍偷拍| 两个人免费观看高清视频| 国产成人精品在线电影| 一区二区三区精品91| 欧美日韩一级在线毛片| 国产男女内射视频| 啦啦啦在线免费观看视频4| 国产日韩欧美视频二区| 国产欧美亚洲国产| 国产一区二区在线观看av| 国产成+人综合+亚洲专区| 国产一区有黄有色的免费视频| 俄罗斯特黄特色一大片| 国产精品欧美亚洲77777| 又紧又爽又黄一区二区| 精品少妇黑人巨大在线播放| 大香蕉久久成人网| 久久久久久久精品精品| 人人妻,人人澡人人爽秒播| 久久中文字幕一级| av一本久久久久| 亚洲午夜精品一区,二区,三区| 啦啦啦免费观看视频1| 99国产精品一区二区蜜桃av | 这个男人来自地球电影免费观看| 亚洲欧美日韩高清在线视频 | 久久亚洲国产成人精品v| 欧美成狂野欧美在线观看| 久久国产亚洲av麻豆专区| 91成人精品电影| 黑丝袜美女国产一区| 亚洲国产欧美在线一区| 亚洲欧美日韩另类电影网站| 午夜福利在线免费观看网站| www.熟女人妻精品国产| 超碰97精品在线观看| 亚洲精品美女久久久久99蜜臀| 久久精品成人免费网站| av国产精品久久久久影院| 悠悠久久av| 国产男女超爽视频在线观看| 欧美日韩亚洲国产一区二区在线观看 | 欧美成人午夜精品| 在线天堂中文资源库| 老司机午夜十八禁免费视频| 国产精品久久久久成人av| 久久精品国产亚洲av香蕉五月 | 黑丝袜美女国产一区| 中文精品一卡2卡3卡4更新| 成年人午夜在线观看视频| h视频一区二区三区| 欧美亚洲 丝袜 人妻 在线| 久久国产精品人妻蜜桃| 精品亚洲成国产av| 嫁个100分男人电影在线观看| 50天的宝宝边吃奶边哭怎么回事| 日本五十路高清| 纵有疾风起免费观看全集完整版| 青青草视频在线视频观看| 婷婷色av中文字幕| 国产亚洲欧美在线一区二区| √禁漫天堂资源中文www| 久久久精品区二区三区| 少妇裸体淫交视频免费看高清 | 韩国精品一区二区三区| a级毛片在线看网站| 亚洲精品久久午夜乱码| 丝袜喷水一区| 久久久国产成人免费| 亚洲中文日韩欧美视频| 欧美日韩亚洲高清精品| 精品人妻在线不人妻| 美女大奶头黄色视频| 亚洲国产中文字幕在线视频| 亚洲专区中文字幕在线| 中文字幕另类日韩欧美亚洲嫩草| 久久ye,这里只有精品| 亚洲精品中文字幕在线视频| 久久国产精品大桥未久av| www.精华液| 蜜桃国产av成人99| 国产福利在线免费观看视频| 免费一级毛片在线播放高清视频 | 成人免费观看视频高清| 91麻豆av在线| 91字幕亚洲| 男女免费视频国产| 美女福利国产在线| 成年动漫av网址| 如日韩欧美国产精品一区二区三区| 国产av国产精品国产| 精品国产乱码久久久久久小说| 午夜福利视频精品| 80岁老熟妇乱子伦牲交| 免费在线观看黄色视频的| 一本大道久久a久久精品| 国产不卡av网站在线观看| 如日韩欧美国产精品一区二区三区| 国产成人免费无遮挡视频| 涩涩av久久男人的天堂| 亚洲精品久久久久久婷婷小说| av在线播放精品| 巨乳人妻的诱惑在线观看| a在线观看视频网站| 国产免费视频播放在线视频| 午夜福利乱码中文字幕| 欧美日韩福利视频一区二区| 一级黄色大片毛片| 91av网站免费观看| 黄片播放在线免费| 老司机靠b影院| 精品少妇黑人巨大在线播放| 亚洲国产av新网站| 在线观看免费高清a一片| 久久人人97超碰香蕉20202| 国产成人精品久久二区二区91| 国产日韩一区二区三区精品不卡| 新久久久久国产一级毛片| 国产精品久久久久久人妻精品电影 | 两个人看的免费小视频| 国产日韩欧美亚洲二区| 亚洲av国产av综合av卡| 亚洲精品国产av成人精品| 欧美黑人精品巨大| 国产麻豆69| 久久人人爽人人片av| 久久精品人人爽人人爽视色| 黄片播放在线免费| 午夜福利一区二区在线看| 日韩有码中文字幕| 黄网站色视频无遮挡免费观看| 美女高潮到喷水免费观看| www.精华液| 韩国精品一区二区三区| 亚洲人成电影观看| 欧美激情久久久久久爽电影 | 欧美亚洲 丝袜 人妻 在线| 精品少妇内射三级| 高清视频免费观看一区二区| 一级片'在线观看视频| 天堂俺去俺来也www色官网| 国产免费现黄频在线看| 青青草视频在线视频观看| 国产日韩欧美亚洲二区| 一区在线观看完整版| 老司机影院成人| 国产一区二区三区综合在线观看| 高清av免费在线| 亚洲精品国产区一区二| 永久免费av网站大全| 日本av手机在线免费观看| 无遮挡黄片免费观看| 国产精品偷伦视频观看了| 黑人操中国人逼视频| 91麻豆精品激情在线观看国产 | 搡老岳熟女国产| 久久久国产成人免费| 精品福利观看| 欧美xxⅹ黑人| 50天的宝宝边吃奶边哭怎么回事| 国产成人系列免费观看| 操美女的视频在线观看| 亚洲 国产 在线| 色精品久久人妻99蜜桃| 免费日韩欧美在线观看| 色94色欧美一区二区| 国产人伦9x9x在线观看| 成年动漫av网址| 老司机午夜十八禁免费视频| 久久久精品区二区三区| 视频在线观看一区二区三区| 色婷婷久久久亚洲欧美| 午夜福利免费观看在线| 中国国产av一级| 国产日韩欧美在线精品| 久久这里只有精品19| 啦啦啦啦在线视频资源| 免费看十八禁软件| 曰老女人黄片| 国产伦人伦偷精品视频| 大片免费播放器 马上看| 久热这里只有精品99| 日本wwww免费看| 国产视频一区二区在线看| 91av网站免费观看| 女人爽到高潮嗷嗷叫在线视频| a级片在线免费高清观看视频| 极品人妻少妇av视频| 蜜桃国产av成人99| 久久九九热精品免费| 欧美国产精品va在线观看不卡| 99国产精品免费福利视频| 日本猛色少妇xxxxx猛交久久| 嫩草影视91久久| 在线观看免费日韩欧美大片| 又紧又爽又黄一区二区| 久久天躁狠狠躁夜夜2o2o| 国产在视频线精品| 啦啦啦 在线观看视频| 成人影院久久| 在线观看免费视频网站a站| 性少妇av在线| 日韩制服丝袜自拍偷拍| 一进一出抽搐动态| 999久久久国产精品视频| 国产成人影院久久av| 国产精品熟女久久久久浪| 国产三级黄色录像| 丰满迷人的少妇在线观看| 亚洲精品久久成人aⅴ小说| 国产区一区二久久| av在线app专区| 两性夫妻黄色片| 亚洲精品一卡2卡三卡4卡5卡 | 窝窝影院91人妻| 1024香蕉在线观看| 伦理电影免费视频| 久久久欧美国产精品| 高清av免费在线| 国产在线一区二区三区精| 嫁个100分男人电影在线观看| av不卡在线播放| 黑丝袜美女国产一区| 一级毛片精品| 亚洲精品国产精品久久久不卡| 久久久精品94久久精品| 一区二区日韩欧美中文字幕| 国产激情久久老熟女| 国产成人系列免费观看| 久久久久国产一级毛片高清牌| 纯流量卡能插随身wifi吗| 超色免费av| 手机成人av网站| 久久人妻福利社区极品人妻图片| 亚洲av成人一区二区三| 一个人免费在线观看的高清视频 | avwww免费| 久久久久视频综合| 狠狠精品人妻久久久久久综合| 久久中文字幕一级| 欧美日韩亚洲国产一区二区在线观看 | 精品国产一区二区三区久久久樱花| 夫妻午夜视频| 搡老岳熟女国产| 日韩大片免费观看网站| 好男人电影高清在线观看| 久久香蕉激情| 亚洲精品美女久久av网站| 国产日韩一区二区三区精品不卡| 国产主播在线观看一区二区| 国产亚洲午夜精品一区二区久久| 在线观看一区二区三区激情| 男人操女人黄网站| 欧美日韩亚洲高清精品| 国产欧美日韩一区二区精品| 欧美亚洲日本最大视频资源| 国产av又大| 亚洲中文字幕日韩| 日本av免费视频播放| 国产亚洲av高清不卡| 亚洲精品中文字幕一二三四区 | 国产无遮挡羞羞视频在线观看| 伊人亚洲综合成人网| 欧美国产精品va在线观看不卡| 日韩视频一区二区在线观看| 久久精品国产综合久久久| 国产亚洲欧美在线一区二区| 大片免费播放器 马上看| 捣出白浆h1v1| 国产1区2区3区精品| www.av在线官网国产| 国产精品一区二区免费欧美 | 一级毛片女人18水好多| 啦啦啦 在线观看视频| av国产精品久久久久影院| 人妻一区二区av| 法律面前人人平等表现在哪些方面 | 成人黄色视频免费在线看| 精品少妇一区二区三区视频日本电影| 一级,二级,三级黄色视频| 九色亚洲精品在线播放| 中国美女看黄片| 一二三四在线观看免费中文在| 久久人妻熟女aⅴ| 日韩免费高清中文字幕av| 久久精品人人爽人人爽视色| 人妻人人澡人人爽人人| 亚洲精品中文字幕一二三四区 | 欧美精品高潮呻吟av久久| 爱豆传媒免费全集在线观看| 亚洲欧美色中文字幕在线| 精品久久蜜臀av无| 日韩精品免费视频一区二区三区| 成人国产av品久久久| 精品国产一区二区三区四区第35| 亚洲国产欧美在线一区| 搡老岳熟女国产| 午夜激情久久久久久久| 国产真人三级小视频在线观看| 18禁裸乳无遮挡动漫免费视频| 久久精品国产a三级三级三级| 精品国产乱码久久久久久小说| 老司机靠b影院| 午夜精品久久久久久毛片777| 一区二区三区激情视频| 最黄视频免费看| 狠狠精品人妻久久久久久综合| 别揉我奶头~嗯~啊~动态视频 | 一进一出抽搐动态| 丁香六月天网| av网站在线播放免费| 美女午夜性视频免费| 一级片'在线观看视频| 国产又爽黄色视频| 成人国语在线视频| 久久久久久人人人人人| a级毛片在线看网站| 一级毛片精品| 免费在线观看完整版高清| 久久久久精品国产欧美久久久 | 女人精品久久久久毛片| 日韩 欧美 亚洲 中文字幕| 亚洲欧洲日产国产| 久久天堂一区二区三区四区| 久久99热这里只频精品6学生| 人人妻人人爽人人添夜夜欢视频| 久久国产精品人妻蜜桃| 大香蕉久久网| 搡老熟女国产l中国老女人| 51午夜福利影视在线观看| a级毛片黄视频| 亚洲精品自拍成人| 欧美成人午夜精品| 免费高清在线观看视频在线观看| 高清欧美精品videossex| 日本黄色日本黄色录像| 亚洲avbb在线观看| 日本五十路高清| 一进一出抽搐动态| 亚洲国产欧美日韩在线播放| 精品福利永久在线观看| 大片免费播放器 马上看| 十八禁网站网址无遮挡| 亚洲中文日韩欧美视频| 中文字幕人妻丝袜制服| 国产精品99久久99久久久不卡| 久久99一区二区三区| 日本欧美视频一区| 日本vs欧美在线观看视频| 一个人免费看片子| 国产亚洲午夜精品一区二区久久| 蜜桃在线观看..| 一级毛片女人18水好多| 国产精品av久久久久免费| 欧美精品啪啪一区二区三区 | 99久久99久久久精品蜜桃| 51午夜福利影视在线观看| 亚洲专区国产一区二区| 久久久国产成人免费| 午夜免费成人在线视频| 欧美中文综合在线视频| 国产精品二区激情视频| 亚洲欧美一区二区三区久久| av网站免费在线观看视频| 97精品久久久久久久久久精品| 免费高清在线观看视频在线观看| 免费在线观看影片大全网站| 国产成人免费观看mmmm| 美女中出高潮动态图| 国产精品一区二区在线不卡| 国产一区二区激情短视频 | 满18在线观看网站| 韩国高清视频一区二区三区| 久久九九热精品免费| 欧美+亚洲+日韩+国产| 18禁裸乳无遮挡动漫免费视频| 久久久久久人人人人人| 欧美人与性动交α欧美精品济南到| 婷婷色av中文字幕| 日本一区二区免费在线视频| 久久久久久人人人人人| 别揉我奶头~嗯~啊~动态视频 | 成人三级做爰电影| 高清视频免费观看一区二区| 午夜福利视频精品| 久久人妻熟女aⅴ| 亚洲精品成人av观看孕妇| 如日韩欧美国产精品一区二区三区| 国产精品国产三级国产专区5o| 亚洲黑人精品在线| 五月开心婷婷网| 亚洲伊人久久精品综合| 美女视频免费永久观看网站| 日本a在线网址| 午夜91福利影院| www.av在线官网国产| 成年女人毛片免费观看观看9 | 蜜桃国产av成人99| 丝袜美腿诱惑在线| 在线观看免费高清a一片| 91麻豆精品激情在线观看国产 | 窝窝影院91人妻| 十八禁网站免费在线| 人人澡人人妻人| a级片在线免费高清观看视频| 2018国产大陆天天弄谢| 新久久久久国产一级毛片| 久久精品久久久久久噜噜老黄| 热re99久久精品国产66热6| 国产精品 国内视频| 久久久久网色| 肉色欧美久久久久久久蜜桃| 老司机在亚洲福利影院| 精品免费久久久久久久清纯 | 欧美精品人与动牲交sv欧美| 欧美精品啪啪一区二区三区 | 法律面前人人平等表现在哪些方面 | 女人久久www免费人成看片| 两个人看的免费小视频| av视频免费观看在线观看| 人人妻人人澡人人看| 亚洲欧洲精品一区二区精品久久久| 久久午夜综合久久蜜桃| 高清视频免费观看一区二区| 50天的宝宝边吃奶边哭怎么回事| 日韩精品免费视频一区二区三区| 18禁裸乳无遮挡动漫免费视频| 亚洲全国av大片| 操美女的视频在线观看| 精品人妻1区二区| 可以免费在线观看a视频的电影网站| 精品人妻1区二区| 十分钟在线观看高清视频www| 乱人伦中国视频| 国产色视频综合| 啦啦啦啦在线视频资源| 一级,二级,三级黄色视频| 国产精品一区二区精品视频观看| 亚洲欧美精品自产自拍| av不卡在线播放| 男女下面插进去视频免费观看| 欧美日韩国产mv在线观看视频| 久久精品人人爽人人爽视色| 亚洲av成人一区二区三| 叶爱在线成人免费视频播放| 美国免费a级毛片| 久久人妻福利社区极品人妻图片| 在线观看舔阴道视频| 中文字幕另类日韩欧美亚洲嫩草| 在线观看舔阴道视频| 久久人妻熟女aⅴ| 亚洲午夜精品一区,二区,三区| 久久精品国产亚洲av高清一级| 亚洲精品国产av成人精品| 精品国产一区二区三区久久久樱花| 丝袜美足系列| 叶爱在线成人免费视频播放| 99国产精品一区二区三区| 日韩视频一区二区在线观看| 国产精品1区2区在线观看. | 日韩一区二区三区影片| 精品国产乱码久久久久久男人| 老司机深夜福利视频在线观看 | 不卡一级毛片| av又黄又爽大尺度在线免费看| 黄片小视频在线播放| 中文字幕av电影在线播放| 90打野战视频偷拍视频| 欧美人与性动交α欧美软件| 麻豆国产av国片精品| 国产真人三级小视频在线观看| 中文精品一卡2卡3卡4更新| 在线观看一区二区三区激情| 欧美另类亚洲清纯唯美| 亚洲av日韩精品久久久久久密| 亚洲精品国产av成人精品| 黄频高清免费视频| 国产精品久久久久久精品古装| av片东京热男人的天堂| 高清视频免费观看一区二区| 亚洲精品国产av成人精品| 亚洲精品日韩在线中文字幕| 老熟女久久久| 久久国产精品人妻蜜桃| 久久人妻福利社区极品人妻图片| 免费av中文字幕在线| 日本精品一区二区三区蜜桃| 精品高清国产在线一区| 18禁观看日本| 久久久久久亚洲精品国产蜜桃av| 悠悠久久av| 操出白浆在线播放| 成人国产一区最新在线观看| 亚洲美女黄色视频免费看| 超碰97精品在线观看| 久久免费观看电影| 日韩一区二区三区影片| 久久久久久久久久久久大奶| 成人av一区二区三区在线看 | 曰老女人黄片| 91国产中文字幕| 国产深夜福利视频在线观看| 国产精品亚洲av一区麻豆| avwww免费| 黄色怎么调成土黄色| 国产一卡二卡三卡精品| 高清在线国产一区| 久久精品亚洲熟妇少妇任你| 叶爱在线成人免费视频播放| 曰老女人黄片| 国产成人影院久久av| 国产无遮挡羞羞视频在线观看| 欧美日韩视频精品一区| 亚洲精品国产av蜜桃| 久久久精品国产亚洲av高清涩受| 亚洲成人免费av在线播放| 男女国产视频网站| 亚洲久久久国产精品| 精品卡一卡二卡四卡免费| 啪啪无遮挡十八禁网站| 丝袜在线中文字幕| 成人国语在线视频| av视频免费观看在线观看| 黄色a级毛片大全视频| 亚洲人成77777在线视频| 在线十欧美十亚洲十日本专区|