萬有引力搜索算法

1、15.2 5.2 比較比較PSOPSO與與RGARGA 我們應(yīng)用GSA到這些最小化函數(shù)，并且比較RGA與PSO所得出的結(jié)果，在這些情況，群體大小設(shè)置為 。維數(shù) ，表1,2的函數(shù)最大迭代次數(shù)為1000，表3為500，在PSO中 ，慣性因素 從0.9到0.2線性下降，在RGA應(yīng)用算法交叉算子，高斯交換以及輪盤算法，交叉和變異的概率分別為0.3和0.1. 在GSA,G用(21)式表示， 為100， 。 是所有迭代次數(shù) (21)50N30n221 cc0G20TTteGtG0)(3近幾年，多種啟發(fā)式優(yōu)化方法得到發(fā)展，這些方法中很多是根據(jù)自然中群體行為得到啟示。本節(jié)課介紹一種基于萬有引力定律和質(zhì)量相互作

2、用的新的優(yōu)化算法引力搜索算法。引力搜索算法在2009年被首次提出，是一種基于萬有引力定律和牛頓第二定律的種群優(yōu)化算法。該算法通過種群的粒子位置移動來尋找最優(yōu)解，即隨著算法的循環(huán)，粒子靠它們之間的萬有引力在搜索空間內(nèi)不斷運動，當(dāng)粒子移動到最優(yōu)位置時，最優(yōu)解便找到了。4. 啟發(fā)式算法回顧. 萬有引力定律. 引力搜索算法（GSA）. 比較研究. 實驗結(jié)果. 引力搜索算法的研究展望5 Heuristic是希臘語，意為“啟發(fā)式”。啟發(fā)式是尋找好的（近似最佳）解的技術(shù)。對于那些受大自然的運行規(guī)律或者面向具體問題的經(jīng)驗、規(guī)則啟發(fā)出來的方法，人們常常稱為啟發(fā)式算法。啟發(fā)式算法是相對于最優(yōu)化算法提出的。很多實際

3、的最優(yōu)化問題的計算是復(fù)雜的。因此，解決這樣問題的實際方法是運用啟發(fā)式算法，這樣可以在合理的計算時間內(nèi)找到一個近似最優(yōu)解。啟發(fā)式算法可以這樣定義：一個基于直觀或經(jīng)驗構(gòu)造的算法，在可接受的花費（計算時間和空間）下給出解決組合優(yōu)化問題每一個實例的一個可行解該可行解與最優(yōu)解的偏離程度一般不能被預(yù)計。（Heuristic algorithms）. 啟發(fā)式算法6啟發(fā)式算法模擬物理或生物過程，例如一些著名的算法，遺傳算法（GA）、模擬退火算法（SA）、蟻群算法（ACO）粒子群優(yōu)化算法（PSO）和細菌覓食算法（BFA）。GA靈感來自于達爾文進化論；SA利用熱力作用設(shè)計；ACO模擬螞蟻覓食行為；BFA來自于搜索

4、和最佳覓食細菌；PSO模擬鳥群的行為。上述提到的啟發(fā)式算法都是隨機行為。然而，F(xiàn)ormato提出了基于引力運動的確定性的啟發(fā)式搜索算法，中心引力優(yōu)化（CFO）。中心引力優(yōu)化算法是根據(jù)物理運動學(xué)的模型建立的一個新型的優(yōu)化算法，通過初始化若干隨機質(zhì)點，進行迭代，直至找到最優(yōu)解。7在一些隨機算法中，像模擬退火算法（SA）搜索開始于一個單一的初始點，并且以一個連續(xù)的方式繼續(xù)。然而，大多數(shù)啟發(fā)式搜索算法用多個初始點以并行方式搜索。例如，群為基礎(chǔ)的算法使用類似于自然的鳥群或者魚群的一系列代理。在一個以群為基礎(chǔ)的算法，每一個體施行一系列的特殊運算，并且分享這些信息給其他個體。這些操作大部分很簡單，然而它們的

5、集體效應(yīng)，稱為群體智能，會產(chǎn)生令人驚訝的結(jié)果。代理之間的局部相互作用提供了一個全局結(jié)果，它允許系統(tǒng)解決問題不需要應(yīng)用任何的中央控制器。這種情況下，個體操作包括隨機搜索、正反饋、負反饋和多元相互作用，進行自組織。群體智能指許多簡單個體通過相互合作產(chǎn)生復(fù)雜智能行為的特性。8我們可以在人群為基礎(chǔ)的啟發(fā)式算法識別兩個常見問題：勘探和開采勘探和開采?？碧接袛U大搜索空間的能力，開采有尋找最佳解決方案能力。在第一次迭代中，啟發(fā)式搜索算法勘探搜索空間尋找新的解。為了避免陷入局部最優(yōu)的陷阱，該算法必須在前幾次迭代中使用勘探。因此，在以人群為基礎(chǔ)的啟發(fā)式算法，勘探是一個重要的問題。通過勘探和開采，算法調(diào)整自己的半

6、最優(yōu)點。要有高性能的搜索，關(guān)鍵點是一個合適的勘探和開采之間的權(quán)衡。然而，所有的以人群為基礎(chǔ)的啟發(fā)式搜索算法采用的勘探和開采方面，他們使用不同的方法和操作。換句話說，所有的搜索算法有一個共同的框架。9從不同的角度來看，一個以群為基礎(chǔ)的搜索算法的個體在每次迭代中通過三個步驟來實現(xiàn)勘探和開采概念:自適應(yīng)，合作和競爭自適應(yīng)，合作和競爭。在自我調(diào)整的步驟，每個個體（代理）提高其性能。在合作中，個體彼此合作形成的信息傳遞。最后，在競爭的一步，個體競爭生存。這些步驟通常隨機形成，可以用不同的方式來實現(xiàn)。這些步驟從自然的啟發(fā)，是以人群為基礎(chǔ)的啟發(fā)式算法的思想。這些概念，引導(dǎo)算法尋找全局最優(yōu)。然而，一個算法在解

7、決一些問題是好的，在解決另外一些問題則不行。因此，提出高性能的新啟發(fā)式算法是非常受歡迎的。我們的目標是建立一個新的考慮到所提到的方面和基于引力規(guī)則的以群為基礎(chǔ)的搜索算法。10. 萬有引力定律萬有引力定律萬有引力定律是Newton于1687年在自然哲學(xué)的數(shù)學(xué)原理上提出的，萬有引力定律解釋物體之間相互作用關(guān)系的定律，是物體間由于它們的引力質(zhì)量而引起的相互吸引力所遵循的規(guī)律。自然界中任何兩個物體都是相互吸引的，萬有引力普遍存在于任意兩個有質(zhì)量的物體之間。萬有引力定律表示如下：自然界中任何兩個物體都是相互吸引的，引力的大小和這兩個物體的質(zhì)量的乘積成正比，和它們之間距離平方成反比。數(shù)學(xué)表達式為：221R

8、MMGF （1）其中，F(xiàn)表示兩個物體間的引力大小，G表示萬有引力常數(shù)，M1，M2分別表示兩個物體的質(zhì)量，R表示兩個物體之間的距離。11牛頓第二定律牛頓第二定律：當(dāng)一個力F作用在一個質(zhì)子上，它的加速度a依賴于力和它的質(zhì)量M：MFa（2）根據(jù)（1）和（2），增加兩個質(zhì)子之間的距離意味著減少他們之間的萬有引力。此外，由于引力減少的影響，引力常數(shù)的實際值依賴于宇宙的實際時間，方程（3）給出了降低引力常數(shù)與時間關(guān)系： 1,00tttGtG（3）其中G(t)是在時間t引力常數(shù)的值，G(t0)是在t0時萬有引力常數(shù)。12理論物理學(xué)中定義三種質(zhì)量：主動引力質(zhì)量，主動引力質(zhì)量，aMpMiM是對物體重力場的強度的

9、測量，小的主動物體的重力場比大的主動引力質(zhì)量的重力場弱。小的被動被動引力質(zhì)量的物體比大的被動引力質(zhì)量的物體改變快。引力質(zhì)量的被動引力質(zhì)量，被動引力質(zhì)量，是對物體重力場中物體相互作用的強度的測量，受到的力小。慣性質(zhì)量慣性質(zhì)量是當(dāng)有一個力作用在物體，改變她位置的移動的力的測量，大的慣性質(zhì)量的物體改變它移動的更慢， 小的慣性13考慮到以上提到的三種質(zhì)量定義，我們重新定義牛頓定律。萬有引力Fij通過物體j作用在物體i，與j 的主動引力質(zhì)量和i 被動引力質(zhì)量乘積成正比，與它們之間距離成反比。ia與Fij成正比，與i 慣性質(zhì)量成反比，則（1）（2）式重新定義如下：2RMMGFpiajijiiijiMFa

10、（4）（5）ajM和piM分別代表質(zhì)子i 的主動引力質(zhì)量，質(zhì)子j 的被動引力質(zhì)量，iiM代表質(zhì)子i 的慣性質(zhì)量。雖然，慣性質(zhì)量，主動引力質(zhì)量，被動引力質(zhì)量有概念上的區(qū)別，但是沒有實驗可以清楚的證明它們之間的不同。追溯到廣義相對論上，假設(shè)慣性質(zhì)量和被動引力質(zhì)量是等價的，這就是所知道的弱等價原理。斯坦-達爾德廣義相對論也假定慣性質(zhì)量和主動引力質(zhì)量是等價的，這個等價有時稱為強等價原理。14在圖1中，F(xiàn)1j是作用在M1到Mj的力，F(xiàn)1是作用在M1和產(chǎn)生加速度 的所有力。1a圖圖115雖然，慣性質(zhì)量，主動引力質(zhì)量，被動引力質(zhì)量有概念上的區(qū)別，但是沒有實驗可以清楚的證明它們之間的不同。追溯到廣義相對論上，

11、假設(shè)慣性質(zhì)量和被動引力質(zhì)量是等價的，這就是所知道的弱等價原理。斯坦-達爾德廣義相對論也假定慣性質(zhì)量和主動引力質(zhì)量是等價的，這個等價有時稱為強等價原理。16.引力搜索算法（GSA）受萬有引力定律啟發(fā)，提出了一種新型群體智能優(yōu)化算法引力搜索算法。引力搜索算法在求解優(yōu)化問題時，搜索個體的位置和問題的解相對應(yīng)，并且還要考慮個體質(zhì)量。個體質(zhì)量用于評價個體的優(yōu)劣，位置越好，質(zhì)量越大。由于引力的作用，個體之間相互吸引并且朝著質(zhì)量較大的個體方向移動，個體運動遵循牛頓第二定律。隨著運動的不斷進行，最終整個群體都會聚集在質(zhì)量最大個體的周圍，從而找到質(zhì)量最大的個體，而質(zhì)量最大個體占據(jù)最優(yōu)位置。因此，算法可以獲得問題

12、的最優(yōu)解。在GSA，每個代理有4個規(guī)格：位置，慣性質(zhì)量，主動引力質(zhì)位置，慣性質(zhì)量，主動引力質(zhì)量和被動引力質(zhì)量量和被動引力質(zhì)量。每個個體的位置對應(yīng)一個問題的解決方法，它們的引力和慣性質(zhì)量確定應(yīng)用的適應(yīng)度函數(shù)。換句話說，每個個體呈現(xiàn)一個解決方法，并且算法通過適當(dāng)?shù)恼{(diào)節(jié)引力和慣性質(zhì)量。17引力搜索算法屬于群體智能優(yōu)化算法，而群體智能優(yōu)化算法最顯著的特點是強調(diào)個體之間的相互作用。這里，相互作用可以是個體間直接或間接的通信。在引力搜索算法中，萬有引力相當(dāng)于是一種信息傳遞的工具，實現(xiàn)個體間的優(yōu)化信息共享，整個群體在引力的作用下進行優(yōu)化搜索。信息的交互過程不僅在群體內(nèi)部傳播了信息，而且群體內(nèi)所有個體都能處理

13、信息，并根據(jù)其所得到的信息改變自身的搜索行為，這樣就能使得整個群體涌現(xiàn)出一些單個個體所不具備的能力和特性，也就是說，在群體中，個體行為雖然簡單，但是個體通過得到的信息相互作用以解決全局目標，信息在整個群體的傳播使得問題能夠比由單個個體求解更加有效的獲得解決。183.1 算法的模型引力搜索算法首先在解空間和速度空間分別對位置和速度進行初始化，其中位置表示問題的解。例如，d維空間中的第i 個搜索個體的位置和速度分別表示為),(1nidiiixxxX),(1nidiiivvvV，dixdiv分別表示個體i 在第d 維的位置分量和速度其中，和分量。通過評價各個個體的目標函數(shù)值，確定每個個體的質(zhì)量和受到

14、的引力，計算加速度，并更新速度和位置。19（1）計算質(zhì)量個體i 的質(zhì)量定義如下：)()()()()(tworsttbesttworsttfittmii（6）NijiitmtmtM1)()()(（7）其中，)(tfiti和)(tMi分別表示在第t 次迭代時第i 個個體的best(t)和worst(t)表示在第t次迭代時所有個體中最優(yōu)適應(yīng)度函數(shù)值和最差適應(yīng)度函數(shù)值，對最小化問題，其定義如下：)(min)(,2,1tfittbestjNj（8）)(max)(, 2, 1tfittworstjNj（9）適應(yīng)度函數(shù)值和質(zhì)量。20（2）計算引力算法源于對萬有引力定律的模擬，但不拘泥于物理學(xué)中的萬有引力公式

15、的精確表達式。在第d 維上，個體j 對個體i 的引力定義如下：)()()()()()()(txtxtRtMtMtGtFdidjijajpidij（10）其中，G(t)表示在t 次迭代時萬有引力常數(shù)的取值， G(t)=G(G0，t),Rij(t)表示個體i和j 之間的歐氏距離，2)(),()(tXtXtRjiij 是一常數(shù)，防止分母為零。21在第d維上，個體i 所受的合力為：NijkbestjdijjditFrandtF,)()(（11）其中，randj表示在0,1之間服從均勻分布的一個隨機變量kbest表示個體質(zhì)量按降序排在前k個的個體，并且k的取值隨迭代次數(shù)線性減小，初值為N，終值為1。22

16、（3）計算加速度根據(jù)牛頓第二定律，個體i 在第d維的加速度方程為：)()()(tMtFtaiididi（12）（4）更新速度和位置)()() 1(tatvrandtvdidiidi（13）) 1()() 1(tvtxtxdididi（14）其中，r表示在0,1之間服從均勻分布的一個隨機變量。23引力搜索算法的目的并不是為了模擬萬有引力定律，而是利用萬有引力定律的特點去解決優(yōu)化問題。算法受萬有引力定律的啟發(fā)，但是不拘泥于萬有引力公式的原始表達式。在算法中引力與兩個個體質(zhì)量的乘積成正比和它們的距離成反比的優(yōu)化效果更好。此外，萬有引力常數(shù)不在固定不變，而是隨著迭代次數(shù)單調(diào)遞減，算法的優(yōu)化效果更好。在

17、計算個體受到的萬有引力合力時，算法只考慮質(zhì)量較大的個體產(chǎn)生的引力。因為在引力搜索算法中，當(dāng)引力較大時，或者有質(zhì)量較大的個體，或者兩個個體間的距離較小。質(zhì)量大的個體占據(jù)較優(yōu)的位置，并且代表較好的解。算法僅考慮來自質(zhì)量較大的個體的引力，可以消除因距離較小而引力較大的影響，引導(dǎo)其他個體向質(zhì)量較大的個體方向移動。在引力的不斷作用下，整個群體逐漸向質(zhì)量較大的個體方向逼近，最終搜索到問題的最優(yōu)解。243.2 算法的流程基本引力搜索算法主要包含三個組成部分：群體初始化、計算個體質(zhì)量和所受的引力以及個體運動更新。算法的主要實現(xiàn)步驟如下：步驟1 隨機初始化群體中各個體的位置，個體的初始速度為零步驟2 個體最適值

18、（個體的適應(yīng)度函數(shù)值）步驟3 更新G(t),best(t),worst(t),Mi(t)步驟4 計算個體所受到的合力步驟5 計算加速度和速度步驟6 更新個體位置步驟7 若滿足終止條件，則輸出當(dāng)前結(jié)果并終止算法，否則轉(zhuǎn)向步驟2.25上述程序的結(jié)構(gòu)流程如圖2所示。圖226對引力搜索算法的特點做如下總結(jié)：（1）每個搜索個體都賦予4個狀態(tài)變量，分別為位置、速度、加速度和質(zhì)量。位置用于表示位置的解，速度用于更新位置，加速度用于更新速度質(zhì)量用于評價個體的優(yōu)劣。（2）整個群體總是尋找質(zhì)量最大的個體，無論是最大值優(yōu)化問題還是最小值優(yōu)化問題，都可以通過質(zhì)量函數(shù)的定義，將優(yōu)化目標轉(zhuǎn)換為搜索質(zhì)量最大的個體。（3）從

19、引力搜索算法設(shè)計的起源來看，算法主要是對萬有引力定律的模擬，是將物理原理和優(yōu)化思想相結(jié)合而產(chǎn)生的。引力搜索算法最顯著的特點是整個群體依靠個體之間的引力作用進行尋優(yōu)，引力相當(dāng)于一種優(yōu)化信息的傳遞工具，根據(jù)算法的設(shè)計，個體的質(zhì)量越大，引力也越大。因此，在引力作用下，整個群體能夠向質(zhì)量最大的個體方向移動，從而能夠搜索到問題的最優(yōu)解。（4）引力搜索算法的流程簡單，參數(shù)設(shè)置少，可以很好的和各種優(yōu)化問題相結(jié)合，易于實現(xiàn)。27除了上述這些特點之外，引力搜索算法也具有智能優(yōu)化算法一些共同的特點。例如，引力搜索算法對目標函數(shù)沒有特別要求，不要求函數(shù)具有連續(xù)和可導(dǎo)等數(shù)學(xué)性質(zhì)，甚至有時連有沒有解析表達式都不做要求，

20、而且對問題中不確定的信息具有一定的適應(yīng)能力，因此，算法的通用性比較強。此外，從算法實現(xiàn)的方法來看，引力搜索算法可以采用串行或者并行的方法實現(xiàn)，可以根據(jù)具體問題，設(shè)計出合理的算法實現(xiàn)方法。284. 比較研究 4.1 粒子群算法（PSO）PSO啟發(fā)于模擬社會行為，這種優(yōu)化方法更新粒子群，通過操作者根據(jù)從環(huán)境中獲得的最適信息，為了群個體可以移向最優(yōu)解。在PSO中， 和 的計算如下：) 1()() 1(tvtxtxdididi（15）)()()()() 1(2211txgbestrctxpbestrctvtwtvdidididiididi（16）dixdivri1,ri2是0,1范圍的隨機變量，c1,

21、c2是位置常數(shù)，w是慣性權(quán)重。pbesti表示第i 個質(zhì)子的最佳位置，gbest表示群中所有質(zhì)子的最佳位置。29從（16）式，我們發(fā)現(xiàn)每個質(zhì)子嘗試更新它的位置（Xi），用當(dāng)前位置和第i 個質(zhì)子最佳位置pbesti之間的距離以及當(dāng)前位置與gbest之間的距離。)()()()() 1(2211txgbestrctxpbestrctvtwtvdidididiididi30 GSA vs PSO GSA和PSO的優(yōu)化在搜索空間通過個體移動獲得，然而移動規(guī)則是不同，一些重要的不同如下：a. PSO代理的方向計算僅用兩個最佳位置pbesti和gbest，GSA的基于整體合力的所有其他代理獲得。b. PSO

22、應(yīng)用一種存儲器來更新速度（pbesti,gbest）,然而，GSA無存儲，在更新中僅需要代理的當(dāng)前位置起作用。c. PSO執(zhí)行更新不需要考慮解之間的距離；GSA的力與解之間距離成反比。d. 兩個算法的搜索理念是不同的，PSO模擬鳥的行為，GSA源于物理現(xiàn)象。314.2 CFO 算法 中心引力優(yōu)化CFO是確定性多維搜索算法。它的模型是穿越搜索空間重力影響下的探針。在開始時，初始探位置用一個確定方式計算。對于初始化，根據(jù)式（17），在零時刻每一個坐標軸的位置矢量排列充滿均勻分布的探針，d=1.n， p=1.nNnNdpi) 1( 1)(1()0(minmaxminnNxxpxxddddi（17），

23、fiti 是探針i 的適應(yīng)度函數(shù)，它必須最大化。n是維數(shù)，N是探針數(shù)量，在CFO，每一次迭代，探針進行評估。M用適用度函數(shù)計算，即式（18）， 是t時刻探針i 的質(zhì)量。)(tiMiifittM)(（18）32隨后，加速度更新應(yīng)用式（19），一個新位置計算應(yīng)用式（20），)(tadi是時間t探針i 的加速度)(txdi是時間t探針i 的位置是兩個常數(shù)。NijjijdidjijijditRtxtxtMtMtMtMUGta, 1)()()()()()()()(（19）)(21)() 1(tatxtxdididi（20）G是引力常數(shù)，Rij(t)是t時刻探針i和j 之間的歐氏距離。U是單位階躍函數(shù)33

24、GSA vs CFO兩者的探索位置和加速度都來源于在引力場中的質(zhì)點運動，但它們使用不同的構(gòu)想（1）其中一個主要的不同是CFO是固有的確定性的并且沒有應(yīng)用任何隨機參數(shù)在它的構(gòu)想，GSA是隨機搜索算法。（2）加速度和運動表現(xiàn)以及群體計算，GSA不同于CFO。（3）CFO初始探針位置分布是系統(tǒng)的（基于確定性的規(guī)則），在算法收斂有很大影響。GSA初始分布是隨機的。（4）CFO的G是常數(shù)，GSA中是控制參數(shù)。34. 實驗結(jié)果實驗結(jié)果 5.1 基準函數(shù) 表1表3中的基準函數(shù)是測試實驗所用到的基準函數(shù)。在這些表中，n代表函數(shù)的維數(shù)，S是 的子集。表1和表2中函數(shù)除了 之外最小值都是0， 的最小值為 并且除了

25、 , 和 以外，它們的最優(yōu)位置 都為 ， 和 的最優(yōu)位置 為 ， 的最優(yōu)位置 為 表1，單峰測試函數(shù)nR8F8Fn982.418-5F8F13FoptX n05F13FoptX n18FoptXn96.42035表2，高維多峰測試函數(shù)表3，多峰低維測試函數(shù)36 三個算法應(yīng)用到基準函數(shù)，結(jié)果如下： （1）單峰高維函數(shù) 到 是單峰高維函數(shù)，這種情況下，因為有其他方法來專門設(shè)計優(yōu)化單峰函數(shù)，因此單峰函數(shù)搜索算法的收斂速度比最終結(jié)果更重要。 在表4中顯示，GSA對所有函數(shù)的運行結(jié)果要比RGA和PSO要好，GSA的收斂速度可由圖3,4得出，根據(jù)這些圖表，GSA比其他算法更快的找到全局最優(yōu)，因此GSA有較

26、高的收斂速度。從表4的結(jié)果顯示，除了F5之外，基于權(quán)值的引力優(yōu)化算法對其他的4個基準函數(shù)的搜索結(jié)果明顯好于引力搜索算法的搜索結(jié)，仿真效果如圖3,4所示 1F7F37表表4 4高維單峰函數(shù)最小值搜索結(jié)果（函數(shù)維數(shù)n為30，最大迭代次數(shù)為1000）38n=30時，GSA,PSO,RGA對最小化優(yōu)化結(jié)果的比較1F圖3，圖4n=30時，GSA,PSO,RGA對最小化優(yōu)化結(jié)果的比較7F39 （2 2）多峰高維函數(shù)）多峰高維函數(shù) 多峰函數(shù)有很多局部極小點并且?guī)缀醵己茈y優(yōu)化，我們進行了在 至 上的局部極小值以指數(shù)形式增加的實驗，函數(shù)維數(shù)設(shè)置為30，平均適應(yīng)度函數(shù)是最佳解的中間值，最后一次迭代函數(shù)在表5中顯示

27、，對 ， 和 ，GSA的表現(xiàn)比其他的算法更好，而對函數(shù) ，GSA無法自身調(diào)整已取得理想的好的結(jié)果。 8F13F10F12F13F8F40表表5 5測試高維多峰函數(shù)最小值搜索結(jié)果（函數(shù)的維數(shù)n為30，最大迭代次數(shù)為1000）41圖5n=30時，GSA,PSO,RGA對最小化優(yōu)化結(jié)果的比較圖6n=30時，GSA,PSO,RGA對最小化優(yōu)化結(jié)果的比較10F12F42（3 3）多峰低維函數(shù)）多峰低維函數(shù) 表6表示了表3的GSA,RGA,PSO的多峰低維函數(shù)基準函數(shù)建的比較，在圖7和圖8，結(jié)果表明RGA,PSO和GSA有相同解并且大部分性能相同表6，表表3 3基準函數(shù)的最小化結(jié)果基準函數(shù)的最小化結(jié)果，（

28、函數(shù)的維數(shù)n為如表1中所示，最大迭代次數(shù)為1000）43圖7n=4時，GSA,PSO,RGA對最小化優(yōu)化結(jié)果的比較15F圖8n=4時，GSA,PSO,RGA對最小化優(yōu)化結(jié)果的比較22F445.3 5.3 與與CFOCFO比較比較 先來比較GSA與CFO，在相同條件下是很難比較區(qū)分這兩種算法。 CFO是一個確定性的算法，有很多性質(zhì)并依賴于初始群的生成和群的規(guī)模，特別是隨著問題的復(fù)雜性與維數(shù)增加,CFO對于群規(guī)模和初始條件更敏感。而且，在這種條件下，CFO必須以一個較大的群規(guī)模以獲得一個可接受的，合理的結(jié)果。這也就意味著，在用CFO解決問題前，我們應(yīng)知道一些先驗信息來建立群數(shù)量或多次嘗試運行算法來

29、獲得合適的群值。 而GSA是一種隨機搜索算法，一種可以廣泛的應(yīng)用于一個固定的小規(guī)模人口問題的優(yōu)化算法。正是由于上述原因，在相同條件下不能比較GSA與CFO對高維函數(shù)的優(yōu)化。因此，我們在低維問題上比較這兩種算法。45 在CFO，在步驟0位置矢量充滿了在，每個坐標軸的探針均勻分布。在GSA，初始值是隨機的，代理數(shù)和迭代次數(shù)的最大值設(shè)置分別為60和500，（因為CFO需要一個大規(guī)模種群，故取代理次數(shù)為60）。GSA的參數(shù)設(shè)置在前一節(jié)已給出，對于CFO， 。注意到，我們在比較GSA和CFO對Formato 提出的函數(shù)時，需要一些適用于CFO的先驗函數(shù)信息來進行函數(shù)的優(yōu)化。2, 2 表表7 7表示，表示

30、，CFOCFO的隨機初始化結(jié)果（不是均勻分布的初始的隨機初始化結(jié)果（不是均勻分布的初始化）平均超過化）平均超過3030分，如表分，如表7 7所示，所示，CFOCFO在所得結(jié)果不能呈現(xiàn)在所得結(jié)果不能呈現(xiàn)隨機初始化時好的結(jié)果，并對初始化條件有重要影響：對隨機初始化時好的結(jié)果，并對初始化條件有重要影響：對于三個函數(shù)每一個算法的性能在圖于三個函數(shù)每一個算法的性能在圖9-119-11中表示，這些結(jié)果中表示，這些結(jié)果顯示，除了函數(shù)顯示，除了函數(shù) ， 外，外，GSAGSA比比CFOCFO有更好的解。有更好的解。5F8F11F46表表7 7，表，表1-31-3的一些函數(shù)的最小化結(jié)果，最大迭代次數(shù)為的一些函數(shù)的

31、最小化結(jié)果，最大迭代次數(shù)為5005001F圖9. n=5時，GSA,PSO,RGA對隨機最小化的比較47圖10，n=5時，GSA,PSO,RGA對隨機最小化的比較圖11，n=2時，GSA,PSO,RGA對隨機最小化的比較5F16F48 對于高維函數(shù)的優(yōu)化，CFO應(yīng)用大量的搜索來運行，由于CFO的確定性，它在最初的幾個迭代收斂。表8總結(jié)概括了對代理次數(shù)N不同的30維的研究結(jié)果，在表4-6，通過比較平均最值的結(jié)果，我們可以看出，除了 ， ， ，GSA提供更好的解法，CFO在前幾次迭代中收斂，在局部最優(yōu)并不能自調(diào)整。5F11F14F表表8 對表1-3的CFO運行優(yōu)化結(jié)果495.4 5.4 結(jié)論結(jié)論 近幾年，已經(jīng)研究了多種啟發(fā)式優(yōu)化方法，有些優(yōu)化方法的靈感來自于自然群集行為，本文中總結(jié)了一種新的優(yōu)化算