3. 扰动与不确定性建模：外部负载扰动、摩擦、参数摄动、未建模动态的数学描述与分类
运动控制系统在实际运行中几乎不可能处于理想的无扰动环境。无论是机床进给轴的切削力突变、机器人关节的负载变化，还是伺服电机转子惯量随负载工装改变，这些因素都会使系统实际输出偏离设计预期。扰动观测器（DOB）的核心价值，正是通过估计并补偿这些扰动来提升系统的鲁棒性。然而，若对扰动本身缺乏清晰的数学描述与分类，观测器的设计便无从谈起。本章将系统梳理运动控制中四类主要扰动的数学本质、建模方法及其在观测器设计中的处理策略，为后续章节的观测器结构设计与参数整定奠定基础。
上图展示了运动控制系统中四类扰动的整体框架及其在等效总扰动中的叠加关系。需要特别强调的是，扰动观测器所估计的“等效扰动”并非某一物理量的直接测量值，而是将参数摄动、未建模动态折算到控制输入端的等效信号。这种折算思想是理解DOB工作原理的关键。下面分别对四类扰动展开深入分析。
3.1 外部负载扰动：时域特征与频域描述
3.1.1 核心原理：负载扰动的物理来源与数学表达
外部负载扰动是指由系统外部施加在运动轴上的力矩或力，其特点是不依赖于系统自身状态（如速度、位置），而是由工作环境决定。在机床进给系统中，切削力的突变是典型的外部负载扰动；在机器人关节中，末端执行器抓取工件瞬间的力矩变化同样属于此类。从数学角度看，外部负载扰动可表示为时间函数 d_load(t)，其拉普拉斯变换为 D_load(s)。在运动控制的标准模型中，被控对象的动力学方程通常写为：
J·dω/dt = T_m - T_load - T_f - B·ω其中 T_m 为电机电磁转矩，T_load 为外部负载转矩，J 为转动惯量，B 为粘性摩擦系数。在DOB设计中，外部负载扰动通常被视为低频信号或阶跃信号，因为大多数实际负载变化（如切削深度变化、工件质量变化）的带宽远低于控制系统的截止频率。
3.1.2 判断方法：负载扰动的辨识与分类
工程实践中，判断外部负载扰动的主要依据是：扰动是否与系统状态无关。具体操作时，可通过以下步骤进行辨识：
- 在系统稳定运行（匀速或静止）时，记录电机转矩指令的稳态值；
- 改变负载条件（如增加切削深度），观察转矩指令的变化量；
- 若转矩指令变化与速度变化无直接关联，且响应时间远大于电气时间常数，则可判定为外部负载扰动。
从频域角度看，外部负载扰动通常集中在低频段（<10 Hz），部分冲击性负载（如急停、碰撞）可能包含较宽的频谱分量，但能量仍集中在低频。因此，DOB中通常采用低通滤波器（Q滤波器）来抑制高频噪声，同时保留低频扰动估计能力。
3.1.3 案例拆解：数控机床切削力突变补偿
某立式加工中心在铣削过程中，当刀具切入工件深度从2mm增至4mm时，主轴负载转矩在约50ms内从5N·m跃升至12N·m。若不进行补偿，速度波动可达±30 r/min，影响加工表面质量。采用DOB后，观测器在10ms内估计出负载变化，并通过前馈通道补偿，速度波动降至±5 r/min以内。关键参数设置如下：
| 参数 | 数值 | 说明 |
|---|---|---|
| Q滤波器截止频率 | 200 Hz | 覆盖负载突变带宽 |
| 观测器增益 | 500 | 保证快速收敛 |
| 补偿限幅 | ±15 N·m | 防止过补偿 |
复盘经验：负载扰动补偿的关键在于Q滤波器的带宽选择。带宽过高会引入测量噪声，过低则无法快速跟踪负载变化。本例中200Hz的截止频率是在噪声水平和响应速度之间的折中。
3.2 摩擦扰动：非线性建模与补偿策略
3.2.1 核心原理：摩擦的物理机制与数学模型
摩擦是运动控制中最常见的非线性扰动，其来源包括轴承滚动摩擦、导轨滑动摩擦、电机换向力矩波动等。摩擦模型经历了从简单到复杂的发展：经典库仑+粘性模型 T_f = T_c·sign(ω) + B·ω 仅适用于中高速段；Stribeck模型进一步描述了低速区的负斜率特性，即摩擦力随速度增加先下降后上升。更精确的LuGre模型则引入了鬃毛形变状态变量，能够描述预滑动位移、摩擦滞后等微观现象。在DOB设计中，摩擦通常被视为状态相关的扰动，即 d_f(ω) = T_f(ω)，其特点是：
- 在速度过零时产生不连续的跳变（库仑摩擦）；
- 在低速区呈现强烈的非线性（Stribeck效应）；
- 在高速区近似线性（粘性摩擦主导）。
3.2.2 判断方法：摩擦参数的实验辨识
摩擦模型的参数辨识通常采用恒速跟踪法：让系统以不同恒定速度运行，记录稳态转矩指令，绘制“速度-转矩”曲线。具体流程为：
- 设定速度序列（如0.1、0.5、1、5、10、50、100 r/min）；
- 对每个速度点，等待系统达到稳态后记录平均转矩；
- 将数据点拟合到Stribeck模型，得到 T_c、T_s、ω_s、B 等参数。
值得注意的是，摩擦参数会随温度、磨损程度变化，因此在线辨识或自适应补偿是工程中的常见需求。DOB可以自动估计摩擦扰动，但若摩擦占主导且变化缓慢，可考虑在DOB前增加前馈摩擦补偿，以减轻观测器负担。
3.2.3 案例拆解：高精度转台低速爬行抑制
某精密转台在0.01 r/min低速运行时出现明显的“爬行”现象，位置波动达±0.002°。经辨识，其库仑摩擦力矩 T_c = 0.8 N·m，Stribeck速度 ω_s = 0.5 r/min。单纯依靠DOB补偿时，由于摩擦在速度过零时的跳变特性，观测器估计值出现振荡。解决方案是采用“前馈+反馈”复合结构：
T_ff = T_c·tanh(ω/ω_s) + B·ω // 前馈摩擦模型
T_dob = Q(s)·(P⁻¹(s)·y - u) // DOB补偿残余扰动实施后，转台位置波动降至±0.0003°，爬行现象完全消除。复盘经验：摩擦前馈模型不需要非常精确，只需补偿主要非线性分量，剩余误差由DOB处理，这种分工可显著提高系统鲁棒性。
3.3 参数摄动：不确定性建模与鲁棒性分析
3.3.1 核心原理：参数摄动的来源与影响
参数摄动是指被控对象模型参数（如转动惯量J、粘性系数B、电气时间常数τ_e）在实际运行中偏离标称值。典型场景包括：机器人负载变化导致关节等效惯量成倍增加；温度变化导致绕组电阻改变；机械磨损导致阻尼系数下降。参数摄动的数学描述通常采用乘性不确定性形式：
P(s) = P₀(s)·(1 + Δ(s)), |Δ(jω)| ≤ W(ω)其中 P₀(s) 为标称模型，Δ(s) 为归一化不确定性，W(ω) 为不确定性上界。在DOB设计中，参数摄动的影响体现在：实际对象与标称模型的逆 P₀⁻¹(s) 不匹配，导致扰动估计出现偏差。若摄动过大，可能引起观测器不稳定。
3.3.2 判断方法：参数摄动的量化与鲁棒性校验
工程中常用以下步骤量化参数摄动：
- 确定参数变化范围：如 J ∈ [J_min, J_max]，B ∈ [B_min, B_max]；
- 计算不确定性上界 W(ω) = max|P(jω)/P₀(jω) - 1|；
- 校验鲁棒稳定性条件：|Q(jω)·W(ω)| < 1 对所有频率成立。
该条件表明，Q滤波器的幅频特性必须低于不确定性上界的倒数。因此，参数摄动越大，Q滤波器的带宽就必须越低，这直接限制了DOB的扰动抑制带宽。这是DOB设计中最基本的权衡关系。
3.3.3 案例拆解：机器人关节惯量变化下的DOB设计
某六轴工业机器人第二关节在空载与满载（抓取10kg工件）时，等效惯量从0.5 kg·m²变化至1.2 kg·m²，变化幅度达140%。若按标称惯量设计DOB，满载时观测器增益过高导致系统振荡。解决方案是采用增益调度策略：
| 负载状态 | 惯量估计值 | Q滤波器带宽 | 观测器增益 |
|---|---|---|---|
| 空载 | 0.5 kg·m² | 300 Hz | 800 |
| 半载 | 0.8 kg·m² | 200 Hz | 500 |
| 满载 | 1.2 kg·m² | 120 Hz | 300 |
通过在线辨识负载惯量（利用加速度与转矩的相关性），自动切换DOB参数，系统在全负载范围内保持稳定且性能一致。复盘经验：参数摄动是DOB设计中必须优先考虑的因素，建议在设计初期就进行不确定性分析，而不是等到调试阶段才处理。
3.4 未建模动态：高频特性与鲁棒性边界
3.4.1 核心原理：未建模动态的本质与影响
未建模动态是指实际系统中存在但未包含在标称模型中的高频特性，主要包括：机械谐振（传动轴柔性、联轴器弹性）、电流环延迟（PWM更新周期、采样延迟）、传感器滤波特性等。这些动态的数学描述通常为乘性不确定性 Δ(s)，其幅值在高频段显著增大。例如，双质量系统的谐振特性可表示为：
P(s) = P₀(s) · (ω_r²)/(s² + 2ζ_r·ω_r·s + ω_r²)其中 ω_r 为谐振频率，ζ_r 为阻尼比。当DOB的Q滤波器带宽接近或超过谐振频率时，观测器可能激励谐振模态，导致系统振荡甚至失稳。
3.4.2 判断方法：未建模动态的识别与边界确定
识别未建模动态的常用方法是频率响应测试：
- 使用扫频信号激励系统，测量输入到输出的Bode图；
- 对比实测频率响应与标称模型的差异，确定谐振峰、延迟时间等；
- 将差异拟合为乘性不确定性 Δ(s)，并计算其上界 W(ω)。
工程经验表明，大多数运动控制系统的未建模动态集中在1 kHz以上，而DOB的Q滤波器截止频率通常设置在100-500 Hz之间。因此，只要Q滤波器设计得当，未建模动态的影响可以被有效抑制。但若系统存在低频谐振（如大型机械臂的柔性模态在20-50 Hz），则必须特别处理。
3.4.3 案例拆解：伺服系统机械谐振的DOB规避
某伺服驱动系统在调试中发现，当DOB的Q滤波器截止频率设为500 Hz时，系统在300 Hz附近出现强烈振荡。经频谱分析，发现机械传动轴存在280 Hz的谐振峰。解决方案是：
1. 将Q滤波器截止频率降至150 Hz，避开谐振区；
2. 在DOB前增加陷波滤波器（Notch Filter），中心频率280 Hz，带宽20 Hz；
3. 重新整定观测器增益，保证相位裕度大于45°。调整后，系统稳定运行，且对200 Hz以下扰动的抑制能力仅下降约15%。复盘经验：未建模动态的处理原则是“避让优先，抑制其次”。若谐振频率与所需扰动抑制带宽重叠，则需考虑机械结构改进或采用更复杂的控制策略（如H∞控制）。
3.5 扰动分类总结与建模流程
综合上述分析，运动控制中的扰动可按以下维度分类：按与状态的关系分为外部扰动（与状态无关）和内部扰动（与状态相关）；按频率特性分为低频扰动（负载变化）、中频扰动（摩擦、参数摄动）和高频扰动（未建模动态）；按确定性分为确定性扰动（可建模）和随机扰动（噪声）。在实际DOB设计中，建议遵循以下建模流程：
| 步骤 | 内容 | 输出 |
|---|---|---|
| 1 | 建立标称模型 P₀(s) | 传递函数、状态空间 |
| 2 | 辨识外部负载扰动特性 | 幅值、带宽、变化率 |
| 3 | 辨识摩擦参数 | T_c、T_s、ω_s、B |
| 4 | 量化参数摄动范围 | ΔJ、ΔB、不确定性上界 |
| 5 | 测量未建模动态 | 谐振频率、延迟时间 |
| 6 | 综合等效扰动模型 | d_total 的时域/频域描述 |
完成上述建模后，设计者即可根据扰动特性选择合适的DOB结构（如线性DOB、非线性DOB、自适应DOB）并确定Q滤波器参数。需要强调的是，扰动建模并非一次性工作，而是贯穿系统全生命周期的持续过程——随着设备老化、工况变化，扰动特性会改变，观测器参数需要相应调整。