如何用MATLAB计算系统稳态误差？

Question

在使用MATLAB分析控制系统时，常需计算系统对阶跃、斜坡或加速度输入的稳态误差。一个常见问题是：**如何通过MATLAB准确计算不同输入类型下系统的稳态误差？** 特别是当系统包含多个环节（如前向通路增益、反馈元件）时，用户难以确定误差传递函数，或混淆开环与闭环结构下的静态误差系数（Kp、Kv、Ka）应用条件。此外，在利用`step()`函数仿真后，如何从响应曲线中精确提取稳态误差值也常引发困惑。例如，为何理论计算的误差与仿真结果存在偏差？这通常涉及系统类型判断错误或未考虑扰动影响。

Answer 1

如何通过MATLAB准确计算控制系统在不同输入下的稳态误差

1. 稳态误差的基本概念与系统类型判定

在控制理论中，稳态误差（Steady-State Error, SSE）是指系统响应进入稳态后，输出与期望输入之间的差值。对于单位反馈系统，其误差信号 e(t) 定义为输入 r(t) 与输出 c(t) 的差：
e(t) = r(t) - c(t)
稳态误差可通过终值定理计算：

$$ e_{ss} = \lim_{s \to 0} sE(s) = \lim_{s \to 0} s \cdot \frac{R(s)}{1 + G(s)H(s)} $$

其中，G(s)H(s) 为开环传递函数。系统根据开环传递函数中积分环节的个数分为0型、I型、II型等，直接影响其对阶跃、斜坡、加速度输入的跟踪能力。

系统类型	Kp	Kv	Ka	阶跃误差	斜坡误差	加速度误差
0型	有限	0	0	1/(1+Kp)	∞	∞
I型	∞	有限	0	0	1/Kv	∞
II型	∞	∞	有限	0	0	1/Ka

2. 静态误差系数的MATLAB计算方法

静态误差系数包括位置常数 Kp、速度常数 Kv 和加速度常数 Ka，分别对应阶跃、斜坡和抛物线输入。在MATLAB中，可通过dcgain()函数直接计算这些系数。

% 定义开环传递函数 G(s)H(s)
s = tf('s');
G = 10 / (s * (s + 2));  % 示例：I型系统
H = 1;

% 计算静态误差系数
Kp = dcgain(G * H);           % 位置常数
Kv = dcgain(s * G * H);       % 速度常数
Ka = dcgain(s^2 * G * H);     % 加速度常数

fprintf('Kp = %.2f, Kv = %.2f, Ka = %.2f\n', Kp, Kv, Ka);

注意：使用s*G*H时需确保系统在原点处无过多极点导致极限发散。若dcgain返回Inf或NaN，应结合系统类型判断是否合理。

3. 构建误差传递函数并仿真验证

当系统结构复杂（如非单位反馈、前向增益可调）时，必须显式构造误差传递函数 E(s)/R(s)：

$$ \frac{E(s)}{R(s)} = \frac{1}{1 + G(s)H(s)} $$

% 构造闭环误差传递函数
Gr = feedback(G*H, 1, 'negative');  % 闭环传递函数 C/R
EoR = 1 - Gr;                      % 误差传递函数 E/R

% 施加阶跃输入并仿真
t = 0:0.01:10;
[y,t] = step(EoR, t);              % e_ss(t) 响应
e_ss_sim = y(end);                 % 提取稳态误差
fprintf('Simulated steady-state error: %.4f\n', e_ss_sim);

4. 使用step()函数提取稳态误差的注意事项

虽然step()能直观展示响应曲线，但直接取末尾值可能因仿真时间不足或数值精度问题产生偏差。建议采用以下策略提高准确性：

• 延长仿真时间以确保真正进入稳态
• 使用final()函数获取稳态值
• 结合lsim()处理任意输入信号（如斜坡）

% 对斜坡输入进行仿真
ramp_input = t;
[y_sys,~] = lsim(tf(Gr), ramp_input, t);
e_ramp = ramp_input - y_sys;
e_ss_ramp = e_ramp(end);

5. 常见误差来源分析与调试流程图

理论与仿真结果不一致的常见原因包括：

1. 误判系统类型（如未识别H(s)中的零点/极点）
2. 忽略了扰动输入的影响
3. 反馈路径非单位增益但未修正误差模型
4. 数值积分误差或采样率不足

graph TD A[开始分析稳态误差] --> B{系统是否为单位反馈?} B -- 是 --> C[构建G(s)H(s)] B -- 否 --> D[修正为一般反馈结构] D --> C C --> E[确定系统类型] E --> F[计算Kp/Kv/Ka] F --> G[理论e_ss预测] G --> H[建立误差传递函数E(s)/R(s)] H --> I[使用step()/lsim()仿真] I --> J[提取末值并对比] J --> K{理论≈仿真?} K -- 否 --> L[检查H(s)、扰动、积分时间] L --> I K -- 是 --> M[完成分析]