差分进化算法
差分进化算法(Differential Evolution Algorithm,DE)是一种高效的全局优化算法。它也是基于群体的启发式搜索算法,群中的每个个体对应一个解向量。
简单来说DE算法可用于求函数的极值点,例如:函数 的最小值点在,最小值为0。
差分进化算法可用于资源调度,笔者在嵌入式容器云系统(私有云)有相关研究。
可以简单地理解DE为球场的勒布朗,三分能力并不是最厉害的,扣篮也不是最厉害的,打球也不是最好看的,身体素质也不是最好的,但是,综合能力却是其他几位顶级无法比拟的。所以DE在很多领域应用很多。
具体的变异因子之类的概念都可以在网上搜到,就不叙述了。
但是需要注意变异因子和交叉因子,两者是影响差分进化的关键因素,并且可以通过多两者的改进来弥补差分进化陷入局部最优解以及收敛速度缓慢的缺点。
以下是关于差分进化的Python源码:
import numpy as np import random class Population: def __init__(self, min_range, max_range, dim, factor, rounds, size, object_func, CR=0.75): self.min_range = min_range self.max_range = max_range self.dimension = dim self.factor = factor self.rounds = rounds self.size = size self.cur_round = 1 self.CR = CR self.get_object_function_value = object_func # 初始化种群 self.individuality = [np.array([random.uniform(self.min_range, self.max_range) for s in range(self.dimension)]) for tmp in range(size)] self.object_function_values = [self.get_object_function_value(v) for v in self.individuality] self.mutant = None def mutate(self): self.mutant = [] for i in range(self.size): r0, r1, r2 = 0, 0, 0 while r0 == r1 or r1 == r2 or r0 == r2 or r0 == i: r0 = random.randint(0, self.size-1) r1 = random.randint(0, self.size-1) r2 = random.randint(0, self.size-1) tmp = self.individuality[r0] + (self.individuality[r1] - self.individuality[r2]) * self.factor for t in range(self.dimension): if tmp[t] > self.max_range or tmp[t] < self.min_range: tmp[t] = random.uniform(self.min_range, self.max_range) self.mutant.append(tmp) def crossover_and_select(self): for i in range(self.size): Jrand = random.randint(0, self.dimension) for j in range(self.dimension): if random.random() > self.CR and j != Jrand: self.mutant[i][j] = self.individuality[i][j] tmp = self.get_object_function_value(self.mutant[i]) if tmp < self.object_function_values[i]: self.individuality[i] = self.mutant[i] self.object_function_values[i] = tmp def print_best(self): m = min(self.object_function_values) i = self.object_function_values.index(m) print("轮数:" + str(self.cur_round)) print("最佳个体:" + str(self.individuality[i])) print("目标函数值:" + str(m)) def evolution(self): while self.cur_round < self.rounds: self.mutate() self.crossover_and_select() self.print_best() self.cur_round = self.cur_round + 1 #测试部分 if __name__ == "__main__": def f(v): return -(v[1]+47)*np.sin(np.sqrt(np.abs(v[1]+(v[0]/2)+47))) - v[0] * np.sin(np.sqrt(np.abs(v[0]-v[1]-47))) p = Population(min_range=-513, max_range=513, dim=2, factor=0.8, rounds=100, size=100, object_func=f) p.evolution()
结果如下图:
。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。