Получение решения ограничений, когда модель дает недопустимое решение

Я кодирую проблему оптимизации. Модель дает недопустимое решение. Я хочу проверить, какое ограничение дает недопустимое решение. Пока что я проверял онлайн, но не смог найти решения проблемы. Может кто-нибудь мне помочь? Например: в приведенном ниже коде модель невозможна из-за ограничения 3. Как определить это из решения? Спасибо

from gurobipy import *

# Create a new model
m = Model("mip1")

# Create variables
x1 = m.addVar(vtype=GRB.INTEGER, name="x1")
x2 = m.addVar(vtype=GRB.INTEGER, name="x2")


# Integrate new variables
m.update()

# Set objective
m.setObjective(7*x1 + 2*x2, GRB.MAXIMIZE)


m.addConstr(-x1 + 2 * x2 <= 4, "constraint-0")


m.addConstr(5*x1 + x2 <= 20, "constraint-1")
m.addConstr(-2*x1 -2*x2 <= -7, "constraint-2")
m.addConstr(x1 <= -2, "constraint-3")
m.addConstr(x2 <= 4, "constraint-4")

m.optimize()

for v in m.getVars():
    print('%s %g' % (v.varName, v.x))

print('Obj: %g' % m.objVal)

python-2.7 optimization gurobi
person Sudipta    schedule 21.11.2016    source источник
comment
Используйте это.   -  person sascha    schedule 22.11.2016
comment
@sudipta, почему x1 должно быть меньше или равно -2? Просто из любопытства   -  person user2567806    schedule 31.10.2017

Ответы (1)


arrow_upward
4
arrow_downward

пример:

from gurobipy import *

# Create a new model
m = Model("mip1")
# Create variables
x1= m.addVar(lb=0,ub=62,vtype=GRB.INTEGER,name="x1")
x2 = m.addVar(lb=0,ub=50, vtype=GRB.INTEGER,name="x2")
m.update()

m.addConstr(-x1 + 2*x2 <= 4, "constraint-0")
m.addConstr(5*x1 + x2 <= 20, "constraint-1")
m.addConstr(-2*x1 -2*x2 <= -25, "constraint-2")
m.addConstr(x1 <= 2, "constraint-3")
#m.addConstr(x2 <= 50, "constraint-4")

m.update()

# Set objective
m.setObjective(7*x1 + 2*x2, GRB.MAXIMIZE)
m.update()

m.optimize()

status = m.status
if status == GRB.Status.OPTIMAL:
    for v in m.getVars():
        print('%s %g' % (v.varName, v.x))
    print('Obj: %g' % m.objVal)
       
elif status == GRB.Status.INFEASIBLE:
    print('Optimization was stopped with status %d' % status)
    # do IIS
    m.computeIIS()
    for c in m.getConstrs():
        if c.IISConstr:
            print('%s' % c.constrName)

person Issouf    schedule 23.11.2016