1. 从“系统会不会振”说起:为什么第五章是自动控制原理的分水岭
学自动控制原理的人,到了第五章通常会经历一次明显的“难度跃迁”。前面几章讲时域分析、根轨迹,好歹还能在纸上画出个大概,到了奈奎斯特稳定判据、相角裕度、幅值裕度、对数频率特性这一块,很多人的第一反应是:这些东西到底在工程上干什么用?
我当年第一次学这一章的时候也有同样的困惑。后来做实际项目才发现,这一章的内容是整个经典控制理论里最“工程化”的部分。它回答的核心问题非常直接:一个闭环系统,在开环频率特性已知的情况下,它到底稳不稳定?如果稳定,它离不稳定还有多远?
这两个问题在工程上极其重要。因为实际系统里,被控对象的参数会漂移,元器件会老化,工作点会变化。你在实验室调好的一个控制器,装到现场可能因为电机温度升高导致参数变化,系统就开始振荡甚至发散。频率法给出的相角裕度和幅值裕度,本质上就是在量化“这个系统有多抗造”。
这一章的知识链条是这样的:先建立频率特性的概念,然后引出奈奎斯特稳定判据作为理论根基,再由判据衍生出相角裕度和幅值裕度这两个工程指标,最后落到对数频率特性(Bode图)这个最常用的工程工具上。整条线索一环扣一环,缺了任何一环,后面的工程应用都会变成死记公式。
这篇文章我打算按自己学习和实际使用的顺序来写,把每个概念的来龙去脉、背后的直觉、手算方法、以及工程上容易踩的坑都讲清楚。适合正在学这门课的学生,也适合工作中需要快速回顾频率法分析的工程师。不需要你数学功底多深,但基本的复数和拉普拉斯变换概念得有。
2. 奈奎斯特稳定判据:从围线到“绕圈数”的直觉建立
2.1 判据到底在说什么
奈奎斯特稳定判据的数学表述很简洁:闭环系统稳定的充要条件是,当频率从负无穷变化到正无穷时,开环传递函数的奈奎斯特曲线绕临界点(-1, j0)逆时针包围的圈数,等于开环传递函数在右半平面的极点数。
公式写出来就是 Z = P - N,其中Z是闭环右半平面极点数,P是开环右半平面极点数,N是奈奎斯特曲线绕(-1, j0)逆时针包围的圈数。系统稳定要求Z=0,也就是N=P。
这个公式本身不难记,难的是理解它为什么成立。我当初学的时候,老师花了两节课讲辐角原理和保角变换,听得云里雾里。后来自己琢磨出一个比较直观的理解方式:奈奎斯特判据本质上是在用开环的频率响应信息,去推断闭环特征方程在右半平面有没有根。它绕过了直接求解闭环特征方程这个难题,因为高阶系统的闭环特征方程往往很难求根。
注意:很多教材在讲这个判据时,默认开环传递函数在虚轴上没有极点。如果开环有积分环节(也就是在原点有极点),奈奎斯特曲线需要做“无穷大半圆”的修正。这个修正步骤是考试和工程计算中最容易出错的地方,后面会专门讲。
2.2 辐角原理的通俗解释
辐角原理是奈奎斯特判据的数学基础。它的核心思想是:如果一个复变函数F(s)在某个闭合围线内部有Z个零点和P个极点,那么当s沿围线走一圈时,F(s)的轨迹绕原点逆时针转的圈数N满足N = P - Z。
你可以这样理解:每个零点让轨迹“多绕一圈”,每个极点让轨迹“少绕一圈”。这个性质对任何解析函数都成立,跟函数的具体形式无关。
奈奎斯特判据的巧妙之处在于,它选择的围线是整个右半平面(由虚轴和无穷大半圆组成),选择的函数是F(s) = 1 + G(s)H(s),也就是闭环特征多项式。这样,F(s)的零点就是闭环极点,F(s)的极点就是开环极点。于是辐角原理就变成了:闭环右半平面极点数 = 开环右半平面极点数 - 绕圈数。
实际使用时,我们不去画F(s) = 1 + G(s)H(s)的轨迹,而是画G(s)H(s)的轨迹,然后看它绕(-1, j0)点的圈数。因为F(s)绕原点的圈数等价于G(s)H(s)绕(-1, j0)的圈数,只是把坐标原点平移了一下。
2.3 手把手画奈奎斯特曲线的步骤
画奈奎斯特曲线是这一章的实操重点。我总结了一套比较顺手的流程,按这个顺序走基本不会乱:
确定开环传递函数的标准形式。把G(s)H(s)写成时间常数形式,也就是每个因式都写成(1 + Ts)或者(1 + 2ζs/ωn + s²/ωn²)的形式。这一步很关键,因为后面的低频段斜率、转折频率都依赖这个形式。
判断开环极点在右半平面的个数P。通常P=0,也就是开环稳定。如果开环不稳定,P就是右半平面极点的个数。
分析低频段行为。令s→0,看G(s)H(s)的幅值和相角。如果有积分环节(1/s^k),低频段幅值趋于无穷,相角趋于-k×90°。
分析高频段行为。令s→∞,看幅值和相角。通常分母阶数高于分子阶数,幅值趋于0,相角趋于-(n-m)×90°,其中n是分母阶数,m是分子阶数。
找与实轴的交点。令虚部为零,解出对应的频率和实部值。这个交点是判断包围圈数的关键。
处理虚轴上的极点。如果有积分环节,需要从低频段开始,沿无穷大半圆逆时针绕过原点,补上这段“修正曲线”。
数包围圈数。从ω=0+开始,沿曲线走到ω=+∞,然后镜像得到ω从-∞到0-的部分,看整条闭合曲线绕(-1, j0)逆时针转了几圈。
这套流程看起来步骤多,但实际操作几次之后就形成条件反射了。我建议初学者至少手画五到六个不同类型的传递函数,把典型环节(惯性、积分、振荡、微分)的奈奎斯特曲线形状记熟。
2.4 积分环节的修正:最容易翻车的地方
开环传递函数里有1/s的时候,奈奎斯特曲线在原点处是断开的。因为s=0时G(s)H(s)趋于无穷,曲线没法直接画出来。这时候需要做一个“无穷大半圆”的修正。
修正的方法是:从ω=0+对应的点开始,以无穷大为半径,逆时针绕过原点,画到ω=0-对应的点。这个半圆扫过的角度是k×180°,其中k是积分环节的个数。
为什么是逆时针?因为围线本身是顺时针绕过右半平面的,映射到G(s)H(s)平面上就变成了逆时针。这个方向问题我见过太多人搞反,包括我自己当年考试的时候。
提示:如果开环传递函数有k个积分环节,低频段相角从-k×90°开始。修正半圆从ω=0+的点开始,逆时针转k×180°,到达ω=0-的点。这个修正半圆不贡献包围圈数,但它决定了曲线的起点和终点位置。
2.5 一个完整的判稳实例
拿一个典型的三阶系统来练手:G(s)H(s) = K / [s(s+1)(s+2)]。
先写成时间常数形式:G(s)H(s) = K / [s × (1+s)(1+0.5s)] × (1/2)。为了简化,令K' = K/2,则G(s)H(s) = K' / [s(1+s)(1+0.5s)]。
开环极点:s=0, s=-1, s=-2,全部在左半平面或原点,所以P=0。系统稳定要求N=0,也就是奈奎斯特曲线不包围(-1, j0)。
低频段:s→0时,幅值→∞,相角→-90°。有一个积分环节,需要做180°的逆时针修正半圆。
高频段:s→∞时,幅值→0,相角→-270°。
找与实轴交点:令G(jω)H(jω)的虚部为零。计算可得ω=√2时,实部为-K'/3。所以曲线与负实轴的交点在(-K'/3, j0)。
要保证不包围(-1, j0),需要-K'/3 > -1,也就是K' < 3,即K < 6。这就是使系统稳定的开环增益范围。
这个例子很经典,考试出现的频率极高。关键点在于:积分环节的修正半圆方向、与实轴交点的计算、以及最终稳定条件的推导。
3. 相角裕度和幅值裕度:量化“稳定程度”的两个标尺
3.1 为什么需要裕度概念
奈奎斯特判据只回答“稳不稳定”这个是非题,但工程上更关心的是“有多稳定”。一个系统即使稳定,如果奈奎斯特曲线紧贴着(-1, j0)点走,那稍微有点参数变化就可能越过临界点,变成不稳定。这种系统在实际中是不能用的。
相角裕度和幅值裕度就是用来量化这种“稳定余量”的。它们的定义都建立在开环频率特性的基础上,物理意义非常明确。
3.2 幅值裕度:增益还能加多少
幅值裕度(Gain Margin,GM)的定义是:在相角为-180°的频率处,开环幅值的倒数。用分贝表示就是GM = -20lg|G(jω)H(jω)|,其中ω是相角穿越频率(也叫相位交界频率),满足∠G(jω)H(jω) = -180°。
这个定义的物理意义是:如果我把开环增益乘以GM倍,系统刚好达到临界稳定状态。所以GM越大,系统对增益变化的容忍度越高。
工程上一般要求GM ≥ 6dB,也就是增益可以增加一倍左右系统仍然稳定。对于要求更高的场合,比如航空航天领域,可能要求GM ≥ 10dB甚至更高。
注意:幅值裕度的定义在不同教材里可能有细微差别。有的用倒数定义,有的直接用|G|的dB值。关键是看最终判断标准:GM > 1(或GM > 0dB)表示稳定。我建议统一用“增益还能放大多少倍”这个物理意义来理解,不容易搞混。
3.3 相角裕度:延迟还能加多少
相角裕度(Phase Margin,PM)的定义是:在开环幅值为1(即0dB)的频率处,开环相角与-180°的差值。用公式表示就是PM = 180° + ∠G(jωc)H(jωc),其中ωc是幅值穿越频率(也叫剪切频率),满足|G(jωc)H(jωc)| = 1。
这个定义的物理意义是:如果系统再增加PM这么多的相位滞后,系统刚好达到临界稳定。相位滞后在实际中可能来自信号传输延迟、滤波器、或者额外的惯性环节。
工程上一般要求PM在30°到60°之间。PM太小,系统响应会剧烈振荡;PM太大,系统响应会过于迟缓。45°左右通常是一个比较均衡的选择。
3.4 两个裕度的计算实例
还是用前面的例子:G(s)H(s) = K' / [s(1+s)(1+0.5s)],取K' = 1来算。
先求幅值穿越频率ωc:|G(jω)H(jω)| = 1。
|G(jω)H(jω)| = K' / [ω × √(1+ω²) × √(1+0.25ω²)]
令K'=1,解方程1 / [ω × √(1+ω²) × √(1+0.25ω²)] = 1。数值求解得ωc ≈ 0.786 rad/s。
在ωc处计算相角:∠G(jωc)H(jωc) = -90° - arctan(0.786) - arctan(0.5×0.786) = -90° - 38.2° - 21.5° = -149.7°。
所以PM = 180° - 149.7° = 30.3°。
再求相角穿越频率ωg:∠G(jω)H(jω) = -180°。
-90° - arctan(ω) - arctan(0.5ω) = -180°,即arctan(ω) + arctan(0.5ω) = 90°。
利用arctan(A) + arctan(B) = 90°当且仅当AB=1,得ω × 0.5ω = 1,即ω² = 2,ωg = √2 ≈ 1.414 rad/s。
在ωg处计算幅值:|G(jωg)H(jωg)| = 1 / [1.414 × √(1+2) × √(1+0.5)] = 1 / [1.414 × 1.732 × 1.225] = 1/3 = 0.333。
所以GM = 1/0.333 = 3,用dB表示就是20lg(3) ≈ 9.5dB。
这个系统PM=30.3°,GM=9.5dB,属于稳定但相角裕度偏小的类型。实际中如果对动态性能要求高,可能需要调整增益或者加校正环节。
3.5 裕度与动态性能的对应关系
相角裕度和幅值裕度不只是稳定性指标,它们和闭环系统的动态性能有很强的相关性。根据经验:
| 相角裕度 | 超调量 | 动态响应特征 |
|---|---|---|
| 20°以下 | 很大 | 剧烈振荡,接近不稳定 |
| 30°左右 | 约30% | 有明显振荡,但可接受 |
| 45°左右 | 约15% | 响应较快,超调适中 |
| 60°以上 | 很小 | 响应迟缓,无超调 |
这个对应关系不是严格的数学定理,但在大量二阶系统和近似二阶系统中验证过,工程上很有参考价值。我实际调控制器的时候,经常先把PM调到45°左右,然后再根据具体需求微调。
幅值裕度主要反映系统对增益变化的鲁棒性。GM小于6dB的系统,如果被控对象增益随工况变化较大,就很容易出问题。比如温度升高导致电机电阻变化,或者负载变化导致增益漂移,这些都会吃掉幅值裕度。
4. 对数频率特性:工程上最顺手的分析工具
4.1 为什么用Bode图而不是奈奎斯特图
奈奎斯特图在理论推导上很漂亮,但实际工程分析中,用得最多的还是对数频率特性,也就是Bode图。原因很简单:Bode图把幅值和相角分开画,横轴用对数坐标,这样可以把复杂的乘除运算变成加减运算。
具体来说,开环传递函数通常是多个典型环节的乘积。在Bode图上,总的对数幅频特性就是各环节对数幅频特性的叠加,总相频特性也是各环节相频特性的叠加。这个性质让手绘Bode图变得非常高效。
另外,Bode图的横轴用对数刻度,可以同时展示很宽的频率范围。从0.01 rad/s到1000 rad/s,一张图就能覆盖五个数量级,这在奈奎斯特图上是做不到的。
4.2 典型环节的Bode图速记
画Bode图的基础是记住几个典型环节的形状。我整理了一个速查表:
| 环节 | 传递函数 | 幅频特性斜率 | 相频特性范围 |
|---|---|---|---|
| 比例 | K | 0 dB/dec | 0° |
| 积分 | 1/s | -20 dB/dec | -90° |
| 微分 | s | +20 dB/dec | +90° |
| 惯性 | 1/(1+Ts) | 低频0,高频-20 dB/dec | 0°到-90° |
| 一阶微分 | 1+Ts | 低频0,高频+20 dB/dec | 0°到+90° |
| 振荡 | 1/(1+2ζs/ωn+s²/ωn²) | 低频0,高频-40 dB/dec | 0°到-180° |
| 二阶微分 | 1+2ζs/ωn+s²/ωn² | 低频0,高频+40 dB/dec | 0°到+180° |
惯性环节的转折频率是1/T,在转折频率处幅值下降3dB,相角为-45°。振荡环节的转折频率是ωn,在转折频率处幅值有一个尖峰,尖峰高度取决于阻尼比ζ。ζ越小,尖峰越高。
提示:画Bode图时,低频段的斜率由积分环节的个数决定。每增加一个积分环节,低频斜率就多-20 dB/dec。高频段的斜率由分母和分子的阶数差决定,等于-(n-m)×20 dB/dec。
4.3 手绘Bode图的完整流程
我平时手绘Bode图有一套固定的流程,分享出来供参考:
把传递函数写成时间常数形式。每个环节都写成(1+Ts)或(1+2ζs/ωn+s²/ωn²)的标准形式。
确定各环节的转折频率,从小到大排列。这些转折频率把频率轴分成若干区间。
计算低频段幅值。取一个参考频率(比如ω=1),计算20lg|G(j1)H(j1)|,确定低频段的起点。
逐段画幅频特性。从低频开始,每经过一个转折频率,斜率就变化一次。惯性环节斜率减20,一阶微分斜率加20,振荡环节斜率减40,以此类推。
画相频特性。每个环节的相频特性单独画,然后叠加。惯性环节在转折频率前后各一个数量级范围内完成-90°的相移,转折频率处为-45°。
检查高频段斜率。最终的高频斜率应该等于-(n-m)×20 dB/dec,可以用来验证画得对不对。
这套流程看起来机械,但熟练之后画一张Bode图也就几分钟的事。关键是转折频率和斜率变化要记准。
4.4 从Bode图读取裕度
Bode图最大的好处就是可以直接读出相角裕度和幅值裕度,不需要额外计算。
读相角裕度:找到幅频特性穿越0dB线的频率ωc,然后看相频特性在ωc处的值,PM = 180° + 相角值。
读幅值裕度:找到相频特性穿越-180°线的频率ωg,然后看幅频特性在ωg处的值,GM = -幅值dB值。
在Bode图上,这两个裕度都表现为“距离”:PM是相频曲线在ωc处距离-180°线的垂直距离,GM是幅频曲线在ωg处距离0dB线的垂直距离。
我实际调控制器的时候,经常在Bode图上直接看这两个距离,然后决定是增大增益还是减小增益。增大增益会让幅频曲线上移,ωc右移,PM减小;减小增益则相反。
4.5 最小相位系统与非最小相位系统
Bode图有一个很重要的性质:对于最小相位系统,幅频特性和相频特性是一一对应的。也就是说,知道了幅频特性,相频特性就唯一确定了。这个性质让Bode图分析变得非常方便。
最小相位系统的定义是:所有零点和极点都在左半平面的系统。如果系统有右半平面零点或者右半平面极点,就是非最小相位系统。
非最小相位系统在Bode图上的表现是:幅频特性相同的两个系统,非最小相位系统的相角滞后更大。这意味着非最小相位系统的相角裕度更小,稳定性更差。
实际中常见的非最小相位环节包括:延迟环节e^(-τs)、右半平面零点(1-Ts)等。延迟环节在低频段近似为1,高频段相角持续下降,对稳定性影响很大。
注意:遇到非最小相位系统时,不能简单地用幅频特性去推断相频特性。必须单独计算相频特性,否则会得出错误的稳定性结论。这是考试和工程中都很容易犯的错误。
5. 常见问题与排查技巧实录
5.1 奈奎斯特判据的包围圈数怎么数
这是问得最多的问题。我的经验是:不要试图一次性看整条曲线,而是分段数。
具体做法:从ω=0+开始,沿着曲线走到ω=+∞,记录曲线从(-1, j0)点的哪一侧穿过负实轴。然后镜像得到ω从-∞到0-的部分。最后看整条闭合曲线绕(-1, j0)逆时针转了几圈。
更简单的方法是用“穿越次数”来判断:如果曲线在(-1, j0)左侧穿过负实轴,算一次穿越。正穿越(从上到下)和负穿越(从下到上)要区分。N = 正穿越次数 - 负穿越次数。
这个方法在考试中特别实用,因为不需要真的画出完整的闭合曲线,只需要看关键穿越点。
5.2 相角裕度为负意味着什么
相角裕度为负,说明在幅值穿越频率处,相角已经低于-180°了。这意味着奈奎斯特曲线已经包围了(-1, j0)点,系统不稳定。
在实际调试中,如果发现PM为负,第一反应应该是减小增益。减小增益会让幅频曲线上移,ωc左移,相角裕度增大。如果减小增益后PM仍然为负,说明系统结构有问题,可能需要加校正环节。
5.3 Bode图手绘时容易犯的错误
我总结了自己和身边人常犯的几个错误:
转折频率处斜率变化搞反。惯性环节是分母,斜率减20;一阶微分是分子,斜率加20。这个方向不能错。
相频特性叠加时忘记乘以环节个数。如果有两个相同的惯性环节,相频特性要叠加两次,总相移是-180°而不是-90°。
低频段幅值计算错误。低频段幅值要把所有环节在低频时的贡献都算进去,包括比例环节的K、积分环节的1/ω等。
忽略积分环节的初始斜率。有积分环节时,低频段斜率不是0而是-20k dB/dec,这个k是积分环节的个数。
5.4 裕度计算中的数值陷阱
计算裕度时,最容易出错的是反三角函数的计算。比如arctan(ω) + arctan(0.5ω) = 90°这个条件,利用的是arctan(A) + arctan(B) = 90°当且仅当AB=1。如果忘了这个性质,硬算会非常麻烦。
另外,在计算幅值穿越频率时,方程通常是高阶的,需要数值求解。我一般用迭代法:先估计一个ωc,代入计算幅值,然后根据结果调整。通常迭代三四次就能得到比较准确的值。
提示:实际工程中,如果手算太麻烦,可以用MATLAB的margin函数直接求裕度。命令是margin(G),会自动画出Bode图并标注PM和GM。但考试时还是得手算,所以基本方法必须掌握。
5.5 常见问题速查表
| 问题现象 | 可能原因 | 排查方法 | 解决措施 |
|---|---|---|---|
| 奈奎斯特曲线包围(-1,j0) | 开环增益过大 | 检查与实轴交点 | 减小增益 |
| PM为负 | 相角滞后过大 | 检查各环节相角贡献 | 减小增益或加超前校正 |
| GM过小 | 高频段幅值衰减不够 | 检查高频斜率 | 增加惯性环节或减小增益 |
| Bode图高频斜率不对 | 分子分母阶数差算错 | 重新数阶数 | 修正传递函数 |
| 相频特性与预期不符 | 有非最小相位环节 | 检查零极点位置 | 单独计算相频特性 |
5.6 几个实用的经验技巧
第一,画Bode图时先用渐近线画幅频特性,然后再在转折频率附近做修正。惯性环节在转折频率处实际值比渐近线低3dB,振荡环节的修正量取决于ζ。这个修正量在精度要求不高时可以忽略。
第二,判断稳定性时,先看PM再看GM。PM通常比GM更敏感,因为大多数系统的相角变化比幅值变化更快。如果PM满足要求,GM一般也不会太差。
第三,实际调控制器时,如果PM在30°到60°之间,GM在6dB到20dB之间,系统基本就是稳定的。如果超出这个范围,就需要仔细检查了。
第四,遇到高阶系统时,不要试图精确计算每个转折频率。抓住主导极点对,忽略远离虚轴的极点和零点,可以大大简化分析。这个近似在工程上完全够用。
第五,奈奎斯特判据和Bode图判据本质上是等价的。Bode图上的PM和GM就是奈奎斯特曲线距离(-1, j0)点的度量。理解了这一点,两种方法就可以互相验证。
6. 从理论到实践:频率法在真实项目中的应用
6.1 一个电机调速系统的实例
我之前做过一个直流电机调速系统的项目,被控对象是电机加负载,传递函数近似为G(s) = K / [(Ts+1)(Tms+1)],其中T是电气时间常数,Tm是机械时间常数。典型参数下,T=0.01s,Tm=0.1s,K=10。
先画Bode图。转折频率分别是1/T=100 rad/s和1/Tm=10 rad/s。低频段幅值为20lg(10)=20dB。经过10 rad/s转折后斜率变为-20 dB/dec,经过100 rad/s转折后斜率变为-40 dB/dec。
求幅值穿越频率:在10 rad/s附近,幅值约为20 - 20lg(ω/1)(近似)。令20 - 20lg(ωc) = 0,得ωc ≈ 10 rad/s。但实际因为有两个转折,需要更精确计算。
在ωc=10 rad/s处,|G| = 10 / [√(1+0.01×100) × √(1+0.1×10)] = 10 / [√2 × √2] = 10/2 = 5。所以ωc应该更大。
继续试ωc=20 rad/s:|G| = 10 / [√(1+0.01×400) × √(1+0.1×20)] = 10 / [√5 × √3] = 10/3.87 = 2.58。
ωc=30 rad/s:|G| = 10 / [√(1+9) × √(1+3)] = 10 / [3.16 × 2] = 1.58。
ωc=40 rad/s:|G| = 10 / [√(1+16) × √(1+4)] = 10 / [4.12 × 2.24] = 1.08。
ωc=42 rad/s:|G| ≈ 10 / [4.2 × 2.3] ≈ 1.03。
所以ωc ≈ 42 rad/s。
计算PM:∠G(j42) = -arctan(0.01×42) - arctan(0.1×42) = -arctan(0.42) - arctan(4.2) = -22.8° - 76.6° = -99.4°。
PM = 180° - 99.4° = 80.6°。
这个PM太大了,系统响应会很慢。实际调试时,我增大了增益K,让ωc右移,PM降到45°左右。具体来说,把K从10增大到30,重新计算ωc和PM,最终得到PM≈47°,系统响应明显加快,超调也在可接受范围内。
6.2 增益调整对裕度的影响规律
从上面的例子可以总结出一个规律:增大开环增益,幅频曲线上移,ωc右移,PM减小,GM也减小。减小增益则相反。
这个规律在调控制器时非常有用。如果系统响应太慢,就增大增益,但要注意PM不能太小。如果系统振荡太厉害,就减小增益,但响应会变慢。实际中需要在快速性和稳定性之间找平衡。
我一般先把增益调到PM≈45°,然后看响应曲线。如果超调还是太大,就稍微减小增益;如果响应太慢,就稍微增大增益。通常调整两三次就能找到合适的值。
6.3 校正环节的引入时机
如果单纯调增益无法同时满足快速性和稳定性的要求,就需要加校正环节。常用的校正方式有超前校正、滞后校正和超前滞后校正。
超前校正的作用是在ωc附近提供正相角,增大PM,改善动态性能。滞后校正的作用是降低高频增益,减小ωc,提高稳态精度。超前滞后校正则结合两者优点。
选择校正方式时,先看系统缺什么。如果PM不够但稳态精度可以,就用超前校正。如果稳态精度不够但PM可以,就用滞后校正。如果两者都不够,就用超前滞后校正。
校正环节的设计方法这里不展开,但核心思想是利用Bode图的叠加原理,把校正环节的频率特性叠加到原系统上,使总的频率特性满足PM和GM的要求。
6.4 数字控制中的频率法
现在很多控制器都是数字实现的,频率法同样适用,但需要注意采样频率的影响。数字控制系统的频率特性在奈奎斯特频率(采样频率的一半)附近会有明显变化,高频段的幅值和相角都会受到采样和保持环节的影响。
实际设计中,采样频率一般取系统带宽的10倍以上,这样在关心的频率范围内,数字控制器的频率特性可以近似为连续控制器。如果采样频率不够高,就需要用z变换或者双线性变换来精确分析。
我在做数字控制器时,通常先用连续域的方法设计,然后用双线性变换离散化,最后再验证离散后的频率特性。如果离散化后PM和GM变化不大,就可以直接用。如果变化较大,就需要重新调整参数。
7. 学习这一章的个人体会
这一章的内容我前前后后学了三四遍,每次都有新的理解。第一遍学的时候,满脑子都是公式和曲线,根本不知道这些东西有什么用。第二遍学的时候,开始能画出Bode图,能算出PM和GM,但还是停留在做题层面。第三遍学的时候,正好在做一个实际项目,突然就明白了这些概念背后的工程意义。
我觉得学这一章最重要的不是记住多少公式,而是建立起“频率域”的思维方式。时域分析看的是信号随时间的变化,频率域分析看的是系统对不同频率信号的响应能力。这两种视角是互补的,很多在时域里看不清楚的问题,换到频率域就一目了然。
另外,Bode图一定要多画。我当年考试前手绘了二十多张Bode图,各种类型的传递函数都画过。画多了之后,看到传递函数就能在脑子里浮现出大致的Bode图形状,这对快速判断稳定性非常有帮助。
最后说一个我踩过的坑:有一次做项目,系统在实验室调试时PM有50°,很稳定。但装到现场后,电机电缆加长了十几米,引入了额外的分布电容和电感,相当于增加了一个小惯性环节。这个环节的转折频率刚好在ωc附近,把PM吃掉了将近20°,系统就开始振荡了。后来在控制器里加了一个超前校正,把PM补回来才解决问题。这个经历让我深刻体会到,裕度不只是纸面上的数字,它是实实在在的工程安全边际。