1. 项目概述：为什么根轨迹是控制工程师的“导航图”？
搞自动控制，无论是做机器人、无人机，还是工业过程控制，最终都要落到一个核心问题上：我这个系统稳不稳？响应快不快？精度够不够？这些问题，本质上都取决于系统闭环传递函数的极点，也就是特征方程的根，在复平面上的位置。但一个稍微复杂点的系统，特征方程可能是个高阶多项式，参数一变，根的位置就跟着变，手动算起来简直是噩梦。这时候，“根轨迹”就登场了。它不是什么高深莫测的理论，而是一套非常直观的图形化工具，专门用来研究当系统某个参数（通常是开环增益K）从0变化到无穷大时，闭环极点在复平面上移动的轨迹。你可以把它想象成一张动态的“地图”，清晰地告诉你：增益调大，系统会变快但可能振荡；增益调小，系统会变稳但可能反应迟钝。这张地图，就是根轨迹。
对于考研党来说，根轨迹是《自动控制原理》的核心章节，无论是胡寿松老师的经典教材，还是卢京潮老师的授课视频，都会花大量篇幅讲解。它不仅是笔试计算题的重头戏，更是理解系统动态性能与参数之间内在联系的关键。对于工程人员，掌握了根轨迹，就等于掌握了系统调试和控制器设计的“预判”能力，不用等到实物做出来再盲目调参，在图纸和仿真阶段就能心中有数。这次，我就结合自己学习和工程应用中的体会，把线性系统根轨迹分析的核心要点、绘制规则、实战应用以及那些容易踩的坑，系统地梳理一遍。咱们不搞纯理论推导，重点放在“怎么用”和“为什么这么用”上。
2. 根轨迹的核心思想与绘制规则全解析
2.1 根轨迹的“灵魂”：幅角条件与幅值条件
根轨迹的所有绘制规则，都源于两个最基础的方程，它们构成了根轨迹理论的基石。假设系统的开环传递函数为 (G(s)H(s))，其中 (K) 为可变的开环增益，那么闭环特征方程为 (1 + KG(s)H(s) = 0)。将其改写为： [ KG(s)H(s) = -1 ] 这个复数方程可以拆解为两个部分：
- 幅角条件：(\angle G(s)H(s) = \pm180^\circ (2k+1), \quad k=0,1,2,...)
- 幅值条件：(|KG(s)H(s)| = 1)，即 (K = \frac{1}{|G(s)H(s)|})
幅角条件是根轨迹的“定性”条件。它决定了复平面上一个点 (s) 是否在根轨迹上。简单说，如果从所有开环零点向该点 (s) 引向量，这些向量的幅角（角度）之和，减去从所有开环极点向该点 (s) 引向量的幅角之和，结果等于 (180^\circ) 的奇数倍，那么这个点 (s) 就一定在根轨迹上。所有满足这个条件的点连成的线，就是根轨迹。这是绘制根轨迹的根本依据，所有后续规则都是为了让手工绘制满足这个条件的曲线变得更简单。
幅值条件是根轨迹的“定量”条件。一旦确定了根轨迹上的某个点 (s)，幅值条件就能告诉你，为了使得闭环极点位于这个点，所需要的开环增益 (K) 是多少。计算方法是 (K) 等于所有开环极点到该点 (s) 的向量长度之积，除以所有开环零点到该点 (s) 的向量长度之积。
注意：很多初学者会混淆开环和闭环。根轨迹画的是闭环极点随开环增益K变化的轨迹，但绘制时用的是开环传递函数(G(s)H(s)) 的零极点。务必时刻清楚你在操作的对象。
2.2 手绘根轨迹的“九阴真经”：八大核心规则
基于幅角条件，可以推导出一套完整的绘图规则。掌握它们，即使不用MATLAB，你也能快速画出根轨迹的草图，这对定性分析和快速判断至关重要。
规则1：起点与终点。根轨迹起始于开环极点（K=0），终止于开环零点（K→∞）。如果开环极点数量 (n) 大于开环零点数量 (m)，则有 (n-m) 条根轨迹趋向于无穷远处。
规则2：分支数、连续性与对称性。根轨迹的分支数等于开环极点数 (n)，且关于实轴对称（因为复极点共轭出现）。
规则3：实轴上的根轨迹。实轴上某一段右侧，开环实数零极点数目之和为奇数时，该段实轴属于根轨迹。这是最常用的规则之一，用于快速确定实轴上的轨迹区间。
规则4：渐近线。当 (s \to \infty) 时，那 (n-m) 条趋向无穷远的根轨迹，会沿着一组渐近线延伸。渐近线与实轴的交点 (\sigma_a) 和夹角 (\phi_a) 由下式决定： [ \sigma_a = \frac{\sum (极点实部) - \sum (零点实部)}{n - m} ] [ \phi_a = \frac{\pm180^\circ(2k+1)}{n - m}, \quad k=0,1,...,n-m-1 ] 这个规则帮助我们确定根轨迹在远离原点时的走向。
规则5：分离点与会合点。两条或更多根轨迹分支在实轴上相遇后又分开的点。它可以通过求解方程 (\frac{dK}{ds} = 0) 得到，其中 (K = -\frac{1}{G(s)H(s)})。更实用的方法是解方程： [ \sum_{i=1}^{n} \frac{1}{d - p_i} = \sum_{j=1}^{m} \frac{1}{d - z_j} ] 其中 (d) 是分离点坐标，(p_i) 和 (z_j) 是开环极点和零点。分离点一定位于实轴上的根轨迹线段上。
规则6：起始角与终止角。根轨迹离开开环复极点时的出射角，和进入开环复零点时的入射角。计算公式源于幅角条件，用于确定复平面根轨迹的初始走向。例如，从某个复极点 (p_k) 出发的出射角 (\theta_{p_k}) 满足： [ \theta_{p_k} = 180^\circ + \sum_{j=1}^{m} \angle(p_k - z_j) - \sum_{\substack{i=1 \ i \neq k}}^{n} \angle(p_k - p_i) ]
规则7：与虚轴的交点。根轨迹穿越虚轴，意味着系统处于临界稳定状态。交点坐标 (j\omega) 和对应的临界增益 (K_c) 可以通过劳斯判据或直接代入 (s = j\omega) 到特征方程求解。
规则8：闭环极点的和与积。对于特征方程 (s^n + a_{n-1}s^{n-1} + ... + a_0 = 0)，当 (n-m \geq 2) 时，闭环极点之和等于开环极点之和，且为常数。这个性质可以用来校验根轨迹的走向或估算其他极点的位置。
实操心得：这八条规则不是孤立的，而是按顺序使用的“流水线”。我的习惯是：1)标零极点；2)画实轴轨迹；3)算渐近线；4)找分离点；5)算出射/入射角；6)求虚轴交点。按这个顺序，一幅完整的根轨迹草图就出来了。考试时，即使时间紧张，完成前四步也能拿到大部分分数。
3. 从理论到实战：典型系统根轨迹绘制与性能分析
3.1 二阶系统：理解根轨迹的“标准模板”
我们从一个最简单的二阶系统开始：开环传递函数 (G(s) = \frac{K}{s(s+2)})。没有开环零点，有两个开环极点：(p_1=0, p_2=-2)。
- 起点与分支：两个起点（0， -2），两条分支。
- 实轴轨迹：实轴上，区间 ([-2, 0]) 右侧零极点总数为1（奇数），故该段是根轨迹。
- 渐近线：(n=2, m=0)，故有2条渐近线。交点 (\sigma_a = (0 + (-2))/2 = -1)，夹角 (\phi_a = \pm90^\circ)。
- 分离点：由公式 (\frac{1}{d} + \frac{1}{d+2} = 0)，解得 (d = -1)。该点在实轴轨迹上，合理。
- 虚轴交点：闭环特征方程 (s^2+2s+K=0)。劳斯表第一列：(s^2: 1, K); (s^1: 2); (s^0: K)。令 (K=0)？不对，令 (s^1) 行全零？这里显然不会。更直接地，令 (s=j\omega) 代入：((j\omega)^2 + 2j\omega + K = -\omega^2 + K + 2j\omega = 0)。解得实部虚部为零：(-\omega^2 + K = 0) 且 (2\omega = 0)。得到 (\omega=0, K=0)。这意味着根轨迹在 (K=0) 时从原点出发，随着K增大沿实轴向-1移动，在 (d=-1) 处分离进入复平面，沿垂直渐近线上下运动。当K继续增大，极点实部恒为-1（因为和不变），虚部增大，系统阻尼比减小，振荡加剧。
这个简单的例子完美展示了根轨迹如何揭示性能：当闭环极点为实数时（K较小时），系统是过阻尼，无超调但响应慢；当极点为复数时（K较大），系统是欠阻尼，响应变快但有超调；极点离虚轴越远（实部更负），响应越快；极点离实轴越远（虚部越大），振荡频率越高。根轨迹上的每一个点，都对应一组确定的时域性能指标。
3.2 带零点的三阶系统：零点的“吸引”效应
考虑 (G(s) = \frac{K(s+1)}{s(s+2)(s+3)})。开环极点：0, -2, -3；开环零点：-1。
- 起点与终点：3个起点，1个终点，有2条分支趋向无穷远。
- 实轴轨迹：分段判断。区间 ((-\infty, -3]) 右侧零极点总数3（奇数），是轨迹；([-3, -2]) 右侧总数2（偶数），不是；([-2, -1]) 右侧总数1（奇数），是；([-1, 0]) 右侧总数2（偶数），不是；([0, +\infty)) 右侧总数3（奇数），是。所以实轴轨迹在 ((-\infty, -3]), ([-2, -1]), ([0, +\infty))。
- 渐近线：(n-m=2)。交点 (\sigma_a = ((0-2-3) - (-1)) / 2 = (-4)/2 = -2)。夹角 (\phi_a = \pm90^\circ)。
- 分离点：计算稍复杂，但通常位于实轴轨迹段 ([-2, -1]) 上。通过解方程或试探法可求得一个分离点。同时，由于有一条从原点出发的轨迹终点在无穷远，它与从-2出发的轨迹可能在复平面某处分离，这需要结合出射角计算。
关键观察：增加了零点 (s=-1) 后，根轨迹被“吸引”向零点的方向弯曲。这使得对于同样的增益K，系统的闭环极点可能比没有零点时具有更小的阻尼（更靠近虚轴），但同时也可能改善稳定性（改变了轨迹走向）。在控制器设计中，我们常常通过添加零点（例如PD控制器）来“塑造”根轨迹，使其通过期望的性能区域。
3.3 条件稳定系统：根轨迹揭示的“稳定域”
考虑一个有趣系统：(G(s) = \frac{K(s+2)}{s(s+1)(s+3)})。绘制其根轨迹会发现，随着K增大，根轨迹会从左侧穿越虚轴进入右半平面（系统不稳定），然后随着K继续增大，又会再次穿越虚轴回到左半平面。这意味着系统只在某个特定的K值区间内是稳定的。这种系统称为“条件稳定系统”。根轨迹清晰地揭示了这一现象，并给出了稳定的K值范围。在工程中，这提醒我们必须谨慎选择工作点，并避免参数漂移导致系统失稳。
4. 基于根轨迹的系统分析与控制器设计思路
4.1 动态性能评估：在轨迹上“读”出指标
根轨迹不仅是一幅图，更是一个性能查询工具。假设我们有一对主导闭环极点 (s = -\sigma \pm j\omega_d)。
- 调节时间 (t_s)：主要由实部 (\sigma) 决定，(t_s \approx \frac{3.5}{\sigma})（按2%误差带）。在根轨迹上，作一条平行于虚轴的直线 (s = -\sigma)，这条线就是“等调节时间线”。轨迹点在这条线左边，调节时间就更短。
- 峰值时间 (t_p) 和振荡频率 (\omega_d)：由虚部 (\omega_d) 决定，(t_p = \frac{\pi}{\omega_d})。在根轨迹上，作一条平行于实轴的直线，其纵坐标即为 (\omega_d)。
- 超调量 (M_p%)：与阻尼比 (\zeta) 相关，(\zeta = \cos \theta)，其中 (\theta) 是极点与负实轴夹角。在根轨迹上，从原点出发的射线与负实轴的夹角，直接决定了阻尼比和超调量。常用的等阻尼比线（射线）有 (\zeta=0.707)（对应约5%超调，最佳阻尼）和 (\zeta=0.5)（对应约16%超调）。
实操方法：在绘制好的根轨迹图上，叠加画出等调节时间线（竖线）、等阻尼比线（射线）和等自然频率线（圆弧，(s = -\zeta\omega_n \pm j\omega_n\sqrt{1-\zeta^2})，其模长为 (\omega_n)）。根轨迹与这些线的交点，就对应了能满足特定性能指标的闭环极点位置及其所需的增益K。
4.2 控制器设计：用零极点“雕刻”轨迹
根轨迹法设计控制器，核心思想就是通过增加新的开环零极点，改变原有根轨迹的形状，使其通过复平面上代表优良性能的区域。
- 比例微分（PD）控制器 / 超前校正：(G_c(s) = K_p + K_d s = K_d(s + \frac{K_p}{K_d}))。它引入了一个零点。这个零点像一个“磁铁”，会把根轨迹向左半平面“吸引”，从而增大系统阻尼，提高稳定性，减少超调，加快响应。通常将零点放置在期望主导极点附近或左侧。
- 比例积分（PI）控制器 / 滞后校正：(G_c(s) = K_p + \frac{K_i}{s} = \frac{K_p(s + \frac{K_i}{K_p})}{s})。它引入了一个极点（原点）和一个零点（靠近原点）。极点会“排斥”根轨迹，使其向右弯曲，可能恶化动态性能。但PI控制器的主要目的是提高系统型别，消除静差。设计时，零极点必须非常接近，形成“偶极子”，这样它们对根轨迹主体形状的影响几乎相互抵消，但积分作用得以保留。
- 比例-积分-微分（PID）控制器：结合了以上两者。设计时，可以先利用PD部分将主导极点“拉”到期望位置，再利用PI部分在不显著改变主导极点位置的前提下提高型别。
注意事项：添加零极点是一把双刃剑。不恰当的零点可能使根轨迹过分弯曲进入右半平面，导致条件稳定甚至不稳定。不恰当的极点（尤其是靠近虚轴的）会严重恶化系统性能。设计时必须结合幅值条件计算实际增益，并在MATLAB中仿真验证。
5. MATLAB辅助分析与工程实践中的常见陷阱
5.1 使用MATLAB高效绘制与分析
虽然手绘草图至关重要，但工程中我们依赖MATLAB进行精确绘制和深入分析。
% 示例：绘制 G(s) = K(s+1)/(s(s+2)(s+3)) 的根轨迹 num = [1 1]; % 分子系数，对应 (s+1) den = conv([1 0], conv([1 2], [1 3])); % 分母系数，对应 s(s+2)(s+3) sys = tf(num, den); % 建立传递函数模型 rlocus(sys); % 绘制根轨迹图 grid on; % 打开网格，便于查看阻尼比和自然频率 sgrid; % 叠加等阻尼比和等自然频率网格 % 交互式获取增益和极点：在图形窗口，点击轨迹上任一点，会显示该点增益、极点位置、阻尼比、超调等信息。rlocus函数是核心。sgrid函数生成的网格线极其有用，能快速评估性能。更高级的用法是结合rlocfind函数用鼠标选取点，或直接用[K, poles] = rlocfind(sys)来编程计算特定极点处的增益。
5.2 工程应用中的典型问题与排查技巧
问题1：理论计算分离点与MATLAB显示不符？
- 原因：分离点方程解出的可能是“会合点”（多条轨迹汇聚到实轴）或“分离点”（从实轴分开进入复平面），也可能存在复平面上的分离点（对于更复杂的系统）。手绘规则主要针对实轴分离点。
- 排查：使用MATLAB精确绘制，观察该点附近轨迹的走向。检查计算过程，确认方程推导无误。对于高阶系统，分离点方程可能有多个解，需根据幅角条件判断哪些是有效的。
问题2：根据根轨迹设计了控制器，但仿真/实际效果很差？
- 原因1：主导极点假设失效。设计时假设了一对共轭复极点主导系统响应，但如果系统存在其他非常靠近虚轴的闭环极点（或零点），它们会显著影响响应，使基于主导极点的预测失效。
- 对策：绘制校正后的系统根轨迹，确认所有闭环极点的位置。检查非主导极点的实部是否至少比主导极点实部大5倍以上。可以使用
pole(feedback(G*Gc, 1))查看所有闭环极点。 - 原因2：未考虑传感器噪声和执行器饱和。根轨迹分析基于线性模型。添加的微分环节（D）可能放大高频噪声；积分环节（I）可能导致积分饱和。在实际系统中，需要加入滤波器（如对D项加低通滤波）和抗饱和措施。
- 对策：在仿真中引入噪声和饱和非线性环节进行测试。考虑使用PID的变体，如带滤波的PDF（伪微分反馈）或抗饱和积分。
问题3：根轨迹显示系统始终稳定，但实际调试中会出现振荡？
- 原因：未建模动态。数学模型（如简单的二阶传递函数）未能完全反映真实系统的所有特性，如高频柔性模态、时滞、非线性摩擦等。这些“隐藏”的动力学可能在某个增益下被激发。
- 对策：进行系统辨识，获得更精确的高阶模型。在根轨迹分析时，使用更复杂的模型。在实际调试中，采用增益调度或自适应控制策略。
问题4：考研真题中遇到非180度根轨迹（0度根轨迹）怎么办？
- 识别：当系统是正反馈结构，或特征方程为 (1 - KG(s)H(s) = 0) 时，幅角条件变为 (\angle G(s)H(s) = \pm360^\circ k)。这就是0度根轨迹。
- 规则变化：绘制规则中，与幅角条件直接相关的规则需要改变：实轴上的根轨迹，其右侧零极点数目之和变为偶数；渐近线夹角公式变为 (\phi_a = \frac{\pm360^\circ k}{n-m})。其他规则（如起点终点、分离点计算）不变。
- 关键：审题时务必看清系统结构图是负反馈（-）还是正反馈（+），或者特征方程的形式。
根轨迹分析，从手绘理解到软件辅助，再到工程实践，是一个层层递进的过程。它提供的是一种图形化的直觉，让你能“看见”参数变化如何影响系统本质。我个人的体会是，初期一定要亲手在纸上多画几个不同零极点配置的根轨迹，感受零极点的“引力”和“斥力”，这比死记硬背规则有效得多。到了工程阶段，MATLAB是你的得力助手，但绝不能替代你的理解和判断。最后，永远对模型和实际系统的差异保持警惕，根轨迹给出的是一种基于模型的理想洞察，而真正的控制智慧，在于如何弥合理想与现实之间的鸿沟。
