GLPK/GMPL：最小化目标函数，但将其保持在 0 >

Question

我正在尝试最小化具有三个参数的目标函数：i，p，j，如下所示：

param mlu{i in I, p in P, j in out[p]} := traffic[i,p]/capacity[j];
minimize MAXLU{i in I, p in P, j in out[p]}: mlu[i,p,j] * x[i,p,j];

但是目标函数必须大于 0，否则它违背了我最小化的目的。

我试图通过在目标函数上添加约束来确保这一点，如下所示：

s.t. constraint1{i in I, p in P, j in out[p]} : MAXLU[i,p,j] != 0;

但是我收到以下错误：

LP.mod:66: invalid reference to status, primal value, or dual value of objective MAXLU above solve statement
Context:  i in I , p in P , j in out [ p ] } : MAXLU [ i , p , j ] !=
glp_mpl_generate: invalid call sequence
Error detected in file glpapi14.c at line 79
Aborted

甚至有可能做到这一点吗？感谢您的任何帮助/建议！

Answer 1

您可以在约束中引用变量和参数，但不能引用目标（MAXLU）。解决此问题的一种简单方法是将constraint1中的MAXLU替换为实际的目标表达式：

s.t. constraint1{i in I, p in P, j in out[p]} : mlu[i,p,j] * x[i,p,j] ...;

此外，您不能使用表达式 ！= 0 定义约束。您可以使用>= epsilon 或 <= epsilon，其中 epsilon 是一个较小的数字：

s.t. constraint1{i in I, p in P, j in out[p]} : mlu[i,p,j] * x[i,p,j] >= epsilon;