奈奎斯特判据的原理:深度解析与稳定性分析指南

奈奎斯特判据的原理:什么是稳定性分析的“黄金标准”?

在自动控制理论中,奈奎斯特判据(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) 是开环传递函数。

  1. 零点关系:F(s) 的零点即为 1 + G(s)H(s) = 0 的根,也就是闭环特征方程的根。因此,F(s) 的零点数 Z 就是闭环系统在右半平面的极点数。
  2. 极点关系:F(s) 的极点即为 G(s)H(s) 的极点。因此,F(s) 的极点数 P 就是开环系统在右半平面的极点数。
  3. 环绕关系:F(s) 环绕原点,等价于 G(s)H(s) 环绕 (-1, j0) 点。

因此,奈奎斯特判据的原理最终简化为公式:

Z = P - N

其中 N 是奈奎斯特曲线 G(jω)H(jω) 逆时针包围 (-1, j0) 点的圈数。若 Z = 0,则闭环系统稳定。

应用步骤:如何绘制与分析奈奎斯特图

掌握奈奎斯特判据的原理后,实际操作分为以下几个关键步骤。为了更清晰地展示,我们使用选项卡来区分不同阶次系统的绘制技巧。

确定开环右极点 P

首先,写出系统的开环传递函数 G(s)H(s)。检查分母多项式,找出所有极点。统计位于右半 s 平面(实部大于0)的极点个数,记为 P。

注意:大多数物理系统是最小相位系统,此时 P=0。但若系统内部存在不稳定环节(如倒立摆的未反馈部分),P 可能大于 0。

绘制奈奎斯特路径与渐近线

奈奎斯特路径是一个包围整个右半 s 平面的半圆。对应的映射曲线由三部分组成:

  1. jω 轴:ω 从 0 变化到 +∞。这是主要的绘图部分。
  2. 无穷大半圆:s = Re^{jθ},R→∞。通常映射到原点或实轴上的某点。
  3. -jω 轴:ω 从 -∞ 变化到 0。它是 jω 轴的镜像(共轭复数)。

对于 I 型系统(含一个积分环节),在 ω=0 处需要补画一个半径无穷大的半圆,顺时针从 -90° 连接到 +90° 方向,以绕过原点处的极点。

计算关键点坐标

为了准确绘图,需计算以下特殊点:

  • 起点 (ω=0):计算 |G(j0)H(j0)| 和 ∠G(j0)H(j0)。
  • 终点 (ω→∞):通常趋向于原点,相角为 -(n-m)×90°。
  • 穿越点:计算曲线与实轴和虚轴的交点。令虚部为0求实轴交点,令实部为0求虚轴交点。

应用 Z = P - 2N 判据

重要修正:由于我们通常只画 ω: 0 → +∞ 的部分,而完整的奈奎斯特图是 ω: -∞ → +∞。因此,实际的包围圈数 N_full 与半平面曲线包围 (-1, j0) 的圈数 N 之间的关系为 N_full = 2N(假设曲线不经过 (-1, j0))。

稳定性判据修正为:

Z = P - 2N

若 Z=0,系统稳定。这意味着:

  • 当 P=0 时,奈奎斯特曲线不能包围 (-1, j0) 点(即 N=0)。
  • 当 P>0 时,曲线必须顺时针包围 (-1, j0) 点 P/2 圈(即 N = -P/2,注意符号定义)。

实例分析:典型系统的奈奎斯特图

通过具体案例来理解奈奎斯特判据的原理是最有效的方法。以下列举三种典型系统。

案例 1:0 型系统(无积分环节)

设开环传递函数 G(s)H(s) = K / ((s+1)(s+2))。

案例 2:I 型系统(含一个积分环节)

设 G(s)H(s) = K / (s(s+1))。

案例 3:III 型系统(含两个积分环节)

设 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° -∞ 需严格计算穿越点

常见问题解答 (FAQ)

Q: 奈奎斯特判据适用于非线性系统吗?

A: 标准的奈奎斯特判据仅适用于线性时不变(LTI)系统。对于非线性系统,可以使用描述函数法(Describing Function Method),它在数学形式上类似于奈奎斯特判据,但引入了非线性环节的平均特性近似,属于频域分析的非线性扩展。

Q: 如果奈奎斯特曲线恰好经过 (-1, j0) 点,系统状态如何?

A: 这表示系统处于临界稳定状态。闭环特征方程有根位于虚轴上,系统输出将产生等幅振荡。在实际工程中,这被视为不稳定,因为任何参数微小变化都会导致发散。

Q: 为什么奈奎斯特判据需要补画无穷大半圆?

A: 这是因为奈奎斯特路径必须包围整个右半 s 平面。当开环传递函数在虚轴上有极点(如积分环节 1/s)时,路径不能直接穿过极点,必须用一个半径无穷小(或无穷大,取决于映射方向)的半圆绕过。这个绕过动作在 G(s)H(s) 的映像平面上表现为一个半径无穷大的圆弧,用于连接 ω=0⁻ 和 ω=0⁺ 的曲线段,保证映射曲线的封闭性。

Q: 如何判断非最小相位系统的稳定性?

A: 非最小相位系统含有右半平面零点或延迟环节。其奈奎斯特图的形状会与最小相位系统不同(通常相角滞后更大)。判定步骤不变:确定 P 值(注意右半平面极点),绘制完整曲线,计算包围圈 N,应用 Z=P-2N。非最小相位系统往往更难稳定,因为相角裕度通常较小。

总结

本文深入探讨了奈奎斯特判据的原理,从复变函数的数学基础到工程应用的具体步骤,再到常见误区与周边热点。奈奎斯特判据不仅是控制理论中的一个定理,更是工程师手中一把衡量系统稳定性的利器。通过理解其背后的映射思想,我们可以更灵活地应对各种复杂的控制系统设计挑战。

掌握奈奎斯特判据,意味着您已经从代数求解的局限中解脱出来,进入了频域分析的广阔天地。建议结合 MATLAB 仿真软件,动手绘制不同参数下的奈奎斯特图,以加深直观理解。

✨ 更多精彩
🏷️ 标签:奈奎斯特判据频域分析法开环传递函数控制工程根轨迹
墨卷研习录
CopyRight@educfa.com.cn