Python GEKKO if3 变量在优化中给出了意想不到的结果

m.if3()函数将默认解算器从IPOPT更改为APOPT，以便可以解决混合整数优化问题.APOPT solver发现了一个不同(更差)的局部解，因此您看到了差异.

此外，随机数生成器每次都需要一个种子值来给出相同的结果，例如np.random.seed(1).与m.if3()函数相关的附加变量和方程也可能需要更多的迭代.try 对用IPOPT解决方案初始化的APOPT使用热启动.

无论m.if3()条语句是否存在，该脚本都给出了相同的答案.

from gekko import Gekko
import numpy as np
m = Gekko()

np.random.seed(1)
electricity_price_in = np.random.uniform(low=0.1, high=1, size=50)
electrictiy_price_out = np.random.uniform(low=0.3, high=3, size=50)

E_battery = m.Var(lb=0, ub=366.2, value=0) #energy in battery at time t, battery size 366 MWh
Pc = m.Var(lb=0, ub=50) #charge power, 50 MW max
Pd = m.Var(lb=0, ub=36.6)  #discharge power, max 36 MW
E_price_in = m.Param(electricity_price_in)
E_price_out = m.Param(electrictiy_price_out)
m.time = np.linspace(0,49, 50)

Revenue =  m.Intermediate(Pd*E_price_out)
Cost = m.Intermediate(Pc*E_price_in)

m.Equation(E_battery.dt() == Pc-Pd)
m.Equation(E_battery >= 0)

m.Maximize(Revenue - Cost)

isCharging = m.if3(-Pc,1,0)
isDischarging = m.if3(-Pd,1,0)

m.options.IMODE = 6

m.options.SOLVER = 3
m.solve(disp=False)
print('Objective (IPOPT)', m.options.objfcnval)

m.options.TIME_SHIFT = 0
m.options.SOLVER = 1
m.solve()
print('Objective (APOPT)', m.options.objfcnval)