本页目录

经典 V · 频域分析与裕度

层次：本科核心｜ 这是工业界最常用的分析语言。频域方法之所以在工程上占主导，有一个极实际的理由：它可以直接从实测数据出发——给系统扫频，测出幅值与相位，就能设计控制器，不需要解析模型。本页讲 Bode 图、Nyquist 判据、稳定裕度，并揭示控制论最深刻的一条限制：水床效应。

学习层：给拾放轴做一次带时滞的开环设计审计

1. 具体工程情境：速度、相位余量和通信延迟谁先吃掉预算？

用一个稳定、可直接计算的对象和超前控制器族：

\[ P(s)=\frac{1}{(1+s/1.2)(1+s/5)}, \qquad C(s)=K\frac{1+s/z}{1+s/p}, \qquad 0<z<p, \]

再把纯延迟 \(e^{-\tau s}\) 放进回路，因此

\[ L(s)=C(s)P(s)e^{-\tau s}, \qquad 1+L(s)=0 \]

是标准负反馈的闭环特征方程。\(P\) 和 \(C\) 的有限极点都在左半平面；延迟不改幅值，却按 \(-\omega\tau\) 弧度侵蚀相位。

预测门，必须先作答：

固定 \(K,z,p\)，增加纯延迟会改变增益交叉频率，还是首先改变交叉处相位？
若 \(|L(j\omega)|\) 从未达到 1，按本页的定义相位裕度应写成多少？
“Nyquist 曲线不包围 \(-1\)”能否脱离开环右半平面极点数 \(P\) 单独使用？

2. 正式公式桥：从 Bode 账本读出交叉与裕度

在相位从低频连续展开、采用负反馈约定时，定义

\[ \omega_{gc}:|L(j\omega_{gc})|=1, \qquad \mathrm{PM}=180^\circ+\angle L(j\omega_{gc}), \]

以及相位交叉

\[ \omega_{pc}:\angle L(j\omega_{pc})=-180^\circ, \qquad \mathrm{GM}=\frac{1}{|L(j\omega_{pc})|}, \qquad \mathrm{GM}_{dB}=-20\log_{10}|L(j\omega_{pc})|. \]

如果扫描区间没有相应交叉，margin 就是未定义，不是 0。若 \(\tau=0\) 且相位只在无穷远渐近 \(-180^\circ\)，有限频率的 GM 也可未定义。对纯延迟，\(|e^{-j\omega\tau}|=1\)，所以幅值曲线不动，交叉处相位变为 \(\angle C(j\omega)+\angle P(j\omega)-\omega\tau\)。

3. 动手实验：Bode、Nyquist 与逐频率设计账

拖动 \(K\)、超前零点 \(z\)、超前极点 \(p\) 和延迟 \(\tau\)。图上同时显示幅值 dB、连续展开相位和 Nyquist 正频支；竖线标出 \(\omega_{gc}\)，水平线标出 \(-180^\circ\)，\(-1\) 点旁给出当前曲线距离。右侧 ledger 明列 \(|P|,|C|,|L|\)、各段相位、延迟相位与总相位，避免把“控制器相位裕度”误读成“闭环极点保证”。

JavaScript 失效时的静态 fallback：实验族固定为 \(P(s)=1/[(1+s/1.2)(1+s/5)]\)、\(C(s)=K(1+s/z)/(1+s/p)\)、\(0<z<p\)，默认 \(K=1,z=0.5,p=4,\tau=0.08\ \mathrm s\)。读账顺序是：

步骤	计算	应检查的边界
1	\(\lvert L\rvert=\lvert P\rvert\lvert C\rvert\)	延迟不改变这一列
2	\(\phi=\tan^{-1}(\omega/z)-\tan^{-1}(\omega/p)-\tan^{-1}(\omega/1.2)-\tan^{-1}(\omega/5)-\omega\tau\)	延迟项单位是弧度，显示时转成度
3	找 \(\lvert L\rvert=1\) 的最低正频率	没有交叉则 PM 未定义
4	找连续相位 \(-180^\circ\) 的最低交叉	没有有限交叉则 GM 未定义
5	代回 PM、GM，并看 Nyquist 到 \(-1\) 的距离	裕度只抽查交叉点，不替代完整稳定性分析

这里的 Nyquist 结论遵守 \(Z=N+P\)：只有在开环极点数、绕行方向和虚轴边界按通常假设处理清楚时，才可把“不包围 \(-1\)”作为开环稳定情形的简化判据。

4. 误区与模型边界

PM/GM 是开环频率量：它们是闭环设计的诊断指标，不是把开环曲线直接当成闭环响应。
没有交叉不是零裕度：它表示当前扫描定义下该裕度未定义；要说明频率范围、相位展开和交叉选择。
延迟不改幅值不等于没有代价：它会在高频持续吃相位，可能让 Nyquist 曲线靠近或穿过 \(-1\)。
Nyquist 简化判据有前提：开环有右半平面极点时必须计入 \(P\)；虚轴极点、重复交叉和非最小相位结构也需要单独处理。
稳定对象族只是实验边界：真实设备的采样保持、饱和、未建模共振和时变延迟可能超出这里的解析族。

迁移问题：远程视觉伺服的延迟从 \(20\) ms 漂到 \(80\) ms。你会优先降低 \(K\)、降低交叉频率，还是调整超前零极点？请同时说明哪一项改变幅值交叉，哪一项改变交叉处相位，并给出一个你会重新测量的闭环指标。

一、频率响应

对 LTI 系统输入正弦 \(\sin(\omega t)\)，稳态输出仍是同频正弦，只是幅值被放大 \(|G(j\omega)|\) 倍、相位滞后 \(\angle G(j\omega)\)。把这两个量对频率作图，就是 Bode 图（幅值用 dB、频率用对数）。

为什么用对数：串联系统的传递函数相乘，取对数后变成相加——于是复杂系统的 Bode 图可以由各环节叠加得到，手工作图成为可能。这是 Bode 在 1930 年代的关键设计。

基本环节的渐近线（务必熟悉）：

环节	幅值斜率	相位
积分 \(1/s\)	−20 dB/dec	−90°
一阶惯性 \(\frac{1}{\tau s+1}\)	转折后 −20 dB/dec	0° → −90°
二阶欠阻尼	转折后 −40 dB/dec，谐振峰	0° → −180°
时滞 \(e^{-\tau s}\)	0 dB（不变）	\(-\omega\tau\)，无限增长

注意最后一行：时滞不衰减幅值，却无限累积相位——这就是第一页说它是"反馈头号敌人"的定量表述。

二、Nyquist 判据

把 \(G(j\omega)\) 在复平面上画成一条曲线（\(\omega\) 从 \(-\infty\) 到 \(+\infty\)），在开环极点没有落在积分路径上、绕行方向和相位展开按通常 Nyquist 约定处理时，Nyquist 判据给出闭环稳定的条件：

\[Z = N + P\]

\(Z\) 是闭环右半平面极点数（要求为 0），\(P\) 是开环右半平面极点数，\(N\) 是曲线顺时针包围 \(-1\) 点的次数。这里的计数依赖上述开环极点与绕行假设，不能把公式脱离约定机械套用。

只有对开环右半平面极点数 \(P=0\)、且没有虚轴极点等特殊边界的通常情形，条件才简化为：Nyquist 曲线不包围 \(-1\) 点。若 \(P\ne0\)，必须按 \(Z=N+P\) 计数；不包围 \(-1\) 本身不能证明闭环稳定。

为什么是 \(-1\)：闭环特征方程是 \(1+L(s)=0\)，即 \(L(s) = -1\)。\(-1\) 点就是"回路增益为 1 且相位为 180°"的那个临界状态——正是第一页说的"负反馈变成正反馈且增益足够"的条件。几何直观与物理直觉在此完全吻合。

三、稳定裕度：离灾难有多远

"不包围 \(-1\)"是二元的，工程需要定量的余量：

增益裕度 GM：相位为 −180° 时，增益还能放大多少倍才到 1。典型要求 > 6 dB（2 倍）；
相位裕度 PM：增益为 1（0 dB）时，相位离 −180° 还差多少。典型要求 45°–60°。

Bode 图与增益/相位裕度 — 图 5-1 开环 Bode 图与稳定裕度。**穿越频率** \(\omega_c\)（幅值过 0 dB 处）大致决定闭环带宽与响应速度；该处相位距 −180° 的余量即**相位裕度**。相位过 −180° 处的幅值余量为**增益裕度**。工程经验：**PM ≈ 60° 对应 \(\zeta\approx0.6\)、超调约 10%**——频域指标与第三页的时域指标由此直接对应。

极其有用的经验关系：

\[\zeta \approx \frac{\mathrm{PM}}{100}\quad(\mathrm{PM} < 70°), \qquad \omega_c \approx \omega_n\]

这条经验把频域设计与时域指标连了起来：客户要"超调小于 10%"，你就知道要设计 PM ≈ 60°；客户要"调节时间 0.5 秒"，你就知道穿越频率大致该放在哪里。这是经典控制最实用的一条转换。

为什么不能只看裕度：GM 与 PM 都只看曲线上的一个点。曲线可能在别处很靠近 \(-1\)（虽然两个裕度都合格）。更严格的指标是向量裕度（曲线到 \(-1\) 的最短距离），它等价于 \(1/\|S\|_\infty\)——这就通向了第十六页的 \(H_\infty\) 理论。

四、回路整形：频域设计的思路

回路整形（loop shaping）是频域设计的核心思想：按频段给回路增益 \(|L(j\omega)|\) 规定形状。

\[\text{低频：}|L|\ \text{大} \quad\Longrightarrow\quad |S| = \frac{1}{|1+L|}\ \text{小} \quad\Longrightarrow\quad \text{跟踪准、抗扰强}\]

\[\text{高频：}|L|\ \text{小} \quad\Longrightarrow\quad |T|\ \text{小} \quad\Longrightarrow\quad \text{抑制噪声、避开未建模动态}\]

\[\text{穿越处：斜率约 -20 dB/dec} \quad\Longrightarrow\quad \text{足够的相位裕度}\]

这是一个清晰、可操作的设计配方：低频抬高（加积分/滞后）、高频压低（加滤波）、穿越区调整斜率与相位（加超前）。工业上大量控制器就是这样设计出来的。

五、水床效应：一条不可逾越的定理

现在讲控制论最深刻、也最少被外行知道的结果。

Bode 灵敏度积分定理：对开环稳定、相对阶 ≥ 2 的系统，

\[\int_0^\infty \ln|S(j\omega)|\,d\omega = 0\]

含义：\(\ln|S|\) 曲线下的正负面积必须相等。你在某个频段把 \(|S|\) 压到 1 以下（抑制干扰），就必然在另一个频段把它抬到 1 以上（放大干扰）。

这就是"水床效应"——按下这头，那头必然翘起来。

它不是设计水平的问题，是数学定理。 工程含义极其实在：

你不可能设计出在所有频段都抑制干扰的控制器；
只能选择把干扰放大放到你不在乎的频段（通常是高频，或没有干扰能量的频段）；
若系统有右半平面极点或零点，情况更糟：积分等式右边变成正数，意味着"必须付出的代价"更大——这是第二页说的"非最小相位系统性能有天花板"的严格版本。

这个定理与其他学科的"不可能性结果"同属一类：🔗 微电子站的 60 mV/dec 玻尔兹曼极限、信息论的香农极限、热力学第二定律。成熟的工程学科都有这样一条告诉你"什么做不到"的定律，而知道边界在哪，正是专业性的体现。

六、从实测数据出发

最后强调频域方法在工业上的决定性优势：

扫频实验（或用伪随机信号、阶跃辨识）直接得到 \(|G(j\omega)|\) 与 \(\angle G(j\omega)\)，无需知道物理模型；
在实测 Bode 图上直接做回路整形；
时滞、柔性模态和选定工作点/激励幅值下的局部小信号效应会反映在实测频响里；强非线性、滞回、饱和及幅值依赖不会被一条 FRF 自动概括，仍需分幅值或时域验证。

这就是为什么现场工程师更常用频域方法：它容忍"我不知道这台设备的精确模型"这个常态。（🔗 第十八页系统辨识会给出更系统的做法。）

七、要点

Bode 图的对数坐标使串联变叠加，是手工设计的基础；
Nyquist 的 \(-1\) 点就是"增益 1、相位 180°"的临界条件；
在低阶、二阶主导且没有强零点等附加动态时，PM ≈ 60°、\(\zeta\) ≈ 0.6 与约 10% 超调可作经验换算；
回路整形给出清晰配方：低频高、高频低、穿越处 −20 dB/dec；
水床效应是定理：干扰抑制的总量守恒，你只能选择把代价放在哪个频段；
频域方法可直接用实测数据，这是它在工业上占主导的根本原因。

下一页：PID。理论讲完，来看工业界真正在用的东西——它占据了实际控制器的绝大多数，而它的工程细节（抗饱和、微分滤波、整定）恰恰是教科书讲得最少的。