相平面分析法在非线性控制系统中的应用与实践

1. 相平面分析的核心价值与工程意义

相平面分析法作为非线性系统研究的经典工具，在电机控制、机器人轨迹规划等工程领域具有不可替代的作用。我最早接触这个方法是在研究生阶段的伺服系统调试中，当时面对一个存在死区特性的直流电机，传统频域分析完全失效，正是相平面图揭示了系统在低速区的极限环振荡现象。

相平面本质上是将二阶系统的时域行为映射到状态变量构成的二维平面上，横轴通常表示位移或位置误差，纵轴对应速度或误差变化率。这种可视化表达比时域波形更直观展现系统的动态特性，特别适合分析存在非线性环节（如饱和、间隙、滞环）的控制系统。

2. 相轨迹绘制方法与关键特征识别

2.1 等倾线法的工程实现步骤

以典型的非线性系统为例：

code复制ẍ + μ(1-x²)ẋ + x = 0 （范德波尔振荡器）

1. 状态空间转换：

   • 令x₁ = x（位移）
   • x₂ = ẋ（速度）
   • 得到状态方程：
     ẋ₁ = x₂
     ẋ₂ = -x₁ + μ(1-x₁²)x₂

2. 等倾线绘制：

   • 令dx₂/dx₁ = k（斜率）
   • 推导等倾线方程：
     k = [-x₁ + μ(1-x₁²)x₂]/x₂
   • 取k=0,±1,±∞等典型值，在x₁-x₂平面绘制曲线簇

实操提示：使用MATLAB的quiver函数可快速生成相轨迹矢量场，配合ode45求解器能验证手工绘制的等倾线。

2.2 奇点类型判据速查表

3. 典型非线性环节的相平面分析

3.1 带死区的继电器特性

考虑位置控制系统中的电机死区模型：

code复制u = { 
  +Umax, e > Δ
  0, |e| ≤ Δ
  -Umax, e < -Δ
}

相轨迹特征：

1. 在|e|≤Δ区域表现为无阻尼自由运动
2. 穿越死区边界时出现轨迹斜率突变
3. 最终形成稳定的极限环振荡

调试经验：死区宽度Δ超过系统自然频率的1/5时，必然出现肉眼可见的稳态抖动。

3.2 饱和非线性系统的相轨迹

以伺服驱动器输出饱和为例：

code复制dx/dt = v
dv/dt = { 
  -k*e, |e|≤emax
  -k*emax, e>emax
  k*emax, e<-emax
}

相平面分区策略：

1. 线性区（|e|≤emax）：椭圆型轨迹
2. 饱和区（|e|>emax）：抛物线轨迹
3. 切换边界处需满足速度连续条件

4. 工程应用中的相平面诊断技巧

4.1 机械臂关节振动分析案例

某SCARA机械臂在拾放作业时出现末端抖动：

1. 采集关节位置误差e与误差导数ė
2. 相图显示存在闭合极限环
3. 结合模型识别出谐波减速器的反向间隙达0.08rad
4. 补偿措施：
   • 在|e|<0.05rad区间增加前馈力矩
   • 调整PID的微分增益抑制穿越速度

4.2 电动汽车再生制动振荡抑制

当电机转矩指令与机械制动存在耦合时：

1. 相轨迹在低速区呈现发散螺旋
2. 根本原因是制动摩擦系数的速度负特性
3. 解决方案：
   • 在相平面划分混合制动模式切换边界
   • 引入ė依赖的转矩补偿项：
     ΔT = K·tanh(ė/ε)

5. 现代控制理论中的相平面延伸

虽然相平面法局限于二阶系统，但其思想可推广到：

1. 高维系统的投影相空间分析
2. 李雅普诺夫函数构造的几何指导
3. 滑模控制切换面的设计依据

在机器人路径规划中，我常用相平面思想处理以下问题：

• 关节空间到操作空间的加速度约束转换
• 基于速度-位置相图的防碰撞策略
• 考虑执行器饱和时的轨迹优化

进阶工具推荐：

• MATLAB的pplane工具箱（教学用）
• Python的PyPhase库（工业级分析）
• ROS中的rqt_phase_plot插件（实时诊断）