Pygad优化：当适应度饱和时动态重置种群以避免局部最优

本文详细介绍了在pygad中实现一种动态种群重置策略的方法。通过在遗传算法运行时监测适应度函数的饱和情况，当算法陷入局部最优时，自动重新初始化种群，从而帮助算法跳出局部极值，提高找到全局最优解的可能性。

在遗传算法（Genetic Algorithm, GA）的优化过程中，一个常见挑战是算法过早收敛到局部最优解，导致无法发现全局最优解。当算法的适应度在连续多代中不再显著提升时，这通常意味着种群多样性下降，算法可能已陷入局部极值。为了克服这一问题，我们可以引入一种动态策略：当检测到适应度饱和时，重新初始化种群，从而为算法注入新的探索能力。

Pygad中的回调机制与动态种群重置

Pygad是一个功能强大的Python遗传算法库，它提供了灵活的回调函数机制，允许开发者在遗传算法的不同阶段执行自定义逻辑。on_generation 回调函数是实现动态种群重置的关键，它会在每一代完成计算后被调用。

通过在 on_generation 回调中实现以下逻辑，我们可以有效地进行种群重置：

1. 监测代数：确保在算法运行足够代数后才开始检查饱和，避免过早重置。
2. 判断适应度饱和：比较当前代与前若干代的最佳适应度值。如果这些值在指定代数内保持不变，则认为适应度已饱和。
3. 重新初始化种群：调用 ga_instance.initialize_population() 方法，使用当前的基因空间参数生成一个新的随机种群。这将替换掉当前陷入局部最优的种群，为算法带来新的搜索方向。

示例代码：实现适应度饱和时的种群重置

以下是一个具体的Pygad示例，展示了如何在适应度饱和10代后，动态重新初始化种群。

燕雀Logo

为用户提供LOGO免费设计在线生成服务

101

查看详情

import pygad

# 定义适应度函数
# 在实际应用中，此函数应根据具体优化问题进行定义
def fitness_func(ga_instance, solution, solution_idx):
    """
    示例适应度函数，返回固定值5。
    在实际场景中，此函数会根据解决方案的质量返回一个评估值。
    """
    return 5

# 定义on_generation回调函数
def on_generation(ga_i):
    """
    每一代完成后执行的回调函数。
    用于检查适应度饱和并重新初始化种群。
    """
    # 确保算法至少运行了10代才开始检查饱和
    if ga_i.generations_completed > 10:
        # 检查过去10代的最佳适应度是否相同，即适应度是否饱和
        # best_solutions_fitness[-1] 是当前代的最佳适应度
        # best_solutions_fitness[-10] 是10代前的最佳适应度
        if ga_i.best_solutions_fitness[-1] == ga_i.best_solutions_fitness[-10]:
            print(f"适应度在第 {ga_i.generations_completed} 代饱和，重新初始化种群...")
            # 重新初始化种群
            # 使用当前GA实例的基因范围和类型参数来生成新种群
            ga_i.initialize_population(low=ga_i.init_range_low,
                                       high=ga_i.init_range_high,
                                       allow_duplicate_genes=ga_i.allow_duplicate_genes,
                                       mutation_by_replacement=True, # 示例中保持与原设定一致
                                       gene_type=ga_i.gene_type)
            # 此时，新的种群已创建并赋值给ga_i.population属性，算法将基于新种群继续迭代。

# 初始化Pygad GA实例
ga_instance = pygad.GA(num_generations=50, # 增加总代数以观察效果
                       sol_per_pop=10,    # 每代解决方案数量
                       num_genes=2,       # 每个解决方案的基因数量
                       num_parents_mating=4, # 参与交配的父代数量
                       fitness_func=fitness_func, # 适应度函数
                       on_generation=on_generation, # 注册回调函数
                       init_range_low=0,  # 基因初始化范围下限
                       init_range_high=10, # 基因初始化范围上限
                       gene_type=float,   # 基因类型
                       mutation_percent_genes=10 # 突变百分比
                      )

# 运行遗传算法
ga_instance.run()

# 打印最终结果（可选）
print("优化完成。")
print(f"最终最佳解决方案: {ga_instance.best_solution()[0]}")
print(f"最终最佳适应度: {ga_instance.best_solution()[1]}")

代码解析:

• fitness_func: 在这个示例中，它简单地返回一个常数 5，以便更容易观察到适应度饱和的情况。在实际应用中，您需要根据具体的优化问题来设计一个能够评估解决方案质量的适应度函数。
• on_generation(ga_i): 这是核心逻辑所在。
  • if ga_i.generations_completed > 10:：设置了一个阈值，表示在算法运行至少10代之后才开始检查适应度饱和。这避免了在算法初期（种群尚未充分演化）就进行不必要的重置。
  • if ga_i.best_solutions_fitness[-1] == ga_i.best_solutions_fitness[-10]:：这一行是判断适应度饱和的关键。它比较了当前代（[-1]）和10代前（[-10]）的最佳适应度值。如果它们相等，则认为适应度在过去10代中没有提升，即已饱和。
  • ga_i.initialize_population(...): 当检测到饱和时，调用此方法重新生成一个全新的种群。参数如 low, high, allow_duplicate_genes, gene_type 等应与GA实例初始化时的设置保持一致，以确保新种群的基因符合问题约束。这个方法会直接更新 ga_i.population 属性，使得下一代算法从新的种群开始演化。

注意事项与最佳实践

1. 饱和判断的阈值: 示例中使用了10代作为判断适应度饱和的阈值。这个值并非固定，应根据具体问题和算法特性进行调整。较小的阈值可能导致频繁重置，增加计算开销；较大的阈值可能使算法在局部最优中停留过久。
2. 重置参数的选择: initialize_population 方法的参数应谨慎选择。通常，我们会沿用GA实例初始化时的基因范围（init_range_low, init_range_high）和基因类型（gene_type），以确保新种群的合法性。
3. 对收敛速度的影响: 动态重置种群可能会打断算法的正常收敛过程，尤其是在算法接近全局最优时。因此，需要权衡其带来的探索能力提升与潜在的收敛效率损失。
4. 与其他策略结合: 动态种群重置可以与自适应变异率、精英保留等其他遗传算法策略结合使用，以进一步提高算法的性能和鲁棒性。
5. 记录与监控: 在 on_generation 回调中添加日志输出（如示例中的 print 语句）可以帮助您监控算法的行为，了解何时发生了种群重置，并据此调整参数。

总结

通过在Pygad中使用 on_generation 回调函数，并结合适应度饱和检测逻辑，我们可以实现一种有效的动态种群重置策略。这种策略能够显著增强遗传算法跳出局部最优的能力，从而提高在复杂优化问题中找到全局最优解的机会。在实际应用中，关键在于根据具体问题调整饱和判断的阈值和重置策略，以达到最佳的优化效果。

以上就是Pygad优化：当适应度饱和时动态重置种群以避免局部最优的详细内容，更多请关注php中文网其它相关文章！