声明:非官方解决方法,不一定正规哈
自定义初始值解决方法
以官方的第一个案例:https://github.com/geatpy-dev/geatpy/tree/master/geatpy/demo/soea_demo/soea_demo1 为例:
- 进入类
MyProblem
里 - 在函数
def __init__(self):
里增加一个第一次运行的判断:self.flag = 0
- 在函数
def aimFunc(self, pop)
里修改,把pop.Phen
改为想要的初始值,比如设置出值全0:if self.flag == 0: pop.Phen = np.ones(shape=(40, 1)) self.flag += 1
修改后的MyProblem如下:
class MyProblem(ea.Problem): # 继承Problem父类
def __init__(self):
name = 'MyProblem' # 初始化name(函数名称,可以随意设置)
M = 1 # 初始化M(目标维数)
maxormins = [-1] # 初始化maxormins(目标最小最大化标记列表,1:最小化该目标;-1:最大化该目标)
Dim = 1 # 初始化Dim(决策变量维数)
varTypes = [0] * Dim # 初始化varTypes(决策变量的类型,元素为0表示对应的变量是连续的;1表示是离散的)
lb = [-100] # 决策变量下界
ub = [2] # 决策变量上界
lbin = [1] * Dim # 决策变量下边界(0表示不包含该变量的下边界,1表示包含)
ubin = [1] * Dim # 决策变量上边界(0表示不包含该变量的上边界,1表示包含)
# 调用父类构造方法完成实例化
self.flag = 0
ea.Problem.__init__(self, name, M, maxormins, Dim, varTypes, lb, ub, lbin, ubin)
def aimFunc(self, pop): # 目标函数
if self.flag == 0: # 只有第一次运行时执行这个:
pop.Phen = np.ones(shape=(40, 1))
self.flag += 1
x = pop.Phen # shape = (40,1)
pop.ObjV = x * 2 + 2.0 # 计算目标函数值,赋值给pop种群对象的ObjV属性
针对多个未知数,需要改为:
pop.Phen = np.array([[-0.00000000, -0.00000000, -4.83461801, -0.00000000, -2.43279756, -1.34510964] for i in range(50)] # 50个种群,6个未知数x
修改前,随机生成的值会从函数值为 -80开始增加:
修改后,直接从4开始: