奈奎斯特图 返回
频率分析

奈奎斯特图·稳定裕量计算工具

选择开环传递函数并调整增益,即可实时绘制奈奎斯特图与波特图。自动计算增益裕量、相位裕量、穿越频率。

开环传递函数
传递函数形式
增益 K
极点 a
极点 b
零点 z
稳定裕量汇总
计算结果
GM (dB)
PM (°)
ωgc (rad/s)
稳定性
Nyquist
Note: 红色×为 (-1,0) 临界点。蓝色轨迹表示 ω>0,虚线表示 ω<0(镜像)。轨迹若逆时针包围 (-1,0),闭环会不稳定。
理论与主要公式

Nyquist稳定判据

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

$N$: Nyquist轨迹逆时针包围(-1,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的滑块,看看图是怎么变化的。
🙋
诶,真的吗?那“包围”到底是什么意思?怎么从图上看出系统稳不稳定?
🎓
问得好!我们有个黄金法则:闭环不稳定极点数 $Z = N + P$。$N$就是奈奎斯特轨迹逆时针绕(-1,0)的圈数,$P$是开环本身不稳定的极点数。如果算出来$Z=0$,系统就稳定。工程现场常见的是,当你把增益K调得太大,这条蓝线就会膨胀,把(-1,0)点包在里面,系统就不稳了。改变参数后你会看到,工具会自动计算并告诉你现在是稳定还是不稳定。
🙋
哦!那旁边的增益裕量和相位裕量又是干嘛的?数字越大越好吗?
🎓
这两个裕量就像是系统的“安全缓冲区”。比如在汽车巡航控制系统中,增益裕量告诉你还能把控制力度加大多少倍而不失稳;相位裕量则反映了系统响应有多“迟钝”,裕量太小,车子的速度就会来回震荡。并不是越大越好,太大意味着系统响应很慢。你可以在模拟器里调整极点a和b的位置,观察这两个裕量如何变化,体验一下在稳定性和快速性之间做权衡的感觉。

物理模型与关键公式

本工具的核心是分析开环传递函数 $G(s)H(s)$ 的频率响应。我们将 $s$ 替换为 $j\omega$,计算其在复平面上的轨迹。

$$G(j\omega)H(j\omega) = K \frac{(j\omega - z)}{(j\omega - a)(j\omega - b)}$$

其中,$K$是系统增益,$z$是零点,$a$和$b$是极点。当$\omega$从$0$变化到$+\infty$时,计算上述复数的实部和Imaginary Part,就得到了奈奎斯特图($\omega<0$的部分与之对称)。

判断系统稳定性的根本依据是奈奎斯特稳定性判据:

$$Z = N + P$$

$Z$:闭环系统在右半平面(RHP)的不稳定极点数,$Z=0$表示稳定。
$N$:奈奎斯特图($\omega$从$-\infty$到$+\infty$的完整闭合曲线)逆时针包围临界点$(-1, 0)$的净圈数。
$P$:开环传递函数 $G(s)H(s)$ 本身在右半平面(RHP)的极点数。

现实世界中的应用

航空航天(飞行控制):在设计飞机或航天器的自动驾驶仪时,必须确保在各种飞行包线内系统绝对稳定。工程师使用奈奎斯特图分析俯仰、滚转通道的控制律,验证增益裕量和相位裕量是否满足严苛的安全标准。

汽车工业(发动机与底盘控制):现代汽车的电子节气门控制、防抱死刹车系统(ABS)和主动悬架都是闭环控制系统。通过奈奎斯特分析,可以优化控制器参数,使发动机响应既迅速又平顺,刹车防抱死动作既有效又不引起车身剧烈抖动。

工业机器人(运动控制):机器人的关节伺服驱动需要高精度定位和快速响应。利用奈奎斯特图,可以设计补偿器(如PID调节),在避免机械臂共振和超调的同时,尽可能提升其运动速度与轨迹跟踪精度。

电力系统(电网稳定):大型电网中,发电机的励磁控制系统对维持电网电压和频率稳定至关重要。奈奎斯特稳定性分析用于评估在负载突变或故障情况下,整个电力系统是否能够保持同步运行,防止发生连锁崩溃事故。

常见误解与注意事项

首先要注意,“奈奎斯特曲线不包围(-1,0)点就绝对稳定”这种说法并不总是成立。因为应用奈奎斯特稳定性判据的大前提是开环系统本身稳定(P=0)。例如,在控制倒立摆这类原本不稳定的被控对象(P>0)时,这一前提就不成立。如果在模拟器中尝试将极点`a`或`b`设为负值(即存在右半平面极点),你会发现情况会发生变化。

其次,不要盲目相信稳定裕度的数值。例如,增益裕度(GM)为10dB、相位裕度(PM)为45度时,教科书上看起来似乎足够。但如果奈奎斯特轨迹曲线形状尖锐,或在(-1,0)点附近擦过,那么模型微小的误差或干扰就很容易使系统进入不稳定区域。请养成在使用此工具调节K值时,同时目视检查裕度的“质量”,即轨迹接近(-1,0)点的方式的习惯。

最后,要理解模拟器的传递函数模型与现实之间的差距。此工具处理的是非常简单的二阶系统模型。实际的控制对象必然包含时滞、更高阶的振动模式以及非线性特性。例如,加入时滞 $e^{-Ls}$ 后相位会进一步滞后,奈奎斯特轨迹将呈螺旋状无限环绕。请始终思考如何将工具中学到的基本原理扩展到更接近现实的复杂模型中。