在自动控制理论中,奈奎斯特判据(Nyquist Stability Criterion)是频域分析法的核心基石。它由哈里·奈奎斯特(Harry Nyquist)于1932年提出,是一种 graphical method(图形法),用于判断线性时不变(LTI)闭环系统的稳定性。与劳斯判据(Routh Criterion)仅适用于代数方程不同,奈奎斯特判据的原理能够处理包含时间延迟(Transport Lag)的系统,这是其在工程实践中极具价值的原因。
简单来说,奈奎斯特判据的原理并不直接求解闭环特征方程的根,而是通过观察开环传递函数的频率响应曲线(即奈奎斯特图)在复平面上的行为,利用复变函数论中的辐角原理,间接推断出闭环系统在右半s平面的极点数量。
无需解高次方程。对于三阶及以上系统,直接求解特征根非常困难,而奈奎斯特判据通过图形即可直观判断。
不仅适用于最小相位系统,也适用于非最小相位系统,特别是含有纯滞后环节 e^{-τs} 的系统。
除了判断绝对稳定性,还能通过穿越点的位置直观看出系统的相位裕度和幅值裕度。
要深刻理解奈奎斯特判据的原理,必须回到复变函数的基础——辐角原理。假设有一个复变函数 F(s),它在 s 平面上有一条封闭曲线 D_s。如果 F(s) 在 D_s 内部有 Z 个零点和 P 个极点,且没有极点在曲线上,那么当 s 沿 D_s 顺时针绕行一周时,F(s) 的映像 F(D_s) 将在复平面上逆时针绕行原点 N 周,且满足以下关系:
N = Z - P
其中:
在控制系统中,我们定义辅助函数 F(s) = 1 + G(s)H(s)。其中 G(s)H(s) 是开环传递函数。
因此,奈奎斯特判据的原理最终简化为公式:
Z = P - N
其中 N 是奈奎斯特曲线 G(jω)H(jω) 逆时针包围 (-1, j0) 点的圈数。若 Z = 0,则闭环系统稳定。
掌握奈奎斯特判据的原理后,实际操作分为以下几个关键步骤。为了更清晰地展示,我们使用选项卡来区分不同阶次系统的绘制技巧。
首先,写出系统的开环传递函数 G(s)H(s)。检查分母多项式,找出所有极点。统计位于右半 s 平面(实部大于0)的极点个数,记为 P。
注意:大多数物理系统是最小相位系统,此时 P=0。但若系统内部存在不稳定环节(如倒立摆的未反馈部分),P 可能大于 0。
奈奎斯特路径是一个包围整个右半 s 平面的半圆。对应的映射曲线由三部分组成:
对于 I 型系统(含一个积分环节),在 ω=0 处需要补画一个半径无穷大的半圆,顺时针从 -90° 连接到 +90° 方向,以绕过原点处的极点。
为了准确绘图,需计算以下特殊点:
重要修正:由于我们通常只画 ω: 0 → +∞ 的部分,而完整的奈奎斯特图是 ω: -∞ → +∞。因此,实际的包围圈数 N_full 与半平面曲线包围 (-1, j0) 的圈数 N 之间的关系为 N_full = 2N(假设曲线不经过 (-1, j0))。
稳定性判据修正为:
Z = P - 2N
若 Z=0,系统稳定。这意味着:
通过具体案例来理解奈奎斯特判据的原理是最有效的方法。以下列举三种典型系统。
设开环传递函数 G(s)H(s) = K / ((s+1)(s+2))。
设 G(s)H(s) = K / (s(s+1))。
设 G(s)H(s) = K / (s²(s+1))。
| 系统类型 | 开环极点 P | 起点相角 (ω→0) | 终点相角 (ω→∞) | 稳定性条件 (P=0时) |
|---|---|---|---|---|
| 0 型 (二阶) | 0 | 0° | -180° | K < 临界增益 |
| I 型 (二阶) | 0 | -90° | -180° | 恒稳定 |
| II 型 (三阶) | 0 | -180° | -270° | 通常不稳定 |
| 含延迟系统 | 0 | 0° | -∞ | 需严格计算穿越点 |
A: 标准的奈奎斯特判据仅适用于线性时不变(LTI)系统。对于非线性系统,可以使用描述函数法(Describing Function Method),它在数学形式上类似于奈奎斯特判据,但引入了非线性环节的平均特性近似,属于频域分析的非线性扩展。
A: 这表示系统处于临界稳定状态。闭环特征方程有根位于虚轴上,系统输出将产生等幅振荡。在实际工程中,这被视为不稳定,因为任何参数微小变化都会导致发散。
A: 这是因为奈奎斯特路径必须包围整个右半 s 平面。当开环传递函数在虚轴上有极点(如积分环节 1/s)时,路径不能直接穿过极点,必须用一个半径无穷小(或无穷大,取决于映射方向)的半圆绕过。这个绕过动作在 G(s)H(s) 的映像平面上表现为一个半径无穷大的圆弧,用于连接 ω=0⁻ 和 ω=0⁺ 的曲线段,保证映射曲线的封闭性。
A: 非最小相位系统含有右半平面零点或延迟环节。其奈奎斯特图的形状会与最小相位系统不同(通常相角滞后更大)。判定步骤不变:确定 P 值(注意右半平面极点),绘制完整曲线,计算包围圈 N,应用 Z=P-2N。非最小相位系统往往更难稳定,因为相角裕度通常较小。
本文深入探讨了奈奎斯特判据的原理,从复变函数的数学基础到工程应用的具体步骤,再到常见误区与周边热点。奈奎斯特判据不仅是控制理论中的一个定理,更是工程师手中一把衡量系统稳定性的利器。通过理解其背后的映射思想,我们可以更灵活地应对各种复杂的控制系统设计挑战。
掌握奈奎斯特判据,意味着您已经从代数求解的局限中解脱出来,进入了频域分析的广阔天地。建议结合 MATLAB 仿真软件,动手绘制不同参数下的奈奎斯特图,以加深直观理解。