提问者:小点点

Pyomo约束中的For循环


我正在努力在Python Pyomo中使用for循环制作多个约束。我的代码概念如下

T = 504
model.times = RangeSet(0, T-1)
model.Machines = Set(initialize = Machines_parameters.keys())

然后我划分了模型。在3个相同长度的集合中执行时间,现在我想为集合的每个部分编写特定的约束。以下是代码的等效简化版本:

for k in range(3): #number of parts in which I've divided the set [(0,167),(168,335),(336,503)]
    for j in range(cluster_index[k][0], cluster_index[k][1]):   #respectively first and last number of the subset
        def logic_constr(model,i):
            if j >= const1[k] and j < const2[k]:
                return model.z[i, j + 1] - model.z[i, j] == model.dsu[i, j + 1] - model.dsd[i, j + 1]
            else j==const2[k]:
                return model.z[i,const2[k]] - model.z[i,j] == model.dsu[i,const1[k]] - model.dsd[i,const1[k]]
        model.logic_constr = Constraint(model.Machines, rule = logic_constr)

我想做的是迭代创建504个不同的约束,每个约束都有自己的规则。你对怎么做有什么建议吗?


共1个答案

匿名用户

您现在制定约束的方式,最终只会在最后有一个约束,因为在每个循环之后,约束都会被覆盖。

因为,正如您所说,根据您的for循环,您也希望每个时间步有一个约束,所以这要简单得多。

首先,您需要将约束定义为:

model.logic_constr = Constraint(model.Machines, model.times, rule = logic_constr)

这意味着约束将应用于集合模型的每个成员。时间,即集合模型的每个元素将有504个约束。机器。

然后,您所有的ifs和for循环都可以进入您对logic_constr函数的定义。