如何为Pyomo变量设置多组边界对？遇元组类型错误

Pyomo的Var组件的bounds参数不支持直接传入多组边界对——它只接受单个(下限, 上限)格式的元组，或者能返回该格式元组的函数，这就是你传入[(0,0.1),(0.2,1)]时触发tuple类型错误的原因。

要实现变量仅在0-0.1或0.2-1区间取值的需求，需要通过添加约束+引入二进制辅助变量的方式来实现，具体操作如下：

1. 定义变量 #

先创建主变量和用来控制区间的二进制辅助变量：

# 主变量，可先设置覆盖所有目标区间的宽泛边界（也可以不设）
model.x = Var(bounds=(0, 1))
# 二进制变量y：y=0时x落在0-0.1区间；y=1时x落在0.2-1区间
model.y = Var(within=Binary)

2. 添加区间控制约束 #

通过两组约束锁定变量的取值范围：

# 当y=0时，x被限制在0-0.1；y=1时，该约束不对x上限做限制
model.constr1 = Constraint(expr=model.x <= 0.1 + 0.9*model.y)
# 当y=1时，x被限制在0.2-1；y=0时，该约束不对x下限做限制
model.constr2 = Constraint(expr=model.x >= 0.2*model.y)

也可以用更直观的分段逻辑写法：

# 约束1：y=0时x≤0.1；y=1时x≤1
model.constr_low_upper = Constraint(expr=model.x <= 0.1 + (1-0.1)*model.y)
# 约束2：y=0时x≥0；y=1时x≥0.2
model.constr_high_lower = Constraint(expr=model.x >= 0.2*model.y)

这样就能确保x只会落在你指定的两个区间内。

注意事项 #

引入二进制变量后，模型会从线性规划（LP）变为混合整数线性规划（MILP），需要使用支持MILP的求解器，比如CBC、Gurobi、CPLEX等。如果是非线性模型，逻辑思路一致，只需调整约束表达式即可。

内容的提问来源于stack exchange，提问作者Yari