控制理论稳定性分析:判据选型与解题实战
1. 先搞清楚“稳定性分析”这道大题在考什么
先说结论:控制理论里稳定性分析这块内容,几乎每年必出,而且出题位置很固定——要么是第3大题(10到15分),要么是最后一题里的一小问(8到10分)。题面形式千变万化,核心考的东西其实就那几样。很多人复习这块最大的问题是判据背了一堆,拿到题目不知道用哪个,或者用了劳斯判据算了一半发现根本算不下去。我见过的绝大多数失分,都不是因为不会算,而是选错了工具。
稳定性这个概念在教材上有三层意思,考试只考其中一层半。第一层是有界输入有界输出的外部稳定性,讲的是传递函数极点的位置;第二层是零输入响应的渐近稳定性,讲的是系统矩阵的特征值;第三层是李雅普诺夫意义下的稳定性,讲的是能量函数。考大题的时候,90%的题目落在第一层,偶尔在状态空间那一章落到第二层,第三层一般只出小题或者干脆不出。你如果一上来就想用李雅普诺夫第二法去解一个传递函数的题,那纯属绕远路。
真正需要你掌握的,是五套工具:劳斯判据、胡尔维茨判据、奈奎斯特判据、伯德图判据、根轨迹法,状态空间那一章再补一个李雅普诺夫直接法。这六样东西看起来五花八门,其实分属三个完全不同的体系——代数判据、频域判据、时域判据。它们之间的换算关系是有路径的,而且路径不止一条。
举个最常见的场景:题目给一个开环传递函数,带一个待定参数K,问“求使闭环系统稳定的K的取值范围”。这种题,你用劳斯判据能做,用奈奎斯特判据也能做,用根轨迹与虚轴交点也能做,三种方法算出来的K范围应该完全一致。考试的时候选哪个?看题目的附加条件。它要是还问了“稳定裕度”或者“穿越频率”,那你必须先走频域;它要是还问了“分离点”或者“主导极点”,那你先走根轨迹;都没问,纯问范围,那就劳斯判据,十分钟内收工。
这里得插一句为什么三个体系的结果必须等价。因为闭环系统的稳定性只由闭环特征方程决定,而开环传递函数和闭环特征方程之间是一条确定的代数关系:1+G(s)H(s)=0。无论是劳斯表、奈氏曲线还是根轨迹,本质上都是在研究这个方程根的位置,只是从不同角度观察同一个东西。这个认知很重要,因为它意味着你可以用两套方法互验,考试的时候省下大量检查时间。我在考场上养成的习惯是:劳斯判据算出临界参数后,立刻用奈奎斯特在脑子里过一遍,看看穿越点是不是在那个位置,两秒钟就能确认,比重新算一遍靠谱得多。
还有一个隐藏的考点,很多人复习的时候会漏掉:离散系统的稳定性。z平面的单位圆和s平面的虚轴之间的映射关系,是这块的命门。考试如果出了离散题,一般会让你判断特征根在不在单位圆内,或者让你做双线性变换把z域问题转成s域再套劳斯判据。这个方法叫w变换,w=(z-1)/(z+1),反变换是z=(1+w)/(1-w)。记住这对公式,离散题的难度会瞬间降一个档次。
2. 工具选型:五套判据怎么选才不浪费时间
2.1 判据功能对照表
这张表我先给你,省得你翻书。你复习到后期,脑子里应该是这张表,而不是六个孤立的公式。
| 判据 | 输入信息 | 输出结论 | 计算量 | 典型适用场景 |
|---|---|---|---|---|
| 劳斯判据 | 闭环特征方程系数 | 右半平面根数、临界参数 | 中(解不等式) | 参数范围题、纯代数题 |
| 胡尔维茨判据 | 特征方程系数 | 稳定性判定 | 中(行列式计算) | 低阶系统(n≤3) |
| 奈奎斯特 | 开环频率特性 | 闭环稳定性、圈数 | 高(画图) | 开环含延迟、非最小相位 |
| 伯德图 | 开环伯德图 | 稳定性、稳定裕度 | 低(对数坐标读图) | 给频率特性的题、裕度题 |
| 根轨迹 | 开环零极点 | 参数变化时根的运动 | 中(画草图) | 分离点、虚轴交点、主导极点 |
| 李雅普诺夫 | 状态方程、能量函数 | 渐近稳定性 | 低(求导判号) | 状态空间章、非线性局部线性化 |
这张表里最关键的一列是“输入信息”。你拿到题目第一件事,不是想“我该背哪个判据”,而是看题目给了什么。给了闭环特征方程,那必然劳斯;给了开环传递函数和K,三条路都通,再看具体问法;给了伯德图或者让你画伯德图,那只能是伯德判据;给了状态方程,李雅普诺夫。
2.2 选型背后的一笔账
有人问为什么不用胡尔维茨判据,它看起来更“正规”。原因很简单:胡尔维茨要算行列式,劳斯表只是做加减乘除。三阶系统两者差不多,四阶以上胡尔维茨的行列式展开就非常痛苦了,而劳斯表每多一行只是多一次除法,复杂度是线性的。考试时间宝贵,劳斯判据是代数判据里的默认选项,胡尔维茨只在一种情况下更优——题目本来就把行列式写出来了,你只要判断符号,那就直接用,别折腾。
奈奎斯特判据和伯德判据的关系也值得说清楚。它们本质上是同一件事,伯德图是奈氏曲线在对数坐标下的展开。区别在于:奈氏图适合看大范围的圈数,尤其是开环右半平面极点数P不为零的场合;伯德图适合看局部细节,因为对数坐标把低频段压得很扁、把高频段拉得很长,穿越点附近的形状一目了然。所以遇到含延迟环节e^(-τs)的系统,我一般直接上伯德图,因为延迟只改相频特性不改幅频特性,在伯德图上就是相角往下多掉一截,判断起来直观得不行;而在奈氏图上,延迟会让曲线变成一个无限延伸的螺旋,画都画不完。
根轨迹的适用边界是个反例。它只能研究单参数变化,一般是增益K从0到无穷。如果题目里有两个参数同时变,根轨迹就废了,必须换频域或者代数判据。这个坑我踩过,某次做一道双参数题,硬画根轨迹画了二十分钟,最后发现根本走不通,用劳斯判据两行就搞定了。所以说工具选型这件事,本质上是对题目结构的一次快速研判,花三十秒做这个研判,能省下二十分钟。
2.3 一个必须刻进肌肉记忆的换算
从开环判断闭环,靠的是这个式子:Z = P - 2N。这里的P是开环传递函数在右半s平面的极点数,N是奈氏曲线包围(-1, j0)点的圈数(逆时针为正,顺时针为负),Z是闭环右半平面的极点数。系统稳定要求Z=0,也就是N = P/2。
大多数题目里开环是稳定的,P=0,这时候退化成一句话:曲线不能包围(-1, j0)。这个简化条件覆盖了八成以上的题目,但它是个陷阱——剩下那一成多,也就是开环本身不稳定(比如含有右半平面的极点,常见于某些反馈结构),你就必须老老实实回去数P。考试里出这种题的意图很明确,就是筛选出那些只会背简化版的人。
3. 劳斯判据:把步骤钉死,特殊情况单独记
3.1 标准流程与一个完整算例
劳斯判据的步骤我按自己的习惯写一遍,你对照着看有没有漏。第一步,写出闭环特征方程,一定要整理成降幂排列、首项系数为正的标准形式。这一步不整理好,后面全错。第二步,如果最高次项系数不是正的,两边同乘-1,注意乘完之后所有不等号要反过来。第三步,如果存在缺项(某次幂系数为0),不用特殊处理,直接在那一位写0往下排。第四步,列劳斯表。第五步,数首列符号变化次数,那就是右半平面的根数。
拿个具体例子走一遍,题目:单位负反馈系统开环传递函数为 G(s) = K / [s(s+1)(s+2)],求闭环稳定的K范围。
闭环特征方程是 1 + G(s) = 0,即 s(s+1)(s+2) + K = 0。展开:s³ + 3s² + 2s + K = 0。首项系数已经是1,可以直接排表。
s³ | 1 2 s² | 3 K s¹ | (6-K)/3 0 s⁰ | K第三行那个 (6-K)/3 是这样来的:用 s² 行的第一项 3 去除交叉相乘的差,(3×2 - 1×K)/3 = (6-K)/3。这是标准的交叉除法,别看错位置。
稳定条件是首列全部同号且不为零。首列是 1、3、(6-K)/3、K。前两个已经为正,所以要求:
(6-K)/3 > 0,即 K < 6; K > 0。
合起来0 < K < 6。
临界值K=6的时候,s¹ 整行为零,说明有根落在虚轴上。这时候用 s² 行构造辅助方程:3s² + 6 = 0,解出 s = ±j√2。也就是说临界状态下系统有一对纯虚根,振荡频率是 √2 rad/s。这个信息在后续问“临界频率”或者“等幅振荡周期”的时候直接拿来用,周期 T = 2π/ω = 2π/√2 ≈ 4.44秒。
这套流程练熟之后,四阶以内的系统三分钟内能出结果。但真正的难点不在流程,在两个特殊情况和参数题的变体上。
3.2 首列出现零怎么办
第一种特殊情况:某一行首列为零,但整行不全为零。这时候交叉除法会除以零,表就排不下去了。标准做法是用一个无穷小正数 ε 代替那个零,继续往下排,最后令 ε→0 取极限判断符号。
但我要说个更实用的技巧:先试试能不能用 (s+a) 因子替换法。具体做法是把原特征方程里的 s 换成 1/x(或者说把多项式系数倒过来写),重新排一次劳斯表。因为“首列为零”意味着有纯虚根或者有对称于原点的根,倒数变换之后这个障碍通常会消失。这个方法比引入 ε 干净,尤其在手算的时候不容易在极限上出错。我用 ε 法错过两次符号,后来就一直用倒数变换。
3.3 整行为零怎么处理
第二种特殊情况:某一行整行全为零。这说明存在关于原点对称的根,可能是一对纯虚根,也可能是一对实根±σ,还可能两者都有。处理方法是用上一行的系数构造辅助方程 A(s)=0,对 s 求导得到 A'(s),用导数的系数替换那一整行,继续排表。
这里有个细节很多人不知道:辅助方程的阶数一定是偶数,而且它的根就是那些对称根。所以做完劳斯表之后,如果题目还问“临界状态的振荡频率”,你直接从辅助方程解就行了,不用再回头验算。辅助方程 2s² + 8 = 0 就意味着 s = ±j2,频率是 2。
还有一句提醒:整行为零的情况下,劳斯判据给出的“首列符号变化次数”只统计非零行,而且结论仍然成立,但你不能直接说“系统稳定”——因为纯虚根的存在意味着临界稳定,不是渐近稳定。考试问“是否稳定”,你必须回答“临界稳定”或者“不稳定”,答“稳定”是要扣分的。
3.4 参数题的三种常见变体
参数题的变体我归了三类。第一类是单参数线性出现,就像上面那个K,直接解不等式组就行。第二类是参数出现在两行以上,这时候要小心,因为不同行给出的不等式可能有交集为空的情况,得逐行列全再取交集,漏一行就前功尽弃。第三类是要求“稳定且满足某项性能指标”,比如“所有闭环极点的实部小于-1”,这时候要做坐标平移,把 s 换成 s-1 再排劳斯表。这个平移操作叫“稳定度”检验,考的频次不高,但一旦出了就是区分度题。
平移的方向别搞反。要求 Re(s) < -1,就是让新变量 s' = s+1 满足 Re(s') < 0,所以代入的是 s = s' - 1。我见过有人代成 s = s' + 1,结论完全反过来。记住一句话:你想让根往左移,就把 s 替换成 s 减掉那个数,因为 s = s' - 1 相当于把原点右移一格,原来的左半平面就变成了更靠左的区域。这句话有点绕,但比死记公式管用。
4. 奈奎斯特与伯德图:频域判据的实操细节
4.1 奈氏判据的计数规则
奈奎斯特判据的公式 Z = P - 2N 我已经给过了,难的是 N 怎么数。N 的定义是奈氏曲线绕(-1, j0)点逆时针转的圈数。手画的时候有几个约定必须遵守。
第一,正穿越和负穿越要分清。曲线从上往下穿过负实轴的 (-∞, -1) 段,叫正穿越,计+1;从下往上穿过,叫负穿越,计-1。半次穿越(曲线正好终止在(-1,j0)上)计±0.5。N 是正穿越减去负穿越。
第二,ω 从 0 到 +∞ 画出来的只是半个奈氏曲线,还得补上 ω 从 -∞ 到 0 的镜像部分,关于实轴对称。很多人只画半条就下结论,圈数直接少一半。
第三,开环含有积分环节(也就是原点处有极点)的时候,要做补圆。补圆的半径是无穷大,方向是顺时针,转过的角度是 ν×180°,ν 是积分环节个数。补圆的方向是最容易搞错的地方,我的记忆方法是:极点用顺时针补,零点用逆时针补,因为极点在分母上,它在无穷远处的行为跟零点相反。
第四,也是最实用的一条:如果开环没有右半平面极点(P=0),那就根本不用数圈,只看曲线有没有把(-1, j0)圈进去。圈进去就不稳定,没圈进去就稳定,穿过去就是临界。这一条能省掉大量纠结。
4.2 伯德图判据的读法
伯德图判据是对数坐标下的奈氏判据,判断规则是:在幅频特性 L(ω) 穿越 0dB 的那个频率(叫穿越频率 ωc)上,看相频特性 φ(ωc) 有没有低于 -180°。低于就说明闭环不稳定,高于就稳定,正好等于就是临界。
这个规则的逻辑其实很直白。0dB 处对应 |G(jω)|=1,-180° 对应相角反相。如果在这个频率上相位已经翻过去了,那闭环就是正反馈,越振越大。伯德图的妙处在于它把这个判断变成了“看两条线在哪个位置交叉”,一眼就能看出来。
我一般会在图上标三个点:L=0dB 处的频率 ωc、φ=-180° 处的频率 ωg、以及两条曲线各自的数值。有了这三个点,稳定裕度就直接读出来了。
4.3 稳定裕度的计算与算例
相角裕度 γ = 180° + φ(ωc),幅值裕度 Kg = -L(ωg),单位是 dB;或者写成线性形式 Kg = 1/|G(jωg)|。两个都要会写,因为题目问法不一样。
拿上面那个例子继续,取 K=2,即 G(s) = 2/[s(s+1)(s+2)]。先求穿越频率 ωc,令 |G(jω)| = 1:
|G| = 2 / (ω·√(ω²+1)·√(ω²+4)) = 1
这个方程直接解不出来,得试算。ω=1 时,分母是 1×√2×√5 = 3.162,|G|=0.632,太小;ω=0.7 时,分母是 0.7×√1.49×√4.49 = 0.7×1.221×2.119 = 1.811,|G|=1.104,偏大;ω=0.75 时,分母是 0.75×1.25×2.136 = 2.002,|G|≈0.999,就它了。所以ωc ≈ 0.75 rad/s。
再算这个频率上的相角:φ(ωc) = -90° - arctan(0.75) - arctan(0.375)。arctan0.75 ≈ 36.9°,arctan0.375 ≈ 20.6°,加起来 147.5°,所以 φ = -147.5°。
相角裕度 γ = 180° - 147.5° = 32.5°。
再求相位穿越频率 ωg,令 arctanω + arctan(ω/2) = 90°。两个反正切加起来等于90°,说明它们的正切互为倒数,即 ω × (ω/2) = 1,解得 ω² = 2,ωg = √2 ≈ 1.414 rad/s。
这时候 |G(jωg)| = 2 / (√2 × √3 × √6) = 2/√36 = 2/6 = 0.333。幅值裕度 Kg = 1/0.333 = 3,换算成 20lg3 ≈ 9.5dB。
这两个数字和前面劳斯判据算出来的临界K=6是对得上的。K=6 的时候,|G(jωg)| = 6/6 = 1,也就是 0dB,幅值裕度刚好归零,系统临界。这个互验关系你可以在考场上用,算完裕度回头看一眼临界参数对不对得上,对不上就是哪一步错了。
4.4 含延迟环节的处理
带 e^(-τs) 的系统是频域判据的主场。延迟环节的幅值是1,不改变幅频特性,只给相频特性加一项 -ωτ(弧度)。所以在伯德图上,幅频曲线完全不动,相频曲线往下多掉一截,而且掉得越来越快。
判断方法非常直接:先按没有延迟的情况算出相角裕度 γ0,再算延迟在 ωc 处造成的额外相移 Δφ = ωc·τ(注意单位,τ 用秒,ωc 用 rad/s,结果是弧度,要换算成角度乘 57.3)。新的裕度是 γ = γ0 - Δφ。裕度变成负值,系统就不稳定。
临界延迟时间也好求,令 γ0 - ωc·τ = 0,得 τ_max = γ0/ωc(γ0 用弧度)。这个公式考试里出过一次,属于送分题,但前提是你得知道延迟只影响相位这件事。
5. 根轨迹:从画图到判稳的一条龙
5.1 八条法则里真正要用的四条
根轨迹的绘制法则教材上列了七八条,考试时间不够,我实际只重点用四条,其余靠草图感觉补。
第一条,分支数与起止点。分支数等于开环极点数 n,从开环极点出发,终止于开环零点或无穷远。这条决定了图形的骨架。
第二条,实轴上的线段。实轴上某一段属于根轨迹的充要条件是:它右侧的实数开环零极点总数为奇数。这条法则判断起来极快,从右往左扫一遍就行,不用背公式。
第三条,渐近线。n-m 条渐近线,角度是 (2k+1)×180°/(n-m),交点是 (Σ极点 - Σ零点)/(n-m)。这两条公式在 n-m≥2 的时候必须用,否则草图会歪得离谱。
第四条,分离点与虚轴交点。分离点解方程 Σ1/(s-pᵢ) = Σ1/(s-zⱼ),虚轴交点要么令 s=jω 代入特征方程分离实虚部,要么直接拿劳斯判据的临界条件。后者更省事,我一般直接用劳斯。
出射角和入射角这两条,只在题目明确问了才用,而且通常是给复数极点的情况,形式是 180° 减去其他零极点造成的相角。考的频次不高,但如果出了,很多人都忘光了。
5.2 用根轨迹判稳的标准动作
题目一般这么问:给开环传递函数,问“K 从 0 变到无穷时,闭环系统的稳定性如何变化”。标准动作是:先写特征方程,用劳斯判据求出临界 K 和对应的虚轴交点坐标,再用分离点判断根怎么跑。
拿前面那个例子继续,G(s)=K/[s(s+1)(s+2)]。开环极点 0、-1、-2,没有零点,n=3,m=0。实轴上 (-1, 0) 和 (-∞, -2) 两段属于根轨迹。渐近线三条,角度分别是 60°、180°、300°,交点在 (0-1-2)/3 = -1。
分离点:1/s + 1/(s+1) + 1/(s+2) = 0。通分整理:
[(s+1)(s+2) + s(s+2) + s(s+1)] / [s(s+1)(s+2)] = 0 分子 = (s²+3s+2) + (s²+2s) + (s²+s) = 3s² + 6s + 2 = 0 s = [-6 ± √(36-24)]/6 = (-6 ± 3.464)/6
得 s₁ = -0.423,s₂ = -1.577。s₂ 不在根轨迹上(因为 (-1.577) 落在 (-2,-1) 区间,右侧零极点数为偶数),舍去。所以分离点是s = -0.423。
代入特征方程求分离点对应的 K:K = -s(s+1)(s+2),s=-0.423 时,K = -(-0.423)(0.577)(1.577) = 0.385。
虚轴交点前面已经算过,K=6,s=±j1.414。
把这三个信息串起来,稳定性结论就出来了:0 < K < 6 时全部根在左半平面,稳定;K=6 临界;K>6 有一对共轭根进入右半平面,不稳定。而且 K>0.385 之后根就从实轴分开了,变成复数,意味着阶跃响应开始出现振荡。这个“从单调到振荡”的转折点,很多题目会专门问,答案就是分离点对应的 K 值。
5.3 主导极点近似的边界
根轨迹经常和“主导极点”一起考,问“当 K 取某值时,系统能否用二阶近似”。判断标准是:一对共轭复极点距离虚轴的距离,要小于其他所有极点到虚轴距离的 1/5。满足这个条件,高阶系统才能降成二阶来估性能指标。
这条规则教材上写得很明白,但用起来有个坑:它是充分不必要条件。不满足的时候,不意味着不能近似,只是误差会大。考试如果问“能否近似”,你按 1/5 判;如果问“近似后的超调量是多少”,那就算完再补一句“实际值会略有偏差”。我吃过这个亏,当时判了不能近似就直接空着,实际上稍微变通一下是可以给出估算的。
6. 离散系统与状态空间:两块容易被忽略的补充
6.1 朱里判据与双线性变换
离散系统的稳定性判据是“特征根全部在单位圆内”。直接解高次方程不现实,考试一般走两条路。
第一条是朱里判据,用特征多项式的系数列一张表,逐行判断。这个表的形式比劳斯表复杂,而且不同教材的写法略有差异,容易记混。我的建议是:能不用就不用。
第二条是双线性变换,把 z 域问题搬到 s 域,然后直接用劳斯判据。变换公式 w = (z-1)/(z+1),反变换 z = (1+w)/(1-w)。这个变换把单位圆内部映射成 w 左半平面,单位圆本身映射成虚轴。所以做完变换之后,判稳就是判 w 特征方程的根是不是全在左半平面。
举个两阶的例子:特征方程 z² + az + b = 0,做变换:
[(1+w)/(1-w)]² + a[(1+w)/(1-w)] + b = 0
两边乘 (1-w)²:
(1+w)² + a(1+w)(1-w) + b(1-w)² = 0
展开整理成 w 的多项式,再排劳斯表。这个过程计算量不小,但步骤机械,不容易出错,比硬记朱里表靠谱。三阶以内我基本都这么做。
有一个细节必须注意:变换后特征方程的系数可能含分母,要先通分整理成标准形式再排表。我见过直接拿没通分的系数去排劳斯表,结果完全错的。
6.2 李雅普诺夫直接法在状态空间题里的用法
状态空间那一章的稳定性题,一般给一个 2×2 的状态矩阵 A,问系统是否渐近稳定。最直接的方法是算特征值,λ² - tr(A)λ + det(A) = 0,二阶系统有现成的判别条件:tr(A) < 0 且 det(A) > 0。这两个条件等价于两个特征值都有负实部,用起来比展开算快得多。
要是题目明确要求“用李雅普诺夫第二法判断”,那就必须构造 V(x)。线性系统的标准套路是设 V(x) = xᵀPx,P 正定,然后解李雅普诺夫方程 AᵀP + PA = -Q,Q 取单位阵最方便。解出 P 之后验证它是不是正定(顺序主子式全大于零)。
二阶的情况可以手算。设 A = [[a, b], [c, d]],P = [[p11, p12], [p12, p22]],代入方程会得到三个线性方程,解出来三个未知数。这个过程大概需要五六分钟,不算快,但拿分稳。
非线性系统那一块,考的是“局部线性化后判断”。步骤是:先找平衡点,在平衡点处求雅可比矩阵,用雅可比矩阵代替 A 判断特征值。注意结论的适用范围——特征值全有负实部,只能说明局部渐近稳定,不能推全局。这个限定条件考试改卷是看关键词的,漏了要扣分。
7. 常见问题与排查技巧实录
7.1 易错点速查表
下面这张表是我自己做题和帮别人看题时总结出来的高频错误,按出现频率从高到低排。
| 错误现象 | 根本原因 | 正确做法 |
|---|---|---|
| 劳斯表首列符号变化数算错 | 忘记把首项系数调成正的 | 排表前先检查最高次项系数 |
| 参数范围开区间闭区间搞混 | 没区分“稳定”和“临界” | 严格不等号对应稳定,等号单独讨论 |
| 奈氏判据圈数少一半 | 只画了 ω>0 的半条曲线 | 补上关于实轴的镜像 |
| 补圆方向反了 | 没搞清极点和零点的区别 | 极点顺时针,零点逆时针 |
| 伯德图裕度正负号弄反 | 幅值裕度定义没记牢 | Kg = -L(ωg),负的L加个负号变成正 |
| 根轨迹分离点解出两个后不筛选 | 忘记验证是否在根轨迹上 | 检查该点右侧零极点总数是否为奇 |
| 相角计算时 arctan 象限错误 | 用了 arctan 而没用 atan2 | 手工判断所在象限后加减180° |
| 离散题变换后没通分 | 直接拿含分母的系数排表 | 先整理成整系数标准形式 |
最后一条错误我想多说两句。相角的象限问题在频域题里非常致命,因为 arctan 的返回值只在 (-90°, 90°) 之间,而真实的相角可能落在任何象限。比如一个环节 1/(1-jω),它的相角是 +arctan(ω),而不是 -arctan(ω)。区别就在分母的实部是正的还是负的。
避免这个错误的方法很简单:每一段频率特性单独判断,先看它是“1+jω”型还是“1-jω”型,再决定符号。我习惯在草稿纸上把每个环节的相角表达式写全,标上正负号,最后再加总,这样即使中途出错也容易查。
7.2 考场时间分配的实战建议
稳定性大题一般 12 到 15 分,我给的时间上限是 18 分钟。分配是这样的:读题和选判据 2 分钟,主计算 10 分钟,验算 3 分钟,写过程 3 分钟。
选判据那 2 分钟千万不能省。我见过太多人上来就排劳斯表,排到一半发现特征方程根本写不出来(因为给的是结构图不是传递函数),只能推倒重来。正确的顺序是:先看清楚题目给的是什么——结构图就先化简求传递函数,传递函数就直接写特征方程,伯德图就直接读数据。
验算的 3 分钟也有讲究。验算不是重新算一遍,而是用另一套方法交叉验证关键结论。比如劳斯算出临界 K=6,那就用奈奎斯特看一眼 ωg 处的幅值是不是 1;算出分离点,就用根轨迹的渐近线检查一下大致位置对不对。整个过程两三个心算就够,比重算快得多。
写过程的 3 分钟同样不能省。这类题是过程分大户,写出特征方程给分,列出劳斯表给分,写出稳定性条件给分,最后结论再给分。我见过结论对但过程跳步被扣一半分的,很亏。反过来,即使最后结论算错了,前面步骤写全也能拿到六成分。
7.3 几条不太见于教材的经验
第一,临界状态一定要单独写一句。题目问“求稳定范围”,你答 (0, 6) 是完整的,但最好补一句“K=6 时系统临界稳定,产生频率为 1.414 rad/s 的等幅振荡”。这句话不占多少篇幅,但展示了你对边界的理解,改卷老师看到会给印象分。
第二,能用表格的地方别用文字。劳斯表、稳定性条件的不等式组、不同判据的对比,写成表格清晰度提升一大截。尤其是多个不等式取交集的时候,画个数轴标出区间,比写一堆“且”字清楚得多。
第三,遇到算不出的方程先估个数。频域题里的穿越频率往往是个超越方程,强行解析求解会耗掉大量时间。我的做法是先估一个近似值(比如 0.75),代进去验证精度,够用就停。考试评分对穿越频率的精度要求一般是两位有效数字,你给 0.75 和给 0.748 得分一样。
第四,题目的最后一句往往是提示。有些题目开头看起来是纯代数题,最后一句却问“系统的相角裕度是多少”,这就明摆着告诉你应该走频域。反过来,如果最后一句问“分离点坐标”,那就走根轨迹。养成先扫一眼问题再动笔的习惯,能省不少弯路。
第五,不要把不同判据的结论混着用。劳斯判据给出的是“右半平面根数”,奈奎斯特给出的是“闭环稳定性”,根轨迹给出的是“根的运动轨迹”。这三样东西描述的是同一件事的不同侧面,但表述方式不能互换。答题的时候用哪个判据,就用那个判据的语言体系来写结论,别串。
我自己在这块栽得最深的一次,是一道四阶系统的参数题。当时图快,直接列了劳斯表,但特征方程是从结构图化简来的,中间有个负反馈的符号我漏了,导致整个方程差一个正负号。排出来的表首列符号变化了两次,我还以为系统不稳定,写了半天分析,最后发现是方程写错了。从那以后我养成了一个习惯:特征方程写完之后,先代 s=0 验证一下常数项,再代 s=1 验证一下各项系数和。这两个数值验算加起来不到二十秒,但能拦住绝大多数低级错误。常数项在 s=0 处等于所有极点的乘积,如果你能从传递函数直接看出极点,就能对得上;各项系数和在 s=1 处等于各极点代入后的乘积,虽然不那么直观,但至少能发现符号错误。
这套方法用到后来,稳定性分析这块基本就成了固定套路:看题、选工具、算临界、写区间、验算、补一句临界状态。整个过程行云流水,剩下的时间可以留给更难的题。要说这块内容有什么真正的门槛,不在计算,在于你得知道每个判据的边界在哪里——什么时候它好用,什么时候它会骗你,这个判断力只能靠做题积累,没有捷径。