遗传算法简介及其实例应用（基于Python）

已于 2024-05-05 22:33:25 修改
阅读量2k 收藏 41
点赞数 52
分类专栏： 算法学习随笔 文章标签： python 算法
于 2024-05-04 10:52:19 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_60699468/article/details/138410349
版权

目录

遗传算法简介

优化问题

编程建模 

初始化种群

适应度计算函数

边界约束破坏处理

选择高适应度的个体

杂交运算

突变运算

主循环

结果分析

约束条件为不等式

检验不等式结果是否正确

约束条件为等式

检验等式结果是否正确

遗传算法简介

遗传算法是一种启发式算法，通过模拟自然界中的生物进化过程来解决优化问题。遗传算法的基本思想源于达尔文的进化论，通过模拟自然选择、交叉、变异等操作，逐步优化解空间中的个体，寻找问题的最优解或者较好的解。具体来说，遗传算法包括以下几个步骤：

  1. 初始化种群：随机生成一定数量的个体作为初始种群。
  2. 选择：根据适应度函数选出优秀个体，通常选择适应度较高的个体用于后续繁殖。
  3. 交叉：通过某种交叉方式，将选出的优秀个体进行交叉繁殖，产生下一代个体。
  4. 变异：对新生代的个体进行一定概率的变异操作，引入新的基因信息，增加种群的多样性。
  5. 竞争：将经过交叉和变异产生的个体与父代进行竞争，从而保留适应度较高的个体作为下一代。

遗传算法在解决复杂优化问题、搜索空间较大或连续的问题时具有一定优势，并在实际应用中得到了广泛的应用

优化问题

数学优化问题如下，这是一个二维非线性优化问题，目标函数和约束条件中都带有非线性约束。

max (x_1-3)^{2}+(x_2-2)^{2}

x_{1}^{2}+x_{2}^{2}- 5\leq 0

x_1\geq 0,x_2\geq 0

x_1+2x_2-4\leq0

编程建模 

初始化种群

常用的编码方式有二进制编码和实数编码，本文采用实数编码方式进行编码。由于这个问题是一个二维非线性优化问题，因此每个个体包含两个基因片段对应x_1和x_2。初始化种群的主要难点在于如何考虑优化问题的约束条件，使得生成的初始种群尽可能可行，减小可行解搜索的空间。

我们观察到约束1x_{1}^{2}+x_{2}^{2}- 5\leq 0，它的可行域区间实际上是以(0, 0)为圆心，以\sqrt{5}为半径的圆，而约束2又要求x_1\geq 0,x_2\geq 0，因此它的可行域区间实际上是第一象限的四分之一圆，并且x_1和x_2还需要满足不等式约束x_1+2x_2-4\leq 0。故可行域为四分之一圆和不等式约束划分的区域内。故生成种群时可以让x_1=random(0, \sqrt{5}),x_2=random(0, 2)，这样生成的初始种群虽然不是都可行，但是缩小了可行解空间。代码如下

def initializing_populations(pop_number):
    """
    生成初始化种群
    :param pop_number: 种群数
    """
    temp_pop_array1 = np.random.uniform(0, math.sqrt(5)+1e-7, pop_number).reshape(-1, 1)
    temp_pop_array2 = np.random.uniform(0, 2, pop_number).reshape(-1, 1)
    pop_array = np.hstack((temp_pop_array1, temp_pop_array2))
    return pop_array

适应度计算函数

这个问题是个最大化问题，我们选用各个体目标函数的计算值作为个体的适应度。

def fitness_calculation_function(pop_array):
    """
    适应度计算函数，即目标函数，求当前约束下目标函数的最大值
    :param pop_array:
    :return: 各个体的目标函数值，每一行表示一个个体目标函数值
    """
    fit_array = (pop_array[:, 0]-3)**2+(pop_array[:, 1]-2)**2
    return fit_array.reshape(-1, 1)

边界约束破坏处理

该优化问题的边界约束都是不等式约束，因此可以采用惩罚函数的方法，对不满足约束条件的个体的适应度进行惩罚，减小其产生下一代的概率。

def penalty_function(pop_array, fit_array):
    """
    对于不满足不等式约束的个体的进行惩罚
    :param fit_array: 适应度函数
    :param pop_array: 种群
    :return: 加入惩罚后的适应度
    """
    penalty_array = 4-(pop_array[:, 0]+2*pop_array[:, 1])
    statis_index = np.where(penalty_array >= 0)  # 满足不等式约束的索引
    penalty_array[statis_index] = 0  # 将满足约束的惩罚值更改为0
    new_fit_array = fit_array+penalty_array.reshape(-1, 1)  # 将惩罚值加入适应度数组，得到新的适应度数组
    zero_index = np.where(new_fit_array < 0)  # 惩罚后的适应度如果小于0，将其更改为0
    new_fit_array[zero_index] = 0
    return new_fit_array

选择高适应度的个体

遗传算法中常用的选择方法有：

  1. 轮盘赌选择：根据每个个体的适应度值比例来选择个体，适应度值高的个体有更大的概率被选中。

  2. 锦标赛选择：随机选取一定数量的个体，根据它们的适应度值选择较优秀的个体作为父代。

  3. 最优选择：直接选择适应度值最高的个体作为父代。

  4. 排序选择：将个体按适应度值从高到低进行排序，根据排名选取个体作为父代。

  5. 随机选择：随机选择个体作为父代，不考虑适应度值

本文选择轮盘赌法，轮盘赌适用于处理适应度为非负数的情况 。根据适应度为每个个体附加相应的杂交概率，适应度越大的个体杂交产生下一代的可能性越大。

def roulette_selection(fit_array):
    """
    轮盘赌选择法，根据适应度的大小，对每个个体附加交叉的概率，适应度越大交叉的可能性越大
    :param fit_array: 适应度数组
    :return: 返回种群内各个个体的交叉概率
    """
    crossover_p_array = np.full((fit_array.size, 1), np.nan, dtype=float)
    fit_sum = np.sum(fit_array)
    if abs(fit_sum) <= 1e-4:
        sys.exit("适应度之和不能为0，请更换选择方法")
    accumulate_p = 0  # 累积概率
    for i in range(fit_array.size):
        crossover_p_array[i, 0] = fit_array[i, 0]/fit_sum+accumulate_p
        accumulate_p = crossover_p_array[i, 0]
    return crossover_p_array

杂交运算

根据输入的杂交概率和种群数量，随机抽取个体进行杂交。

def crossover_operation(pop_number, crossover_p, parent_pop, crossover_p_array):
    """
    交叉运算，高适应度的个体交叉的概率大，低适应度的个体交叉的概率小
    :param crossover_p_array: 种群的交叉概率数组
    :param pop_number: 种群数量
    :param parent_pop: 父代
    :param crossover_p: 交叉概率
    :return: 返回交叉后的种群
    """
    new_parent_pop = copy(parent_pop)
    for i in range(int(crossover_p*pop_number)):  # 总交叉次数
        temp_crossover_p1 = random.uniform(0, 1)  # 随机选择父代概率1
        temp_crossover_p2 = random.uniform(0, 1)  # 随机选择父代概率2
        parent_index1 = np.where(temp_crossover_p1 <= crossover_p_array)[0][0]  # 父代1行索引
        parent_index2 = np.where(temp_crossover_p2 <= crossover_p_array)[0][0]  # 父代2行索引
        random_int = np.random.randint(2)  # 生成0-1之间的随机整数确定交叉的位置
        # 根据父代进行交叉得到子代
        new_parent_pop[parent_index1, random_int] = parent_pop[parent_index2, random_int]
        new_parent_pop[parent_index2, random_int] = parent_pop[parent_index2, random_int]
    return new_parent_pop

突变运算

种群中所有的个体突变的概率都相同。

def mutation_operation(pop_number, mutation_p, parent_pop):
    """
    突变函数，对种群进行突变运算
    :param pop_number:
    :param mutation_p: 突变概率
    :param parent_pop:
    :return:
    """
    new_parent_pop = parent_pop
    for i in range(int(pop_number*mutation_p)):
        random_int1 = np.random.randint(pop_number)  # 突变个体
        random_int2 = np.random.randint(2)  # 突变位置
        if random_int2 == 0:
            new_parent_pop[random_int1, random_int2] = random.uniform(0, math.sqrt(5) + 1e-7)
        else:
            new_parent_pop[random_int1, random_int2] = random.uniform(0, 2)
    return new_parent_pop

主循环

根据输入的最大繁殖代数进行循环迭代。

def main_loop(pop_number, reproduction_number, crossover_p, mutation_p):
    """
    遗传算法实现的主循环
    :param pop_number: 种群数量
    :param reproduction_number: 最大的繁殖代数
    :param crossover_p: 杂交概率
    :param mutation_p: 突变概率
    """
    init_pop_array = initializing_populations(pop_number)  # 初始化种群
    for i in range(reproduction_number):
        if i == 0:
            parent_pop = init_pop_array  # 父代种群
        else:
            parent_pop = child_pop_array
        fit_array = penalty_function(parent_pop, fitness_calculation_function(parent_pop))  # 获取父本中个体的适应度
        crossover_p_array = roulette_selection(fit_array)  # 获取各个体的杂交概率
        temp_pop_array = crossover_operation(pop_number, crossover_p, parent_pop, crossover_p_array)  # 进行交叉运算，生成中间结果的种群
        temp_pop_array = mutation_operation(pop_number, mutation_p, temp_pop_array)  # 进行突变运算
        temp_fit_array = penalty_function(temp_pop_array, fitness_calculation_function(temp_pop_array))  # 临时种群适应度
        child_pop_array = competitiveness_operation(pop_number, temp_pop_array,temp_fit_array, parent_pop, fit_array)

结果分析

约束条件为不等式x_1+2x_2-4\leq 0

输入参数：种群数20，最大繁殖代数1000，杂交概率0.9，突变概率0.3

if __name__ == "__main__":
    main_loop(20, 1000, 0.9, 0.3)

 结果表示适应度值,x_1,x_2的值。

12.990864951472929 [0.00074035 0.00117371]  # 第一次运行结果
12.983539152035894 [0.00057131 0.00326098]  #  第二次运行结果
12.986867969022803 [0.00157115 0.00092712]  # 第三次运行结果

检验不等式结果是否正确

 首先上述结果x_1,x_2都满足约束条件，我们采用开源非线性求解器ipopt来检验结果是否最优。Ipopt（Interior Point OPTimizer）是一种开源软件包，用于求解非线性连续优化问题。它采用内点方法来解决大规模非线性优化问题，可以处理包括等式约束、不等式约束和边界约束等多种类型的约束条件。Ipopt能够高效地求解复杂的非线性优化问题，通常用于在数学建模、工程优化、经济学等领域中的应用。Ipopt在连续优化领域有着广泛的应用，且性能表现良好，是许多优化问题的首选求解器之一。（本文重点不在此，用法不做介绍）

    from pyomo.environ import *

    model = ConcreteModel()

    model.x1 = Var(within=NonNegativeReals)
    model.x2 = Var(within=NonNegativeReals)
    model.obj_value = Var()

    def rule1(m):
        return m.x1**2+m.x2**2-5 <= 0
    model.constraint1 = Constraint(rule=rule1)

    def rule2(m):
        return m.x1+2*m.x2-4 <= 0
    model.constraint2 = Constraint(rule=rule2)

    def get_obj_value(m):
        return model.obj_value == (m.x1 - 3) ** 2 + (m.x2 - 2) ** 2
    model.get_obj_value_con = Constraint(rule=get_obj_value)

    def obj_rule(m):
        return model.obj_value
    model.objective_constraint = Objective(rule=obj_rule, sense=maximize)

    opt = SolverFactory('ipopt')
    opt.solve(model, tee=True)
    print(value(model.obj_value), value(model.x1), value(model.x2))

所得结果如下：

13.0 0.0 0.0

和遗传算法求解的一致，说明执行的算法是有效的。

约束条件为等式x_1+2x_2-4= 0

当我们将约束条件修改为等式之后，采用惩罚方法还能得到稳定的结果吗？

def penalty_function2(pop_array, fit_array):
    """
    将约束更改为等式约束
    :param fit_array: 适应度函数
    :param pop_array: 种群
    :return: 加入惩罚后的适应度
    """
    penalty_array = 4-(pop_array[:, 0]+2*pop_array[:, 1])
    index_lower = np.where(penalty_array >= 0)
    index_upper = np.where(penalty_array <= 1e-6)

    # 取交集，即既满足大于等于0，又满足小于等于1e-6的索引
    statis_index = np.intersect1d(index_lower, index_upper)
    penalty_array[statis_index] = 0  # 将满足约束的惩罚值更改为0
    new_fit_array = fit_array-np.abs(penalty_array.reshape(-1, 1))  # 将惩罚值加入适应度数组，得到新的适应度数组
    zero_index = np.where(new_fit_array < 0)
    new_fit_array[zero_index] = 0
    return new_fit_array

其余参数和条件不变，运行几次输出结果

8.987463656361355 [1.87469013e-03 1.99841554e+00]  # 第一次运行结果
8.998749946120093 [2.44311335e-04 1.42785488e-05]  # 第二次运行结果
8.993601963219662 [0.00088544 0.00098629]  # 第三次运行结果

从上面可以看出输入条件相同，但是所得结果却相差较大，并且不满足等式约束，无法像不等式约束一样获得较稳定的最优解，那么如何解决这个问题呢？利用等式约束增加一个修正，修正个体的基因编码。

def genetic_modification(pop_array):
    """
    对个体的基因型进行修正，根据等式约束
    :param pop_array:
    :return:
    """
    new_pop_array = copy(pop_array)
    new_pop_array[:, 1] = 2-0.5*new_pop_array[:, 0]
    return new_pop_array

运行三次结果如下：

8.990181777990122 [1.63692857e-03 1.99918154e+00]  # 第一次运行结果
8.991384439727147 [1.43635653e-03 1.99928182e+00]  # 第二次运行结果
8.998369968364907 [2.71687317e-04 1.99986416e+00]  # 第三次运行结果

检验等式结果是否正确

model = ConcreteModel()

model.x1 = Var(within=NonNegativeReals)
model.x2 = Var(within=NonNegativeReals)

model.obj_value = Var()

def rule1(m):
    return m.x1 ** 2 + m.x2 ** 2 - 5 <= 0

model.constraint1 = Constraint(rule=rule1)

def rule2(m):
    return m.x1 + 2 * m.x2 - 4 == 0

model.constraint2 = Constraint(rule=rule2)

def get_obj_value(m):
    return model.obj_value == (m.x1 - 3) ** 2 + (m.x2 - 2) ** 2

model.get_obj_value_con = Constraint(rule=get_obj_value)

def obj_rule(m):
    return model.obj_value

model.objective_constraint = Objective(rule=obj_rule, sense=maximize)

opt = SolverFactory('ipopt')
opt.solve(model, tee=True)
print(value(model.obj_value), value(model.x1), value(model.x2))

结果如下：

9.0 0.0 2.0

完整代码：

import sys
import random
from copy import copy

import numpy as np
import math


# 初始化种群
def initializing_populations(pop_number):
    """
    生成初始化种群
    :param gen_number: 每个个体的基因数量
    :param pop_number: 种群数
    """
    temp_pop_array1 = np.random.uniform(0, math.sqrt(5)+1e-7, pop_number).reshape(-1, 1)
    temp_pop_array2 = np.random.uniform(0, 2, pop_number).reshape(-1, 1)
    pop_array = np.hstack((temp_pop_array1, temp_pop_array2))
    return pop_array

def fitness_calculation_function(pop_array):
    """
    适应度计算函数，即目标函数，求当前约束下目标函数的最大值
    :param pop_array:
    :return: 各个体的目标函数值，每一行表示一个个体目标函数值
    """
    fit_array = (pop_array[:, 0]-3)**2+(pop_array[:, 1]-2)**2
    return fit_array.reshape(-1, 1)

def roulette_selection(fit_array):
    """
    轮盘赌选择法，根据适应度的大小，对每个个体附加交叉的概率，适应度越大交叉的可能性越大
    :param fit_array: 适应度数组
    :return: 返回种群内各个个体的交叉概率
    """
    crossover_p_array = np.full((fit_array.size, 1), np.nan, dtype=float)
    fit_sum = np.sum(fit_array)
    if abs(fit_sum) <= 1e-4:
        sys.exit("适应度之和不能为0，请更换选择方法")
    accumulate_p = 0  # 累积概率
    for i in range(fit_array.size):
        crossover_p_array[i, 0] = fit_array[i, 0]/fit_sum+accumulate_p
        accumulate_p = crossover_p_array[i, 0]
    return crossover_p_array

def penalty_function(pop_array, fit_array):
    """
    对于不满足等式约束的个体的进行惩罚
    :param fit_array: 适应度函数
    :param pop_array: 种群
    :return: 加入惩罚后的适应度
    """
    penalty_array = 4-(pop_array[:, 0]+2*pop_array[:, 1])
    statis_index = np.where(penalty_array >= 0)  # 满足不等式约束的索引
    penalty_array[statis_index] = 0  # 将满足约束的惩罚值更改为0
    new_fit_array = fit_array+penalty_array.reshape(-1, 1)  # 将惩罚值加入适应度数组，得到新的适应度数组
    zero_index = np.where(new_fit_array < 0)
    new_fit_array[zero_index] = 0
    return new_fit_array


def crossover_operation(pop_number, crossover_p, parent_pop, crossover_p_array):
    """
    交叉运算，高适应度的个体交叉的概率大，低适应度的个体交叉的概率小
    :param crossover_p_array: 种群的交叉概率数组
    :param pop_number: 种群数量
    :param parent_pop: 父代
    :param crossover_p: 交叉概率
    :return: 返回交叉后的种群
    """
    new_parent_pop = copy(parent_pop)
    for i in range(int(crossover_p*pop_number)):  # 总交叉次数
        temp_crossover_p1 = random.uniform(0, 1)  # 随机选择父代概率1
        temp_crossover_p2 = random.uniform(0, 1)  # 随机选择父代概率2
        parent_index1 = np.where(temp_crossover_p1 <= crossover_p_array)[0][0]  # 父代1行索引
        parent_index2 = np.where(temp_crossover_p2 <= crossover_p_array)[0][0]  # 父代2行索引
        random_int = np.random.randint(2)  # 生成0-1之间的随机整数确定交叉的位置
        # 根据父代进行交叉得到子代
        new_parent_pop[parent_index1, random_int] = parent_pop[parent_index2, random_int]
        new_parent_pop[parent_index2, random_int] = parent_pop[parent_index2, random_int]
    return new_parent_pop


def mutation_operation(pop_number, mutation_p, parent_pop):
    """
    突变函数，对种群进行突变运算
    :param pop_number:
    :param mutation_p: 突变概率
    :param parent_pop:
    :return:
    """
    new_parent_pop = parent_pop
    for i in range(int(pop_number*mutation_p)):
        random_int1 = np.random.randint(pop_number)  # 突变个体
        random_int2 = np.random.randint(2)  # 突变位置
        if random_int2 == 0:
            new_parent_pop[random_int1, random_int2] = random.uniform(0, math.sqrt(5) + 1e-7)
        else:
            new_parent_pop[random_int1, random_int2] = random.uniform(0, 2)
    return new_parent_pop


def competitiveness_operation(pop_number, temp_pop_array, temp_fit_array,  parent_pop, parent_fit_array):
    """
    令生成的temp_pop_array与父代进行竞争， 降序排列，选取前pop_number作为子代
    :param parent_fit_array: 父代适应度
    :param pop_number: 种群数
    :param temp_fit_array: 临时种群适应度
    :param temp_pop_array: 临时种群
    :param parent_pop: 父代
    :return: 子代
    """
    fit_array = np.vstack((parent_fit_array, temp_fit_array))
    sort_fit_array = sorted(list(fit_array[:, 0]), reverse=True)
    child_generation = copy(parent_pop)
    for i in range(pop_number):
        child_index = np.where(sort_fit_array[i] == fit_array)[0][0]
        if child_index > pop_number-1:
            child_generation[i, :] = temp_pop_array[child_index-pop_number, :]
        else:
            child_generation[i, :] = parent_pop[child_index, :]
    print(sort_fit_array[0], child_generation[0, :])
    return child_generation

def main_loop(pop_number, reproduction_number, crossover_p, mutation_p):
    """
    遗传算法实现的主循环，当约束为不等式
    :param pop_number: 种群数量
    :param reproduction_number: 最大的繁殖代数
    :param crossover_p: 杂交概率
    :param mutation_p: 突变概率
    """
    init_pop_array = initializing_populations(pop_number)  # 初始化种群
    for i in range(reproduction_number):
        if i == 0:
            parent_pop = init_pop_array  # 父代种群
        else:
            parent_pop = child_pop_array
        fit_array = penalty_function(parent_pop, fitness_calculation_function(parent_pop))  # 获取父本中个体的适应度
        crossover_p_array = roulette_selection(fit_array)  # 获取各个体的杂交概率
        temp_pop_array = crossover_operation(pop_number, crossover_p, parent_pop, crossover_p_array)  # 进行交叉运算，生成中间结果的种群
        temp_pop_array = mutation_operation(pop_number, mutation_p, temp_pop_array)  # 进行突变运算
        temp_fit_array = penalty_function(temp_pop_array, fitness_calculation_function(temp_pop_array))  # 临时种群适应度
        child_pop_array = competitiveness_operation(pop_number, temp_pop_array,temp_fit_array, parent_pop, fit_array)

def penalty_function2(pop_array, fit_array):
    """
    将约束更改为等式约束
    :param fit_array: 适应度函数
    :param pop_array: 种群
    :return: 加入惩罚后的适应度
    """
    penalty_array = 4-(pop_array[:, 0]+2*pop_array[:, 1])
    index_lower = np.where(penalty_array >= 0)
    index_upper = np.where(penalty_array <= 1e-6)

    # 取交集，即既满足大于等于0，又满足小于等于1e-6的索引
    statis_index = np.intersect1d(index_lower, index_upper)
    penalty_array[statis_index] = 0  # 将满足约束的惩罚值更改为0
    new_fit_array = fit_array-np.abs(penalty_array.reshape(-1, 1))  # 将惩罚值加入适应度数组，得到新的适应度数组
    zero_index = np.where(new_fit_array < 0)
    new_fit_array[zero_index] = 0
    return new_fit_array

def genetic_modification(pop_array):
    """
    对个体的基因型进行修正，根据等式约束
    :param pop_array:
    :return:
    """
    new_pop_array = copy(pop_array)
    new_pop_array[:, 1] = 2-0.5*new_pop_array[:, 0]
    return new_pop_array


def main_loop2(pop_number, reproduction_number, crossover_p, mutation_p):
    """
    遗传算法实现的主循环，当约束为等式
    :param pop_number: 种群数量
    :param reproduction_number: 最大的繁殖代数
    :param crossover_p: 杂交概率
    :param mutation_p: 突变概率
    """
    init_pop_array = initializing_populations(pop_number)  # 初始化种群
    for i in range(reproduction_number):
        if i == 0:
            parent_pop = init_pop_array  # 父代种群
        else:
            parent_pop = child_pop_array
        parent_pop = genetic_modification(parent_pop)
        fit_array = fitness_calculation_function(parent_pop) # 获取父本中个体的适应度
        crossover_p_array = roulette_selection(fit_array)  # 获取各个体的杂交概率
        temp_pop_array = crossover_operation(pop_number, crossover_p, parent_pop, crossover_p_array)  # 进行交叉运算，生成中间结果的种群
        temp_pop_array = mutation_operation(pop_number, mutation_p, temp_pop_array)  # 进行突变运算
        temp_pop_array = genetic_modification(temp_pop_array)
        temp_fit_array = fitness_calculation_function(temp_pop_array)  # 临时种群适应度
        child_pop_array = competitiveness_operation(pop_number, temp_pop_array, temp_fit_array, parent_pop, fit_array)

from pyomo.environ import *
# 检验不等式结果是否成立
def check_result1():

    model = ConcreteModel()

    model.x1 = Var(within=NonNegativeReals)
    model.x2 = Var(within=NonNegativeReals)

    model.obj_value = Var()

    def rule1(m):
        return m.x1**2+m.x2**2-5 <= 0
    model.constraint1 = Constraint(rule=rule1)

    def rule2(m):
        return m.x1+2*m.x2-4 <= 0
    model.constraint2 = Constraint(rule=rule2)

    def get_obj_value(m):
        return model.obj_value == (m.x1 - 3) ** 2 + (m.x2 - 2) ** 2
    model.get_obj_value_con = Constraint(rule=get_obj_value)

    def obj_rule(m):
        return model.obj_value
    model.objective_constraint = Objective(rule=obj_rule, sense=maximize)

    opt = SolverFactory('ipopt')
    opt.solve(model, tee=True)
    print(value(model.obj_value), value(model.x1), value(model.x2))


# 检验等式结果是否成立
def check_result2():
    model = ConcreteModel()

    model.x1 = Var(within=NonNegativeReals)
    model.x2 = Var(within=NonNegativeReals)

    model.obj_value = Var()

    def rule1(m):
        return m.x1 ** 2 + m.x2 ** 2 - 5 <= 0

    model.constraint1 = Constraint(rule=rule1)

    def rule2(m):
        return m.x1 + 2 * m.x2 - 4 == 0

    model.constraint2 = Constraint(rule=rule2)

    def get_obj_value(m):
        return model.obj_value == (m.x1 - 3) ** 2 + (m.x2 - 2) ** 2

    model.get_obj_value_con = Constraint(rule=get_obj_value)

    def obj_rule(m):
        return model.obj_value

    model.objective_constraint = Objective(rule=obj_rule, sense=maximize)

    opt = SolverFactory('ipopt')
    opt.solve(model, tee=True)
    print(value(model.obj_value), value(model.x1), value(model.x2))

if __name__ == "__main__":
    main_loop2(20, 1000, 0.9, 0.3)
    check_result2()

知情人士黄某
关注
  • 52
    点赞
  • 踩
  • 41
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 2
    评论
专栏目录
遗传算法应用实例
qq_43390943的博客
03-30 2万+
h
遗传算法及python实现
qq_54543084的博客
12-18 6982
遗传算法及python实现
遗传算法的各种应用示例（含源码）
12-27
遗传算法的各种应用示例（含源码） 问题1 求f(x)=x 10*sin(5x) 7*cos(4x)的最大值，其中0<=x<=9 问题2 求下列函数的最大值 % f(x)=10*sin(5x)+7*cos(4x) x∈[0,10] % 问题3 已知n个城市之间的相互距离，现有一个推销员必须遍访这n个城市，并且每个城市 只能访问一次，最后又必须返回出发城市。如何安排他对这些城市的访问次序，可使其 旅行路线的总长度最短？ 用图论的术语来说，假设有一个图g=(v,e)，其中v是顶点集，e是边集，设d=(dij) 是由顶点i和顶点j之间的距离所组成的距离矩阵，旅行商问题就是求出一条通过所有顶 点且每个顶点只通过一次的具有最短距离的回路。
现代优化算法 （二）： 遗传算法 及应用举例
冷月无声的博客
04-29 2万+
组合优化算法系列： 现代优化算法 （一）：模拟退火算法 及应用举例 现代优化算法 （二）： 遗传算法 及应用举例 现代优化算法(三）：禁忌搜索算法 现代优化算法（四）：改进的遗传算法 现代优化算法（五）： 蚁群算法 目录 遗传算法简介 2 模型及算法 遗传算法简介 遗传算法（Genetic Algorithms，简称 GA）是一种基于自然选择原理和自然遗传机制的搜索（寻优...
遗传算法小结及算法实例（附Matlab代码）
热门推荐
LXzerorb的博客
04-21 5万+
目录 1、遗传算法流程 2、关键参数说明 （1）群体规模 \(NP\) （2）交叉概率 \(P_c\) （3）变异概率 \(P_m\) （4）进化代数 \(G\) 3、MATLAB仿真实例 3.1 遗传算法求解一元函数的极值 3.2 遗传算法求解旅行商问题（TSP） 4、遗传算法的特点 1、遗传算法流程 遗传算法的运算流程如下图所示： 具体步骤如下： （1）初始化。设置进化代数计数器 \(g=0\)，设置最大进化代数 \(G\)，随机生成 \(NP\) 个个体作为初始群体..
一个简单的遗传算法的应用实例（java代码）
但行好事，莫问前程
09-05 7799
求下述二元函数的最大值： 下参考Mitchell《机器学习》中的遗传算法解决上面问题的java代码如下： import java.lang.Math; import java.util.ArrayList; import java.util.List; /** * 采用遗传算法求下列二元函数的最大值 * max f(x1,x2)=x1^2+x2^2; * s.t.x1∈{1,2,
遗传算法概念、步骤、应用解析(案例直白--黄含驰）
bob595078694的博客
01-16 7785
遗传算法将“优胜劣汰，适者生存”的生物进化原理引入优化参数形成的编码串联群体中
基于python的遗传算法设计与实现
04-29
《基于Python的遗传算法设计与实现》 遗传算法（Genetic Algorithm, GA）是一种借鉴生物进化理论中的自然选择、遗传和突变等机制的优化算法，广泛应用于解决复杂问题的求解。Python作为一门易读性强、库丰富的编程...
基于遗传算法的图像分割python实现.rar
04-02
《基于遗传算法的图像分割Python实现》 在计算机视觉领域，图像分割是一项核心任务，它旨在将图像划分为多个有意义的区域，以便更好地理解和分析图像内容。本项目着重探讨了一种利用遗传算法进行图像分割的方法，并...
基于python实现的遗传算法.rar
04-24
本项目是基于Python编程语言实现的遗传算法，旨在帮助用户理解和应用这种算法到实际问题中。以下是对这个压缩包文件内容的详细解读。 首先，让我们了解遗传算法的基本原理。遗传算法受到生物进化论的启发，包括选择...
Python遗传算法求一元函数最大值
12-21
遗传算法是一种基于生物进化原理的全局优化方法，它模拟了自然选择、遗传和突变等生物进化过程来寻找问题的最优解。在这个案例中，我们将使用Python实现遗传算法来求解一元函数的最大值。目标函数是2*sin(x) + cos(x...
基于遗传算法的CVRP建模求解-Python代码
03-10
**基于遗传算法的CVRP建模求解-Python代码** 遗传算法（Genetic Algorithm, GA）是一种模拟生物进化过程的优化方法，广泛应用于解决复杂问题的全局优化。在这个场景中，我们关注的是车辆路径问题（Vehicle Routing ...
遗传算法及实例
Zhoudaxia2000的博客
03-07 3330
遗传算法是模拟生物在自然环境中的遗传和进化的过程而形成的自适应全局优化搜索算法
遗传算法应用案例
他们都叫我神奇的博客
08-13 8128
&lt;&lt;MATLAB在数学建模中的应用&gt;&gt; 案例一.无约束目标函数最大值遗传算法求解策略 求解问题            %主程序:用遗传算法求解y=200*exp(-0.05*x).*sin(x)在[-2,2]上的最大值 clc; clear all; close all; global BitLength global boundsbegin global boun...
遗传算法应用（实例详细演示最小生成树的prufer编码和Cayley定理）
qq_43391414的博客
04-05 2471
prufer编码：用n-2位自然数唯一的表达出一棵n个节点的生成树。而且两者相互可逆，即给定一颗生成树的连接方式，可以唯一确定这棵树的编码。 Cayley定理：n个顶点的完全图中有nn−2n^{n-2}nn−2棵不同的生成树。 显然这两个描述具有很强的联系，n个顶点编号为1,2,…n。有n-2个位置，随便填入1,2,…n中的一个数，这有nn−2n^{n-2}nn−2种可能。 所以我们要证明Cayley定理，只需要证明prufer编码。 prufer编码 1.设有如下生成树（n=6)，根据prufer编码，我
遗传算法及其应用
我一点也不皮的博客
06-27 1万+
遗传算法定义： 遗传算法（Genetic Algorithm，GA）是计算数学中用于解决最优化的搜索算法，是进化算法的一种。进化算法最初是借鉴了进化生物学中的一些现象而发展起来的，这些现象包括遗传、突变、自然选择以及杂交等。
遗传算法（附简单案例及matlab详细代码）
最新发布
你的问题在于，读书太少而想得太多。
06-04 3462
10分钟彻底学会遗传算法！（附典型案例及matlab详细代码）
遗传算法入门
晨晨的博客
07-22 1070
遗传算法理论 遗传算法的原理就和达尔文进化论一样，只有适应环境的个体才能存活并更有可能产生更多后代。分为以下几个主要步骤： 1、初始化，即设置初始群体大小 2、设定适应度，即按照实际需求设置区分个体好坏的标准 3、选择，即把适应度高的个体筛选出来，并用他们通过交叉产生后代，也能通过基因突变的方法产生新后代 4、替代、即用产生的优良后代替换之前适应度低的亲本 具体流程图所示如下： 背包问题解释遗传...
"基于遗传算法的旅行商问题优化求解（python）开发平台与算法原理详解
遗传算法在群智能优化算法中是应用最广泛、最成功、最具代表性的智能优化方法之一。它基于达尔文的生物进化论和孟德尔的遗传变异理论，模拟生物界的进化过程，从而对种群进行优化搜索。遗传算法的设计和应用，是为了...
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

知情人士黄某

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值