【问题标题】:AMPL finds optimal solution zero, all variables are set to zeroAMPL 找到最优解为零,所有变量都设置为零
【发布时间】:2016-04-10 06:43:34
【问题描述】:

这是我第一次使用 AMPL,我对它并不是很熟悉。 我必须优化我的模型,但我认为我的代码是错误的,我一直得到最佳解决方案和“NOS 5.51:忽略 20 个变量的完整性”。 我的模型是 最小(废物) 受制于 I attached a picture of constraints

这是我的代码:

`set DAY;     #the day we buy ingredient
 set INGRED;   #fresh ingredients
 set TIME;     #day that keeps track of the inventory
 param M default 0;
 param cost{INGRED} > 0; #cost for each ingredient (per pound)
 param demand{TIME, INGRED} >= 0; #the expected demand for each ingredient fo each day t
 param min_pur_req > 0;  #minimum total cost of order required to get a delivery
 param expiry{INGRED} ; #shelf life of each ingredient
 var amount{t in TIME, i in DAY, j in INGRED} >=0; #is defined only for t = i
 var is_bought {t in TIME,i in DAY,j in  INGRED : i <= t }binary; #only for t =i
 var inventory {t in TIME,i in DAY, j in INGRED : i<=t <= i+expiry[j] } >= 0; # t >= i
 var used{t in TIME,i in DAY, j in INGRED :  i <= t <= i+expiry[j]} >= 0; 
 var waste{t in TIME,i in DAY, j in INGRED : t <= i+expiry[j] } >= 0;`
 minimize Total_waste: 
        sum{t in TIME,i in DAY, j in INGRED : i <= t <= i+expiry[j]} waste[t,i,j];
subject to invalid_amount{t in TIME, i in DAY, j in INGRED: i <= t <= i+expiry[j]}: 
      if t != i then amount[t,i,j] = 0;
subject to invalid_binary{t in TIME, i in DAY, j in INGRED: i <= t <= i+expiry[j]}:
     if t != i then is_bought[t,i,j] = 0;
subject to usage{t in TIME, j in INGRED}:
    sum {i in DAY: i <= t <= i + expiry[j]} used[t,i,j] = demand[t, j];
subject to inventory_formula {t in TIME, i in DAY, j in INGRED: i+1 <= t <= i+expiry[j]}: 
    waste[t,i,j] = amount[t,i,j] + inventory[t-1,i,j] - used[t,i,j] - inventory[t,i,j] ;
subject to cost_check {t in TIME, i in DAY: t = i}:
    sum {j in INGRED} cost[j]*amount[t,i,j] >=  sum{j in INGRED: t <= i+expiry[j]} min_pur_req*is_bought[t,i,j];
subject to blah{t in TIME,i in DAY, j in INGRED: i <= t <= i+expiry[j]}:
        amount[t,i,j] <= M*is_bought[t,i,j];
subject to  invalid_inventory{t in TIME, i in DAY, j in INGRED: i+1 <= t <= i+expiry[j]}: 
 #i+1 because for i = 1 i get error in invetory[0,1,j] 
      if t=i then inventory[t-1,i,j] = 0;

这就是我得到的:

ampl: model waste.mod;
ampl: data waste.dat;
ampl: solve;
MINOS 5.51: ignoring integrality of 20 variables
MINOS 5.51: optimal solution found.
23 iterations, objective 0

【问题讨论】:

    标签: binary continuous ampl


    【解决方案1】:

    MINOS 求解器不支持二进制变量,仅支持连续变量。因此,它忽略(放松)完整性并通过消息“忽略 20 个变量的完整性”来警告它。您应该使用混合整数规划 AMPL 求解器,例如 CPLEXGurobi 来获得整数解。

    【讨论】:

      猜你喜欢
      • 1970-01-01
      • 1970-01-01
      • 1970-01-01
      • 2012-11-11
      • 1970-01-01
      • 2012-04-18
      • 1970-01-01
      • 1970-01-01
      • 2012-05-09
      相关资源
      最近更新 更多