一、Python中傳遞函數

def Test(a,b):

print(a,b)

class TEST:

def __init__(self, FitnessFunction):

self.FitnessFunction = FitnessFunction

def initial(self):

self.FitnessFunction(1,2)

if __name__ == '__main__':

fa = TEST(Test)

fa.initial()

二、常用優化算法測試函數及實現

這里學到的兩個小技巧是：

(1)np.linalg.norm()來求一范式和二范式

(2)累加、累成的便捷寫法：reduce(lambda x,y: x+y, x_)，當x_為[x1, x2, x3, x4]時，reduce這個函數相當于執行f(f(f(x1,x2),x3),x4)，其中f(a, b)是由lambda創建的。

參考資料：https://zhidao.baidu.com/question/397090143281977845.html

def f1(x_): #Sphere

return np.linalg.norm(x_) ** 2

def f2(x_): #Schwefel 2.22

xa = reduce(lambda x,y: x*y, abs(x_))

xb = np.linalg.norm(x_, ord=1)

return xa + xb

def f3(x_): #Rosenbrock

xa = 100 * (x_[1:] - x_[:-1] * x_[:-1])

xb = (x_[:-1] - 1) * (x_[:-1] - 1)

return reduce(lambda x,y: x+y, xa+xb)

def f4(x_): #Step

return np.linalg.norm(np.floor(x_ + 0.5))

def f5(x_): #Schwefel 2.26

x_new = x_ * np.sin(np.sqrt(abs(x_)))

return reduce(lambda x,y: x+y, x_new)

def f6(x_): #Rastrigin

x_new = (x_ * x_) - (10 * np.cos(2 * np.pi * x_)) + 10

return reduce(lambda x,y: x+y, x_new)

def f7(x_): #Ackley

D = len(x_)

x_a = 20 * np.exp(-0.2 * np.sqrt(np.linalg.norm(x_)) / np.sqrt(D))

x_b = np.exp(reduce(lambda x,y: x+y, np.cos(2 * np.pi * x_)) / D)

return (20 - x_a) + (np.e - x_b)

def f8(x_): #Girewank

x_a = (np.linalg.norm(x_) ** 2)/4000

x_b = reduce(lambda x,y: x*y,

np.cos(x_ / np.array(range(1,len(x_)+1))))

return x_a - x_b + 1

這8個函數是目前主流的優化算法測試函數。

參考文獻：[1]張哲辰,劉三陽.基于拓撲改進與交叉策略的螢火蟲算法[J].計算機工程與應用,2019,55(07):1-8.

三、前人提出的螢火蟲改進算法

目前知網最新的對螢火蟲算法從拓撲結構上改進的文章是《基于拓撲改進與交叉策略的螢火蟲算法》，它提出了以馮諾依曼拓撲結構來組織螢火蟲，并以此確定更新路徑時的對象，但是我覺得描述得比較模糊，沒辦法代碼實現，所以這里不考慮。

主要參考這兩篇：

[1]??Wang H，Wang W，Sun H，et al.Firefly algorithm with random attraction[J].International Journal of Bio-Inspired Computation，2016，?8 (1)：33-41.

[2] Wang H，Wang W，Zhou X，et al.Firefly algorithm with neighborhood attraction[J].Information Sciences，2017， 382/383：374-387.

總結

以上是生活随笔為你收集整理的python萤火虫算法_FA：萤火虫算法的测试及Python实现的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。