从根轨迹到系统响应：5个典型阻尼比案例的MATLAB仿真对比

从根轨迹到系统响应：5个典型阻尼比案例的MATLAB仿真对比

如果你刚开始学习自动控制原理，面对课本上抽象的根轨迹图和一堆关于阻尼比ζ的公式，是不是感觉有点云里雾里？根轨迹上那些弯弯曲曲的线，到底和系统实际输出有什么关系？为什么老师说极点靠近虚轴系统就会“振”起来？这些问题，光靠理论推导和想象确实有点困难。

我刚开始学的时候也卡在这里，直到后来在MATLAB里亲手把不同参数的根轨迹画出来，再同时跑一下系统的阶跃响应，那种“原来如此”的顿悟感才真正出现。你会发现，根轨迹不是一幅静态的、孤立的数学图画，而是一个动态系统的“基因图谱”，上面每一个点的位置，都直接决定了系统“性格”的方方面面——是反应敏捷但容易“过冲”，还是稳重迟缓但绝不出错。

这篇文章，我们就用最直观的案例教学法，抛开复杂的公式堆砌，直接动手。我会带你用MATLAB搭建一个经典二阶系统，然后通过调整开环增益，让它的闭环极点沿着根轨迹“跑”起来，分别停留在阻尼比ζ为0.2、0.5、0.7、1.0和1.5这五个典型位置。每一个位置，我们都会同时观察根轨迹图、极点坐标和对应的阶跃响应曲线。你看完就能明白：极点实部如何决定响应速度，虚部如何决定振荡频率，而阻尼角这个几何概念又是如何把两者统一起来，直观地告诉你系统“振不振”。

文末我会提供完整的、可下载的.m文件源码，你可以在自己的MATLAB上原样复现，甚至修改参数进行探索。我们的目标是：让根轨迹和阻尼比，从课本上的考点，变成你设计系统时工具箱里一件趁手的可视化武器。

1. 搭建实验环境：一个经典系统的根轨迹

在开始对比五个案例之前，我们得先有一个可供观察的“舞台”。这个舞台就是一个具体的控制系统，它的根轨迹要足够典型，能清晰展示极点从欠阻尼到过阻尼的完整迁移路径。

我选择了一个在教材和工程中都非常常见的单位负反馈系统。它的开环传递函数是： [ G(s) = \frac{K}{s(s+2)} ] 这里的 (K) 就是我们将要调整的开环增益。为什么选这个？因为它简单，其闭环特征方程是二阶的（(s^2 + 2s + K = 0)），根轨迹清晰可预测，非常适合教学。同时，它又包含了所有阻尼比状态的可能性。

1.1 在MATLAB中创建系统模型

打开你的MATLAB，我们首先定义这个系统。最直接的方法是使用 tf 函数创建传递函数对象。

% 定义开环传递函数 G(s) = K / (s*(s+2))
% 先固定K=1，用于绘制根轨迹的基础模型
num = [1]; % 分子多项式系数，对应 K，我们先设为1
den = conv([1, 0], [1, 2]); % 分母多项式系数，s*(s+2) 的乘积
G = tf(num, den); % 创建传递函数对象

% 查看系统信息
disp('开环传递函数 G(s):');
G

运行这段代码，命令行会输出 G(s) = 1 / (s^2 + 2s)，确认系统模型创建正确。

接下来，我们绘制这个系统的根轨迹。根轨迹展示的是闭环极点随 (K) 从 (0) 变化到 (+\infty) 时的运动轨迹。

% 绘制根轨迹图
figure('Position', [100, 100, 800, 600]); % 设置图形窗口大小
rlocus(G);
grid on;
title('系统 G(s) = K / [s(s+2)] 的根轨迹');
xlabel('实轴 (Real Axis)');
ylabel('虚轴 (Imaginary Axis)');
sgrid; % 添加等阻尼比线和等自然频率线网格，便于观察

运行后，你会看到一条从开环极点（s=0 和 s=-2）出发，在实轴上某点汇合后垂直向上、下延伸至无穷远的根轨迹。sgrid 命令生成的扇形网格线就是等阻尼比线，每条射线与负实轴的夹角 (\beta) 满足 (\cos\beta = \zeta)。这是我们后续定位特定阻尼比极点的关键工具。

1.2 理解根轨迹与阻尼比的几何关系

在我们得到的根轨迹图上，你可以清晰地看到两条分支。当 (K) 很小时，两个闭环极点都是负实数（位于实轴上的根轨迹部分）。随着 (K) 增大，两个极点沿着实轴相向移动。当 (K=1) 时，两个极点在 (s=-1) 处重合，此时对应临界阻尼 ((\zeta=1))。继续增大 (K)，极点离开实轴，成为一对共轭复根，并沿着垂直方向上下运动。

注意：sgrid 默认生成的阻尼比线是 0.1:0.1:0.9，你可能需要更密的网格来精确定位。可以使用 sgrid(z, wn) 自定义，例如 sgrid(0:0.05:1, []) 来生成从0到1，步长0.05的阻尼比线。

此时，根轨迹上任意一点（即某个 (K) 值下的闭环极点）的阻尼比，就等于该点与原点连线（即向径）和负实轴夹角的余弦值。这个夹角就是阻尼角。

• 夹角为0度：极点在负实轴上，(\cos0=1)，(\zeta=1) 或 >1（过阻尼）。
• 夹角为60度：(\cos60^\circ=0.5)，(\zeta=0.5)（欠阻尼）。
• 夹角接近90度：(\cos\beta) 接近0，(\zeta) 接近0，振荡非常剧烈。

这个几何关系是连接抽象的根轨迹图与具体的系统响应的桥梁。下面，我们就让极点在这条“桥梁”上选取五个有代表性的位置停下来，看看系统的“表情”各是什么样。

2. 案例一：欠阻尼振荡 ((\zeta = 0.2))

首先，我们来看一个阻尼很弱的情况：(\zeta = 0.2)。根据几何关系，阻尼角 (\beta = \arccos(0.2) \approx 78.5^\circ)。我们需要在根轨迹图上，找到位于这条78.5度射线上、且处于左半平面的那个交点。

2.1 确定极点位置与增益K

对于我们的系统 (G(s) = K/[s(s+2)])，其闭环特征方程为 (s^2 + 2s + K = 0)。闭环极点为： [ s = -1 \pm \sqrt{1 - K} ] 当 (K > 1) 时，极点为共轭复根：(s = -1 \pm j\sqrt{K-1})。 此时，极点的实部 (\sigma = -1)，虚部 (\omega_d = \sqrt{K-1})。阻尼比 (\zeta)