pyomo - 获取 Pyomo MIP 中约束的边际值（对偶）

Question

我想使用 Pyomo 在 python 中开发的 MIP 问题访问对偶变量。据我了解，对偶不是为 MIP 问题创建的，但我认为应该有一个解决方法。

这应该可以用作一个最小的工作示例，我自己正在使用 Gurobi。

我可以看到两种可能的解决方案；1. 修复二进制/整数变量并解析为 LP 并重新创建对偶。2. 仅在必要的约束条件下获得对偶。

我一直无法想出尝试第二种方法的方法，首先我做了类似的事情：

m.solve()                      #solves the MIP problem
m.instance.x.fix()             #fixing the binary variables (only have one in actual model)
m.instance.y.fix()
m.instance.z.fix()
m.instance.preprocess()   
m.instance.dual = Suffix(direction=Suffix.IMPORT) 
m.solve()                      #hopefully solving LP problem with dual

如果有任何不清楚或需要更多信息，请告诉我，我们将不胜感激。