4. 多变量系统特性分析
多变量过程控制的核心难点,并不在于回路数量多,而在于回路之间通过同一物理过程相互影响。一个操纵变量变化,往往同时改变多个被控变量;一个被控变量的偏差,也可能由多个操纵变量共同引起。这种“你中有我、我中有你”的特性就是耦合。若忽视耦合,单回路整定得再好,也可能在投运后出现振荡、反向调节甚至失稳。因此,多变量系统分析的第一步,是把耦合“看清楚”,再判断它是否严重、是否需要解耦、以及解耦后系统是否仍然可控可观测。
一、耦合现象:从物理通道到传递矩阵
以典型的双输入双输出过程为例,两个操纵变量 u₁、u₂ 分别通过四条通道影响两个被控变量 y₁、y₂。其传递关系可写为:
y1(s) = G11(s)·u1(s) + G12(s)·u2(s)
y2(s) = G21(s)·u1(s) + G22(s)·u2(s)
矩阵形式:
[ y1 ] [ G11 G12 ] [ u1 ]
[ y2 ] = [ G21 G22 ] [ u2 ]
其中 G11、G22 是主通道,G12、G21 是耦合通道。若 G12 与 G21 近似为零,系统可视为两个独立单回路;若耦合通道增益与主通道相当,甚至符号相反,则必须进行解耦设计。工程上常见的耦合来源包括:同一加热器同时影响出口温度和压力、精馏塔回流比同时影响塔顶与塔底组成、反应器进料流量同时影响液位与浓度等。
二、相对增益阵列 RGA:配对与耦合程度的量化
相对增益阵列(Relative Gain Array, RGA)是判断回路配对最常用的工具。其元素 λij 定义为:其他回路全部开环时 uj 对 yi 的增益,与其他回路全部闭环时 uj 对 yi 的增益之比。对于 2×2 系统,只需计算一个元素:
λ11 = 1 / (1 - (G12·G21)/(G11·G22))
RGA = [ λ11 1-λ11 ]
[ 1-λ11 λ11 ]
RGA 的工程含义非常明确:λij 越接近 1,说明该通道受其他回路影响越小,适合配对;λij 接近 0,说明该通道几乎被其他回路“抢占”,不宜配对;λij 为负值,说明闭环后增益符号反转,极易引起不稳定。对于 2×2 系统,通常选择 λ11 接近 1 的配对方案;若 λ11 在 0.5 附近,说明两个配对方案几乎等价,耦合严重,需考虑解耦。
| RGA 元素 λij | 含义 | 工程建议 |
|---|---|---|
| λ ≈ 1 | 该通道几乎不受其他回路影响 | 优先配对,耦合弱 |
| 0.5 < λ < 1 | 存在一定耦合,但可接受 | 可先按单回路整定,观察效果 |
| λ ≈ 0.5 | 两个配对方案等价,耦合严重 | 建议解耦或采用多变量控制器 |
| λ < 0 | 闭环增益符号反转 | 禁止直接配对,必须解耦或重构 |
三、奇异值分解与条件数:方向敏感性的度量
RGA 主要反映稳态配对关系,而奇异值分解(SVD)从“方向增益”角度揭示系统在不同输入方向上的放大能力。对传递矩阵 G 做 SVD:
G = U · Σ · V^T
Σ = diag(σ1, σ2, ..., σn),σ1 ≥ σ2 ≥ ... ≥ σn ≥ 0
条件数 κ = σ1 / σn
最大奇异值 σ1 表示系统在最敏感方向上的增益,最小奇异值 σn 表示最不敏感方向上的增益。条件数 κ 越大,说明系统在不同方向上的增益差异越大,即“方向性”越强。高条件数系统对输入方向非常敏感:某些方向的扰动会被大幅放大,而另一些方向的控制作用却非常微弱。工程上,κ 接近 1 表示系统各方向增益均衡,易于控制;κ 大于 10 时需警惕;κ 大于 100 时通常认为系统病态,常规解耦和单回路控制难以奏效。
四、可控性与可观测性:系统“能不能控、能不能看”
可控性回答的是:能否通过操纵变量在有限时间内把系统状态驱动到任意目标状态;可观测性回答的是:能否通过被控变量测量值重构出系统全部内部状态。对于线性状态空间模型:
ẋ = A·x + B·u
y = C·x + D·u
可控性矩阵: Co = [B, AB, A²B, ..., A^(n-1)B]
可观测性矩阵:Ob = [C; CA; CA²; ...; CA^(n-1)]
若 rank(Co) = n,则系统完全可控;
若 rank(Ob) = n,则系统完全可观测。
在多变量过程控制中,可控性与可观测性直接决定解耦方案是否可行。若某个被控变量对任何操纵变量都不敏感,则它不可控;若某个关键状态无法由测量值推断,则它不可观测。此时,再精巧的解耦器也无法弥补信息缺失。工程实践中,应优先保证关键被控变量可控、关键状态可观测,再讨论耦合与解耦。
五、案例与复盘:精馏塔双回路耦合分析
某精馏塔以回流量控制塔顶温度、以再沸器蒸汽量控制塔底温度。通过阶跃测试得到传递矩阵:
G = [ 1.2 -0.8 ]
[ 0.9 1.1 ]
计算 RGA:
λ11 = 1 / (1 - (-0.8×0.9)/(1.2×1.1))
= 1 / (1 + 0.72/1.32)
≈ 0.65
RGA = [ 0.65 0.35 ]
[ 0.35 0.65 ]
λ11 约为 0.65,说明主配对方案(回流控顶温、蒸汽控底温)可行,但耦合不可忽略。进一步做 SVD 分析,得到 σ1 = 1.85、σ2 = 0.42,条件数 κ ≈ 4.4,属于中等偏大。复盘结论是:该塔不宜直接投用两个独立 PID,应先做静态解耦,再整定控制器;同时应检查塔顶温度与塔底温度的测量滞后是否匹配,避免耦合通道动态差异引发振荡。
六、常见误区
- 只看稳态增益,忽略动态差异。耦合通道响应比主通道快时,静态 RGA 可能显示耦合不严重,但实际过渡过程振荡明显。
- 把 RGA 当作唯一判据。RGA 只反映稳态配对,不能替代 SVD 对方向增益和条件数的评估。
- 条件数大就强行解耦。若系统本身不可控或不可观测,解耦器会放大噪声和模型误差,反而恶化控制品质。
- 忽略变量重构。当 RGA 出现负值或条件数极高时,重新选择操纵变量与被控变量,往往比设计复杂解耦器更有效。
综上,多变量系统特性分析不是单一指标的判断,而是一条从耦合识别、RGA 配对、SVD 方向性评估到可控可观测性验证的完整链条。只有把这条链条走通,后续的解耦设计才有可靠依据。