【正文】 于可重慶大學(xué)本科學(xué)生畢業(yè)設(shè)計 （論文） 3 PSO 算法的實現(xiàn) 11 視化顯示的目的是用 此效果來反應(yīng)粒子的運動情況所以具體實現(xiàn)這里就不作具體說明了。 int iResult = 0。 } return iResult。 i()()。 j++) { SwarmAdjust[i].iPbestPosition[j] = SwarmAdjust[i].iPos[j]。 for (i=0。 } 4．方法名稱： AdjustBirdPosition 方法功能： 調(diào)整粒子群中各個粒子位置 int SwarmBird::AdjustBirdPosition(void) { int iResult = 0。 //對微粒位置進行調(diào)整 for (i=0。 kDIM。 //最大速度限制 if (SwarmAdjust[i].iVec[k] MAX_V) { SwarmAdjust[i].iVec[k] = MAX_V。 } 重慶大學(xué)本科學(xué)生畢業(yè)設(shè)計 （論文） 4 基于 PSO 算法的函數(shù)優(yōu)化 15 4 基于 PSO 算法的函數(shù)優(yōu)化 引 言 隨著現(xiàn)代科技和生產(chǎn)的發(fā)展，人們在解決一個工程應(yīng)用或經(jīng)濟應(yīng)用問題時，總是希望得到一個最優(yōu)的方案。在第二次世界大戰(zhàn)中，由于軍事上的需要產(chǎn)生了運籌學(xué)，提出了大量不能用古典方法解決的最優(yōu)化問題，從而產(chǎn)生了如線性規(guī)劃、非線性規(guī)劃、動態(tài)規(guī)劃、網(wǎng)絡(luò)流規(guī)劃等新的方法。因此，要尋求最優(yōu)的設(shè)計與決策，以獲得最好的經(jīng)濟或技術(shù)效果，這就促使最優(yōu)化技術(shù)迅速發(fā)展，為最優(yōu)化技術(shù)的發(fā)展提供了有力的工具 [18]。特別，若約束集為 nDR? ，則最優(yōu)化問題稱為無約束最優(yōu)化問題 重慶大學(xué)本科學(xué)生畢業(yè)設(shè)計 （論文） 4 基于 PSO 算法的函數(shù)優(yōu)化 16 min ( )nxRfx? (4. 2) 約束 最優(yōu)化問題通常寫為 m in ( ). . ( ) 0 1 , 2 , ,( ) 0 1 , 2 , ,ijfxs t g x i mh x j l???? (4. 3) 當(dāng)目標函數(shù)和約束函數(shù)中到少有一個是 x 的非線性函數(shù)，式 和式 稱為非線性規(guī)劃。求解線性規(guī)劃的方法中的單純形算法是一種通用的有效算法，于 1947 年由 首先提出， 50 多年的計算實踐表明，它不僅是求解線性規(guī)劃的基本方法而且是整數(shù)規(guī)劃和非線性規(guī)劃某些算法的基礎(chǔ)。 PSO 算法求進行函數(shù) 優(yōu) 化問題 PSO 算法在求解復(fù)雜的函數(shù)優(yōu)化 問題時對函數(shù)的要求相當(dāng)寬，并不要求函數(shù)有解析解、可微，甚至連續(xù)。三個權(quán)重因子，主要影響收斂速度方面。為了在窗口中方便顯示數(shù)據(jù)，把式 修改為 21( ) ( 2 0 0 )niif x x???? (4. 6) 我們先對 參數(shù) 進行測試， 群體規(guī)模 25，迭代次 為 數(shù) 300。 由于式 我們可以看出，其最小值點為（ 200， 200）。 其中各個色素塊代表一個粒子，粒子的初始位置為圖 中選取的位置，每條線代表一個粒子的運動軌跡?，F(xiàn)我們選取幾個典型非線性函數(shù)來運用算法求其解。 圖 Rastrigrin 函數(shù)測試一 需要指出的是， 在粒子數(shù)量不多的時候，位置會稍微有些偏差， 一般取 2040后效果就會很好。 2． Rosenbrock 函數(shù)，也叫香蕉函數(shù)（ Banana ） 2211( ) (1 0 0 ( ) ( 1 ) )ni i iif x x x x??? ? ? ?? (4. 9) 與對 Rastrigrin 函數(shù)一樣，對 Rosenbrock 作平移處理，得式 。其最小值 st val ue ? ， 最小值 位置 gb e st po st i on =( 50 2, 30 2, 10 3) 3． Schaffe 函數(shù) 2 2 22 2 2( s in ) 0 . 5( ) 0 . 5 (1 . 0 0 . 0 0 1 ( ) )xyfx xy???? ?? (4. 11) 平移后的函數(shù)。 PSO 算法與演化規(guī)劃 (EP)很相似。 PSO 基本都能更快地達到全局優(yōu)值， PSO 基本不 受問題峰數(shù)增加的影響，受問題維數(shù)的影 響也很小。對三個測試函數(shù)分別用程序進行測試，結(jié)果表明微粒群算法的正確性和它在函數(shù)優(yōu)化方面優(yōu)點。 感謝陳世林同學(xué)提供電腦，讓我能順利完成論文的書寫。 胡老師 淵博的學(xué)識、嚴謹?shù)闹螌W(xué)態(tài)度、敏銳的學(xué)術(shù)思想、高尚的為師品德無時無刻不在影響著我，使我受益匪淺。該算法為人們提供了如下一種思路 ：使智慧出現(xiàn)而不是努力強迫它；模擬自然而不是力圖控制它；尋求使事情簡單化而不是讓它復(fù)雜。也就是說， PSO 算法執(zhí)行 一種有“意識 (Conscience)”的變異。 重慶大學(xué)本科學(xué)生畢業(yè)設(shè)計 （論文） 4 基于 PSO 算法的函數(shù)優(yōu)化 24 圖 Schaffe 函數(shù) 優(yōu)化 情況 小 結(jié) PSO 最直接的應(yīng)用就是多元函數(shù)的優(yōu)化問題，包括帶約束的優(yōu)化問題。 圖 二維 Rosenbrock 函數(shù) 優(yōu)化 效果圖 三維情況如圖 所示。 020040060080002004006008000200400600800xR a s t r i g r i n 函數(shù)三維情況時粒子位置圖yz 圖 Rastrigrin 函數(shù)三維情況 重慶大學(xué)本科學(xué)生畢業(yè)設(shè)計 （論文） 4 基于 PSO 算法的函數(shù)優(yōu)化 22 其中 每 條連續(xù)的線表示一個粒子在三維空間中的一條運動軌跡。 1( ) ( ( 2 0 0 ) 1 0 c o s ( 2 ( 2 0 0 ) ) 1 0 )niiif x x x??? ? ? ? ?? (4. 8) 重慶大學(xué)本科學(xué)生畢業(yè)設(shè)計 （論文） 4 基于 PSO 算法的函數(shù)優(yōu)化 20 用 我 們的程序?qū)κ? 進行求解，如圖 所示， 從 圖 中可知 在二維情況下，最小值是 m in ( ) 9. 83 39 96fx ?? ，位置是（ 200, 200）；重新選取初始化粒子的測試結(jié)果如圖 所示。 圖 4. 3 程序運行結(jié)果 其中 速度 參數(shù)為： m ax m in5, 5。 圖 程序初始化注極值點 然后用鼠標在程序窗口中隨機點取點， 作 為初始化時粒子的位置 ， 如圖 所示 。 表 4. 1: 參數(shù)測試表 ? 1c 2c 全局最佳值 全局最佳位置 維數(shù) （ 201， 200） 2 （ 200， 200） 2 （ 200， 200） 2 （ 201， 201， 201） 3 （ 200， 200， 201） 3 （ 200， 201， 200） 3 從表 的測試結(jié)果，我們可以看出，當(dāng)權(quán)重系數(shù)取后兩組時得出較好的結(jié)果 。 我們只對無約束情況進行討論。 要想達到較好的測試效果，在測試時，還可能 需要 對粒子最大 速度進行適重慶大學(xué)本科學(xué)生畢業(yè)設(shè)計 （論文） 4 基于 PSO 算法的函數(shù)優(yōu)化 17 當(dāng)?shù)恼{(diào)整。 常用的求解無約束非 線性規(guī)劃的方法有：一維搜索法、梯度法、 NEWTON 法、擬 NEWTON 法、共軛方向法和 POWELL 方法等；求解約束非線性規(guī)劃的方法有：罰函數(shù)法與乘子法、可行方向法、二次 規(guī)劃法等。 此外，根據(jù)決策變量、目標函數(shù)和要求不同，最優(yōu)化分類為整數(shù)規(guī)劃、動態(tài)規(guī)劃、網(wǎng)絡(luò)規(guī)劃、幾何規(guī)劃、多目標規(guī)劃等若干分支。所謂的最佳表現(xiàn)為一個目標函數(shù)在滿足一定約束條件下的極大或極小。特別是 20 世紀 50 年代以來，由于 近代科技與生產(chǎn)發(fā)展的需要和計算機技術(shù)的飛速發(fā)展，使得最優(yōu)化技術(shù)得到了迅速發(fā)展，現(xiàn)最優(yōu)化技術(shù)已成為一門得到廣泛應(yīng)用的新興學(xué)科。隨著計算機科學(xué)的發(fā)展和應(yīng)用，應(yīng)用最優(yōu)化方法最解決問題領(lǐng)域不斷擴大，最優(yōu)化的理論和方法不斷發(fā)展。 } //位置調(diào)整 SwarmAdjust[i].iPos[k] = SwarmAdjust[i].iPos[k] + SwarmAdjust[i].iVec[k]。 fRand2 = (double)rand() / (double)RAND_MAX。 i++) { int k =0。 double fRand2 = 。 i++) { if (SwarmAdjust[i].dPbestValue SwarmAdjust[gbest].dPbestValue) { gbest = i。 } 3．方法名稱： GetGobleParticleValue 方法功能： 取得全局最優(yōu)值，并且把它放到 gbest 中 重慶大學(xué)本科學(xué)生畢業(yè)設(shè)計 （論文） 3 PSO 算法的實現(xiàn) 13 int SwarmBird::GetGobleParticleValue(void) { int iResult = 0。 //給局部最佳位置賦值 for (int j=0。 int iResult = 0。 i()()。 d l g )4 : E v a l u a t e P a r t i cl e V a l u e ( v o i d )6 : G e t P a r t P a r t i cl e V a l u e ( v o i d )7 : G e t G o b l e P a r t i cl e V a l u e ( v o i d )8 : A d j u st B i r d P o si t i o n ( v o i d )5 : E v a l u a t e P a r t i cl e V a l u e ( v o i d )2 : I n i t ( ) 圖 微粒群系統(tǒng)時序圖 具體實現(xiàn) 現(xiàn)在將 SwarmBird 中的關(guān)鍵方法實現(xiàn)代碼做簡要說明。在實現(xiàn) PSO 算法時， 其時序圖見圖 所示。dlg) ：設(shè)置權(quán)重參數(shù)值 EvaluateParticleValue(void) ：計算每個粒子的評價值 GetPartParticleValue(void) ：設(shè)置每個粒子的局部最優(yōu)值 GetGobleParticleValue(void) ：設(shè)置微粒群的全局最優(yōu)值索引 AdjustBirdPosition(void) ：調(diào)整微粒群中的微粒的位置 Clear(void) ：清除數(shù)據(jù) 算法實現(xiàn)流程 在面向?qū)ο缶幊虝r，我們用 UML(統(tǒng)一 建模語 言 )來控制程序的開發(fā)過程。其類圖如圖 所示。 重慶大學(xué)本科學(xué)生畢業(yè)設(shè)計 （論文） 3 PSO 算法的實現(xiàn) 9 3 PSO 算法的實現(xiàn) 算法實現(xiàn) 分析 PSO 算法 的提出是 美國社會心理學(xué) 家 James Kennedy 和電氣工程師 Russell Eberhar 受他們早期對許多鳥類的群體行為 進行建模與仿真研究結(jié)果 的啟發(fā)。 2． 權(quán)重因子 在 PSO算法中有 3個權(quán)重 因子：慣性權(quán)重 ? ，加速常數(shù) 1c 和 2c 。 關(guān)于 參數(shù)總結(jié) 如下： 1． 最大速度 maxv ，最小速度 minv maxv 決定當(dāng)前位置與最好位置之間的區(qū)域的分辨率 (或精度 )。 文獻中 提出 PSO 的異步模式 ， 采用 Java 的多線程技術(shù)實現(xiàn) ， 把每個粒子的行為看成為一個獨立的線程或進程 ， 運行中的粒子充分表現(xiàn)出高度的獨立性 ， 而在種群層次上表現(xiàn)為異步性 。 Juang 提出將 GA 與