Python Gekko中MPC问题的分段式PID设定点实现问询

完全可以实现你需要的分段式设定点，核心是通过约束多个时间步的设定点等于同一个FV变量，让MPC在15分钟步长下优化2小时间隔的设定值。下面是具体实现步骤和错误原因分析：

1. 核心思路 #

• 保留97个15分钟的时间步长，将PID设定点（SP）拆分为12个分段（每8步为一段，最后一段覆盖剩余1步），每个分段对应一个FV变量
• 逐个时间步建立约束，让同一段内的所有SP值等于对应的FV
• 状态变量T₁/T₂是随时间步变化的数组，必须索引到具体步长来计算PID误差

2. 具体代码实现 #

2.1 划分分段时间步索引 #

先明确每个FV对应的时间步范围：

dt = 15  # 分钟步长
n_steps = 97  # 总时间步数
segment_len = 8  # 每段包含的步数
n_segments = (n_steps + segment_len - 1) // segment_len  # 向上取整得到12个分段

# 生成每个分段对应的时间步索引列表
segment_indices = []
for i in range(n_segments):
    start = i * segment_len
    end = min(start + segment_len, n_steps)
    segment_indices.append(range(start, end))

2.2 定义分段FV变量 #

创建两个区域的分段设定点FV，允许MPC优化这些值：

# 区域1、2的分段设定点FV（共12个）
SP1_seg = m.Array(m.FV, n_segments)
SP2_seg = m.Array(m.FV, n_segments)

# 初始化并开启优化状态
for sp in SP1_seg:
    sp.value = 20.0  # 初始设定20℃
    sp.STATUS = 1    # 允许MPC优化该值
    sp.LOWER = 18    # 最低设定温度
    sp.UPPER = 24    # 最高设定温度
for sp in SP2_seg:
    sp.value = 20.0
    sp.STATUS = 1
    sp.LOWER = 18
    sp.UPPER = 24

2.3 绑定FV到每个时间步的SP #

创建与时间步等长的SP变量，通过等式约束让同一段内的SP等于对应的FV：

# 每个时间步的设定点Var数组
SP1 = m.Array(m.Var, n_steps)
SP2 = m.Array(m.Var, n_steps)

# 绑定分段FV到对应时间步的SP
for seg_idx, time_indices in enumerate(segment_indices):
    for t_idx in time_indices:
        m.Equation(SP1[t_idx] == SP1_seg[seg_idx])
        m.Equation(SP2[t_idx] == SP2_seg[seg_idx])

2.4 修正PID误差与模型方程 #

现在可以逐个时间步计算PID误差，避免之前的数组运算错误：

# 假设T1、T2是已定义的状态变量数组（长度n_steps）
err1 = m.Array(m.Intermediate, n_steps)
err2 = m.Array(m.Intermediate, n_steps)

# 逐时间步计算误差
for t_idx in range(n_steps):
    err1[t_idx] = SP1[t_idx] - T1[t_idx]
    err2[t_idx] = SP2[t_idx] - T2[t_idx]

# 后续基于err1、err2构建PID逻辑，以及RC模型微分方程
# 示例RC模型方程（以T1为例）
# m.Equation(T1.dt() == (Q1 - (T1-T2)/R12 - (T1-T_out)/R1o) / C1)

2.5 目标函数与MPC配置 #

保持原有的能耗+最大功率优化目标：

# 计算总能耗（转换为kWh）
total_energy = m.Intermediate(m.sum(Q1 * dt/60) + m.sum(Q2 * dt/60))
# 计算最大加热功率
max_power = m.Intermediate(m.max3(*Q1, *Q2))

# 加权组合目标函数（权重可根据需求调整）
m.Obj(total_energy + 0.1 * max_power)

# MPC模式配置
m.options.IMODE = 6
m.options.NODES = 2
m.options.SOLVER = 3

3. 之前的错误原因 #

你遇到的Intermediate Definition Error是因为直接用数组相减（SP_1-T_1），Gekko无法自动识别数组的逐元素运算，必须显式索引到每个时间步来定义Intermediate变量。

4. 额外优化建议 #

• 如果不需要单独的SP Var数组，可直接在误差计算时引用FV，比如err1[t_idx] = SP1_seg[seg_idx] - T1[t_idx]，减少变量数量
• 可以给FV添加DCOST参数，让设定点的变化更平滑，避免大幅跳变：SP1_seg[i].DCOST = 0.1

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