奈奎斯特图 返回
分析工具

奈奎斯特图·稳定裕度模拟器

选择开环传递函数并调整增益,即可实时绘制奈奎斯特图和伯德图。自动计算增益裕度、相位裕度和交叉频率。

开环传递函数
传递函数形式
增益 K
极点 a
极点 b
零点 z
稳定裕度概览

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

计算结果
增益裕度 (dB)
相位裕度 (°)
ωgc (rad/s)
0
N(环绕−1次数)
稳定性 Z=N+P
奈奎斯特图
注:红色×表示(-1, 0)临界点。蓝色轨迹为ω>0部分,虚线为ω<0部分(镜像)。轨迹反时针包围(-1,0)点时闭环系统不稳定。
理论·主要公式

奈奎斯特稳定判别法则

闭环不稳定极数 $Z = P - N$

$N$: 奈奎斯特轨迹 $L(j\omega)$($\omega:-\infty\to+\infty$)包围(-1,0)点的净圈数(逆时针为正)。对开环稳定($P=0$)系统,顺时针净包围(-1,0)($N<0$)则闭环不稳定。

$P$: 开环传递函数右半平面(RHP)的极数

$Z = 0$ → 闭环系统稳定

增益裕度: $GM = 20\log_{10}\frac{1}{|G(j\omega_{pc})|}$ (dB)

相位裕度: $PM = 180° + \angle G(j\omega_{gc})$

奈奎斯特图与稳定裕度概述

🙋
奈奎斯特图就是在复平面上画一条曲线对吧?怎么用它来判断"稳定"和"不稳定"呢?
🎓
简单来说,这条曲线是否绕过"临界点"(-1, 0)就决定了系统的稳定性。这个模拟器中,你可以拖动上面的增益K滑块,看到轨迹逐渐靠近(-1, 0)点,最后可能包围它。一旦包围就说明系统不稳定。这就是奈奎斯特稳定判别法的核心。
🙋
明白了。那为什么还要显示"增益裕度"和"相位裕度"呢?这跟奈奎斯特图有什么关系?
🎓
很好的问题。即使轨迹没有包围(-1, 0)点,如果离得太近,系统也会很脆弱。只要有一点干扰或者参数变化,就可能掉进不稳定区间。所以我们定义了"增益裕度"和"相位裕度"来量化这个安全距离。增益裕度就是说还能增大多少倍增益而不失稳定;相位裕度就是相位还能再滞后多少度。通过调节参数a、b、z,你可以看到这些余裕值怎么变化,体验设计稳定系统的过程。
🙋
那右边同时显示的伯德图是什么意思?它和奈奎斯特图是什么关系?
🎓
两个图其实是同一个系统的两种表示方式。伯德图显示频率响应的增益和相位随频率的变化;奈奎斯特图把同样的信息表示为复平面上的一条曲线。在伯德图上,增益达到0dB的频率(增益交叉频率)时的相位值,就对应了相位裕度。而相位达到-180°的频率(相位交叉频率)时的增益,就对应了增益裕度。两个图可以互相验证和补充理解,我建议你边调参数边观察两个图的变化。

物理模型与主要公式

本模拟器处理的标准开环传递函数 $L(s)$ 如下。屏幕上的滑块可以改变参数K、a、b和z。

$$L(s) = G(s)H(s) = K \frac{(s + z)}{(s + a)(s + b)}$$

$K$:增益常数,$z$:零点值,$a, b$:极点值。将 $s = j\omega$ (虚数单位×角频率)代入,得到频率响应 $L(j\omega)$,由此可以绘制奈奎斯特图。

奈奎斯特稳定判别法则:仅从开环系统信息即可判断闭环系统稳定性的强大定理。

$$Z = N + P$$

$P$:开环传递函数 $L(s)$ 右半平面(实部为正)的极数。
$N$:奈奎斯特轨迹 $L(j\omega)$ ($\omega: -\infty \to +\infty$) 反时针包围点 $(-1, 0)$ 的圈数。
$Z$:闭环系统不稳定极的数目。当 $Z=0$ 时,闭环系统稳定。

常见问题

增益裕度是相位为-180°时增益[dB]的绝对值;相位裕度是增益为0dB时相位与-180°的差值。在伯德图上找到对应频率并读取数值,或者参考屏幕右侧的自动计算结果。正值表示稳定,负值表示不稳定。
松开滑块或在数值输入框按Enter键后会实时更新。检查JavaScript是否启用,必要时刷新页面或清除缓存后重试。
减小增益K可以增大增益裕度,提高稳定性。增加零点z(向左半平面)会产生相位超前效果,改善相位裕度。减小极点a、b会降低系统速度但增加稳定性。调整参数使奈奎斯特轨迹在(-1,0)点的左侧。
P是开环传递函数L(s)右半平面(实部>0)的极数。本工具的传递函数形式为 K(s+z)/((s+a)(s+b)),极点位于 s=−a、−b;a、b>0 时极点都在左半平面,P=0。在滑块范围(a,b≥0.1)内恒有 P=0。

实际应用

飞行器与无人机姿态控制:飞行中的机体不断受风等外部扰动。通过奈奎斯特图确保充足的稳定裕度(增益裕度、相位裕度),能使控制系统在外扰下也保持稳定,避免危险的振荡(荷兰滚现象)。

汽车定速巡航系统:在坡度变化下保持设定速度需要动力系统和制动器的精确控制。设计时通过奈奎斯特分析,使系统既能快速响应又不会超调,平衡了跟踪性能和稳定性。

机器人臂位置伺服:生产线上的精密组装机器人需要以设定位置停止,且不能振荡。相位裕度(PM)与过渡过程的振荡特性直接关联,是关键的设计指标。

电力系统稳定度评估:大型输电网中发电机的同步关系影响系统整体稳定性。奈奎斯特判别法在这样复杂的高维系统中仍然是评估稳定性的古典且有效的工具。

常见误区与注意事项

首先注意,"奈奎斯特轨迹不包围(-1,0)就一定稳定"这个说法有前提条件。这个结论仅在开环系统本身稳定(P=0)的情况下才成立。对于倒立摆这样原本就不稳定的对象(P>0),情况会更复杂。注意:本工具的极点 a、b 只能取正值(极点 s=−a、−b 恒在左半平面,P=0),因此该前提在本工具中始终成立。

其次,不能只看裕度数字而忽视曲线形状。即使增益裕度显示10dB、相位裕度45°,看起来充足,但如果轨迹在(-1,0)附近急速转向或尖锐弯曲,系统对模型误差和扰动的抵抗力仍然很弱。养成既看数字也看图形,观察轨迹距离(-1,0)有多"舒缓"的习惯。

最后,要认识到这个工具的模型与现实的差距。这里用的都是简单的2阶系统。实际控制对象往往有时间延迟(死时间),更高阶的振动模态,以及非线性特性。加入死时间 $e^{-Ls}$ 后相位会进一步滞后,奈奎斯特轨迹会像螺旋一样反复绕圈。要把这里学到的原理迁移到更现实、更复杂的模型中去。

使用指南

  1. 选择传递函数形式,设置增益K、极点a、b和零点z
  2. 用增益调整滑块将K在0.1~20范围内改变,观察奈奎斯特图上轨迹与(-1, 0)点的接近程度
  3. 增益裕度GM(dB)、相位裕度PM(°)、增益交叉频率ωgc(rad/s)实时自动输出,判断稳定性

具体计算示例

选择 K/(s(s+a)) 形式,取 K=2、a=1,即开环传递函数 G(s)=2/(s(s+1)):相位从-90°渐近趋向-180°而不会穿越,因此增益裕度 GM=∞;增益交叉频率 ωgc≈1.25rad/s 处的相位裕度 PM≈38.7°,闭环系统判定为稳定。增大增益K会使 ωgc 向高频移动、PM 减小(如 K=10 时 PM≈18.0°),可观察奈奎斯特轨迹逐渐逼近(-1, 0)点。

工程实施注意

  1. 存在控制延迟(死时间)时,相位随频率线性下降,如水泵控制中50ms延迟造成15°/rad·s的相位损失,在ωgc附近PM急剧减小
  2. 含有积分环节(1/s)的系统,ω→0时增益趋于∞,奈奎斯特轨迹沿y轴伸展,低频稳定性验证必不可少
  3. 液压伺服和BLDC电动机因频率响应非线性,计算值与实测值通常有10~15%的误差,必须进行验证试验

依据标准与假设

依据/公式: 奈奎斯特稳定判别法 \(Z = P - N\)(\(N\):开环传递函数 \(L(j\omega)\) 包围(-1,0)的净圈数,逆时针为正,\(\omega:-\infty\to+\infty\);\(P\):开环右半平面极数)。稳定 \(\Leftrightarrow Z=0\)。裕度 \(GM=20\log_{10}(1/|L(j\omega_{pc})|)\)、\(PM=180^\circ+\angle L(j\omega_{gc})\)(Ogata/Nise)。

模型假设: 一阶/二阶有理传递函数。本工具的极点位于 \(s=-a,-b\)(\(a,b>0\)),恒在左半平面,故 \(P=0\)。稳定性由 \(L(j\omega)\) 频率扫描得到的增益/相位裕度判定。

适用范围与局限: 学习经典控制稳定裕度的教学工具。未计入纯滞后 \(e^{-Ls}\)、高阶模态与非线性,存在时轨迹会更复杂。设计目标参考:\(GM\ge 6\,\mathrm{dB}\)、\(PM\ge 30^\circ\)。