av一区二区在线观看_亚洲男人的天堂网站_日韩亚洲视频_在线成人免费_欧美日韩精品免费观看视频_久草视

您的位置:首頁技術文章
文章詳情頁

Python實現(xiàn)粒子群算法的示例

瀏覽:3日期:2022-06-27 18:55:53

粒子群算法是一種基于鳥類覓食開發(fā)出來的優(yōu)化算法,它是從隨機解出發(fā),通過迭代尋找最優(yōu)解,通過適應度來評價解的品質(zhì)。

Python實現(xiàn)粒子群算法的示例

PSO算法的搜索性能取決于其全局探索和局部細化的平衡,這在很大程度上依賴于算法的控制參數(shù),包括粒子群初始化、慣性因子w、最大飛翔速度和加速常數(shù)與等。

PSO算法具有以下優(yōu)點:

不依賴于問題信息,采用實數(shù)求解,算法通用性強。

需要調(diào)整的參數(shù)少,原理簡單,容易實現(xiàn),這是PSO算法的最大優(yōu)點。

協(xié)同搜索,同時利用個體局部信息和群體全局信息指導搜索。

收斂速度快, 算法對計算機內(nèi)存和CPU要求不高。

更容易飛越局部最優(yōu)信息。對于目標函數(shù)僅能提供極少搜索最優(yōu)值的信息,在其他算法無法辨別搜索方向的情況下,PSO算法的粒子具有飛越性的特點使其能夠跨過搜索平面上信息嚴重不足的障礙,飛抵全局最優(yōu)目標值。比如Generalized Rosenbrock函數(shù)全局最小值在原占附近.但是此函數(shù)全局最優(yōu)值與可到達的局部最優(yōu)值之間右一條獨長的山路,曲面山谷中點的最速下降方向幾乎與到函數(shù)最小值的最佳方向垂直,找到全局最小值的可能性微乎其微, 但是PSO算法完全有可能找到全局最優(yōu)值。

同時, PSO算法的缺點也是顯而易見的:

算法局部搜索能力較差,搜索精度不夠高。

算法不能絕對保證搜索到全局最優(yōu)解。

PSO算法設計的具體步驟如下: 初始化粒子群(速度和位置)、慣性因子、加速常數(shù)、最大迭代次數(shù)、算法終止的最小允許誤差。 評價每個粒子的初始適應值。 將初始適應值作為當前每個粒子的局部最優(yōu)值,并將各適應值對應的位置作為每個粒子的局部最優(yōu)值所在的位置。 將最佳初始適應值作為當前全局最優(yōu)值,并將最佳適應值對應的位置作為全局最優(yōu)值所在的位置。 依據(jù)公式更新每個粒子當前的飛翔速度。 對每個粒子的飛翔速度進行限幅處理,使之不能超過設定的最大飛翔速度。 依據(jù)公式更新每個粒子當前所在的位置。 比較當前每個粒子的適應值是否比歷史局部最優(yōu)值好,如果好,則將當前粒子適應值作為粒子的局部最優(yōu)值,其對應的位置作為每個粒子的局部最優(yōu)值所在的位置。 在當前群中找出全局最優(yōu)值,并將當前全局最優(yōu)值對應的位置作為粒子群的全局最優(yōu)值所在的位置。 重復步驟(5)~(9),直到滿足設定的最小誤差或最大迭代次數(shù) 輸出粒子群的全局最優(yōu)值和其對應的位置以及每個粒子的局部最優(yōu)值和其對應的位置。

本文中我們假設要求解一個維度為10的向量,這里的適應度函數(shù)采用簡單的線性誤差求和。

#基本粒子群算法#vi+1 = w*vi+c1*r1*(pi-xi)+c2*r2*(pg-xi) 速度更新公式#xi+1 = xi + a*vi+1 位置更新公式(一般a=1)#w = wmax -(wmax-wmin)*iter/Iter 權重更新公式#iter當前迭代次數(shù) Iter最大迭代次數(shù) c1、c2學習因子 r1、r2隨機數(shù) pi粒子當前最優(yōu)位置 pg粒子群全局最優(yōu)#初始化 wmax=0.9 wmin=0.4 通常c1=c2=2 Iter對于小規(guī)模問題(10,20)對于大規(guī)模(100,200)#算法優(yōu)劣取決于w、c1和c2,迭代結束的條件是適應度函數(shù)的值符合具體問題的要求#初始化粒子群,包括尺寸、速度和位置#本算法假設想要的輸出是長度為10的矩陣,y=[1.7]*10,適應度函數(shù)f(x)= |x-y| <=0.001符合要求import numpy as npswarmsize = 500partlen = 10wmax,wmin = 0.9,0.4c1 = c2 = 2Iter = 400def getwgh(iter): w = wmax - (wmax-wmin)*iter/Iter return wdef getrange(): randompv = (np.random.rand()-0.5)*2 return randompvdef initswarm(): vswarm,pswarm = np.zeros((swarmsize,partlen)),np.zeros((swarmsize,partlen)) for i in range(swarmsize): for j in range(partlen): vswarm[i][j] = getrange() pswarm[i][j] = getrange() return vswarm,pswarmdef getfitness(pswarm): pbest = np.zeros(partlen) fitness = np.zeros(swarmsize) for i in range(partlen): pbest[i] = 1.7 for i in range(swarmsize): yloss = pswarm[i] - pbest for j in range(partlen): fitness[i] += abs(yloss[j]) return fitnessdef getpgfit(fitness,pswarm): pgfitness = fitness.min() pg = pswarm[fitness.argmin()].copy() return pg,pgfitnessvswarm,pswarm = initswarm()fitness = getfitness(pswarm)pg,pgfit = getpgfit(fitness,pswarm)pi,pifit = pswarm.copy(),fitness.copy()for iter in range(Iter): if pgfit <= 0.001: break #更新速度和位置 weight = getwgh(iter) for i in range(swarmsize): for j in range(partlen): vswarm[i][j] = weight*vswarm[i][j] + c1*np.random.rand()*(pi[i][j]-pswarm[i][j]) + c2*np.random.rand()*(pg[j]-pswarm[i][j]) pswarm[i][j] = pswarm[i][j] + vswarm[i][j] #更新適應值 fitness = getfitness(pswarm) #更新全局最優(yōu)粒子 pg,pgfit = getpgfit(fitness,pswarm) #更新局部最優(yōu)粒子 for i in range(swarmsize): if fitness[i] < pifit[i]: pifit[i] = fitness[i].copy() pi[i] = pswarm[i].copy()for j in range(swarmsize): if pifit[j] < pgfit: pgfit = pifit[j].copy() pg = pi[j].copy()print(pg)print(pgfit)

下面的結果分別是迭代300次和400次的結果。

Python實現(xiàn)粒子群算法的示例

可以看到400次迭代雖然適應度沒有達到預期,得到的向量已經(jīng)很接近期望的結果了。

寫在最后:粒子群算法最重要的參數(shù)就是慣性權重和學習因子,針對這兩個參數(shù)有了新的優(yōu)化粒子群算法(IPSO)。還有初始化粒子群時速度和位置范圍的確定,包括種群的大小和迭代次數(shù)的選擇,這些都是‘摸著石頭過河’,沒有標準答案。

以上就是Python實現(xiàn)粒子群算法的示例的詳細內(nèi)容,更多關于Python 粒子群算法的資料請關注好吧啦網(wǎng)其它相關文章!

標簽: Python 編程
相關文章:
主站蜘蛛池模板: 亚洲成人免费观看 | 日韩精品一区二区三区免费观看 | 神马久久久久久久久久 | 国产成人一区 | 欧洲毛片| 日本在线看片 | 国产99在线 | 欧美 | 自拍偷拍亚洲欧美 | 日日噜噜夜夜爽爽狠狠 | 羞羞视频一区二区 | 亚洲在线久久 | 日本精品一区二区三区视频 | 美女黄色在线观看 | 亚洲二区在线 | 亚洲黄色成人网 | 精品1区2区| 狠狠躁天天躁夜夜躁婷婷老牛影视 | 天天操操 | 男女搞网站 | 午夜专区 | 欧美日韩综合 | 偷拍自拍在线观看 | 九九99九九精彩46 | 日韩成人在线免费观看 | 午夜伦4480yy私人影院 | 国产一级一级国产 | 日韩欧美电影在线 | 亚洲一级毛片 | 亚洲成人日韩 | 男女爱爱网站 | 欧美色性 | 91亚洲精品在线 | 免费永久av| 中文字幕国产在线 | 激情毛片 | 精品久久久久一区二区国产 | 羞羞在线视频 | 精品欧美二区 | 日韩精品中文字幕一区二区三区 | 精国产品一区二区三区四季综 | 亚洲一区在线日韩在线深爱 |