控制系统性能分析 2 种视角：从时域超调量到频域稳定裕度的 5 步映射

控制系统性能分析 2 种视角：从时域超调量到频域稳定裕度的 5 步映射

在工程实践中，控制系统的性能评估往往需要同时考虑时域和频域两个维度的指标。时域指标如超调量、上升时间直观反映了系统的动态响应特性，而频域指标如相位裕度、增益裕度则揭示了系统的稳定边界。本文将构建这两种视角之间的桥梁，通过具体案例演示如何实现指标间的相互转换与映射。

1. 时域与频域性能指标的核心关联

控制系统设计本质上是在响应速度与稳定性之间寻找平衡。时域分析通过阶跃响应曲线直接呈现这一矛盾关系：

• 超调量（Overshoot）：反映系统阻尼特性，与频域中的相位裕度密切相关
• 上升时间（Rise Time）：表征系统快速性，对应频域带宽指标
• 稳态误差（Steady-state Error）：体现系统精度，与开环增益直接相关

典型二阶系统的时频域指标存在精确的数学对应关系。对于传递函数为 $G(s)=\frac{\omega_n^2}{s^2+2\zeta\omega_n s+\omega_n^2}$ 的系统：

提示：上述关系式在ζ=0.4~0.7范围内具有较高精度，对于高阶系统需进行适当修正

2. 从时域响应推导频域指标的5步流程

2.1 采集时域响应数据

通过实验或仿真获取系统的阶跃响应曲线，测量以下关键参数：

% MATLAB示例：获取阶跃响应数据 sys = tf([25],[1 6 25]); % 二阶系统示例 [y,t] = step(sys); stepinfo(y,t) % 自动计算时域指标

2.2 计算阻尼比与自然频率

利用超调量σ%反推阻尼比ζ： $$ \zeta = \frac{-\ln(\sigma/100)}{\sqrt{\pi^2+\ln^2(\sigma/100)}} $$

通过峰值时间Tp计算自然频率ωn： $$ ω_n = \frac{\pi}{T_p\sqrt{1-\zeta^2}} $$

2.3 构建传递函数模型

根据得到的ζ和ωn建立标准二阶系统模型： $$ G(s) = \frac{ω_n^2}{s^2+2\zetaω_n s+ω_n^2} $$

2.4 转换频域特征参数

计算关键频域指标：

• 相位裕度：$PM = \tan^{-1}\left(\frac{2\zeta}{\sqrt{\sqrt{1+4\zeta^4}-2\zeta^2}}\right)$
• 带宽频率：$ω_b = ω_n\sqrt{1-2\zeta^2+\sqrt{2-4\zeta^2+4\zeta^4}}$
• 谐振峰值：$M_r = \frac{1}{2\zeta\sqrt{1-\zeta^2}}$ (当ζ<0.707)

2.5 验证与修正

将理论计算结果与实际频域测量值对比，典型验证方法包括：

1. 伯德图法：绘制理论传递函数的频响曲线
2. 正弦扫频实验：实测系统在不同频率下的响应
3. 参数辨识：使用系统辨识工具包验证模型准确性

3. 频域到时域的逆向映射方法

3.1 从伯德图提取关键参数

通过开环频率特性曲线获取：

• 增益穿越频率ωc：幅值曲线穿越0dB时的频率
• 相位裕度PM：ωc处相位与-180°的差值
• 增益裕度GM：相位达到-180°时的增益余量

3.2 估算时域性能指标

对于相位裕度PM在30°~60°范围内的系统：

1. 超调量估算： $$ \sigma% \approx 0.01PM^2 - 0.9PM + 20 \quad (\text{经验公式}) $$

2. 调节时间预测： $$ T_s \approx \frac{4}{ω_c}\left(\frac{100}{PM}\right)^{1.3} $$

3. 上升时间关联： $$ T_r \approx \frac{1.8}{ω_c} $$

3.3 高阶系统的近似处理

对于非标准二阶系统，可采用以下策略：

• 主导极点法：保留影响最大的极点对
• 模型降阶：使用平衡截断等降阶技术
• 数值优化：基于频域数据拟合等效二阶模型

4. 工程应用中的决策指南

4.1 优先使用时域分析的场景

• 需要直观评估动态响应特性的场合
• 系统阶数较低（n≤3）且线性度良好
• 关注瞬态过程的具体形态（如机器人轨迹跟踪）

4.2 优先采用频域分析的场景

• 系统存在明显噪声或高频扰动
• 需要分析稳定裕度的设计问题
• 多变量系统的解耦设计（如飞行器控制）

4.3 混合分析方法的最佳实践

1. 初始设计阶段：用时域指标确定基本需求
2. 详细设计阶段：通过频域分析优化稳定裕度
3. 验证阶段：时域仿真检验最终性能

5. 典型应用案例：伺服系统调参

某位置伺服系统的设计要求：

• 时域要求：超调量<15%，上升时间<0.2s
• 频域要求：相位裕度>45°，带宽>10Hz

调参步骤：

1. 根据时域要求推算目标阻尼比： $$ \zeta = \frac{-\ln(0.15)}{\sqrt{\pi^2+\ln^2(0.15)}} \approx 0.52 $$

2. 计算需要的自然频率： $$ ω_n = \frac{1.8}{T_r\sqrt{1-\zeta^2}} = \frac{1.8}{0.2×0.854} \approx 10.5 \text{rad/s} $$

3. 设计PID参数满足频域指标：

   # Python控制库示例 import control as ct sys = ct.tf([25], [1, 6, 25]) # 初始系统 pid = ct.tf([0.5, 2.5], [1, 0]) # PID控制器 ol_sys = pid * sys # 开环系统 ct.bode_plot(ol_sys) # 验证相位裕度

4. 时频域联合验证：

   • 时域：阶跃响应超调量12.8%，上升时间0.18s
   • 频域：相位裕度48°，带宽10.3Hz

实际调试中发现，当负载惯量增加30%时，需将微分增益提高15%以维持原有相位裕度。这种时频域参数的动态调整正是工程实践中的关键技巧。