1. 伺服系统低速摩擦问题本质解析
当伺服电机运行在每分钟几转的低速工况时，工程师们常会观察到速度曲线出现明显的"平台期"——就像汽车在雪地打滑，油门踩下去却不见车速提升。这种现象的本质是静摩擦（Stiction）与库伦摩擦（Coulomb friction）在速度过零点附近形成的非线性死区。传统PID控制在这个区域会表现出明显的极限环振荡，就像新手司机在拥堵路段不断踩放油门导致的车辆顿挫。
摩擦力的非线性特性在数学上可表示为分段函数：
code复制F_friction = {
F_c * sign(v) + σ2*v, |v| > v_threshold (运动状态)
min(F_ext, F_s) * sign(F_ext), |v| = 0 (静止状态)
}
其中F_c是库伦摩擦系数，F_s是最大静摩擦力（通常F_s≈1.5F_c），σ2是粘滞摩擦系数，v_threshold是静动摩擦转换的临界速度（通常0.1-1rpm）。这个数学模型揭示了为什么普通线性控制器难以应对——在速度过零区域，系统动态特性发生了突变。
2. LuGre模型深度拆解
相比经典摩擦模型，LuGre模型的先进性在于引入了虚拟鬃毛变形量z这个内部状态变量，其微分方程为：
code复制dz/dt = v - σ0*|v|/g(v)*z
g(v) = F_c + (F_s - F_c)*exp(-(v/v_s)^2)
摩擦力 = σ0*z + σ1*dz/dt + σ2*v
这个模型巧妙地将微观接触面的弹性形变（σ0项）、形变阻尼（σ1项）与宏观粘滞摩擦（σ2项）统一起来。参数v_s控制静动摩擦转换的平滑度，实测表明v_s=0.1mm/s时对多数金属接触面拟合效果最佳。
在Matlab中实现时需特别注意数值稳定性问题：
matlab复制function [tau_f, z_new] = LuGre(v, z, F_c, F_s, sigma, v_s, Ts)
g = F_c + (F_s - F_c)*exp(-(v/v_s)^2);
dz = v - sigma(1)*abs(v)/g*z;
z_new = z + dz*Ts; % 欧拉离散化
tau_f = sigma