别再死记硬背公式了！用MATLAB手把手教你玩转根轨迹，分析系统稳定性
用MATLAB实战根轨迹分析：从图形读懂系统稳定性
打开MATLAB，输入几行代码，你就能看到抽象的控制理论在屏幕上"活"过来——这就是根轨迹法的魅力。作为自动控制原理中的核心分析方法，根轨迹不仅能帮你避开繁琐的数学推导，更能通过直观的图形展示系统参数变化时极点的移动轨迹。本文将带你用MATLAB完成一次完整的根轨迹分析实战，从基础绘图到稳定性判断，再到性能优化，让你真正掌握这个强大的工具。
1. 准备工作：认识根轨迹与MATLAB基础
根轨迹法是W.R. Evans于1948年提出的一种图形化分析方法，它描绘了系统开环增益K从0变化到∞时，闭环系统极点在复平面上的运动轨迹。这种方法之所以经久不衰，正是因为它将复杂的数学问题转化为了直观的几何问题。
在开始前，确保你的MATLAB安装了Control System Toolbox。我们将使用以下传递函数作为示例：
num = [1 3]; % 分子多项式系数 (s+3) den = conv([1 1], conv([1 2], [1 1 1.25])); % 分母多项式系数 (s+1)(s+2)(s²+s+1.25) G = tf(num, den); % 创建传递函数对象提示：conv函数用于多项式乘法，可以方便地展开分母中的乘积项。
理解这个传递函数的组成很重要：
- 零点：使分子为零的s值（本例中为s=-3）
- 极点：使分母为零的s值（本例中为s=-1, s=-2, s=-0.5±j）
2. 绘制根轨迹：从基础到进阶
2.1 基础绘图与解读
绘制根轨迹的最简单命令是：
rlocus(G) grid on执行后你将看到一幅曲线图，其中：
- 红色×表示开环极点
- 蓝色○表示开环零点
- 曲线表示K增大时闭环极点的移动路径
关键观察点：
- 根轨迹起始于开环极点（K=0）
- 终止于开环零点或无穷远处（K→∞）
- 实轴上的根轨迹遵循"奇偶规则"
- 渐近线角度和中心点可以通过公式计算
2.2 获取特定增益下的极点位置
要查看K=1时的极点位置：
[K, poles] = rlocfind(G)在图形窗口点击你感兴趣的位置，MATLAB会返回该点的增益值和对应的极点坐标。这对于确定系统稳定性边界非常有用。
2.3 可视化增强技巧
为了让图形更易读，可以添加以下代码：
rlocus(G) sgrid % 添加阻尼比和自然频率网格 title('系统根轨迹图') xlabel('实轴') ylabel('虚轴') axis equal % 保持纵横比一致3. 稳定性分析：从图形到数值
3.1 图形法判断稳定性
根轨迹与虚轴的交点是判断系统稳定性的关键：
- 所有极点位于左半平面→ 系统稳定
- 任何极点位于右半平面→ 系统不稳定
- 极点在虚轴上→ 临界稳定
在我们的示例中，可以看到一对共轭极点随着K增大逐渐向右移动，最终会穿越虚轴进入右半平面。
3.2 计算临界增益值
除了目测，我们可以精确计算临界稳定时的K值：
[K, poles] = rlocfind(G) % 点击根轨迹与虚轴的交点MATLAB会返回临界增益K≈1.94，此时极点位于±1.39j处。这意味着：
- 当K<1.94时，系统稳定
- 当K>1.94时，系统不稳定
- 当K=1.94时，系统处于临界稳定状态
3.3 验证计算结果的正确性
为了验证图形法的准确性，我们可以使用Routh判据进行理论计算：
- 构造闭环特征方程
- 建立Routh表
- 求解使第一列出现零行的K值
你会发现两种方法得到的结果一致，但图形法明显更直观高效。
4. 系统性能优化：超越稳定性
4.1 利用阻尼比设计系统
根轨迹上的每一点都对应特定的阻尼比(ζ)和自然频率(ωₙ)。在MATLAB中，可以使用sgrid命令叠加等阻尼比和等自然频率线：
rlocus(G) sgrid(0.6, []) % 绘制ζ=0.6的等阻尼比线点击根轨迹与这条线的交点，MATLAB会给出对应的K值。选择这个K值，系统将具有约0.6的阻尼比，通常能提供良好的动态响应。
4.2 主导极点与降阶
当系统存在主导极点时（距离虚轴最近且附近无零点），可以忽略其他极点的影响。在我们的例子中：
- 当K较小时，实极点(-1,-2)可能主导系统响应
- 当K增大时，复极点(-0.5±j)可能成为主导
通过比较极点实部的绝对值，可以判断哪些极点可以忽略。
4.3 添加零极点改善性能
有时我们需要修改系统结构而不仅仅是调整增益。例如，添加一个零点可以"吸引"根轨迹向左移动：
G_new = tf([1 5],1)*G; % 在s=-5处添加一个零点 rlocus(G_new)比较新旧根轨迹，你会发现新增零点确实改善了系统的稳定性裕度。
5. 实战技巧与常见问题
5.1 提高绘图精度
默认的根轨迹可能在某些区域分辨率不足。可以通过指定K的范围来提高精度：
K = logspace(-2,2,1000); % 从10^-2到10^2的1000个对数间隔点 rlocus(G,K)5.2 处理复杂系统
对于高阶系统，根轨迹可能非常复杂。这时可以：
- 使用zoom命令放大感兴趣的区域
- 结合使用rlocus和rlocfind精确定位关键点
- 分多次绘制不同K范围的根轨迹
5.3 根轨迹与频域分析的结合
根轨迹展示了极点位置随K的变化，而Bode图展示了频响特性。结合两者可以全面评估系统性能：
K = 1.5; % 选择一个稳定增益 T = feedback(K*G,1); % 构造闭环系统 bode(T)5.4 常见错误与排查
- 根轨迹不显示：检查传递函数是否正确输入，确保Control System Toolbox已安装
- 结果不符合预期：确认分子分母多项式的顺序是否正确（降幂排列）
- 临界增益计算不准确：尝试缩小K的范围，提高计算精度
6. 扩展应用：从理论到工程实践
6.1 参数根轨迹
除了增益K，我们还可以研究其他参数变化的影响。例如，研究时间常数τ的变化：
s = tf('s'); tau = 1; G_param = 1/(s*(tau*s+1)); rlocus(G_param)6.2 多变量系统分析
对于MIMO系统，可以使用sigmaplot分析奇异值轨迹，这是根轨迹在多变量系统中的推广。
6.3 实时交互工具
MATLAB的Control System Designer提供了图形化界面，可以交互式地调整控制器参数并实时观察根轨迹变化：
controlSystemDesigner(G)6.4 与Simulink的集成
在Simulink模型中，你可以：
- 使用LTI System块导入MATLAB工作空间中的传递函数
- 在仿真过程中实时调整参数
- 使用Linear Analysis Tool自动生成根轨迹
7. 从根轨迹到控制器设计
理解了根轨迹后，你可以开始设计控制器了。以PID控制器为例：
Kp = 1; Ki = 0.5; Kd = 0.1; C = pid(Kp,Ki,Kd); rlocus(C*G) % 观察控制器如何改变根轨迹设计时考虑：
- 根轨迹是否通过期望的性能区域
- 稳态误差要求
- 抗干扰能力
8. 实际工程案例：位置控制系统
考虑一个直流电机位置控制系统，其开环传递函数为：
J = 0.01; b = 0.1; K = 1; R = 1; L = 0.5; s = tf('s'); P_motor = K/((J*s+b)*(L*s+R)+K^2); rlocus(P_motor)通过根轨迹分析，我们可以：
- 确定使系统稳定的最大增益
- 选择合适的增益满足超调量要求
- 评估增加速度反馈的影响
9. 高级技巧：自定义根轨迹分析
对于特殊需求，你可以直接访问根轨迹数据：
[K, poles] = rlocus(G); % 手动处理数据 plot(real(poles), imag(poles), 'b-')这允许你：
- 创建自定义可视化
- 与其他分析工具集成
- 实现自动化的参数优化
10. 资源与进一步学习
要深入掌握根轨迹法，建议：
- 实践不同的系统模型
- 比较根轨迹法与频域方法的结果
- 阅读MATLAB文档中的相关案例
- 尝试在实时硬件上验证你的设计
记住，熟练使用根轨迹法不仅能帮你通过考试，更能让你在实际工程中快速评估和设计控制系统。下次当你面对一个复杂的控制系统问题时，不妨先画出它的根轨迹——图形化的视角往往会带给你意想不到的洞见。