Python linprog 最大化目标函数

Question

我已经有一段时间没有这样做了，所以我有点生疏，但是等式是：

max t(C)*x
s.t. Ax <=b

我有我的 A 约束矩阵，即 (1448x1359)：

[[ 1.  1.  0. ...,  0.  0.  0.]

 ..., 


 [ 0.  0.  0. ...,  1.  1.  1.]]

然后我有我的绑定 b (1448x1)：

[ 1.  1.  7. ...,  2.  1.  2.]

我要最大化的目标函数是一个向量 (1359,1)。

现在在其他包中，我的最大化目标函数是 841，但是使用 linprog：

res = linprog(c=OBJ_N, A_ub=A, b_ub=b, options={"disp": True})

它成功地优化到了​​ -0.0，所以我想知道我是否在 python 中使用了正确的命令并且我的约束是正确的？

编辑：好吧，这是有道理的，它试图最小化。我现在已经重写了（交换了 c 和 b 并转置了 A 以最小化）。

# (max t(C)*x s.t. Ax <=b) = min t(b)*x s.t. ATy = c, y ≥ 0
# (i): minimise number of shops no bounds
ID = np.ones(len(w[0]))
print(ID)
print(ID.shape)  #1359

At = A.transpose()

need_divest = (A.dot(ID)) - 1
print(need_divest)
print(need_divest.shape)  #1448

res = linprog(c=need_divest, A_eq=At, b_eq=ID, options={"disp": True})
print(res)

但是，我收到“消息：'优化失败。无法找到可行的起点。'”

Answer 1

我猜你可能是minimizing 而不是maximizing 你的目标函数。 试试这个（在你的目标函数系数前面插入一个 -）：

res = linprog(c=-OBJ_N, A_ub=A, b_ub=b, options={"disp": True})

你的结果应该是 -841。

这很有效，因为：

min(f(x))=-max(-f(x))