1. 项目概述
在工业控制领域，PI控制器是最基础也最经典的控制算法之一。不同于教科书上理想化的理论推导，实际工程中基于传递函数设计PI控制器时，需要综合考虑系统动态特性、稳定性裕度和实际工况限制。这次我将分享一个从传递函数出发，完整设计并实现PI控制器的实战案例，包含频域分析、参数整定和实际调试的全过程。
这个项目特别适合两类读者：一是刚接触自动控制的学生，可以通过实际案例理解抽象理论；二是需要快速实现控制方案的工程师，文中的参数计算方法和调试技巧都能直接套用。我们将使用MATLAB/Simulink作为主要工具，但所有原理和步骤同样适用于其他编程环境。
2. 核心原理与传递函数分析
2.1 被控对象建模
我们以一个直流电机转速控制系统为例，其传递函数可表示为：
code复制G(s) = K / (τs + 1)
其中K=0.8为静态增益，τ=0.5秒为时间常数。在MATLAB中建立该模型的代码如下：
matlab复制s = tf('s');
G = 0.8 / (0.5*s + 1);
bode(G); grid on
关键提示：实际系统中，传递函数参数需要通过阶跃响应实验或频域辨识获得。建议至少采集3组不同工况的数据取平均值。
2.2 PI控制器结构
标准PI控制器的传递函数为：
code复制C(s) = Kp + Ki/s = Kp(1 + 1/(Ti*s))
包含两个可调参数：
- 比例增益Kp：决定系统响应速度
- 积分时间Ti：消除稳态误差（Ki=Kp/Ti）
其频域特性表现为：
- 低频段：积分作用主导，提供高增益保证跟踪精度
- 中频段：比例作用主导，影响相位裕度
- 高频段：增益衰减，抑制噪声
3. 参数整定方法与实现
3.1 频域设计法
我们采用相位裕度法进行设计，目标相位裕度设为60°：
- 绘制未补偿系统的Bode图：
matlab复制margin(G)
- 找到相位穿越频率ωgc（相位为-120°时的频率）
- 计算需要提升的幅值：|G(jωgc)| ≈ 0.45 → Kp=1/0.45≈2.2
- 设置积分时间Ti=10/ωgc≈6.3s
最终控制器：
matlab复制Kp = 2.2; Ti = 6.3;
C = Kp * (1 + 1/(Ti*s))
3.2 时域验证
闭环系统阶跃响应仿真：
matlab复制sys_cl = feedback(C*G, 1);
step(sys_cl); grid on
关键指标应满足：
- 超调量<10%
- 调节时间<3τ=1.5秒
- 稳态误差<1%
实测技巧：如果响应振荡剧烈，先降低Kp；如果有稳态误差，减小Ti（但不能小于系统主要时间常数的1/5）
4. 实际调试中的问题解决
4.1 积分饱和现象
当执行机构达到限幅值时，积分项会持续累积导致控制失效。解决方法：
- 抗饱和处理：检测到饱和时暂停积分
- 积分分离：误差较大时仅用比例控制
matlab复制% 伪代码实现
if abs(error) > threshold
u = Kp * error;
else
u = Kp*error + Ki*integral(error);
end
4.2 噪声敏感问题
高频测量噪声会被积分环节放大，导致控制量抖动。对策：
- 在反馈通道增加一阶低通滤波器：
matlab复制F = 1/(0.1*s + 1); - 改用不完全积分（伪微分反馈）：
matlab复制C = Kp*(1 + 1/(Ti*s)) / (Tf*s + 1); % Tf通常取Ti的1/10~1/5
4.3 参数自整定技巧
对于时变系统，可采用在线调整策略：
- 继电反馈法：自动获取临界增益和周期
- 梯度下降法：根据性能指标实时优化
- 模糊规则：经验丰富的工程师可设置20-30条调整规则
5. 进阶优化方向
5.1 两自由度结构
将命令跟踪与干扰抑制分开设计：
matlab复制% 前馈控制器
Cf = (Ti*s + 1)/((Ti/Kp)*s + 1);
% 反馈控制器仍用PI
5.2 增益调度控制
针对不同工作点设计多组参数：
matlab复制% 根据工作点选择参数
if operating_mode == 1
Kp = 2.2; Ti = 6.3;
else
Kp = 1.8; Ti = 8.0;
end
5.3 数字实现要点
离散化时推荐采用双线性变换：
matlab复制C_d = c2d(C, Ts, 'tustin');
采样周期Ts应满足：
- 小于系统最小时间常数的1/10
- 大于计算周期（考虑实时性）
我在实际项目中验证过，当Ts=0.01秒时，数字控制效果几乎与连续系统一致。但要注意避免整定参数时出现数值溢出问题，特别是在低精度嵌入式设备上。