1. 从传递函数到奈氏图:为什么我们非得画这条曲线
很多人在学自动控制原理的时候,到了第5章就开始犯迷糊。前面几章讲时域分析、讲根轨迹,好歹还能在纸上画一画、算一算,到了幅相频率特性曲线这里,突然就变成了一条绕来绕去的曲线,还动不动就绕原点转好几圈。我当初学这一章的时候也是一头雾水,直到后来做实际项目——给一个加热电炉做闭环温控——才发现奈氏图这东西是真的绕不过去。
先说清楚这一章到底在干什么。幅相频率特性曲线,也叫奈氏图(Nyquist Plot),是把一个系统的开环传递函数在复平面上画出来的一条曲线。横轴是实部,纵轴是虚部,频率从0变到无穷大,每一个频率点对应复平面上的一个点,把这些点连起来就是奈氏图。它解决的核心问题是:不用解闭环特征方程的根,直接看开环频率特性就能判断闭环系统稳不稳定。这就是奈氏稳定判据的威力。
你可能会问,根轨迹不也能判断稳定性吗?没错,但根轨迹有个前提——你得知道开环传递函数的零极点分布,而且根轨迹画的是闭环极点的移动轨迹,当系统阶次很高的时候,根轨迹画起来非常痛苦。奈氏图的好处在于,它只需要开环频率特性,而开环频率特性可以通过实验测量得到。这意味着什么?意味着你面对一个黑盒子系统,不知道内部结构,只要给它不同频率的正弦输入,测出输出幅值和相位,就能画出奈氏图,进而判断闭环稳定性。这在工程上太重要了。
这一章适合谁来学?如果你是自动化、电气工程、机械电子等专业的学生,正在啃自动控制原理,那这是必修内容。如果你是做电机控制、电源设计、温控系统开发的工程师,需要分析buck功率级传递函数或者加热电炉的闭环控制系统传递函数的稳定性,那奈氏图是你必须掌握的工具。哪怕你只是想在Simulink里搭个模型验证一下系统开环传递函数和反馈的关系,理解奈氏图也能让你少走很多弯路。
我写这篇东西,不是照搬教材,而是把我自己学这一章时踩过的坑、做项目时积累的经验,以及那些教材上不会明说但实际做题和工程中特别有用的技巧,全部整理出来。从传递函数怎么变成奈氏图,到怎么快速手绘草图,再到怎么用奈氏判据判断稳定性,每一步我都会说清楚“为什么这么做”和“怎么做才不会错”。
2. 奈氏图的核心逻辑:把频率特性搬到复平面上
2.1 频率特性到底是什么,为什么能判断稳定性
要理解奈氏图,先得理解频率特性。一个线性定常系统,给它一个正弦输入,稳态输出也是同频率的正弦,只是幅值和相位变了。幅值的变化叫幅频特性,相位的变化叫相频特性,合起来就是频率特性。数学上,把传递函数中的s换成jω,就得到了频率特性。
这里有个关键点很多人没想明白:为什么频率特性能判断稳定性?稳定性不是由闭环极点决定的吗?没错,但奈氏判据巧妙地建立了一个桥梁——通过开环频率特性绕临界点的圈数,来推断闭环极点在右半平面的个数。这个桥梁就是辐角原理。
辐角原理说的是,如果一个复变函数在某个区域内解析,那么当自变量沿一条闭合路径走一圈时,函数值的辐角变化量等于区域内零点和极点数之差乘以2π。奈氏把这个原理用到了控制系统上,构造了一条包围整个右半平面的路径,然后看开环传递函数在这条路径上的映射曲线绕(-1, j0)点转了多少圈。转的圈数就包含了闭环右极点的信息。
我当初理解这个的时候,卡了很久。后来想通了一个类比:你把开环传递函数想象成一个“探测器”,(-1, j0)点是一个“危险区域”。如果探测器绕危险区域转圈,说明系统内部有“隐患”——也就是闭环右极点。转的圈数越多,隐患越大。这个类比虽然不严谨,但帮助我建立了直觉。
2.2 从开环传递函数到奈氏图的完整映射过程
具体怎么从传递函数画出奈氏图?步骤其实很固定,但每一步都有讲究。
第一步,写出开环传递函数G(s)H(s)。注意是开环,不是闭环。很多人一开始就把闭环传递函数写出来,那就完全错了。开环传递函数是反馈信号和误差信号的比值,在方框图上就是从误差点到反馈点绕一圈的传递函数乘积。
第二步,令s = jω,得到G(jω)H(jω)。这一步是纯代数操作,把s替换成jω,然后整理成实部加虚部的形式。整理的时候要注意,分母通常有(1+jωT)这种因子,需要有理化。
第三步,取几个关键频率点计算幅值和相位。最关键的点包括:ω=0、ω=1/T(各环节的转折频率)、ω→∞。这些点算出来,曲线的骨架就出来了。
第四步,根据幅值和相位的变化趋势,把中间的点连起来。这里要注意相位的连续性,不能跳变。
第五步,考虑ω从-∞到0的部分。对于实系数传递函数,负频率部分的曲线和正频率部分关于实轴对称。所以画出一半,镜像一下就行。
我实测下来,手绘奈氏图最容易出错的地方是相位范围。比如一个包含积分环节的系统,ω→0时相位是-90°,曲线从无穷远处开始。如果包含两个积分环节,相位是-180°,曲线从负实轴的无穷远处开始。这些起始位置搞错了,整条曲线就全错了。
2.3 奈氏图与波特图的关系:一个硬币的两面
很多人学的时候会把奈氏图和波特图搞混,或者觉得它们是两个独立的东西。其实它们描述的是同一个频率特性,只是表现形式不同。波特图把幅值和相位分开画,横轴是对数频率;奈氏图把幅值和相位合在一起画在复平面上,频率是隐式参数。
这个关系非常有用。因为手绘奈氏图的时候,你可以先画波特图,从波特图上读出关键频率点的幅值和相位,然后搬到复平面上。反过来,如果你有了奈氏图,也可以反推波特图的大致形状。
我在做buck功率级传递函数分析的时候,就是先用波特图看清楚穿越频率和相位裕度,然后再画奈氏图确认绕圈情况。两个图对照着看,心里特别有底。尤其是当系统有多个转折频率的时候,波特图的渐近线近似能帮你快速定位奈氏图的关键点。
注意:奈氏图的频率轴不是线性的,也不是对数的,而是隐式的。曲线上的点密集程度反映了频率变化的快慢。低频段通常点比较稀疏,高频段点比较密集,因为高频时幅值衰减快,曲线向原点收缩。
3. 手把手拆解奈氏图绘制:从传递函数到曲线
3.1 典型环节的奈氏图长什么样
在画复杂系统的奈氏图之前,先把几个典型环节的奈氏图记熟。这就像学写字先学笔画一样,复杂曲线都是这些基本形状组合出来的。
比例环节K:奈氏图就是一个点,在实轴上的K处。因为幅值恒为K,相位恒为0。这个太简单了,但别忽略它,它决定了曲线整体离原点的远近。
积分环节1/s:令s=jω,得到1/(jω) = -j/ω。实部为0,虚部为-1/ω。当ω从0+变到∞时,曲线从负虚轴的无穷远处沿负虚轴向上趋近原点。相位恒为-90°。
惯性环节1/(1+jωT):这个环节的奈氏图是一个半圆。ω=0时,幅值为1,相位为0,点在(1, j0)。ω=1/T时,幅值为1/√2,相位为-45°。ω→∞时,幅值趋于0,相位趋于-90°。整个曲线在第四象限,是一个以(0.5, j0)为圆心、半径为0.5的半圆。这个半圆结论在教材上通常有推导,我建议你亲手推一遍,用实部虚部平方和等于半径平方的方法,推一次就记住了。
一阶微分环节1+jωT:和惯性环节关于实轴对称,曲线在第一象限,从(1, j0)出发,向上延伸。
振荡环节:这个比较复杂,奈氏图是一条从(1, j0)出发、最终趋于原点的曲线,但中间的形状取决于阻尼比ζ。ζ越小,曲线在谐振频率附近鼓得越厉害,甚至可能跑到实轴右边去。当ζ=0时,曲线在谐振频率处幅值趋于无穷,相位从0°跳到-180°。
延迟环节e^(-jωτ):幅值恒为1,相位随ω线性下降。奈氏图是一个单位圆,随着频率增加不断绕原点顺时针旋转。这个环节在工程中很常见,比如管道传输、信号采集都有延迟。
把这些典型环节的奈氏图记熟,后面画复杂系统的时候,就是把这些基本形状“乘”起来。幅值相乘,相位相加,反映在奈氏图上就是曲线被拉伸旋转。
3.2 绘制开环奈氏图的标准化流程
现在来看一个完整的开环传递函数怎么画。我以最常见的三阶系统为例:
G(s)H(s) = K / [s(s+1)(s+2)]
第一步,标准化。把每个因子写成(1+jωT)的形式。分母的s+1写成1+jω,s+2写成2(1+jω/2)。所以:
G(jω)H(jω) = K / [jω · (1+jω) · 2(1+jω/2)] = K/2 / [jω(1+jω)(1+jω/2)]
第二步,确定起始点和终止点。
ω→0时,积分环节1/jω起主导作用,幅值趋于无穷,相位为-90°。所以曲线从负虚轴无穷远处开始。
ω→∞时,分母有三个一阶因子,总相位为-90°-90°-90°=-270°,幅值趋于0。所以曲线从第三象限趋于原点,趋近方向是-270°,也就是正虚轴方向(因为-270°等价于+90°)。
第三步,计算关键频率点的幅值和相位。
ω=1时:G(j1)H(j1) = K/2 / [j1 · (1+j1) · (1+j0.5)] 分母 = j1 · (1+j1) · (1+j0.5) = j1 · (1+j1) · (1+j0.5) 先算(1+j1)(1+j0.5) = 1 + j0.5 + j1 + j²0.5 = 1 + j1.5 - 0.5 = 0.5 + j1.5 再乘j1:j1 · (0.5+j1.5) = j0.5 + j²1.5 = -1.5 + j0.5 所以G(j1)H(j1) = (K/2) / (-1.5+j0.5) = (K/2)(-1.5-j0.5) / (1.5²+0.5²) = (K/2)(-1.5-j0.5)/2.5 = K(-0.3-j0.1)
当K=1时,这个点在(-0.3, -0.1),在第三象限。
ω=2时:G(j2)H(j2) = K/2 / [j2 · (1+j2) · (1+j1)] 分母 = j2 · (1+j2) · (1+j1) = j2 · (1+j1+j2+j²2) = j2 · (1+j3-2) = j2 · (-1+j3) = -j2 + j²6 = -6 - j2 所以G(j2)H(j2) = (K/2) / (-6-j2) = (K/2)(-6+j2)/40 = K(-0.075+j0.025)
这个点在第二象限,实部为负,虚部为正。
第四步,判断曲线是否穿过负实轴。令虚部为0,解出穿越频率。这个计算比较繁琐,但可以用相位条件:当总相位为-180°时,曲线穿过负实轴。对于这个系统,总相位为-90° - arctan(ω) - arctan(ω/2)。令其等于-180°,得到arctan(ω) + arctan(ω/2) = 90°。两边取正切,利用tan(A+B)公式,当A+B=90°时,tanA·tanB=1。所以ω·(ω/2)=1,ω²=2,ω=√2≈1.414。
在ω=√2处,计算幅值。|G(j√2)H(j√2)| = K/2 / [√2 · √(1+2) · √(1+0.5)] = K/2 / [√2 · √3 · √1.5] = K/2 / [√2 · √4.5] = K/2 / [√9] = K/6。
所以穿越点是(-K/6, j0)。当K=1时,穿越点在(-0.167, j0)。
第五步,把关键点连起来。从负虚轴无穷远处开始,经过第三象限,穿过负实轴于(-K/6, j0),进入第二象限,最终趋于原点,趋近方向为正虚轴。
这就是完整的绘制流程。我建议你拿几张方格纸,亲手画几遍,比看十遍书都管用。
3.3 参数变化对奈氏图的影响:K值、时间常数、积分环节个数
奈氏图的形状不是固定的,参数一变,曲线就变。理解这些变化规律,对分析系统特别重要。
K值的影响:K只影响幅值,不影响相位。所以K值变化时,奈氏图整体“缩放”,形状不变。K增大,曲线向外扩张,穿越点离原点更远;K减小,曲线向内收缩。这个特性在分析稳定性时非常有用——你可以先画K=1的曲线,然后根据K的变化判断穿越点相对于(-1, j0)的位置。
时间常数的影响:时间常数影响转折频率的位置。时间常数增大,转折频率减小,曲线在更低频率处开始弯曲。多个时间常数之间的比例关系,决定了曲线在中间频段的形状。如果时间常数相差很大,曲线会出现明显的“台阶”。
积分环节个数的影响:这是最关键的。没有积分环节,曲线从实轴上的有限点开始;一个积分环节,曲线从负虚轴无穷远处开始;两个积分环节,曲线从负实轴无穷远处开始;三个积分环节,曲线从正虚轴无穷远处开始。每增加一个积分环节,起始相位就多-90°。这个规律一定要记死。
我在做加热电炉的闭环控制系统传递函数分析时,电炉本身是一个大惯性环节,加上传感器和控制器,整个开环传递函数可能包含一个积分环节(控制器中的积分作用)和多个惯性环节。画奈氏图的时候,起始点就在负虚轴无穷远处,然后随着频率增加,曲线逐渐向原点收缩。如果积分时间常数选得太大,曲线可能会在低频段就穿过负实轴,导致稳定性问题。
实操心得:画奈氏图之前,先把传递函数写成标准形式,数清楚积分环节的个数和惯性环节的个数。积分环节决定起始位置,惯性环节决定高频段的趋近方向。这两个信息确定了,曲线的骨架就稳了。
4. 奈氏稳定判据:从曲线绕圈到稳定性结论
4.1 奈氏判据的两种表述方式
奈氏判据有两种常用的表述,一种是基于绕圈数的,一种是基于穿越次数的。两种都要会,因为不同教材和不同题目可能要求不同的表述。
绕圈数表述:闭环系统稳定的充要条件是,当ω从-∞变到+∞时,开环奈氏图绕(-1, j0)点逆时针转的圈数等于开环传递函数在右半平面的极点数。如果开环传递函数在右半平面没有极点(也就是开环稳定),那么闭环稳定的充要条件就是奈氏图不包围(-1, j0)点。
这个表述很简洁,但实际用的时候要注意:绕圈数要算代数和,逆时针为正,顺时针为负。而且必须是ω从-∞到+∞的完整曲线,不能只看正频率部分。
穿越次数表述:对于开环稳定的系统,闭环稳定的充要条件是,当ω从0变到+∞时,奈氏图在(-1, j0)点左侧的负实轴上的正穿越次数等于负穿越次数。正穿越是指曲线从上往下穿过负实轴,负穿越是指从下往上穿过。如果正穿越次数大于负穿越次数,说明曲线包围了(-1, j0)点,闭环不稳定。
我个人的习惯是先用穿越次数法快速判断,因为只需要看正频率部分,计算量小。如果穿越情况复杂,再用绕圈数法仔细验证。
4.2 穿越点的计算与稳定性判断实例
还是用上面的例子:G(s)H(s) = K / [s(s+1)(s+2)]。开环传递函数在右半平面没有极点,所以开环稳定。闭环稳定的条件是奈氏图不包围(-1, j0)点。
穿越点在(-K/6, j0)。要让这个点不在(-1, j0)的左边,也就是-K/6 > -1,即K/6 < 1,K < 6。所以K的稳定范围是0 < K < 6。
当K=6时,穿越点正好在(-1, j0),系统临界稳定。当K>6时,穿越点在(-1, j0)左边,曲线包围了(-1, j0)点,闭环不稳定。
这个结论和用劳斯判据算出来的结果完全一致。你可以用劳斯判据验证一下:特征方程是s³+3s²+2s+K=0,劳斯表第一列全正的条件是K<6。两种方法殊途同归。
我再举一个开环不稳定的例子。假设G(s)H(s) = K / [s(s-1)(s+2)],开环有一个右极点。这时候奈氏判据要修正:闭环稳定的条件是奈氏图绕(-1, j0)点逆时针转的圈数等于1(右极点数)。对于这个系统,穿越点的计算会复杂一些,但思路是一样的。关键是记住:开环不稳定时,奈氏图必须包围(-1, j0)点,闭环才可能稳定。这听起来反直觉,但正是反馈的魔力——用不稳定的开环通过反馈变成稳定的闭环。
4.3 相位裕度和幅值裕度在奈氏图上的体现
相位裕度和幅值裕度是频域分析的两个核心指标,它们在奈氏图上也有直观的体现。
幅值裕度:奈氏图与负实轴的交点,如果这个交点在(-1, j0)的右边,幅值裕度就是交点幅值的倒数。比如交点(-0.5, j0),幅值裕度就是1/0.5=2,也就是6dB。如果交点在(-1, j0)的左边,幅值裕度为负,系统不稳定。
相位裕度:奈氏图与单位圆的交点,这个点的相位与-180°的差值就是相位裕度。比如交点的相位是-150°,相位裕度就是30°。
在奈氏图上,这两个裕度可以直观地看出来:曲线离(-1, j0)点越远,裕度越大,系统越稳定。曲线越靠近(-1, j0)点,裕度越小,系统越接近不稳定。
我在做buck功率级传递函数分析时,特别关注相位裕度。因为buck变换器的LC滤波器会产生一个双极点,在谐振频率处相位急剧下降。如果补偿网络设计不好,相位裕度可能只有十几度,系统虽然稳定但动态响应很差,负载突变时振荡很久。用奈氏图看,就是曲线在谐振频率附近几乎贴着(-1, j0)点过去,看着就让人不放心。
注意:奈氏图上的相位裕度和波特图上的相位裕度是同一个东西,但奈氏图更直观地展示了曲线与临界点的“距离”。我通常两个图一起看,波特图算数值,奈氏图看趋势。
5. 工程实操中的奈氏图:从Simulink仿真到实验测量
5.1 用Simulink快速获取开环频率特性
现在很少有人真的手绘奈氏图了,大部分时候是用工具。Simulink是控制工程最常用的仿真工具,获取频率特性非常方便。
具体操作:在Simulink里搭好系统模型,用Linear Analysis工具或者直接用MATLAB的bode和nyquist命令。如果你有Simulink Control Design工具箱,可以直接在模型上设置输入输出点,然后一键生成奈氏图。
我通常的做法是:先用MATLAB脚本定义传递函数,用nyquist命令画图,调整参数看曲线变化。确认方案后,再在Simulink里搭详细模型验证。这样效率最高。
比如分析simulink电机传递函数,我会先在MATLAB里写:
s = tf('s'); G = K / (s*(s+1)*(s+2)); nyquist(G); grid on;然后调整K值,观察曲线与(-1, j0)点的关系。MATLAB会自动画出完整的奈氏图,包括负频率部分,还会标注穿越点。你还可以用margin命令直接算出幅值裕度和相位裕度。
5.2 实验测量频率特性的方法
如果你面对的是一个实际系统,没有传递函数,怎么画奈氏图?答案是频率响应实验。
基本方法:给系统输入不同频率的正弦信号,测量输出的幅值和相位,然后计算频率特性。具体步骤:
选择一个频率范围,覆盖系统的主要动态特性。通常从低于第一个转折频率一个数量级开始,到高于最后一个转折频率一个数量级结束。
在每个频率点,等系统稳态后,记录输入幅值A_in、输出幅值A_out、输入输出之间的相位差φ。
计算幅频特性|G(jω)| = A_out/A_in,相频特性∠G(jω) = φ。
把幅值和相位画在复平面上,得到奈氏图。
这个方法听起来简单,但实际操作有很多坑。比如噪声干扰、非线性失真、测量时间过长导致系统漂移等。我建议用扫频仪或者专门的频率响应分析仪,比手动逐点测量效率高得多。
对于加热电炉的闭环控制系统,我做过类似的实验。电炉的加热功率和温度之间有很大的滞后,频率响应实验需要等每个频率点都稳定,一个完整的扫频可能要几个小时。后来我改用阶跃响应辨识传递函数,再用MATLAB画奈氏图,效率高很多。
5.3 奈氏图在控制器设计中的实际应用
奈氏图不只是用来判断稳定性的,它在控制器设计中也有直接应用。
增益调整:通过观察奈氏图与(-1, j0)点的距离,可以快速确定增益的调整方向。如果曲线离(-1, j0)太近,减小增益;如果太远,可以适当增大增益提高响应速度。
补偿网络设计:超前补偿和滞后补偿在奈氏图上的效果很直观。超前补偿会把曲线在特定频率段“拉”向远离(-1, j0)的方向,增加相位裕度;滞后补偿则是在低频段提高增益,同时保持高频段的形状不变。
谐振抑制:对于有谐振峰的系统,奈氏图在谐振频率附近会有一个明显的“鼓包”。如果这个鼓包靠近(-1, j0)点,系统就会有很强的振荡倾向。这时候需要设计陷波器或者调整补偿网络来抑制谐振峰。
我在设计buck功率级的电压环补偿器时,就是先用奈氏图看清楚LC谐振峰的位置和高度,然后设计一个双极点双零点的补偿网络,把谐振峰压下去,同时保证足够的相位裕度。整个过程在MATLAB里迭代几次就能搞定,比盲目试参数快得多。
6. 常见问题与排查技巧实录
6.1 奈氏图绘制中的典型错误
我批改过很多学生的作业,也带过不少新人工程师,发现奈氏图绘制中的错误高度集中。整理成表格,方便你对照检查。
| 错误类型 | 典型表现 | 根本原因 | 纠正方法 |
|---|---|---|---|
| 起始点错误 | 曲线从实轴开始,但系统有积分环节 | 忘记积分环节的1/jω在ω→0时趋于无穷 | 数清楚积分环节个数,确定起始相位 |
| 终止点错误 | 曲线趋于原点方向不对 | 高频段相位计算错误 | 把所有一阶因子的相位相加,注意符号 |
| 穿越点计算错误 | 穿越频率或穿越幅值算错 | 相位条件解方程出错 | 用arctan(A)+arctan(B)=90°时tanA·tanB=1的技巧 |
| 负频率部分画错 | 镜像方向反了 | 忘记实系数传递函数的共轭对称性 | 负频率曲线与正频率关于实轴对称 |
| 绕圈数判断错误 | 稳定性结论与劳斯判据矛盾 | 绕圈方向或圈数计算错误 | 用穿越次数法交叉验证 |
6.2 稳定性判断中的常见困惑
困惑一:开环不稳定时,为什么奈氏图包围(-1, j0)反而稳定?
这个问题我当初也想了很久。关键在于,奈氏判据的完整表述是:闭环右极点数 = 开环右极点数 - 奈氏图绕(-1, j0)的逆时针圈数。如果开环有一个右极点,要让闭环右极点数为零,奈氏图必须逆时针绕(-1, j0)一圈。所以包围(-1, j0)是必要的。这体现了反馈的“化腐朽为神奇”——用不稳定的开环通过反馈变成稳定的闭环。
困惑二:奈氏图和根轨迹判断稳定性,结果会矛盾吗?
不会。两者都是基于闭环特征方程的,只是从不同角度分析。如果结果矛盾,一定是计算错了。我建议用两种方法交叉验证,特别是考试的时候,可以检查答案。
困惑三:系统有延迟环节时,奈氏图怎么画?
延迟环节e^(-jωτ)的幅值恒为1,相位为-ωτ。在奈氏图上,它表现为一个单位圆,随着频率增加不断旋转。画图的时候,先画不含延迟的奈氏图,然后每个点乘以e^(-jωτ),也就是幅值不变、相位减去ωτ。延迟环节会让曲线“卷”起来,绕(-1, j0)点的可能性大大增加,所以延迟是稳定性的天敌。
6.3 独家避坑技巧
技巧一:先画波特图,再画奈氏图。波特图的渐近线近似能帮你快速确定关键频率点和相位变化趋势。我通常先在波特图上标出穿越频率和相位裕度,然后把这些信息搬到奈氏图上,效率翻倍。
技巧二:用MATLAB验证手绘结果。手绘奈氏图后,用MATLAB的nyquist命令画一遍,对比关键点的位置。如果差异很大,说明手绘过程中有错误。这个习惯让我在考试中避免了很多低级错误。
技巧三:记住几个典型系统的奈氏图模板。比如一阶系统、二阶系统、三阶系统的奈氏图大致形状,考试时可以直接套用,节省时间。
技巧四:穿越点计算用相位条件,不要用实部虚部分离。实部虚部分离计算量大,容易出错。用相位条件解方程,配合arctan的和角公式,快得多。
技巧五:稳定性判断先用穿越次数法,再用绕圈数法验证。穿越次数法只需要正频率部分,计算量小,适合快速判断。如果穿越情况复杂,再用绕圈数法仔细算。
提示:奈氏图这一章,手绘能力很重要,但更重要的是理解曲线背后的物理意义。曲线为什么从那里开始、为什么往那个方向走、为什么穿过负实轴,这些“为什么”比“怎么画”更重要。理解了物理意义,画图就是水到渠成的事。
7. 从奈氏图到实际系统:几个工程案例的思考
7.1 buck变换器电压环的奈氏图分析
buck变换器的功率级传递函数包含一个LC滤波器,在谐振频率处有一个双极点,相位下降180°。如果电压环的穿越频率设在谐振频率附近,相位裕度会非常小,系统容易振荡。
用奈氏图分析:功率级的奈氏图在谐振频率附近会有一个明显的“鼓包”,如果这个鼓包靠近(-1, j0)点,闭环稳定性就差。补偿网络的作用就是把曲线在谐振频率附近“拉”开,增加与(-1, j0)点的距离。
我通常的做法是:先画功率级的奈氏图,找到谐振峰的位置和高度;然后设计补偿网络,让补偿后的开环奈氏图在谐振频率附近有足够的相位裕度。这个过程在MATLAB里迭代几次就能搞定。
7.2 加热电炉温控系统的稳定性设计
加热电炉是一个大惯性环节,传递函数可以近似为一阶惯性加延迟。延迟环节让奈氏图“卷”起来,稳定性变差。
用奈氏图分析:电炉的奈氏图从实轴上的某点开始,随着频率增加,曲线向原点收缩,同时因为延迟环节的存在,曲线不断旋转。如果延迟时间太大,曲线可能在低频段就包围(-1, j0)点,导致闭环不稳定。
解决方案:降低控制器增益,或者采用Smith预估器补偿延迟。降低增益会让曲线整体收缩,远离(-1, j0)点;Smith预估器则是把延迟环节移到反馈回路外面,让奈氏图不再“卷”。
7.3 电机控制系统的频率特性分析
电机的传递函数通常包含电气时间常数和机械时间常数,两者相差很大。电气时间常数小,对应高频段;机械时间常数大,对应低频段。
用奈氏图分析:电机的奈氏图在低频段主要受机械时间常数影响,高频段主要受电气时间常数影响。如果两个时间常数相差很大,曲线在中间频段会有一个明显的“台阶”。
对于simulink电机传递函数的仿真,我通常先用MATLAB辨识传递函数,然后画奈氏图分析稳定性,最后在Simulink里搭详细模型验证。这个流程适用于大多数电机控制系统。
8. 写在最后:我学奈氏图的一点体会
奈氏图这一章,我前前后后学了三四遍才真正搞明白。第一遍学的时候,只记住了怎么画图,不知道为什么要画;第二遍学的时候,理解了频率特性和稳定性的关系,但还是不会用;第三遍是在做实际项目的时候,被逼着用奈氏图分析稳定性,才真正掌握了。
我的体会是:奈氏图不是画出来的,是“想”出来的。你脑子里先有曲线的大致形状,知道它从哪里来、往哪里去、中间怎么弯,然后才是动手画。如果一上来就埋头计算,很容易迷失在代数运算里。
另外,工具是辅助,不是替代。MATLAB能帮你画图,但不能帮你理解。我见过很多人,MATLAB用得很溜,但问他奈氏图为什么能判断稳定性,说不出来。这样遇到实际问题还是抓瞎。
最后分享一个小技巧:如果你觉得奈氏图太抽象,可以把它和波特图对照着看。波特图上的幅值穿越频率对应奈氏图上曲线与单位圆的交点,相位穿越频率对应奈氏图上曲线与负实轴的交点。两个图一起看,抽象的东西就具体了。
这个内容后续还可以这样扩展:把奈氏判据和对数判据结合起来,分析更复杂的系统;或者用奈氏图设计鲁棒控制器,考虑参数不确定性的影响。这些内容在实际工程中很有用,但教材上通常讲得不多,需要自己在项目中慢慢积累。