尼科尔斯线图模拟器 返回
控制工程

尼科尔斯线图模拟器

开环传递函数的频率响应,用相位(横轴)和增益(纵轴)合并到一个图上绘制的工具。改变增益、时常数、无差拍时间,增益裕度·相位裕度·共振峰会实时显示,可以直观地研究闭环系统有多稳定。

参数设置
开环增益 K
传递函数的比例增益。越大速度越快但裕度越小
时常数 τ₁
s
主导一阶滞后的时常数
时常数 τ₂
s
第二个一阶滞后的时常数(高频侧)
无差拍时间 T_d
s
传输延迟·计算延迟。单向延迟相位,削弱裕度
计算结果
增益裕度 GM (dB)
相位裕度 PM (°)
增益交叉频率 ωgc (rad/s)
相位交叉频率 ωpc (rad/s)
共振峰 Mr (dB)
稳定性判定
尼科尔斯线图 — 开环轨迹

横轴为开环相位(度),纵轴为开环增益(dB)。十字为临界点(−180°, 0 dB)。标记随频率 ω 增加在轨迹上移动。

开环 波德图(增益)
闭环频率响应 |T(jω)|
理论·主要公式

$$L(j\omega)=\frac{K\,e^{-j\omega T_d}}{j\omega\,(j\omega\tau_1+1)(j\omega\tau_2+1)}$$

分析对象的开环传递函数(I型系统)。K:开环增益,τ₁·τ₂:时常数,T_d:无差拍时间。包含一个积分器 1/(jω)。

$$\text{PM}=180^\circ+\angle L(\omega_{gc}),\qquad \text{GM}=-20\log_{10}|L(\omega_{pc})|$$

相位裕度 PM 是增益交叉频率 ωgc(|L|=1)处的相位,增益裕度 GM 是相位交叉频率 ωpc(∠L=−180°)的增益。

$$T(j\omega)=\frac{L(j\omega)}{1+L(j\omega)},\qquad M_r=\max_\omega\,20\log_{10}|T(j\omega)|$$

闭环传递函数 T 和共振峰 Mr。临界点(−180°, 0 dB)是闭环稳定性的边界,轨迹离它有多远就是裕度的表现。

尼科尔斯线图简介

🙋
我第一次听到"尼科尔斯线图"。它与波德图和奈奎斯特图有什么不同?
🎓
这三个都是用来看"开环传递函数 L(jω) 的频率响应"的图。波德图把增益和相位分成2个图。奈奎斯特图在复平面上用实部和虚部绘制。尼科尔斯线图介于两者之间,横轴是相位(度),纵轴是增益(dB),用1个图统合。可以理解为把波德的两个图合成一个。当频率 ω 变化时,那一点就在图上形成一条曲线。
🙋
合成成一个图有什么好处呢?
🎓
最大的优点是可以"用眼睛读出"增益裕度和相位裕度。看图上有没有十字标记的点?那是临界点(−180°, 0 dB),闭环稳定的分界点。轨迹越接近这个点越危险。相位裕度是在0 dB高度时从轨迹到临界点的水平间隙,增益裕度是在−180°竖线上的竖直间隙。左边的"开环增益 K"试着提高一下。轨迹会向上移动,靠近临界点,你就能看到裕度在减少。
🙋
我真的把 K 提到14左右,判定就变成"不稳定"了。发生了什么?
🎓
提高增益会让轨迹整体上移。相位不变但高度(dB)增加了。这样的话,穿过0 dB的点(增益交叉频率 ωgc)处的相位越来越接近 −180°,最后超过了。相位裕度 PM=180°+∠L(ωgc) 变成负数的时刻,那个闭环系统就发散了。实际的话,马达会尖叫,振动停不下来。所以"增益提高就会变快"无脑提高的话,一定在某个地方会不稳定。这是反馈制御的宿命啊。
🙋
增加无差拍时间 T_d 时,为什么裕度一下子减少这么多?我没有提高增益啊。
🎓
好问题。无差拍时间 e^(−jωT_d) 不改变增益(大小),但"单向地"将相位延迟 −ωT_d。而且频率 ω 越高延迟越大。在尼科尔斯线图上看,轨迹在高频侧会被一直向左拉,逼近临界点。所以加一点无差拍时间就会大幅削弱相位裕度。运输延迟的过程、网络的制御、采样粗糙的数字制御中,无差拍时间是威胁稳定性的最大敌人。看下面的闭环响应图,增加 T_d 时共振峰 Mr 也会增加,可以确认。
🙋
共振峰 Mr 最后怎么读呢?
🎓
Mr 是闭环频率响应 T(jω)=L/(1+L) 的山的高度。某个频率处输入增幅多大,即闭环响应有多振荡,超调有多大,Mr 就表示这个。粗略来说,Mr 在 0〜3 dB 很好,超过 5 dB 超调明显。Mr 与相位裕度有个反向的关系,PM 小的话 Mr 就大。所以设计时的目标是"相位裕度 45°左右,Mr 3 dB以下",这是定石。尼科尔斯线图还可以重叠 M等高线,那样的话看轨迹通过等高线的地方就能从图上直接读出共振峰。

常见问题

尼科尔斯线图是一个图表,其中开环传递函数 L(jω) 的"相位(横轴·度)"和"增益(纵轴·dB)"被合并到一个图上。当频率 ω 变化时,绘制的轨迹将波德图的增益图和相位图集成为一条曲线。最大的优势是增益裕度和相位裕度作为"临界点(−180°, 0 dB)的距离"一目了然。相位裕度表现为0 dB高度到临界点的水平间隙,增益裕度表现为−180°线上的竖直间隙。
首先对频率 ω 进行对数扫描,计算每个点的开环 L(jω) 的大小和相位。增益交叉频率 ωgc 是 |L|=1(0 dB)的频率,那里的相位为 ∠L(ωgc),相位裕度为 PM=180°+∠L(ωgc)。相位交叉频率 ωpc 是 ∠L 为 −180° 的频率,相位裕度从那里的增益计算为 GM=−20log₁₀|L(ωpc)| [dB]。本工具通过插值交叉点前后,实时显示这些值。
一般控制系统的目标是相位裕度 30〜60°,增益裕度 6 dB(约2倍)以上。相位裕度太小时,闭环响应会产生振荡,超调和建立时间增加。相反,过大则响应变慢。相位裕度 45〜60° 时,速度性和稳定性的平衡最好。本工具在 PM < 30° 或 GM < 6 dB 时警告"裕度较小",在 PM≤0 或 GM≤0 时判定"不稳定"。
共振峰 Mr 是闭环频率响应 T(jω)=L/(1+L) 的最大值(dB表示)。Mr 越大,某个频率处的输入增幅越强,即闭环响应越振荡,超调越大。经验上,Mr 在 0〜3 dB 时良好,超过 5 dB 时振荡明显,对应相位裕度较小的设计。Mr 与相位裕度密切相关,PM 越小 Mr 越大。

实际应用

伺服·运动制御:工作机械的进给轴、机器人关节、硬盘磁头定位等伺服系统,需要保证增益裕度和相位裕度,同时用尽可能高的增益提高速度性。尼科尔斯线图一眼就看出增益提高时轨迹如何接近临界点,可以直接用来决定增益的上限。如果叠加 M等高线,就能从图上直接读出保证目标共振峰以下的增益。

过程制御(化工·厂房):温度·流量·压力的控制回路包含配管的运输延迟和测量响应延迟,有大量无差拍时间。无差拍时间只延迟相位削弱相位裕度,所以尼科尔斯线图上高频侧的轨迹向左伸出。PID 增益的设定时,要确认包含无差拍时间的相位裕度是否充分。

电源·电力电子:开关电源和逆变器的电压·电流控制回路,设计最终确认也用频率响应看稳定裕度。数字制御的计算延迟·采样延迟作为无差拍时间有效,尼科尔斯线图或波德图会检查相位裕度 45°以上、增益裕度 6 dB以上的基准是否满足。

制御教育和回路整形:尼科尔斯线图作为回路整形(loop shaping)的教材很优秀。加入补偿器时轨迹如何动、进补偿加上去相位裕度如何回复,这些都可以视觉上理解。本工具中改变增益·时常数·无差拍时间,体感轨迹和裕度的关系,是培养制御工程感觉的很好的入口。

常见误解和注意事项

首先一个常见误解是,"增益裕度和相位裕度都充分就一定稳定性能好"。GM 和 PM 是只看一个交叉点的"局部"指标。复杂的传递函数中,增益曲线可能多次穿过0 dB,相位多次通过 −180°,会出现多个交叉点。那样的话不看最近的临界点距离(向量裕度或尼科尔斯线图上的最近距离)而判断,会见落危险的设计。本工具用的是I型单纯系统所以交叉点只有一对,但实机一定要确认整个轨迹的形状。

其次,"无差拍时间很小可以忽视"这个想法是错的。无差拍时间 e^(−jωT_d) 维持增益为1,只将相位延迟 −ωT_d。低频延迟很小,但增益交叉频率 ωgc 越高 ωgc·T_d 越大,相位裕度会大幅削弱。"为了提高速度提高增益"导致 ωgc 上升,同样的无差拍时间相位损失就增加,形成负的连锁。采样周期的一半程度的等效无差拍时间,传感器·执行器的延迟也要记得纳入模型。

最后,"尼科尔斯线图是开环的图,与闭环性能无关"这个误解也是错的。尼科尔斯线图是开环 L(jω) 的轨迹,但本来的强处是可以叠加闭环增益 |T| 的等高线(M等高线)和闭环相位的等高线(N等高线)。看轨迹通过哪里的 M等高线就能直接知道闭环共振峰 Mr 和带宽。也就是说,尼科尔斯线图是"开环的操作"和"闭环的性能"在一个图上连接的,不只是开环。本工具也是为此把开环轨迹和闭环响应 |T(jω)| 与 Mr 一起显示的。

使用指南

  1. 输入增益系数 k (0.5~5.0)、时常数 T1 (0.01~1.0s)、T2 (0.001~0.5s)、无差拍时间 Td (0~0.2s) 定义开环传递函数 G(s)=k/((1+T1·s)(1+T2·s))·e^(-Td·s)
  2. 在频率范围 0.01~100 rad/s 内自动计算,在尼科尔斯线图上绘制开环频率响应
  3. 读取增益裕度 (GM)、相位裕度 (PM)、共振峰 (Mr),判定闭环稳定性

具体计算例

工业过程制御系统:k=2.0、T1=0.5s、T2=0.1s、Td=0.05s 时,增益交叉频率 ωgc≈3.2 rad/s 处增益裕度 GM≈8.5 dB,相位交叉频率 ωpc≈1.8 rad/s 处相位裕度 PM≈35°,共振峰 Mr≈2.8 dB 判定"稳定"。无差拍时间增加时 PM 下降,Td=0.15s 时 PM≈12°。

实务中的注意

  1. 马达驱动系中无差拍时间 Td>0.2s 时,相位裕度下降到10°以下,外扰响应恶化,需要减少增益或设计补偿器
  2. 共振峰 Mr 超过 6 dB 时,闭环瞬态响应会出现超调,需要重新调整制御参数
  3. 化学反应器等一阶滞后系中,要正确同定 T1, T2,用频率响应测量值和模拟结果的差异进行验证