H∞控制模拟器 返回
控制工程模拟器

H∞控制模拟器 — 混合灵敏度设计和H∞范数

对于过程 G_p=1/(T_p s+1) 的PI控制,可视化灵敏度函数S和互补灵敏度T、权重W_S和W_T的混合灵敏度峰值(H∞范数)。学习扰动抑制与噪声抑制的设计权衡。

参数设置
PI 比例增益 K_c
灵敏度峰值约束 M_s
衰减频率 ω_BW
rad/s
互补灵敏度衰减 ω_BT
rad/s

过程为 T_p=2s 的一阶延迟,控制器为PI(积分时间 T_i=2s)。权重假设 A_s=0.01、A_t=0.01、M_t=2.0。

暂停时,拖动滑块即可即时更新结果。

计算结果
H∞范数 ‖S‖∞
峰值频率 ω_p
‖W_S·S‖∞
鲁棒性能 ‖W_S·S‖∞<1
Bode 图与 H∞ 范数实时扫描

红=|S|, 蓝=|T|, 绿虚线=|W_S|, 紫虚线=|W_T|, 灰虚线=M_s 上限。标记=‖S‖∞ 峰值(绿=鲁棒, 黄=未达)

理论与主要公式

过程 $G_p(s)=\dfrac{1}{T_p s+1}$ 和 PI 控制器 $C(s)=K_c\left(1+\dfrac{1}{T_i s}\right)$ 组成开环 $L=CG_p$。

灵敏度函数 $S$ 和互补灵敏度 $T$(代数约束 $S+T=1$):

$$S(s)=\frac{1}{1+L(s)},\qquad T(s)=\frac{L(s)}{1+L(s)}$$

灵敏度权 $W_S$(低频大)和互补灵敏度权 $W_T$(高频大):

$$W_S(s)=\frac{s/M_s+\omega_{BW}}{s+\omega_{BW} A_s},\qquad W_T(s)=\frac{s+\omega_{BT}/M_t}{A_t\,s+\omega_{BT}}$$

混合灵敏度的H∞范数(设计评价指标):

$$\|T_{zw}\|_\infty=\sup_{\omega}\sqrt{|W_S(j\omega)S(j\omega)|^2+|W_T(j\omega)T(j\omega)|^2}$$

设计条件是 $\|W_S S\|_\infty<1$ 且 $\|W_T T\|_\infty<1$。同时满足两者时,扰动抑制和噪声抑制、模型误差容限均得到保证。

H∞控制模拟器说明

🙋
我听说过「H∞控制」这个名字,但不清楚它是什么样的设计方法。与普通的PID有什么不同呢?
🎓
简单来说,H∞控制是一种直接降低闭环系统频率增益最大值的设计方法。PID是从时间域角度看超调、稳定时间等响应指标,但H∞是从频率域确保"坏频率上的增益不会暴增"。在上面的模拟器中,试试增大K_c。你会看到红色的|S|曲线峰值(Ms)在变化。这就是鲁棒性的直接指标。
🙋
|S|和|T|在低频和高频交换的现象是什么意思呢?
🎓
这就是H∞设计的核心。S表示"扰动在输出中剩余的程度",T表示"目标值传递到输出的程度、传感器噪声通过的程度"。有一个数学约束叫S+T=1,所以两个都不能同时为零。做法是:低频下降低S来消除扰动,高频下降低T来切断噪声——这个分界点由ω_BW和ω_BT决定,这就是加权函数设计。
🙋
绿色和紫色虚线是「权重」吧。它们的上方不能越过S和T吗?
🎓
准确地说,我们要使|W_S·S|和|W_T·T|都小于1。这叫「混合灵敏度问题」。从图上看,红色的|S|不能超过绿色的|W_S|^(-1),蓝色的|T|不能超过紫色的|W_T|^(-1)。模拟器中的「||W_S·S||_∞」数值如果小于1,说明扰动抑制的要求满足了。如果超过1,说明要求没达到。
🙋
当把M_s从1.4变到5时,图形差异很大。在实际应用中怎么决定M_s呢?
🎓
实务中的标准范围是1.4~2.0。1.4比较保守,鲁棒性高但追随性能差。超过2.0时,对模型误差的容限(增益余量、相位余量)就会变小,实际设备中只要摩擦力变了一点、管路长度有偏差,就容易振荡。建议以M_s=2.0为基准,通过调整K_c来让|S|峰值恰好接触M_s线。

常见问题

H∞控制用最大奇异值范数对不确定性进行统一评价,因此对具有块对角结构的结构化不确定性会比较保守。μ-综合(结构化奇异值设计)利用块结构获得更精确的鲁棒性能指标,但需要通过D-K迭代计算,成本高且容易陷入局部最优。实际做法是先从H∞开始,必要时再考虑μ设计。
鲁棒控制是「对模型不确定性保证稳定性和性能」的总体概念。H∞控制是其中一个代表性设计方法,用传递函数H∞范数量化不确定性。其他鲁棒方法还有LQG/LTR、QFT、滑模控制、μ-综合等,H∞在频率域严格保证方面的优势最明显。
在线性时不变的稳定系统中,log|S(jω)|在全频率上的积分等于一个非负的定值(与不稳定极点和有关)。这意味着,如果在某个频带降低S(扰动抑制),必定在另一个频带增大S(扰动放大)。比喻来说就像水面,一端被压下去,另一端必然抬起来。H∞设计就是在这个约束下,权衡哪些频率需要控制。
标准的做法是:W_S在低频增益大(扰动抑制),在ω_BW处通过0 dB,高频逼近1/M_s(允许峰值)的水平渐近线。W_T则相反,高频增益大(噪声抑制),在ω_BT处通过0 dB。在ω_BW和ω_BT之间留出一定裕度会让数值求解更容易收敛。如果硬要在同一频率同时压低两个权重,会与S+T=1约束矛盾,导致无解。

实际应用

航空器与火箭飞行控制:机体参数随飞行速度、高度、载重变化很大,鲁棒性至关重要。H∞控制在F-16CAS、各类无人机姿态控制、火箭推力矢量控制等领域广泛应用,在整个飞行包线上保证稳定余裕。

硬盘寻头定位、半导体曝光机:需要同时抑制外部振动/冲击和实现高速追踪,这类机电一体化系统标准采用H∞控制。柔性结构的共振峰可以通过权重内的陷波滤波器自动处理,这是H∞的独特优势。

工业过程控制(化工厂、电力系统):多变量耦合、死时间大、模型不确定性严重的场景,H∞和混合灵敏度设计是标配。蒸馏塔组分控制、火电机组锅炉-汽轮机控制、电网稳定装置都有实例。

汽车底盘控制:主动悬架、电动助力转向、ESP等系统会因路面和乘员数变化,特性偏离很大。乘坐舒适性与路面跟踪的权衡——正好是H∞混合灵敏度的典型应用场景。

常见误解与注意事项

最常见的误解是「H∞范数越小,实机性能就自动变好」。H∞是频率域峰值增益的最小化,并不能直接最小化时间响应的稳定时间、超调。在模拟器中增大K_c来降低||W_S·S||_∞,可能导致过渡响应变振荡。H∞设计的目标是「保证鲁棒性」,时间响应品质需要单独验证和调整。

另一个常见错误是「权重函数设置得越严格越好」。如果W_S在低频过度放大,或W_T在高频过度放大,会导致问题无解(infeasible),或输出极高阶的控制器。试试在模拟器中把ω_BW改得接近ω_BT——你会看到权重曲线重叠,Ms峰值突然激增。这说明设计要求已接近物理极限。实务中要注意Bode感度积分(waterbed效应)的约束,分出需求的轻重。

最后要说明,本模拟器只是在固定PI结构下评估混合灵敏度指标的可视化工具。真正的H∞控制器综合是通过求解Riccati方程或LMI来进行的。这里展示的是「改变PI增益K_c时混合灵敏度如何变化」的直观教学演示。完整的H∞控制器设计应该用MATLAB鲁棒控制工具箱的hinfsyn函数、Python控制库,或SLICOT库。得到的控制器通常阶数高于PI(等于过程阶数加权重阶数)。

使用指南

  1. 调整灵敏度权重函数W_S的增益(slMsVal)在0.5~2.0范围内,设定低频扰动抑制强度
  2. 设置互补灵敏度权重函数W_T的增益(slBTVal)在0.1~1.0,控制高频噪声灵敏度,确保|T|max不超过1.9
  3. 同步调整PI控制增益Kc(slKcVal)和积分时间τ_i(slBWVal),使闭环带宽ω_c收敛到0.5~5.0 rad/s目标值,确认||W_S·S||_∞小于1.0

具体计算示例

对一阶延迟过程(τ=2s, K_p=1)采用Kc=1.2、τ_i=3.0s的PI控制时,灵敏度函数S(jω)峰值|S|max≈1.35,互补灵敏度函数T(jω)峰值|T|max≈1.18。设W_S=0.8(低频优先)、W_T=0.3时,||W_S·S||_∞≈0.95,闭环带宽ω_c≈1.2rad/s,满足混合灵敏度条件||W_S·S||_∞ + ||W_T·T||_∞ < 2.0

实务中的注意事项