二维线性动力系统

来自集智百科 - 复杂系统|人工智能|复杂科学|复杂网络|自组织
跳到导航 跳到搜索

二维线性动力系统是线性动力系统在二维状态空间中的基本形式，其系统状态随时间的演化规律由线性动力学方程决定。状态可表示为[math]\displaystyle{ \mathbf{x}(t)=[x_1(t),x_2(t)]^T }[/math]，其中[math]\displaystyle{ x_1(t) }[/math]和[math]\displaystyle{ x_2(t) }[/math]为两个状态变量，变量会随时间演化。连续时间二维线性动力系统通常写作[math]\displaystyle{ \dot{\mathbf{x}}=A\mathbf{x} }[/math]，离散时间形式则写作[math]\displaystyle{ \mathbf{x}_{k+1}=A\mathbf{x}_k }[/math]，其中 [math]\displaystyle{ A\in\mathbb{R}^{2\times 2} }[/math]为系统矩阵，用于描述两个状态变量之间的线性耦合关系。该类系统常用于分析平面相图、振荡行为和稳定性问题[1][2]。

性质[2]

二维向量场与相平面

对于二维线性动力系统[math]\displaystyle{ \dot{\mathbf{x}}=A\mathbf{x} }[/math]，状态向量可写为[math]\displaystyle{ \mathbf{x}=(x_1,x_2)^T }[/math]。在相平面中，每一个状态点[math]\displaystyle{ ((x_1,x_2) }[/math]都通过矩阵[math]\displaystyle{ A }[/math]对应一个状态变化率向量[math]\displaystyle{ (\dot{x}_1,\dot{x}_2)^T }[/math]，表示系统在该点处的瞬时变化方向和变化速率。所有状态点上的这些向量共同构成系统的向量场，直观展示了系统在二维状态空间中的流动趋势，是分析相轨迹、平衡点和稳定性的基础。

相轨迹的几何形态

在二维线性动力系统中，系统状态可写为 [math]\displaystyle{ \mathbf{x}(t)=(x_1(t),x_2(t))^T }[/math]，对应的二维状态空间称为相平面。给定初始状态[math]\displaystyle{ \mathbf{x}(0)=(x_{1}(0),x_{2}(0))^T }[/math]后，状态点会随时间 t 在相平面中运动，其描绘出的曲线称为相轨迹。相轨迹用于直观表示系统从初始状态出发后的演化路径，并可用于判断系统是否趋于平衡、发散或呈现振荡行为。

相图的定性分析

在二维线性动力系统中，研究重点不仅在于求出解析解[math]\displaystyle{ \mathbf{x}(t) }[/math]，也在于通过相平面中的整体轨迹分布判断系统的定性行为。通过观察相轨迹的几何形状、流动方向以及其与平衡点的关系，可以在不显式求解微分方程的情况下，分析系统的稳定性、振荡性和长期演化趋势。例如，无阻尼弹簧振子在相平面中的轨迹通常表现为围绕平衡点的闭合曲线，反映了系统能量守恒和周期振荡的特征。

平衡点附近的相平面结构

在二维线性动力系统中，满足[math]\displaystyle{ A\mathbf{x}=\mathbf{0} }[/math]的状态称为系统的平衡点或不动点。对于连续时间系统[math]\displaystyle{ \dot{\mathbf{x}}=A\mathbf{x} }[/math]，若系统矩阵 A 非奇异，则原点 (0,0) 是唯一平衡点；若 A 奇异，则可能存在一条由平衡点组成的直线，而当[math]\displaystyle{ A=0 }[/math] 时，整个相平面上的点均为平衡点。平衡点附近相轨迹的走向反映了系统的局部稳定性：若附近轨迹随时间趋向平衡点，则该平衡点是渐近稳定的；若轨迹沿某些方向远离平衡点，则系统是不稳定的；若轨迹围绕平衡点运动，则系统可能表现为振荡行为，例如中心或焦点型相图。

基于迹—行列式平面的分类

对于二维线性系统，平衡点附近的局部几何结构主要由矩阵[math]\displaystyle{ \displaystyle{ A } }[/math]的特征值完全刻画。设矩阵[math]\displaystyle{ \displaystyle{ A } }[/math]的两个特征值为[math]\displaystyle{ \displaystyle{ \lambda_1,\lambda_2 } }[/math]，其迹（Trace）与行列式（Determinant）分别为：

[math]\displaystyle{ \displaystyle{ T=\lambda_1+\lambda_2,\quad \Delta=\lambda_1\lambda_2 } }[/math]

根据特征方程[math]\displaystyle{ \displaystyle{ \lambda^2-T\lambda+\Delta=0 } }[/math]，由求根公式可得：

[math]\displaystyle{ \displaystyle{ \lambda_{1,2}=\frac{T\pm\sqrt{T^2-4\Delta}}{2} } }[/math]

因此，可以利用[math]\displaystyle{ \displaystyle{ T } }[/math]、[math]\displaystyle{ \displaystyle{ \Delta } }[/math]以及判别式[math]\displaystyle{ \displaystyle{ T^2-4\Delta } }[/math]对二维线性动力系统的平衡点类型进行分类。对于重复特征值或零特征值等边界情形，还需要结合特征向量及矩阵结构进一步判断。各类典型情况可概括如下：

条件 特征值类型 相图类型 稳定性
[math]\displaystyle{ \Delta\lt 0 }[/math] 一正一负实特征值 鞍点 不稳定
[math]\displaystyle{ \Delta\gt 0,\ T^2-4\Delta\gt 0,\ T\lt 0 }[/math] 两个负实特征值 稳定结点 渐近稳定
[math]\displaystyle{ \Delta\gt 0,\ T^2-4\Delta\gt 0,\ T\gt 0 }[/math] 两个正实特征值 不稳定结点 不稳定
[math]\displaystyle{ \Delta\gt 0,\ T^2-4\Delta\lt 0,\ T\lt 0 }[/math] 共轭复特征值，实部为负 稳定焦点 渐近稳定
[math]\displaystyle{ \Delta\gt 0,\ T^2-4\Delta\lt 0,\ T\gt 0 }[/math] 共轭复特征值，实部为正 不稳定焦点 不稳定
[math]\displaystyle{ \Delta\gt 0,\ T=0 }[/math] 一对纯虚特征值 中心 稳定但非渐近稳定
[math]\displaystyle{ T^2-4\Delta=0,\ \Delta\gt 0 }[/math] 重复实特征值 星形结点或退化结点 由 [math]\displaystyle{ T }[/math]的符号决定
[math]\displaystyle{ \Delta=0 }[/math] 至少一个零特征值 非孤立平衡点 需进一步判断

如下图 Fig.1 所示，以行列式[math]\displaystyle{ \Delta }[/math]为横轴、迹[math]\displaystyle{ T }[/math]为纵轴，可以在迹—行列式平面上对上述平衡点类型及其稳定性进行直观的几何划分。

Fig.1 二维线性系统不动点的迹-行列式（Trace-Determinant）分类平面图
  • 鞍点（Saddle Point）： 当 [math]\Delta < 0[/math]（即左半平面）时，特征值为一正一负的实数。 此时鞍点在两个特征方向上表现出不同的动力学行为。负特征值对应'''稳定方向'''，该方向上的轨迹随时间逐渐趋近平衡点；正特征值对应'''不稳定方向'''，轨迹则随时间远离平衡 点。由于只有严格位于稳定方向上的轨迹才能收敛到平衡点，而附近的微小扰动产生的不稳定方向分量会随时间增长，因此鞍点整体上是不稳定的。例如，系统 [math]\dot{x}=x,\ \dot{y}=-y[/math] 的两个特征值分别为 [math]1[/math] 和 [math]-1[/math]，其中 [math]y[/math] 轴为稳定方向，[math]x[/math] 轴为不稳定方向，其余轨迹呈双曲线形状，由此形成如图Fig.2所示的典型鞍点相图
    Fig.2 鞍点相图
  • 结点（Node）： 当 [math]\Delta>0[/math] 且 [math]T^2-4\Delta>0[/math] 时，系统具有两个同号且互异的实特征值，对应于结点（Node）。由于两个特征值均为实数，相轨迹通常沿特征方向趋近或远离平衡点，而不表现出围绕平衡点的旋转行为。 当 [math]T<0[/math] 时，两个特征值均为负，系统形成稳定结点（Stable Node）。如 Fig.3(a) 所示，相轨迹从不同方向逐渐趋向平衡点，表明该平衡点为渐近稳定。 当 [math]T>0[/math] 时，两个特征值均为正，系统形成不稳定结点（Unstable Node）。如 Fig.3(b) 所示，相轨迹从平衡点向外发散，附近任意微小扰动都会随时间远离平衡点，因此该平衡点是不稳定的。
  • 焦点（Focus）： 当 [math]\Delta>0[/math] 且 [math]T^2-4\Delta<0[/math] 时，系统具有一对共轭复特征值，对应于焦点（Focus）。此时特征值可写为 [math]\lambda_{1,2}=T/2\pm i\omega[/math]，其中虚部使相轨迹围绕平衡点旋转，而实部 [math]T/2[/math] 决定轨迹是逐渐收缩还是向外发散。 当 [math]T<0[/math] 时，共轭复特征值的实部为负，系统形成稳定焦点（Stable Focus）。如 Fig.3(c) 所示，相轨迹以螺旋形式逐渐趋近平衡点，因此该平衡点为渐近稳定。 当 [math]T>0[/math] 时，共轭复特征值的实部为正，系统形成不稳定焦点（Unstable Focus）。如 Fig.3(d) 所示，相轨迹以螺旋形式逐渐远离平衡点，因此该平衡点是不稳定的。
Fig.3 结点与焦点的典型相图
  • 中心（Center）： 当 [math]\Delta>0[/math] 且 [math]T=0[/math] 时，系统具有一对纯虚共轭特征值，可写为 [math]\lambda_{1,2}=\pm i\sqrt{\Delta}[/math]。由于特征值的实部为零，轨迹既不会随时间趋近平衡点，也不会向外发散，而是围绕平衡点形成闭合轨道。 如 Fig.4 所示，相轨迹在平衡点周围持续作周期运动，因此该平衡点称为中心（Center）。此时系统保持周期性振荡，但轨迹不会收敛到平衡点。
    • Fig.4 中心相图
  • 星形结点（Star Node）： 当 [math]T^2-4\Delta=0[/math] 且 [math]\Delta>0[/math] 时，系统具有重复实特征值 [math]\lambda_1=\lambda_2=T/2[/math]。若系统矩阵 [math]A[/math] 具有两个线性无关的特征向量，则形成星形结点（Star Node）。此时相轨迹沿不同方向直接趋近或远离平衡点，如 Fig.5(a) 所示。当 [math]T<0[/math] 时，重复特征值为负，轨迹随时间趋近平衡点，形成稳定星形结点；当 [math]T>0[/math] 时，重复特征值为正，轨迹从平衡点向外发散，形成不稳定星形结点。
  • 退化结点（Improper Node）： 当 [math]T^2-4\Delta=0[/math] 且 [math]\Delta>0[/math] 时，若系统矩阵 [math]A[/math] 只有一个线性无关的特征向量，则形成退化结点（Improper Node）。与星形结点不同，退化结点附近的相轨迹通常呈弯曲形态，并在接近平衡点时逐渐与唯一的特征方向相切，如 Fig.5(b) 所示。当 [math]T<0[/math] 时，轨迹趋近平衡点，形成稳定退化结点；当 [math]T>0[/math] 时，轨迹远离平衡点，形成不稳定退化结点。
  • 非孤立平衡点（Non-isolated Equilibrium Points）： 当 [math]\Delta=0[/math] 时，系统至少存在一个零特征值，因此系统矩阵 [math]A[/math] 为奇异矩阵，平衡点不一定是孤立的。若 [math]T\neq0[/math]，两个特征值分别为 [math]0[/math] 和 [math]T[/math]，矩阵的零空间通常形成一条由平衡点组成的直线。当 [math]T<0[/math] 时，附近轨迹趋向这条平衡点直线，如 Fig.5(c) 所示；当 [math]T>0[/math] 时，附近轨迹则远离该平衡点直线。 当 [math]T=0[/math] 且 [math]\Delta=0[/math] 时，两个特征值均为零。若 [math]A=0[/math]，则整个相平面上的点均为平衡点；若 [math]A\neq0[/math]，则系统存在一条平衡点直线，其余轨迹沿与该直线平行的方向发生线性漂移。
Fig.5 退化情况: 星形结点、退化结点和非孤立平衡点

经典实例：二维谐振子系统[2]

为了更直观地理解线性系统的降维与重构，可以以经典的无阻尼弹簧振子为例[3]。设一个质量为[math]\displaystyle{ \displaystyle{ m } }[/math]的物体连接在劲度系数为[math]\displaystyle{ \displaystyle{ k } }[/math]的弹簧上，其偏离平衡位置的位移记为[math]\displaystyle{ \displaystyle{ x } }[/math]。根据牛顿第二定律与胡克定律，该系统满足以下二阶常微分方程：

[math]\displaystyle{ \displaystyle{ m\ddot{x}+kx=0 } }[/math]

为了将其转化为动力系统分析中标准的一阶方程组形式，引入速度变量[math]\displaystyle{ \displaystyle{ v=\dot{x} } }[/math]。此时，原方程可降阶并改写为：

[math]\displaystyle{ \displaystyle{ \dot{v}=-\frac{k}{m}x } }[/math]

联立上述两式，即可得到该二维线性动力系统的矩阵表达形式：

[math]\displaystyle{ \displaystyle{ \begin{pmatrix} \dot{x}\\ \dot{v} \end{pmatrix}=\begin{pmatrix} 0 & 1\\ -\frac{k}{m} & 0 \end{pmatrix}\begin{pmatrix} x\\ v \end{pmatrix} } }[/math]

从该演化方程可以看出，要完全确定系统未来的状态，仅仅知道位置[math]\displaystyle{ \displaystyle{ x } }[/math]这一物理量是不够的，还必须掌握系统当下的速度[math]\displaystyle{ \displaystyle{ v } }[/math]。因此，描述该系统的自然状态空间不再是一维的直线，而是由位置与速度构成的二维相平面（Phase plane），即[math]\displaystyle{ \displaystyle{ (x,v) } }[/math]空间。如 Fig.6 所示，弹簧振子的实际运动与相平面中的状态演化具有一一对应关系。

Fig.6 无阻尼弹簧振子的运动与相轨迹对应关系。 左侧展示质量块在弹簧作用下的周期往复运动，右侧展示对应状态点 ((x(t),v(t))) 在相平面中的演化。随着质量块的位置和速度周期性变化，状态点沿闭合椭圆轨迹运动。由于系统矩阵的迹 (T=0)、行列式 (\Delta=k/m>0)，其原点对应中心（Center）型平衡点。

















跨学科案例：罗密欧与朱丽叶的恋爱动力学[2]

二维线性系统不仅可以描述机械振动、电路演化等物理过程，也可以作为抽象模型刻画个体之间的相互作用。为了说明线性系统分类在不同学科中的普适性，非线性动力学研究者 Steven Strogatz 引入了一个经典的“恋爱动力学”（Love Dynamics）模型，用以描述两个人情感状态之间的线性反馈关系。

设[math]\displaystyle{ \displaystyle{ R(t) } }[/math]表示时刻[math]\displaystyle{ \displaystyle{ t } }[/math]罗密欧对朱丽叶的情感强度，[math]\displaystyle{ \displaystyle{ J(t) } }[/math]表示朱丽叶对罗密欧的情感强度，其中正值表示爱慕，负值表示厌恶。则两人的情感互动可表示为如下二维线性系统：

[math]\displaystyle{ \displaystyle{ \dot{R}=aR+bJ } }[/math]

[math]\displaystyle{ \displaystyle{ \dot{J}=cR+dJ } }[/math]

写成矩阵形式即：

[math]\displaystyle{ \displaystyle{ \begin{pmatrix} \dot{R}\\ \dot{J} \end{pmatrix}=\begin{pmatrix} a & b\\ c & d \end{pmatrix}\begin{pmatrix} R\\ J \end{pmatrix} } }[/math]

其中，参数 [math]\displaystyle{ \displaystyle{ a } }[/math]和[math]\displaystyle{ \displaystyle{ d } }[/math]描述个体对自身当前情感状态的内在反应；参数[math]\displaystyle{ \displaystyle{ b } }[/math]和[math]\displaystyle{ \displaystyle{ c } }[/math]描述双方对彼此情感变化的外部响应方式。不同参数符号和大小的组合，对应于不同的互动模式和长期演化结果。这一模型说明，即使系统变量来自社会或心理意义上的抽象量，其数学分析仍然可以采用与物理系统相同的线性动力学框架。

特例分析：相爱相杀的循环

作为一个典型情形，考虑如下系统：

[math]\displaystyle{ \displaystyle{ \dot{R}=bJ,\quad \dot{J}=cR } }[/math]

其中，[math]\displaystyle{ \displaystyle{ b\gt 0 } }[/math]，[math]\displaystyle{ \displaystyle{ c\lt 0 } }[/math]。这表示罗密欧主要依据朱丽叶的态度来调整自己：朱丽叶越爱他，他越热情；而朱丽叶则对罗密欧的接近表现出回避倾向：罗密欧越热情，她越倾向于后退。此时系统矩阵为：

[math]\displaystyle{ \displaystyle{ A=\begin{pmatrix} 0 & b\\ c & 0 \end{pmatrix} } }[/math]

计算其迹与行列式可得：

[math]\displaystyle{ \displaystyle{ T=0,\quad \Delta=-bc\gt 0 } }[/math]

根据二维线性系统的不动点分类，该系统的特征值为一对纯虚数，原点对应于一个中心（Center）。在相平面上，系统轨迹表现为围绕原点的闭合曲线。这意味着两人的情感状态既不会收敛到静止平衡，也不会无限发散，而是呈现出周期性的循环变化。

Fig.7 罗密欧与朱丽叶恋爱动力学的情感状态循环。 左侧展示两人情感状态的周期性交替，右侧展示对应状态点沿闭合椭圆相轨迹的演化。四种情感组合依次出现，体现了中心型平衡点的周期性动力学特征。

















如 Fig.7 所示，随着状态点沿闭合轨迹运动，两人的情感关系依次经历以下四种典型状态：

  • 彼此爱慕： 当 (R>0,\ J>0) 时，两人都处于积极的情感状态。此时罗密欧的情感继续增强，而朱丽叶的情感开始减弱。
  • 罗密欧爱慕，朱丽叶回避： 当 (R>0,\ J<0) 时，罗密欧仍然保持爱慕，而朱丽叶表现出回避或消极态度。此时两人的情感强度均在下降。
  • 彼此冷淡： 当 (R<0,\ J<0) 时，双方均处于消极的情感状态。罗密欧的情感继续下降，而朱丽叶的情感开始回升。
  • 朱丽叶爱慕，罗密欧冷淡： 当 (R<0,\ J>0) 时，朱丽叶重新表现出积极情感，而罗密欧仍处于消极状态。此时两人的情感均开始上升，系统逐渐回到彼此爱慕的阶段。

这一循环不断重复，使两人的情感状态既不会收敛到静止平衡，也不会无限发散，而是在积极与消极状态之间周期性交替。对应的相轨迹为闭合椭圆，体现了中心型平衡点所具有的周期性动力学特征。


这一案例展示了动力系统分析中的结构相似性：二维线性动力系统的不动点分类不仅可以用于解释弹簧振子等经典物理系统中的周期振动，也可以用于刻画社会或心理模型中的相互作用关系。在数学上，这类问题都可以归结为对系统矩阵的特征值、迹与行列式的分析；在几何上，它们都可以通过相平面中的轨迹结构进行直观理解。

  1. ↑ Alan V. Oppenheim, Alan S. Willsky, and S. Hamid Nawab, Signals and Systems, 2nd ed., Prentice Hall, 1997.
  2. ↑ 2.0 2.1 2.2 2.3 S. H. Strogatz, 非线性动力学与混沌, 孙梅, 译. 北京, 中国: 机械工业出版社, 2016.
  3. ↑ J. R. Taylor, Classical Mechanics. Sausalito, CA, USA: University Science Books, 2005.
取自“https://wiki.swarma.org/index.php?title=二维线性动力系统&oldid=51548”