python使用gurobi用法解析和案例

1. Gurobi Python接口的基本使用

在Python中使用Gurobi进行优化，通常需要按以下步骤操作：

1. 导入Gurobi包 ：
   使用import gurobipy as gp导入Gurobi Python接口包，并通常使用缩写gp。

import gurobipy as gp  # 导入Gurobi的Python接口，并缩写为gp
from gurobipy import GRB  # 导入Gurobi中的常量，通常用于设置变量类型和求解状态

2. 创建模型 ：
   通过gp.Model()创建一个新的优化模型。

model = gp.Model("my_model")  # 创建一个新的Gurobi模型，并命名为"my_model"

3. 添加变量 ： 详情见 pycharm中gurobi使用指南-变量类型及声明
   使用model.addVar()或model.addVars()方法添加决策变量。变量的类型可以是连续型、整数型或二元型。

x = model.addVar(vtype=GRB.CONTINUOUS, name="x")  # 添加一个连续型变量x
y = model.addVar(vtype=GRB.BINARY, name="y")  # 添加一个二元变量y

4. 设置目标函数 ：
   使用model.setObjective()方法设置优化的目标函数，可以是最大化或最小化。

model.setObjective(2 * x + 3 * y, GRB.MAXIMIZE)  # 设置目标函数为2x + 3y，并指定最大化

5. 添加约束 ：详情见gurobi约束条件使用大全(model.addConstr()添加单个约束和model.addConstrs()添加多个约束和强不等式约束）
   通过model.addConstr()或model.addConstrs()方法添加约束条件。可以是等式约束、不等式约束（≥，≤ )，遇到强不等式约束解决方法见 python中使用gurobi遇到强不等式约束(只有大于或者小于而不是大于等于或者小于等于的形式)的解决办法

model.addConstr(x + y <= 10, "c1")  # 添加约束条件x + y <= 10，并命名为"c1"

6. 求解模型 ：
   使用model.optimize()方法求解优化问题。

model.optimize()  # 求解模型

7. 获取结果 ：
   优化完成后，可以通过var.x属性获取变量的最优解，通过model.objVal获取最优目标值。

print(f"Optimal value of x: {x.x}")  # 输出变量x的最优解
print(f"Optimal value of the objective function: {model.objVal}")  # 输出目标函数的最优值

---

2. 变量类型

Gurobi支持多种变量类型，常见的包括：

• GRB.CONTINUOUS : 连续变量，可以取任意实数值。
• GRB.BINARY : 二元变量，仅取值0或1。
• GRB.INTEGER : 整数变量，仅取整数值。
• GRB.SEMICONT : 半连续变量，非零时必须大于等于给定的下界。
• GRB.SEMIINT : 半整数变量，非零时必须大于等于给定的下界并取整数值。

例子：

x = model.addVar(vtype=GRB.CONTINUOUS, lb=0, ub=10, name="x")  # 添加一个取值范围在[0, 10]的连续变量
y = model.addVar(vtype=GRB.INTEGER, name="y")  # 添加一个整数变量
z = model.addVar(vtype=GRB.SEMICONT, lb=1, name="z")  # 添加一个下界为1的半连续变量

• lb和ub参数分别表示变量的下界和上界。
• name参数为变量命名，便于调试和结果分析。
• 总结：主要包括变量类型，变量上下界和变量命名

---

3. 目标函数

Gurobi支持线性和二次目标函数，分别通过GRB.MINIMIZE和GRB.MAXIMIZE指定最小化或最大化目标。

例子：

model.setObjective(x + 2 * y, GRB.MINIMIZE)  # 设置目标函数为x + 2y，并指定最小化
model.setObjective(0.5 * x * x + y * y, GRB.MINIMIZE)  # 设置二次目标函数

---

4. 约束条件

约束可以是不等式或等式约束。Gurobi支持添加单个约束和多重约束。

例子：

model.addConstr(x + y == 5, "eq_constr")  # 添加等式约束x + y = 5
model.addConstr(x - 2 * y >= 3, "ineq_constr")  # 添加不等式约束x - 2y >= 3

5. 模型求解和结果分析

求解模型后，必须检查求解状态以确定是否找到最优解。

例子：

model.optimize()  # 求解模型

if model.status == GRB.OPTIMAL:
    print("Optimal solution found")  # 如果找到最优解，输出结果
else:
    print("No optimal solution found")  # 如果未找到最优解，输出相应信息

---

6. 常见注意事项

1. 模型求解的时间和内存 ：
   对于较大的模型，求解时间和内存消耗可能会非常高。建议通过设置TimeLimit、MIPGap等参数控制求解过程。

2. 变量的松弛和固定 ：
   可以使用model.addVars()和model.addConstrs()结合松弛变量或拉格朗日松弛方法对难解问题进行松弛处理。

3. 模型的更新与写入 ：
   在大规模模型中，频繁的变量添加和约束添加后，需要调用model.update()来刷新模型。可以通过model.write("model.lp")导出模型文件。

---

7. gurobi代码示例

以下是一个综合使用上述方法的完整代码示例，演示了从模型创建到求解的整个过程：

import gurobipy as gp  # 导入Gurobi的Python接口，并缩写为gp
from gurobipy import GRB  # 导入Gurobi中的常量

# 创建模型
model = gp.Model("example_model")

# 添加变量
x = model.addVar(vtype=GRB.CONTINUOUS, lb=0, name="x")  # 添加连续变量x
y = model.addVar(vtype=GRB.INTEGER, name="y")  # 添加整数变量y

# 设置目标函数
model.setObjective(2 * x + 3 * y, GRB.MAXIMIZE)  # 设置目标函数为2x + 3y，最大化

# 添加约束条件
model.addConstr(x + y <= 10, "c1")  # 添加约束x + y <= 10
model.addConstr(x - y >= 3, "c2")  # 添加约束x - y >= 3

# 求解模型
model.optimize()

# 获取并输出结果
if model.status == GRB.OPTIMAL:
    print(f"Optimal value of x: {x.x}")  # 输出x的最优解
    print(f"Optimal value of y: {y.x}")  # 输出y的最优解
    print(f"Optimal value of the objective: {model.objVal}")  # 输出目标函数的最优值
else:
    print("Optimal solution not found")  # 如果未找到最优解，输出信息

---