GEKKO实时MPC时间窗移位失效问题求助（非零起始被重置）

核心原因 #

GEKKO内部强制时间轴起始点为0，所以你直接把m.time设成非零开头时，求解器会自动把时间轴拉回从0开始，这就是那个警告的由来。直接扩展时间窗会让问题规模越来越大，计算自然越来越慢，这路子不对。

正确实现方式 #

保持m.time为固定长度的相对时间窗（比如预测步长20，就用np.linspace(0,20,21)），通过以下步骤实现滚动优化：

1. 固定时间窗初始化
   一开始就定义好固定长度的预测时间窗，后续循环里不用改：

   prediction_horizon = 20
   m.time = np.linspace(0, prediction_horizon, prediction_horizon + 1)
   

2. 循环中更新初始状态
   每次循环开始，把当前实际系统的状态值赋值给GEKKO的状态变量，作为下一次优化的起点：

   # 假设x1、x2是你的状态变量
   x1.value = current_system_x1
   x2.value = current_system_x2
   

   如果是动态系统，也可以用m.fix_initial()固定初始状态，确保每次优化都从当前真实状态开始。

3. 处理时间相关扰动/参数
   如果有随时间变化的参考轨迹、外界干扰这类参数，别去扩展数组长度，直接循环移位数组：把新时刻的数据放到数组末尾，同时删掉最旧的第一个数据：

   # 假设ref是和m.time长度一致的参考轨迹数组
   # 先获取当前时刻的新参考值new_ref_val
   ref = np.roll(ref, -1)
   ref[-1] = new_ref_val
   # 再把更新后的ref赋值给GEKKO的参数变量
   m.Param(value=ref)
   

4. 求解并应用第一个控制动作
   每次求解完成后，只取优化结果里的第一个控制值（对应当前时间步的动作），给到实际系统，然后进入下一轮循环：

   m.solve(disp=False)
   # 假设u是控制变量
   current_u = u.value[0]
   # 把current_u发送给实际系统执行
   

完整示例片段 #

import numpy as np
from gekko import GEKKO

# 初始化GEKKO模型
m = GEKKO(remote=False)
prediction_horizon = 20
# 固定时间窗，后续不再修改
m.time = np.linspace(0, prediction_horizon, prediction_horizon + 1)

# 定义状态、控制变量
x = m.Var(value=0)
u = m.MV(value=0, lb=-1, ub=1)
u.STATUS = 1

# 示例一阶系统动态方程
m.Equation(x.dt() == -x + u)

# 目标函数：跟踪参考轨迹+控制量惩罚
ref = np.zeros(prediction_horizon + 1)
ref_param = m.Param(value=ref)
m.Obj((x - ref_param)**2 + 0.1*u**2)

# 实时循环模拟
current_time = 0
current_x = 0  # 这里模拟实际系统的状态，真实场景替换为读取硬件数据
while True:
    # 更新优化初始状态
    x.value = current_x
    
    # 更新参考轨迹（示例：随时间变化的正弦参考）
    new_ref = np.sin(current_time / 10)
    ref = np.roll(ref, -1)
    ref[-1] = new_ref
    ref_param.value = ref
    
    # 求解优化问题
    m.solve(disp=False)
    
    # 获取当前控制动作，模拟系统状态更新
    current_u = u.value[0]
    current_x += (-current_x + current_u) * 1  # 步长为1，真实场景替换为硬件反馈
    
    current_time += 1
    # 可添加延时模拟真实实时环境
    # import time
    # time.sleep(1)

关键说明 #

• 这种方式保持优化问题的规模始终固定，不会随着运行时间推移变慢。
• GEKKO的时间轴始终是相对的（从0到预测步长），通过更新初始状态和参数数组来模拟时间滚动，完全符合实时MPC的逻辑。
• 既避开了修改m.time导致的自动移位问题，也不用扩展数组长度。

内容的提问来源于stack exchange，提问作者Mohammed Uzair Khan