1. 项目概述：扰动观测器（DOB）在抗摩擦控制中的应用
在工业控制领域，摩擦力的非线性特性一直是精密运动控制的难点。传统PID控制器面对这种非线性扰动时往往表现不佳，而扰动观测器（Disturbance Observer, DOB）通过实时估计并补偿系统扰动，能显著提升控制性能。Simulink作为MATLAB的图形化建模环境，为DOB的实现提供了直观的仿真验证平台。
这个示例将展示如何从零开始构建一个完整的抗摩擦控制系统。我们会先建立包含摩擦模型的被控对象，然后设计DOB结构，最后通过Simulink模块搭建闭环系统。整个过程不仅适用于旋转运动系统，其方法论也可推广到直线电机、机械臂等高精度控制场景。
2. 系统建模与问题定义
2.1 摩擦模型建立
库仑摩擦+粘滞摩擦的组合模型最能反映大多数机械系统的特性：
code复制F_friction = Fc * sign(v) + Fv * v
其中Fc为库仑摩擦系数，Fv为粘滞摩擦系数，v为运动速度。在Simulink中可用Sign模块和Gain模块组合实现。
注意：实际系统中还存在静摩擦（Stiction）现象，当速度过零时需要特殊处理。建议采用连续近似函数避免仿真发散。
2.2 被控对象建模
以典型的直流电机为例，其传递函数可表示为：
code复制P(s) = 1/(Js + B)
J为转动惯量，B为阻尼系数。在Simulink中可直接用Transfer Fcn模块实现。将摩擦模型作为外部力矩扰动接入系统。
2.3 控制难点分析
摩擦力的主要挑战在于：
- 速度过零时的非线性突变
- 参数随温度、磨损等时变
- 传统线性控制器无法完全补偿
通过DOB实时估计这些扰动并前馈补偿，可大幅降低对反馈控制器的性能要求。
3. 扰动观测器设计与实现
3.1 DOB基本原理
DOB的核心思想是通过逆模型Q滤波器分离出系统扰动：
code复制d_hat = Q(s)P^-1(s)(u - y)
其中Q(s)通常选择低通滤波器，截止频率决定观测带宽。
3.2 Simulink实现步骤
-
构建名义模型：
在Library Browser中添加Transfer Fcn模块，参数设为1/(J_nom*s + B_nom) -
设计Q滤波器：
使用Analog Filter Design模块，选择Butterworth低通型matlab复制[b,a] = butter(2, 100/(2*pi)); % 二阶100Hz截止 -
搭建观测回路：
- 将控制输入u和实际输出y接入减法器
- 通过名义逆模型后接Q滤波器
- 输出扰动估计d_hat
-
前馈补偿：
将d_hat直接叠加到控制信号u上
3.3 参数调试技巧
- Q滤波器带宽应覆盖扰动主要频率，但不超过系统带宽的1/3
- 名义模型参数误差控制在20%以内
- 初始调试时可暂时关闭主控制器，单独测试DOB补偿效果
4. 完整控制系统集成
4.1 闭环结构设计
code复制[Reference] --> [PID Controller] --+--> [Plant] --> [Output]
^ |
| v
[DOB] <--[Q Filter]
4.2 Simulink模块选型建议
-
执行器饱和处理：
在PID输出端添加Saturation模块，设置合理的电压限制 -
抗积分饱和：
配置PID模块的Anti-windup模式为"back-calculation" -
信号监测：
使用To Workspace模块记录关键信号便于后期分析
4.3 仿真参数设置
matlab复制simout = sim('DOB_Control.slx', 'StopTime', '10', ...
'Solver', 'ode45', 'MaxStep', '0.001');
5. 典型问题排查与优化
5.1 高频振荡问题
现象：系统出现未预期的抖动
解决方法：
- 降低Q滤波器带宽
- 检查执行器饱和情况
- 在DOB输出端加入小的死区
5.2 稳态误差问题
现象：最终存在固定偏差
排查步骤：
- 确认库仑摩擦系数Fc设置准确
- 检查DOB输出是否达到限幅值
- 适当提高PID的积分增益
5.3 实时性优化
对于需要代码生成的实时系统：
- 将Q滤波器转换为离散形式
matlab复制[bz,az] = bilinear(b,a,1/Ts); - 使用Coder工具生成优化代码
- 在Speedgoat等实时目标机上验证
6. 进阶应用与扩展
6.1 自适应DOB设计
当系统参数变化较大时，可在线更新名义模型：
matlab复制function UpdateModel(J_new, B_new)
set_param('DOB_Control/Plant','Denominator',['[',num2str(J_new),' ',num2str(B_new),']']);
end
6.2 多自由度系统扩展
对于机械臂等多变量系统：
- 设计对角化的Q滤波器矩阵
- 考虑各轴间的耦合扰动
- 使用MATLAB的Robust Control Toolbox辅助设计
6.3 实验验证建议
- 先进行开环阶跃测试辨识真实摩擦参数
- 逐步增加DOB带宽观察效果
- 记录不同工况下的跟踪误差统计量
我在实际项目中发现，DOB配合前馈控制在纳米定位平台上可将跟踪误差降低60%以上。关键是要确保名义模型的结构正确，参数精度反而可以放宽。一个实用的技巧是在系统启动时自动进行参数辨识，定期更新模型参数以适应设备老化。