📚 课程目录第3章多变量系统建模

3. 多变量系统建模

多变量系统建模是解耦控制与多变量预测控制的第一道工程关口。单回路控制中，我们习惯用“一个操纵变量对应一个被控变量”的直觉来整定 PID；但在多输入多输出（MIMO）过程中，任何一个操纵变量的变化都会通过过程内部的耦合通道影响多个被控变量。若建模阶段忽略这种耦合结构，后续无论采用静态解耦、动态解耦还是模型预测控制，都会因模型失配而出现振荡、超调甚至失稳。因此，本节从过程描述、传递函数矩阵、状态空间模型三条主线出发，说明机理建模与辨识建模各自适用的场景与衔接方式。

一、多输入多输出过程描述：从“通道”到“矩阵”

一个典型的多变量过程可写成 y(s) = G(s)u(s) + G_d(s)d(s)。其中 u 为操纵变量向量，y 为被控变量向量，d 为可测或不可测扰动向量。G(s) 是 n×m 的传递函数矩阵，其第 i 行第 j 列元素 g_ij(s) 表示第 j 个操纵变量到第 i 个被控变量的通道动态。对角元素是“主通道”，非对角元素是“耦合通道”。

工程上判断一个过程是否“强耦合”，不能只看非对角元素是否为零，而要看相对增益矩阵（RGA）。RGA 由 Λ = G(0) ⊗ (G(0)^{-1})^T 计算，元素 λ_ij 表示在其他回路开环时，第 j 个操纵量对第 i 个被控量的相对影响。若 λ_ij 接近 1，说明配对合理；接近 0 表示该通道几乎无控制作用；若出现负值或远大于 1，则说明耦合严重，需要重新配对或引入解耦器。

建模的第一项任务不是拟合参数，而是确定“谁控制谁”的配对结构。配对错误会让再精确的模型也无法转化为稳定的控制效果。

二、传递函数矩阵：频域描述与耦合可视化

传递函数矩阵是频域建模的核心表达。对于 2×2 过程，可写为：

G(s) = [ g11(s)  g12(s) ]
       [ g21(s)  g22(s) ]

g11(s) = K11 / (T11 s + 1) · exp(-τ11 s)
g12(s) = K12 / (T12 s + 1) · exp(-τ12 s)
...

每个通道通常由增益 K、时间常数 T、纯滞后 τ 描述。耦合通道的增益符号可能与主通道相反，这会直接影响解耦器的稳定性。例如在精馏塔中，回流量增加会使塔顶温度下降，而再沸器蒸汽量增加会使塔底温度上升，但两者通过塔压和组分分布相互影响，形成典型的 2×2 耦合结构。

传递函数矩阵的优点是物理意义清晰、便于频域分析和解耦器设计；缺点是当过程非线性强、工况变化大时，单一线性矩阵难以覆盖全操作范围。此时需要按工况分段建模，或转向状态空间与辨识建模。

三、状态空间模型：内部机理的统一表达

状态空间模型把过程内部储能、物料平衡和能量平衡统一写成：

dx/dt = A x + B u + E d
y     = C x + D u + F d

其中 x 为状态向量，A 描述状态自身演化，B 描述操纵变量对状态的作用，C 描述状态到输出的映射。与传递函数矩阵相比，状态空间模型能自然处理多变量、非线性、初始条件和内部不可测状态，也便于扩展为扩展卡尔曼滤波、模型预测控制等现代控制结构。

在机理建模中，状态空间形式尤其适合化工过程：物料平衡给出 dx/dt，能量平衡给出温度状态方程，反应动力学给出组分变化率。通过线性化工作点，可得到局部线性状态空间模型，再转换为传递函数矩阵用于解耦设计。

四、机理建模与辨识建模：两条路径的衔接

机理建模从守恒定律出发，参数具有物理意义，外推能力较强，但建模周期长、部分参数难以获取。辨识建模通过阶跃响应、伪随机二进制序列（PRBS）或正弦扫频实验，从输入输出数据中拟合模型参数，适合复杂工业过程，但依赖激励信号设计和数据质量。

对比维度	机理建模	辨识建模
基础	质量、能量、动量守恒	输入输出实验数据
参数意义	物理量明确	拟合参数，需校验
外推能力	较强	受激励范围限制
典型方法	状态空间、微分方程	阶跃响应、ARX、子空间辨识
适用场景	新工艺、安全关键过程	已有装置、数据充足过程

实际工程中常采用“机理定结构、辨识定参数”的混合策略：先用机理分析确定模型阶次、耦合通道和主要时间常数，再用现场阶跃或 PRBS 数据辨识增益和滞后。这样既保留物理可解释性，又提高模型对实际装置的贴合度。

建模精度不必追求“全频段完美”。对解耦控制而言，关键频段（穿越频率附近）的幅值与相位误差才是决定闭环稳定性的核心。

五、案例：2×2 耦合过程的建模与校验

某换热网络中被控变量为出口温度 T1、T2，操纵变量为冷剂流量 F1、F2。通过机理分析得到主通道为一阶加纯滞后，耦合通道为反向增益。现场做阶跃实验后，用最小二乘拟合得到：

G(0) = [ 1.20   -0.35 ]
       [ 0.28    0.95 ]

RGA = [ 1.09   -0.09 ]
      [-0.09    1.09 ]

RGA 对角元素接近 1，说明原配对基本合理，但耦合通道不可忽略。若直接投用两个独立 PID，出口温度会出现明显交互振荡；引入静态解耦器 D = G(0)^{-1} 后，主通道等效为近似对角系统，PID 整定难度显著下降。随后再用状态空间模型校验动态解耦效果，确认穿越频率附近相位裕度满足要求。

六、常见误区与复盘

误区一：只建主通道，忽略耦合通道。 耦合通道的增益和滞后往往决定解耦器是否稳定，不能省略。
误区二：用单回路阶跃数据直接拼成 MIMO 模型。 多变量实验必须考虑其他回路的保持状态，否则通道间会相互污染。
误区三：辨识数据未充分激励耦合方向。 若操纵变量变化幅度过小，耦合通道信噪比不足，辨识结果不可信。
误区四：模型校验只看拟合度，不看闭环频段。 应重点检查穿越频率附近的幅相误差和 RGA 随频率的变化。

复盘时建议按“配对结构 → 通道参数 → 频域校验 → 闭环仿真”四步走：先用 RGA 确认配对，再用阶跃或 PRBS 辨识各通道参数，接着在 Bode 图上检查关键频段误差，最后用闭环仿真验证解耦器与控制器联合投用时的稳定性。只有模型在关键频段可信，后续解耦控制才有工程意义。