4. 根轨迹与虚轴的交点：临界稳定点的确定
做控制系统设计这么多年，我有个很深的体会——系统稳不稳定，往往就在虚轴那根线上。根轨迹一旦跨过虚轴进入右半平面，系统就失控了。所以找到根轨迹与虚轴的交点，说白了就是找到系统从稳定到不稳定的那个临界点。
这一节，咱们就专门聊聊怎么确定这个临界稳定点。我习惯用两种方法：一种是直接求解，另一种是用劳斯判据辅助判断。两种方法各有千秋，咱们一个一个说。
4.1 临界稳定点的物理意义
先想一个问题：根轨迹上的点，什么时候会落在虚轴上？
嗯，其实很简单。当系统的闭环极点落在虚轴上时，系统处于等幅振荡状态。这时候既不发散也不收敛，就是临界稳定。你想想看，这对工程应用来说有多重要——我们既希望系统有足够的稳定裕度，又想知道它能承受的最大增益是多少。
我在做电机伺服系统时遇到过这种情况：增益调得太高，电机就开始抖。一测，根轨迹正好穿过了虚轴。所以找到这个交点，就能确定增益的极限值。
4.2 方法一：直接代入法求解
这个方法最直接。既然根轨迹上的点满足特征方程，那虚轴上的点有什么特点？
虚轴上的点，实部为0。所以我们可以设 s = jω，代入特征方程，然后分别令实部和虚部等于0，解出ω和K。
来看个例子。假设开环传递函数为：
G(s)H(s) = K / [s(s+1)(s+2)]闭环特征方程是：
1 + G(s)H(s) = 0
s(s+1)(s+2) + K = 0
s³ + 3s² + 2s + K = 0令 s = jω：
(jω)³ + 3(jω)² + 2(jω) + K = 0
-jω³ - 3ω² + 2jω + K = 0整理实部和虚部：
实部：K - 3ω² = 0
虚部：-ω³ + 2ω = 0从虚部方程解出：ω(-ω² + 2) = 0，ω = 0 或 ω = ±√2
ω=0对应根轨迹起点，我们关心的是ω=√2的情况。代入实部：
K = 3ω² = 3 × 2 = 6所以交点坐标是 s = j√2，临界增益 Kc = 6。
4.3 方法二：劳斯判据辅助法
直接代入法虽然直观，但遇到高阶系统时计算量很大。这时候我更喜欢用劳斯判据。
劳斯判据能告诉我们系统稳定时K的取值范围。当系统处于临界稳定时，劳斯表会出现全零行。利用这个特性，我们可以反推出临界增益和振荡频率。
还是用刚才的例子：
特征方程：s³ + 3s² + 2s + K = 0列劳斯表：
| s³ | 1 | 2 |
|---|---|---|
| s² | 3 | K |
| s¹ | (6-K)/3 | 0 |
| s⁰ | K |
系统稳定的条件是第一列元素全部大于0：
3 > 0（成立）
(6-K)/3 > 0 → K < 6
K > 0所以稳定范围是 0 < K < 6。临界增益 Kc = 6。
那振荡频率怎么求？当K=6时，s²行对应的辅助方程为：
3s² + K = 0 → 3s² + 6 = 0 → s² = -2 → s = ±j√2看，和直接代入法结果完全一致。
4.4 两种方法的对比
我平时怎么选？给你个参考：
- 系统阶次低（3阶以下）：直接代入法更快，计算量不大
- 系统阶次高（4阶及以上）：用劳斯判据，避免解高次方程
- 需要同时求多个参数：劳斯判据更系统化
不过说实话，我更喜欢两种方法互相验证。有一次做项目，直接代入法算出来一个值，劳斯判据算出来另一个值，一检查发现是代入时符号写错了。所以交叉验证是个好习惯。
4.5 知识体系总览
下面这张图把本节的核心逻辑串起来了，你可以对照着看：
4.6 避坑指南
最后分享几个我踩过的坑：
- 符号别搞反：代入 jω 时，注意 j² = -1，j³ = -j。我见过有人把符号弄反，结果算出来的K是负的。
- ω=0的解要排除：虚部方程解出来通常有ω=0，这是根轨迹的起点，不是我们要的交点。
- 劳斯表要仔细：列劳斯表时系数别抄错，尤其是高阶系统。我曾经因为漏了一个系数，多花了半天时间排查。
- 辅助方程别忘：用劳斯判据求出K后，一定要用辅助方程求ω。很多人只求了K就以为完事了。
好了，关于根轨迹与虚轴的交点，核心就是这些。记住：临界稳定点不是设计目标，而是设计边界。知道边界在哪，你才能放心地在安全区域内调参数。
公众号:蓝海资料掘金营，微信deep3321