GEKKO迭代仿真控制中约束违例错误的原因与优化方案问询

问题背景 #

我实例化了两个GEKKO对象：Simulation（IMODE=4）和Controller（IMODE=6），通过循环迭代运行：控制器设置仿真器的独立变量，下一周期读取仿真器的状态变量，目的是评估各类模型误差下的控制器响应，后续还要演示Estimator对象在自适应控制中的应用。

控制器正常运行37次迭代后，触发错误：

Restoration phase is called a point that is almost feasible, with constraint violation 8.305823e-009. Abort.

尝试设置SCALING=2后，可运行至42次迭代，但仍出现相同错误（约束违例数值不同）；调整MV的DMAX等参数仅能延后错误出现时间。

模型包含两个等式约束：

m.Equation(m.TankLevel.dt()==m.Balance)  
m.Equation(m.Flare==m.Rundown-(m.Unit1_Feed+m.Unit2_Feed+m.Fuel))

判断为数值精度问题后，将第二个等式约束改为不等式：

m.Equation(m.Flare<=m.Rundown-(m.Unit1_Feed+m.Unit2_Feed+m.Fuel)+10e-8)

问题得以解决。现询问该错误的根本原因及更优解决方式。

版本信息：GEKKO 1.0.5, APMonitor 1.0.0

根本原因 #

这个错误的核心是循环迭代中的数值精度累积导致等式约束微小违例：

• GEKKO依赖APMonitor求解器处理优化问题，求解器默认会以严格的可行性标准（通常是1e-8量级的绝对/相对误差）检查约束满足情况。
• 在动态循环迭代场景中，每一步的状态变量传递、浮点运算都会产生微小误差，这些误差累积到一定程度后，会让第二个物料平衡等式Flare = Rundown - (Unit1_Feed+Unit2_Feed+Fuel)的两边出现极微小的偏差（如8.3e-9），刚好超出求解器的可行阈值，触发恢复阶段；若恢复阶段无法将偏差拉回容忍范围内，就会抛出中止错误。

更优解决方式 #

1. 放宽求解器精度容忍度 #

直接调整求解器的相对容忍度（RTOL）和绝对容忍度（ATOL），允许微小的约束偏差，无需修改约束结构：

# 示例：将容忍度从默认1e-8放宽到1e-7
m.options.RTOL = 1e-7
m.options.ATOL = 1e-7

这种方式不改变模型的物理意义，是最直接的解决方案。

2. 引入松弛变量软化等式约束 #

如果必须保留物料平衡的逻辑等式，可通过松弛变量允许微小偏差，并通过目标函数惩罚偏差，兼顾约束刚性和数值稳定性：

# 定义松弛变量，限制其最大偏差范围
slack = m.Var(lb=-1e-6, ub=1e-6)
# 将等式约束转为带松弛变量的形式
m.Equation(m.Flare == m.Rundown - (m.Unit1_Feed+m.Unit2_Feed+m.Fuel) + slack)
# 添加惩罚项，让求解器尽量缩小偏差
m.Obj(1e6 * slack**2)

3. 优化变量缩放 #

除了自动缩放（SCALING=2），手动统一变量量级能减少浮点运算的精度损失：

• 检查所有变量的量级，若存在量级差异过大的变量（比如Rundown是1e4、Flare是1e2），可对变量进行缩放（如Rundown_scaled = m.Param(Rundown/100)），让所有变量的取值范围尽量集中在1~100之间，提升求解器的计算稳定性。

4. 控制状态变量传递精度 #

在循环中传递仿真器的状态变量给控制器时，对变量进行适当舍入（如保留7~8位有效数字），避免累积过多的微小浮点误差：

# 示例：读取仿真器状态后舍入到8位有效数字
tank_level = round(m.TankLevel.value[0], 8)
controller.TankLevel.MEAS = tank_level

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