Gekko求解器无法找到x2全局最大值的技术咨询

我正在学习使用Gekko作为动力学数据的优化工具，基于Gekko官网示例实现了一个简单的ODE系统。尝试通过m.Maximize(x2)寻找组分x2的全局最大值，但求解器始终无法找到，仅当强制dx2/dt>0时才能得到正确的最大值。我想知道是否存在设置错误，或是我的解决方案存在问题？

代码如下：

from gekko import GEKKO 
import numpy as np
import matplotlib.pyplot as plt

m = GEKKO()  # initialize gekko

nt=200
m.time = np.linspace(0,1,nt)
# Variables
x1 = m.Var(value=20,lb=0,)
x2 = m.Var(value=0,lb=0,)
x3 = m.Var(value=0,lb=0,)

p = np.zeros(nt)  # mark final time point
p[-1] = 1.0

# optimize final time
tf = m.FV(value=0.7,lb=0.001,ub=100.0)
tf.STATUS = 1

# Equations
m.Equation(x1.dt() == -5*x1*tf)
m.Equation(x2.dt() == (5* x1 - 0.5*x2)*tf)
m.Equation(x3.dt() == 0.5*x2*tf)


#Minifunc
#m.Minimize(tf)
m.Maximize(x2)

#options
m.options.IMODE = 6
m.options.SOLVER = 3 #IPOPT optimizer
m.options.RTOL = 1E-8
m.options.OTOL = 1E-8
m.options.NODES = 5
m.solve(disp=True,debug=1)

print('Final Time: ' + str(tf.value[0]))
#Plotting stuff 
tm = np.linspace(0,tf.value[0],nt)

plt.figure(1)  # plot results
plt.plot(tm, x1.value, "k-", label=r"$x_1$")
plt.plot(tm, x2.value, "b-", label=r"$x_2$")
plt.plot(tm, x3.value, "r--", label=r"$x_3$")

plt.legend(loc="best")
plt.xlabel("Time")
plt.ylabel("Value")
plt.show()

结果显示求解器明显未找到最大值。

核心问题 #

直接调用m.Maximize(x2)时，Gekko默认只会最大化最终时刻的x2值，而非整个时间区间内的全局最大值。这就是为什么你只有强制dx2/dt>0时才能得到预期结果——此时求解器会让x2持续增长到约束边界，但这并非真正的全局最大值寻优逻辑。

解决方法 #

要寻找整个时间区间内x2的全局最大值，需要通过以下方式调整：

1. 定义全局最大值变量：新增一个变量x2_max，约束它大于等于所有时刻的x2值，再对这个变量进行最大化操作。
2. 添加逐时间点约束：通过循环为每个时间点设置x2_max >= x2[t]的约束，确保x2_max能覆盖x2的所有峰值。

修改后的关键代码片段：

# 新增变量存储x2的全局最大值
x2_max = m.Var(value=0, lb=0)
# 约束x2_max >= 所有时刻的x2值
for t in range(nt):
    m.Equation(x2_max >= x2[t])
# 最大化全局最大值变量
m.Maximize(x2_max)

额外优化建议 #

• 保留tf作为优化变量时，求解器会自动调整时间尺度来找到最优的x2_max，无需额外添加dx2/dt>0的约束。
• 若求解出现收敛问题，可切换到SOLVER=1（APOPT求解器），它在处理非凸优化场景时稳定性更强；也可尝试将NODES参数调整为3，在保证精度的同时减少计算量。

内容的提问来源于stack exchange，提问作者toomany questions