对角解码协方差矩阵适应 (DDCMA)

class pypop7.optimizers.es.ddcma.DDCMA(problem, options)[source]

对角解码协方差矩阵适应(DDCMA)。

注意

DDCMA 是设计精良的 CMA-ES 算法的最新改进版本,它通过自适应对角解码,同时享受了 SEP-CMA-ES(在几乎可分离问题上更快的适应)和 CMA-ES(在病态非可分离问题上更稳健的适应)两个世界的优势。由于它具有二次时间复杂度(相对于每次采样),因此强烈建议首先尝试其他ES变体(例如,LMCMALMMAES)进行大规模黑盒优化。

Parameters:
  • 问题 (字典) –

    问题参数包含以下常见设置 ():
    • ’fitness_function’ - 需要最小化的目标函数 (函数),

    • ’ndim_problem’ - 维度数 (整数),

    • ’upper_boundary’ - 搜索范围的上限 (类数组),

    • ’lower_boundary’ - 搜索范围的下限 (类数组).

  • options (dict) –

    具有以下常见设置的优化器选项 (keys):
    • ’max_function_evaluations’ - 函数评估的最大次数 (int, 默认: np.inf),

    • ’max_runtime’ - 允许的最大运行时间 (float, 默认: np.inf),

    • ’seed_rng’ - 随机数生成的种子,需要显式设置 (int);

    以及以下特定设置 (keys):
    • ’sigma’ - 初始全局步长,也称为变异强度 (float),

    • ’mean’ - 初始(起始)点,也称为高斯搜索分布的均值 (array_like),

      • 如果未给出,将从均匀分布中随机抽取一个样本,其搜索范围由 problem[‘lower_boundary’]problem[‘upper_boundary’] 限定。

    • ’n_individuals’ - 后代数量,也称为后代种群大小 (int, 默认: 4 + int(3*np.log(problem[‘ndim_problem’]))).

示例

使用黑盒优化器 DDCMA 来最小化著名的测试函数 Rosenbrock:

 1>>> import numpy  # engine for numerical computing
 2>>> from pypop7.benchmarks.base_functions import rosenbrock  # function to be minimized
 3>>> from pypop7.optimizers.es.ddcma import DDCMA
 4>>> problem = {'fitness_function': rosenbrock,  # to define problem arguments
 5...            'ndim_problem': 2,
 6...            'lower_boundary': -5.0*numpy.ones((2,)),
 7...            'upper_boundary': 5.0*numpy.ones((2,))}
 8>>> options = {'max_function_evaluations': 5000,  # to set optimizer options
 9...            'seed_rng': 2022,
10...            'mean': 3.0*numpy.ones((2,)),
11...            'sigma': 3.0}  # global step-size may need to be fine-tuned for better performance
12>>> ddcma = DDCMA(problem, options)  # to initialize the optimizer class
13>>> results = ddcma.optimize()  # to run the optimization/evolution process
14>>> print(f"DDCMA: {results['n_function_evaluations']}, {results['best_so_far_y']}")
15DDCMA: 5000, 3.0714e-19

关于Python代码的正确性检查,请参考此基于代码的可重复性报告以获取所有详细信息。对于基于pytest的自动化测试,请参见test_ddcma.py

mean

初始(起始)点,也称为高斯搜索分布的平均值。

Type:

array_like

n_individuals

后代数量,也称为后代种群大小。

Type:

int

sigma

最终的全局步长,也称为变异强度。

Type:

float

参考文献

Akimoto, Y. 和 Hansen, N., 2020. Diagonal acceleration for covariance matrix adaptation evolution strategies. 进化计算, 28(3), 第405-435页.

请参考Akimoto教授的官方 Python实现: https://gist.github.com/youheiakimoto/1180b67b5a0b1265c204cba991fa8518