4. H∞范数与系统增益：H∞范数的定义、物理意义、计算方法

好，咱们今天聊聊H∞范数。说实话，我刚接触鲁棒控制那会儿，看到这个符号就头大——一个无穷符号，加个H，到底想表达什么？后来做项目做多了，才慢慢品出它的味道。

H∞范数，说白了就是衡量一个系统“放大能力”的指标。你给它一个输入，它给你一个输出，这个输出能有多大？H∞范数就是回答这个问题的。

4.1 数学定义：先别怕，其实很简单

先看定义。对于一个稳定的线性时不变系统，传递函数为G(s)，它的H∞范数定义为：

||G(s)||∞ = sup_ω σ_max[G(jω)]

其中σ_max表示最大奇异值，sup表示上确界，也就是在所有频率ω上取最大值。

嗯，这里要注意：对于SISO系统，奇异值就是幅值，所以定义简化为：

||G(s)||∞ = sup_ω |G(jω)|

说白了，就是Bode图上幅频曲线的最高点。你想想看，一个系统的频率响应，在不同频率下增益不一样，H∞范数就是那个最坏情况下的增益。

4.2 物理意义：为什么它这么重要？

我在做伺服系统项目时，遇到过一个问题：电机在某个频率下会产生共振，导致定位精度下降。当时我就在想，怎么量化这个“最坏情况”？后来发现，H∞范数就是干这个的。

它的物理意义可以从三个角度理解：

• 能量增益角度： 输入信号的能量经过系统后，输出能量与输入能量的最大比值。说白了就是系统对能量的放大倍数。
• 最坏情况分析： 系统对哪个频率最敏感？H∞范数告诉你答案。它代表系统在最不利频率下的响应强度。
• 鲁棒稳定性： 在鲁棒控制中，H∞范数直接关联到系统对不确定性的容忍度。小范数意味着系统更“抗造”。

4.3 计算方法：手算与工具

计算H∞范数，有几种路子。我按从简单到复杂给你捋一捋。

4.3.1 从Bode图直接读

对于SISO系统，最直观的方法就是画Bode图。幅频曲线的最高点对应的dB值，就是H∞范数。

举个例子：G(s) = 10/(s+1)
Bode图上，低频增益20dB，高频滚降。
H∞范数 = 20dB = 10（线性值）

这个方法我经常用，尤其是在现场调试时，用示波器扫个频率响应，一眼就能看出问题。

4.3.2 状态空间法

对于MIMO系统，或者需要精确计算时，用状态空间法。我记得有一次做飞行器控制，系统有6个输入4个输出，Bode图根本看不清楚，只能用这个方法。

基本思路是解一个Riccati方程：

给定状态空间模型 (A, B, C, D)
H∞范数 < γ 的充要条件是：
存在对称正定矩阵X，满足：
A^T X + X A + C^T C + (X B + C^T D)(γ^2 I - D^T D)^{-1}(B^T X + D^T C) < 0

嗯，这个公式看着吓人，但实际用MATLAB算起来很快。

4.3.3 数值迭代法

工程上最常用的是二分法加γ-迭代。我给大家看一段MATLAB代码，这是我做项目时常用的套路：

% 计算H∞范数的二分法
sys = ss(A, B, C, D);
gamma_low = 0;
gamma_high = 1e6;

for iter = 1:50
    gamma = (gamma_low + gamma_high) / 2;
    % 检查是否满足H∞条件
    if hinfnorm_check(sys, gamma)
        gamma_high = gamma;
    else
        gamma_low = gamma;
    end
end

% 最终结果
hinf_norm = gamma_high;

4.4 知识体系：一张图看懂

下面这张图是我自己整理的，把H∞范数的核心逻辑串起来了：

4.5 实际应用中的注意事项

最后，我分享几个实战经验：

1. 别只看数值： H∞范数小不代表系统好。有时候为了压低某个频率的峰值，会牺牲其他性能。我见过有人把范数压到1以下，结果系统响应慢得像蜗牛。
2. 注意频率范围： 实际系统只在有限频率范围内有意义。高频段模型通常不准，算出来的范数可能虚高。
3. MIMO系统要小心： 奇异值分解的结果要结合物理意义看。哪个通道的增益大？哪个方向最敏感？这些信息比单纯一个范数数值更有价值。

---