Модель Гуроби становится невозможной после введения набора переменных - PullRequest
Новая
       10

Модель Гуроби становится невозможной после введения набора переменных

0 голосов
/ 30 июня 2019

Я пытаюсь смоделировать модель оптимизации, используя gurobi в python.Модель становится невозможной, когда я добавляю в нее набор переменных, даже если нет ограничений, связанных с этими переменными.Проблемный код: 'queue_length = model.addVars (number_of_movememts, total_time + 1, name =' queue_length ')'. 

from gurobipy import *

order_of_miu = 3
movement_to_phase_dict = {0:[0],1:[1]}
saturated_flow_rate = [1, 2/3]
number_of_phase = 2
number_of_movememts = 2
C_min = 60
C_max = 180
lost_time_per_phase = 4
lost_time_per_cycle = lost_time_per_phase * number_of_phase
total_time = 1500


def OptimizationModel():
    model = Model('Signal Optimization')

    gr_min = (C_min - lost_time_per_cycle) / C_min
    gr_max = (C_max - lost_time_per_cycle) / C_max


    gr_coeff = model.addVars(number_of_phase,order_of_miu + 1, lb = -1 * GRB.INFINITY, name = 'gr_coeff')

    queue_length = model.addVars(number_of_movememts,total_time + 1, name = 'queue_length')

    model.addConstrs((sum([gr_coeff[p,k]*t**k for p in range(number_of_phase) 
        for k in range(order_of_miu + 1)]) >= gr_min for t in range(total_time + 1)), "C7-1")

    model.addConstrs((sum([gr_coeff[p,k]*t**k for p in range(number_of_phase) 
        for k in range(order_of_miu + 1)]) <= gr_max for t in range(total_time + 1)), "C7-2")

    model.optimize()


    if model.status == GRB.Status.OPTIMAL:
        solution = model.getAttr('x',gr_coeff)
        coeffs_of_all_phase = [[solution[p,k] for k in range(order_of_miu + 1)] for p in range(number_of_phase)]

        return coeffs_of_all_phase

    else:
        return -1


if __name__ == "__main__":
    coeffs_of_all_phase = OptimizationModel()
...