欠阻尼二阶系统单位阶跃响应与动态性能指标详解

1. 项目概述：为什么"欠阻尼"值得单独拿出来分析

如果你学过自控原理，最开始接触二阶系统时，看到一堆参数、曲线和公式，很容易一头雾水。尤其是"欠阻尼"这个词，听着很玄乎，但它其实是工程实际里最常见、也最值得研究的一种状态。这里的核心关键词是：欠阻尼、二阶系统、单位阶跃响应、动态性能指标。我在做控制系统设计与调试时，几乎每个项目都要跟它打交道，因为现实中的机械系统、电路系统、伺服驱动系统，很多都属于欠阻尼类型。分析清楚它的单位阶跃响应，你才能真正理解"超调量""调节时间""峰值时间"这些指标的物理含义，也才能做出合理的控制器参数设计。这篇文章适合正在学习自控原理的本科生、刚入门做嵌入式控制或者伺服调参的工程师，以及对系统动态响应有困惑的开发者。我会从数学模型、动态指标推导、仿真实操、常见坑点这几个角度把欠阻尼二阶系统讲透。

2. 欠阻尼二阶系统建模与特征参数

2.1 标准二阶系统数学描述与共振直觉

二阶系统最标准的闭环传递函数长这样：

Φ(s) = ωn² / (s² + 2ζωn·s + ωn²)

这个公式里，ωn 是自然频率，ζ 是阻尼比。阻尼比 ζ 的取值直接决定了系统的响应形态：

• ζ < 0：负阻尼，发散，工程上要避免（除非故意用正反馈做振荡器）
• ζ = 0：无阻尼，等幅振荡，永远不衰减
• 0 < ζ < 1：欠阻尼，响应振荡衰减，最终趋于稳定
• ζ = 1：临界阻尼，最快不振荡地回到稳态
• ζ > 1：过阻尼，慢吞吞地逼近稳态，无振荡

为什么欠阻尼最常见？因为很多物理系统本身的固有阻尼很小。比如一个弹簧-质量-阻尼系统，理想化情况下阻尼系数可能为零，但实际总有摩擦、空气阻力等带来少量阻尼，这就形成了欠阻尼。再比如直流电机带负载，转动惯量和粘性摩擦系数往往使系统的阻尼比落在0.1~0.8之间。所以研究欠阻尼段，比研究临界阻尼和过阻尼更有实战价值。

从根轨迹角度看，欠阻尼系统有一对共轭复极点：

s1,2 = -ζωn ± j·ωn·√(1-ζ²)

其中 σ = ζωn 叫衰减系数，ωd = ωn·√(1-ζ²) 叫阻尼振荡频率。这对极点位于左半平面，但不在实轴上，所以响应既有指数衰减项，又有正弦振荡项。

2.2 单位阶跃输入的响应推导过程

给系统输入 r(t) = 1(t)，即单位阶跃信号，其拉普拉斯变换 R(s) = 1/s。输出为：

C(s) = Φ(s)·R(s) = ωn² / [s·(s² + 2ζωn·s + ωn²)]

对欠阻尼情况（0 < ζ < 1），部分分式展开后做拉普拉斯反变换，得到时域表达式：

c(t) = 1 - (e^(-ζωn·t) / √(1-ζ²)) · sin(ωd·t + β)

其中 β = arctan(√(1-ζ²) / ζ)，也就是极点在复平面上的相角。

这个公式看着有点吓人，但拆开看很清晰：稳态值为1；动态部分是一个振幅按 e^(-ζωn·t) 衰减的正弦波。振幅衰减的速度由 ζωn 决定，振荡的频率由 ωd 决定。初看这条曲线就像一条弹簧在阻尼作用下的往复运动，最后停在目标位置。

我在实际推导时有个小习惯：不要只记最终公式，自己动手从部分分式展开走一遍。因为只有你自己展开过，才知道分母的配方法、实部虚部怎么对应，后面遇到高阶系统才能迁移。

2.3 极点在复平面上的位置决定一切

欠阻尼系统的极点具有统一的几何表达：位于以原点为圆心、半径为 ωn 的圆上，并落在与负实轴夹角为 β 的直线上。

ζ 越小，极点越靠近虚轴，β 越接近90°，振荡越剧烈，衰减越慢。ζ 越大，极点越靠近负实轴，β 越小。当 ζ → 1 时，两个极点聚到实轴上，进入临界阻尼。

这个几何关系很有用。比如你想设计一个系统，要求超调量不超过20%，可以直接反推出 ζ ≥ 0.456。反过来，给定 ζ = 0.707（工程上常说的最佳阻尼比），极点就在与负实轴成45°的直线上，超调量大约只有4.3%，调节时间也较短。很多控制器设计都拿 ζ = 0.707 作为单目标调试的起点，因为此时系统性能非常均衡。

3. 动态性能指标的硬核计算与工程解读

3.1 峰值时间、超调量、调节时间的完整推导

欠阻尼系统有几个核心指标，考试和工程都要用。我直接给结论，并结合推导思路。

上升时间 tr：响应从稳态值的10%上升到90%所需的时间（也有定义为0~100%的，看教材）。对欠阻尼系统可以近似为：

tr ≈ (π - β) / ωd

峰值时间 tp：响应达到第一个峰值的时间。对欠阻尼阶跃响应求导并令其为零，可得：

tp = π / ωd

这个结果很优雅：峰值时间只与阻尼振荡频率有关，与 ζ 无关（当然 ωd 本身含有 ζ）。物理意义是振荡波形的第一个半周期结束时刻。

超调量 Mp：峰值减去稳态值再除以稳态值的百分比：

Mp = e^(-ζπ / √(1-ζ²)) × 100%

这个公式只需要 ζ 一个参数。我在做调参时经常用这个公式反向算：如果工艺要求超调不能超过10%，则要求 ζ ≥ 0.591。这项指标对系统安全很重要，比如电机转速环如果超调过大，可能触发过流报警。

调节时间 ts：响应进入并保持在一个误差带内所需时间，工程上常用 ±5% 或 ±2% 误差带。对欠阻尼系统，由包络线近似：

ts = -ln(Δ·√(1-ζ²)) / (ζωn) （Δ为误差带宽度，如0.05或0.02）

工程简化计算公式，5%误差带常用：

ts ≈ 3 / (ζωn)

2%误差带常用：

ts ≈ 4 / (ζωn)

注意，这些公式是基于指数包络线的近似，实际曲线可能因为正弦项的零点与包络交点有些微出入，但工程上足够准确。唯一要注意的是：调节时间的计算必须从响应进入误差带且不再越界算起，而不是从首次进入误差带算起。首次进入后可能会再次超出去，这些细节在数值仿真的指标提取中很关键。

3.2 指标之间的耦合关系与实际调参直觉

把上面四个指标放在一起看，会发现一个矛盾：想要超调小，就要增大 ζ；但 ζ 增大后 ωd 减小，峰值时间和上升时间都会变长。同时调节时间 ts 与 ζωn 成反比，增大 ζ 可以减小 ts，但 ζ 过大后系统变成过阻尼，调节时间反而变长。所以欠阻尼区间内存在一个最优折中。

有一个经典结论：在 ζ ≈ 0.707 附近，系统的调节时间最短，同时超调量也不大（约4.3%），因此常作为"最佳阻尼比"。这不是绝对的，不同应用场景有不同最优值。比如：

• 数控机床进给系统：要求超调极低，甚至无超调，以免撞刀，所以 ζ 往往取0.8~1.0
• 雷达伺服系统：允许一定超调，但要求响应快，ζ 可能取0.5~0.6
• 电源变换器的电压环：超调会造成过压，ζ 取0.7以上更安全
• 音频系统或测量仪表：更关心无超调和无振铃，往往按临界阻尼或过阻尼去设计

理解了指标耦合，你就不会盲目追求某一个指标，而是结合工况做权衡。

3.3 一个完整算例：给定 ζ 和 ωn 计算全部指标

假设系统参数 ζ = 0.5，ωn = 10 rad/s，我们来算一遍：

• 阻尼振荡频率：ωd = 10 × √(1 - 0.25) = 10 × 0.866 = 8.66 rad/s
• 峰值时间：tp = π / 8.66 = 0.363 s
• 相角：β = arctan(8.66 / 5) = arctan(1.732) = 60°
• 上升时间：tr = (π - 1.047) / 8.66 = 0.242 s
• 超调量：Mp = e^(-0.5π / 0.866) = e^(-1.814) = 16.3%
• 衰减系数：σ = 5
• 5%调节时间：ts ≈ 3 / 5 = 0.6 s
• 2%调节时间：ts ≈ 4 / 5 = 0.8 s

然后你可以在这组参数基础上，试着把 ζ 提高到0.7，超调量会降到4.6%，但 ¥调节时间变短了吗？算一下：ωd = 10 × √(1-0.49)=7.14，σ=7，ts(2%)=4/7=0.571s，确实更短。这就说明了为什么0.707附近是最佳区域。多算几组，系统的参数变化手感就建立起来了。

4. 仿真验证：从传递函数到响应曲线

4.1 MATLAB/Simulink 中的快速仿真方法

仿真不是目的，是用来验证理论推导是否正确、观察参数变化趋势的手段。最直接的方法是用 MATLAB 的step函数。

% 欠阻尼二阶系统阶跃响应仿真 zeta = 0.5; wn = 10; sys = tf([wn^2], [1, 2*zeta*wn, wn^2]); % 获取阶跃响应数据 [y, t] = step(sys); % 计算指标 % 超调量 [peak, peak_idx] = max(y); Mp = (peak - 1) * 100; % 稳态值为1 % 峰值时间 tp = t(peak_idx); % 2% 调节时间（从响应进入误差带后不再越界算起） idx = find(abs(y - 1) <= 0.02*1, 1, 'first'); % 注意：这里首次进入误差带可能不准确，需要往后找最长连续满足条件的起点 ts = t(idx); % 这里简化，实际工程更严谨 fprintf('超调量: %.2f%%\n', Mp); fprintf('峰值时间: %.4f s\n', tp); fprintf('调节时间: %.4f s\n', ts); plot(t, y, 'b', 'LineWidth', 1.5); grid on; xlabel('时间 (s)'); ylabel('输出量'); title(['欠阻尼二阶系统阶跃响应 ζ=', num2str(zeta), ', ωn=', num2str(wn)]);

上面的代码会输出响应曲线，并打印出超调量和峰值时间。一个容易踩的坑是find找到的 idx 只是首次进入误差带的时间点，如果曲线在此之前振荡穿越了误差带，那就不准确。更严谨的方法是从后往前扫描，找到最后一个越界点的位置再加一小段。

4.2 Python 仿真与自定义绘图

如果你没有 MATLAB，用 Python 的 control 库也很方便。一个示例：

import numpy as np import matplotlib.pyplot as plt from control import tf, step_response zeta = 0.5 wn = 10.0 num = [wn**2] den = [1, 2*zeta*wn, wn**2] sys = tf(num, den) T, yout = step_response(sys, T=np.linspace(0, 2, 1000)) # 绘制响应曲线 plt.figure(figsize=(10, 6)) plt.plot(T, yout, 'b', linewidth=1.5) plt.axhline(1.0, color='gray', linestyle='--', linewidth=1) plt.axhline(1.0 + 0.02, color='r', linestyle=':', linewidth=1) plt.axhline(1.0 - 0.02, color='r', linestyle=':', linewidth=1) plt.xlabel('Time (s)') plt.ylabel('Output') plt.title(f'Unit Step Response (zeta={zeta}, wn={wn})') plt.grid(True) plt.xlim([0, 2]) plt.show()

运行这段代码，你会看到一条在稳态值1附近振荡衰减的曲线，峰值大约1.163，与理论计算的16.3%超调吻合。我建议你把 zeta 设成0.2、0.5、0.707、1.0、1.5分别跑一遍，做成一张对比图，这样对"阻尼比如何影响响应"会有非常直观的记忆。

4.3 仿真中不能忽视的数值陷阱

第一个陷阱是仿真时间的范围。如果你用的是固定步长求解器，步长过大可能导致峰值点采样不到，算出的超调量偏小。我在调试一段仿真脚本时遇到过：步长0.05s，ωn=20 rad/s，系统自振周期0.314s，一个周期只采6个点，峰值几乎总是落在采样间隔之间，超调量计算明显失真。解决方法是把仿真步长缩短到自振周期的1/20以下，或者使用自适应步长的求解器。

第二个陷阱是稳态值的归一化处理。有的传递函数增益不为1，阶跃响应稳态值就不是1。实际被控对象经常有静态增益 K，你需要先把传递函数归一化，或者用dcgain(sys)求出稳态值后再算超调量。

第三个陷阱是二阶系统模型中的零点。标准二阶系统没有零点，但很多实际系统化简后带有零点，比如 PID 控制器作用下的闭环系统常常包含零点。零点会改变响应曲线的形状，尤其是上升阶段的初始斜率，导致峰值时间和超调量明显偏离无零点近似。分析时如果发现理论计算与实际仿真偏差大，先检查闭环传递函数是否有零点，再用完整的传递函数做仿真。

5. 工程场景下的欠阻尼现象与应对措施

5.1 机械伺服系统中的欠阻尼表现与调试实例

我在现场调试某个三轴机械臂的伺服驱动时，发现 Z 轴在定位停止后出现来回小幅振荡，持续了将近两秒才停稳。示波器看位置反馈曲线，就是典型的欠阻尼二阶响应：启动阶段先冲到目标位置上方，然后衰减振荡。查参数时发现速度环增益调得偏高，等效阻尼比只有0.3左右，超调量超过30%。

处理方案不是单纯降低增益，因为增益低了响应就慢。我采取了两步：第一步降低速度环比例增益，把系统阻尼比提到0.5附近，超调降到了16%，但定位时间还是偏长。第二步在前向通道中加入一个低通滤波器，滤除高频扰动信号，同时稍微增加微分项的阻尼作用，最终等效阻尼比到了0.65，超调量降到8%，调节时间缩短了30%。这种"既要又要"的矛盾就是欠阻尼分析的核心价值：你得知道当前阻尼比是多少，才知道往哪个方向调。

5.2 电气系统中的欠阻尼振荡实例

电源转换器的输出 LC 滤波器也是一个典型二阶系统。当负载突变时，输出电压会振荡。如果阻尼太小，输出过冲会打坏后级芯片。很多电源控制芯片内部会刻意补偿阻尼，或者在环路中引入"斜坡补偿"，本质就是增大等效阻尼比，降低振荡幅度。

我见过一个案例：某板卡 DC-DC 电路在满载切换时输出电压波动达到500mV，而芯片要求不超过200mV。检查控制环参数后发现，补偿网络零点频率设置过高，导致环路的等效阻尼比只有0.2。把零点频率下移后，阻尼比提高到0.4，波动降到150mV以下。这类问题的根源就是没有对闭环二阶动态特性做定量的欠阻尼指标分析。

5.3 欠阻尼条件下控制器的补偿思路

针对欠阻尼被控对象，补偿手段有几类：

• 增大系统阻尼：在控制律中加入微分（D）或速度反馈，相当于在原二阶系统中引入额外阻尼项
• 改变系统增益：适当降低开环增益，通常能增大阻尼，但也会降低响应速度，需要折中
• 前置滤波器：在输入通道加低通或平滑滤波器，减少参考突变给对象带来的冲击，但前提是对象本身具有欠阻尼特性
• 状态反馈设计：现代控制理论里的极点配置法，把极点放到你期望的位置，直接确定 ζ 和 ωn。这是最直接的"定制"方法

如果项目用单片机做控制，比例项、积分项、微分项的整定本质上就是在一个未知二阶模型里寻找合适的阻尼比和自然频率。你每次调 Kp 和 Kd，等于在设计极点的位置。理解了欠阻尼分析后，可以从"试凑"升级为"有计划的调整"。

6. 常见问题与排查技巧实录

6.1 超调量比理论值高很多，怎么回事

现象：用标准传递函数算超调量16%，实际系统测出来超调30%。

排查思路：

1. 检查模型中是否存在零点。闭环系统中，控制器零点会改变响应。带零点的二阶系统超调量通常比无零点时大，可以提高等效阻尼比来压制。
2. 检查执行机构是否存在饱和非线性。如果控制器输出在响应初期就进入饱和，等效输入不再是单位阶跃的理想线性叠加，响应会出现更大的超调。
3. 检查信号采集与滤波延迟。反馈通道的低通滤波器会引入滞后，相当于增加了系统相角延迟，使有效阻尼比降低。
4. 检查是否存在齿隙、间隙、摩擦等非线性环节，这些机械间隙会在换向时造成额外的冲击，使超调量异常增大。

6.2 调节时间算不准，误差带判断逻辑要注意

现象：脚本里find(abs(y-1)<=0.02)找出的调节时间只有实际曲线的一半，明显不对。

原因：首次进入2%误差带后，系统还可能在后面再一次越界。正确的调节时间定义是"以后不再越界"的时刻。实现方法建议从数据末尾往前扫描：

% 从后往前扫描，找到最后一个越界点 ts = 0; for ii = length(t):-1:1 if abs(y(ii) - 1) > 0.02 ts = t(ii); break; end end % 最终调节时间为 ts 之后再经过一个交点 % 但工程上可以近似取 ts + 一个采样周期

如果是 Python 也可以用np.where(np.abs(yout-1) > 0.02)[0][-1]来定位最后一次越界点，取它之后的时间作为调节时间。不要只盯着第一次进入误差带的点。

6.3 仿真曲线发散，但模型明明是稳定的

现象：自己搭的 Simulink 模型，积分步长固定为0.01，仿真时长10秒，后半段曲线开始发散。

原因：很可能是你用了连续传递函数模块，但求解器选取了非刚性固定步长，比如 ode1 或 ode3，遇到欠阻尼系统的高频振荡分量会因为截断误差累积而数值发散。对策是改用ode45或者减小步长到振荡周期的1/50以下。另外检查模型是否不小心形成了代数环，代数环也会造成数值不稳定。

6.4 快速调试速查表

现象	可能原因	快速排查方法	处理建议
响应振荡不衰减	ζ≈0 或系统有正反馈	看闭环极点实部是否为负	增加阻尼项或检查反馈极性
超调大、上升快	ζ过小	计算/测量超调量反推ζ	增大微分作用或降低增益
响应慢、无超调	ζ过大	看阶跃曲线是否逼近一阶惯性	减小阻尼、提高增益
稳态误差不为0	系统不是1型或增益有限	求系统型别和开环增益	添加积分环节
高频振铃	存在零点或噪声	用正弦扫描确认振荡频率	加前置滤波或陷波器

7. 从理论到实践的个人体会

我在做实际的控制器整定项目时，越来越觉得欠阻尼二阶分析不是一个学院派概念，而是一把尺子。每次看到响应曲线，我心里都会快速估算：振荡频率是多少、衰减快慢如何、超调多少，然后反推系统的阻尼比大概在什么范围。这个过程一旦熟练，调试效率会明显提高。

分享一个个人习惯：我在调电机速度环时，会先给一个小幅阶跃，比如额定转速的5%，观察响应曲线。如果超调在10%到20%之间，并且振荡1到2个周期就衰减完，说明系统处于健康的欠阻尼状态。如果曲线猛冲上去还来回荡好几下才稳定，那就该降低增益或增加阻尼。如果曲线慢悠悠爬上去一点超调都没有，可能在安全和精度要求高的场合可以接受，但在快速响应场合就太慢了。

另一个建议：不要只分析一个固定参数。我们可以在纸面上先画出 ζ 从0.3到0.9变化时超调量和调节时间的曲线，找到系统性能的交汇点。比如某系统要求超调≤10%，同时调节时间尽量短，那么 ζ 取0.65左右通常是最优的。然后针对这个 ζ 反推需要配置的极点和控制器参数，再用仿真确认一次，基本就八九不离十了。

最后想强调一个容易被忽视的细节：工程系统的阻尼比不是固定不变的，它随负载、温度、老化而变化。你今天把阻尼比调到0.707，明天负载加了，摩擦变了，实际阻尼比可能变成0.6。所以设计时要留出裕量，比如指标要求超调20%，你最好按超调10%去设计，用那部分裕量去容纳系统参数漂移。这也是为什么我在调试时总喜欢把曲线保存下来，隔一段时间再测一遍，看看参数漂移了多少。欠阻尼二阶系统分析教给我们的不只是那几个公式，更是一种"用振荡的视角去看待系统"的思维习惯。只要反馈控制存在，振荡就一定存在，我们要做的不是消除振荡，而是理解它的规律，把它控制在可以接受的范围里。