Я написал код оптимизации Gorubi, но из-за некоторых проблем мне нужно преобразовать его в код Scipy. Все еще есть трудности, чтобы преобразовать это. Вот часть кода, связанная с Gorubi:
m = Model()
#x is charging, discharging variable
x = m.addVars(n,lb=-1.5,ub=1.5,vtype=GRB.INTEGER, name="x")
#Y is SOC variable
Y = m.addVars(n+1,lb=0,ub=100,vtype=GRB.CONTINUOUS, name="Y")
# Add constraint: SOC[start]=50, initial SOC
m.addConstr(Y[0]==initialsoc,name='c1')
#Final targeted SOC
m.addConstr(Y[n]>=65,name='c2')
m.addConstrs((Y[i+1]-Y[i] == 3.75*x[i] for i in range(n)), name='c0')
#Objective function. 6 comes from capacity of inverter.
obj1=sum(((load[i+1]-(6*x[i]))*(load[i+1]-(6*x[i])) for i in range (n)))
m.setObjective(obj1,GRB.MINIMIZE)
m.optimize()
Мое ограничение x может иметь только значения -1, 0 или 1. Другим ограничением является Y, где на каждом шаге y (i + 1) -y (i) равно 3,75 * X (i).
Возможно ли преобразовать этот код в Scipy? Или вы рекомендуете другие библиотеки?