4. 线性系统的相平面分析:一阶线性系统的相轨迹、二阶线性系统的相轨迹、奇点类型与系统稳定性

各位好,我是老张。今天咱们聊聊线性系统的相平面分析。说实话,这块内容我当年学的时候也觉得有点抽象,直到后来做项目时被一个振荡问题折磨了三天,才真正体会到相平面的威力。你想想看,一个系统到底稳不稳定?响应是快是慢?有没有振荡?这些问题,相平面都能给你一个直观的答案。

4.1 一阶线性系统的相轨迹

先看最简单的。一阶线性系统,微分方程长这样:

dx/dt = ax

其中 a 是常数。这个方程的解,大家都会求:x(t) = x₀·e^(at)。

那相轨迹呢?说白了,就是状态变量 x 和它的导数 dx/dt 之间的关系。对于一阶系统,相平面就是一条直线——横轴是 x,纵轴是 dx/dt。

从方程直接可得:

dx/dt = a·x

这是一条过原点的直线。斜率就是 a。

我给大家总结一下规律:

  • a < 0:斜率为负。相轨迹从远处指向原点。系统稳定。
  • a > 0:斜率为正。相轨迹从原点向外发散。系统不稳定。
  • a = 0:斜率为零。相轨迹是水平线。系统临界稳定。

核心要点:一阶系统的相轨迹只有一种形状——直线。奇点就是原点,它决定了系统的全部命运。

我在项目中遇到过一个小问题:有个温度控制回路,一阶模型,按理说 a 应该是负的。但实际调试时发现温度一直往上飘。画了相轨迹一看,斜率是正的!后来查出来是传感器接反了,反馈变成了正反馈。嗯,相平面一眼就揪出了这个低级错误。

4.2 二阶线性系统的相轨迹

二阶系统就丰富多了。标准形式:

d²x/dt² + 2ζωₙ·dx/dt + ωₙ²·x = 0

写成状态方程:

dx₁/dt = x₂
dx₂/dt = -ωₙ²·x₁ - 2ζωₙ·x₂

其中 x₁ = x,x₂ = dx/dt。

相轨迹的形状取决于阻尼比 ζ。我习惯把 ζ 分成几个区间来看:

阻尼比 ζ 特征根 相轨迹形状 稳定性
ζ < 0 实部为正 螺旋向外发散 不稳定
ζ = 0 纯虚根 同心椭圆 临界稳定
0 < ζ < 1 共轭复根(负实部) 螺旋向内收敛 稳定
ζ = 1 相等负实根 趋向原点的抛物线 稳定
ζ > 1 两个负实根 趋向原点的曲线 稳定

为什么会这样?说白了,特征根决定了相轨迹的「走向」。负实部让轨迹往原点跑,虚部让轨迹转圈圈。

我的经验:实际项目中,0.4 < ζ < 0.8 是最常用的区间。太小的 ζ 振荡太多,太大的 ζ 响应太慢。我曾经调一个伺服电机,ζ 设到 0.3,结果位置超调了 30%,差点撞到限位开关。后来改成 0.6,稳得很。

4.3 奇点类型与系统稳定性

奇点,就是相平面上 dx₁/dt = 0 且 dx₂/dt = 0 的点。对于线性系统,原点 (0,0) 是唯一的奇点。

奇点的类型,完全由系统矩阵 A 的特征值决定。我给大家画个分类图:

奇点类型与稳定性分类 奇点 (0,0) 实特征值 复特征值 零特征值 同号 异号 均为负→稳定节点 均为正→不稳定节点 鞍点(不稳定) 实部≠0 实部=0 实部负→稳定焦点 实部正→不稳定焦点 中心点(临界稳定) 退化情况 非孤立奇点 核心规律:特征值实部决定稳定性,虚部决定振荡特性

这张图我花了不少心思。你顺着看就明白了:先看特征值是实是复,再往下细分。每个叶子节点对应一种奇点类型。

稳定性判据:所有特征值实部 < 0 → 系统稳定;只要有一个实部 > 0 → 系统不稳定;有实部 = 0 且其余实部 < 0 → 临界稳定。

我给大家几个具体的例子:

例1:dx/dt = [-2 0; 0 -3]·x
特征值:λ₁ = -2, λ₂ = -3
→ 稳定节点,相轨迹直接指向原点

例2:dx/dt = [0 1; -4 0]·x
特征值:λ₁ = 2j, λ₂ = -2j
→ 中心点,相轨迹是椭圆,系统临界稳定

例3:dx/dt = [1 1; -2 -1]·x
特征值:λ₁ = j, λ₂ = -j
→ 中心点,但注意!实际系统有扰动就可能发散

避坑指南:我曾经在设计一个飞行器姿态控制系统时,算出来的特征值实部都是负的,理论上稳定。但实际试飞时发现姿态一直在小幅振荡。后来用相平面一画,发现有个共轭复根的实部虽然为负,但非常接近零,阻尼比只有 0.05。这种「弱稳定」系统在实际中很容易被噪声激励出振荡。所以我的建议是:不要只看正负,还要看实部的绝对值大小。

4.4 奇点类型与系统稳定性的对应关系

咱们把前面的内容串起来。奇点类型和稳定性之间,其实是一一对应的:

  • 稳定节点:两个负实特征值。系统单调收敛,没有振荡。
  • 稳定焦点:共轭复根,实部为负。系统振荡收敛。
  • 中心点:纯虚根。系统等幅振荡,临界稳定。
  • 不稳定焦点:共轭复根,实部为正。系统振荡发散。
  • 不稳定节点:两个正实特征值。系统单调发散。
  • 鞍点:一正一负实特征值。系统从一个方向收敛,另一个方向发散。

你想想看,这其实很直观。特征值实部就是系统的「阻尼」,虚部就是「刚度」。阻尼够大,系统就稳;刚度够大,系统就爱振荡。

实用技巧:我在做控制器设计时,习惯先画出系统的相平面草图。不用精确计算,大致判断奇点类型和轨迹走向就行。这能帮我快速排除一些明显不合理的参数组合。比如看到相轨迹往外跑,那肯定是有个特征值实部大于零,得赶紧改参数。

好了,这一章的内容就到这儿。相平面分析说白了就是「看图说话」——看轨迹的形状,就知道系统稳不稳、怎么稳。下一章咱们会聊非线性系统的相平面分析,那才是真正考验功力的时候。


公众号:蓝海资料掘金营,微信deep3321