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 (2k1), \quad k0,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起点与终点。根轨迹起始于开环极点K0终止于开环零点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(2k1)}{n - m}, \quad k0,1,...,n-m-1 ] 这个规则帮助我们确定根轨迹在远离原点时的走向。规则5分离点与会合点。两条或更多根轨迹分支在实轴上相遇后又分开的点。它可以通过求解方程 (\frac{dK}{ds} 0) 得到其中 (K -\frac{1}{G(s)H(s)})。更实用的方法是解方程 [ \sum_{i1}^{n} \frac{1}{d - p_i} \sum_{j1}^{m} \frac{1}{d - z_j} ] 其中 (d) 是分离点坐标(p_i) 和 (z_j) 是开环极点和零点。分离点一定位于实轴上的根轨迹线段上。规则6起始角与终止角。根轨迹离开开环复极点时的出射角和进入开环复零点时的入射角。计算公式源于幅角条件用于确定复平面根轨迹的初始走向。例如从某个复极点 (p_k) 出发的出射角 (\theta_{p_k}) 满足 [ \theta_{p_k} 180^\circ \sum_{j1}^{m} \angle(p_k - z_j) - \sum_{\substack{i1 \ 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(s2)})。没有开环零点有两个开环极点(p_10, p_2-2)。起点与分支两个起点0 -2两条分支。实轴轨迹实轴上区间 ([-2, 0]) 右侧零极点总数为1奇数故该段是根轨迹。渐近线(n2, m0)故有2条渐近线。交点 (\sigma_a (0 (-2))/2 -1)夹角 (\phi_a \pm90^\circ)。分离点由公式 (\frac{1}{d} \frac{1}{d2} 0)解得 (d -1)。该点在实轴轨迹上合理。虚轴交点闭环特征方程 (s^22sK0)。劳斯表第一列(s^2: 1, K); (s^1: 2); (s^0: K)。令 (K0)不对令 (s^1) 行全零这里显然不会。更直接地令 (sj\omega) 代入((j\omega)^2 2j\omega K -\omega^2 K 2j\omega 0)。解得实部虚部为零(-\omega^2 K 0) 且 (2\omega 0)。得到 (\omega0, K0)。这意味着根轨迹在 (K0) 时从原点出发随着K增大沿实轴向-1移动在 (d-1) 处分离进入复平面沿垂直渐近线上下运动。当K继续增大极点实部恒为-1因为和不变虚部增大系统阻尼比减小振荡加剧。这个简单的例子完美展示了根轨迹如何揭示性能当闭环极点为实数时K较小时系统是过阻尼无超调但响应慢当极点为复数时K较大系统是欠阻尼响应变快但有超调极点离虚轴越远实部更负响应越快极点离实轴越远虚部越大振荡频率越高。根轨迹上的每一个点都对应一组确定的时域性能指标。3.2 带零点的三阶系统零点的“吸引”效应考虑 (G(s) \frac{K(s1)}{s(s2)(s3)})。开环极点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-m2)。交点 (\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(s2)}{s(s1)(s3)})。绘制其根轨迹会发现随着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) 是极点与负实轴夹角。在根轨迹上从原点出发的射线与负实轴的夹角直接决定了阻尼比和超调量。常用的等阻尼比线射线有 (\zeta0.707)对应约5%超调最佳阻尼和 (\zeta0.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(s1)/(s(s2)(s3)) 的根轨迹 num [1 1]; % 分子系数对应 (s1) den conv([1 0], conv([1 2], [1 3])); % 分母系数对应 s(s2)(s3) 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是你的得力助手但绝不能替代你的理解和判断。最后永远对模型和实际系统的差异保持警惕根轨迹给出的是一种基于模型的理想洞察而真正的控制智慧在于如何弥合理想与现实之间的鸿沟。