二维线性动力系统

2026-06-25T21:00:21Z

张新宇：

'''二维线性动力系统'''是线性动力系统在二维状态空间中的基本形式，且系统状态随时间的演化规律由线性方程决定的动力系统。状态可表示为<math>\mathbf{x}(t)=[x_1(t),x_2(t)]^T</math>，其中<math>x_1(t)</math>和<math>x_2(t)</math>为两个状态变量。连续时间二维线性动力系统通常写作<math>\dot{\mathbf{x}}=A\mathbf{x}</math>，离散时间形式则写作<math>\mathbf{x}_{k+1}=A\mathbf{x}_k</math>，其中 <math>A\in\mathbb{R}^{2\times 2}</math>为系统矩阵，用于描述两个状态变量之间的线性耦合关系。该类系统常用于分析平面相图、振荡行为和稳定性问题<ref>Alan V. Oppenheim, Alan S. Willsky, and S. Hamid Nawab, ''Signals and Systems'', 2nd ed., Prentice Hall, 1997.</ref><ref name=":4">S. H. Strogatz, 非线性动力学与混沌, 孙梅, 译. 北京, 中国: 机械工业出版社, 2016.</ref>。
=='''性质<ref name=":4" />'''==
===二维向量场与相平面===
对于二维线性动力系统<math>\dot{\mathbf{x}}=A\mathbf{x}</math>，状态向量可写为<math>\mathbf{x}=(x,y)^T</math>。在相平面中，每一个状态点<math>(x,y)</math>都通过矩阵 A 对应一个速度向量<math>(\dot{x},\dot{y})^T</math>，表示系统在该点处的瞬时变化方向和变化速率。所有状态点上的速度向量共同构成系统的'''向量场'''。向量场直观展示了系统在二维状态空间中的流动趋势，是分析相轨迹、平衡点和稳定性的基础。
===相轨迹的几何形态===
在二维线性动力系统中，若两个状态变量分别表示位置和速度，则系统状态可写为<math>\mathbf{x}(t)=(x(t),v(t))^T</math>，对应的二维状态空间称为'''相平面'''。给定初始状态<math>\mathbf{x}_0=(x_0,v_0)^T</math>后，状态点会随时间 t 在相平面中运动，其描绘出的曲线称为'''相轨迹'''。相轨迹用于直观表示系统从初始状态出发后的运动路径，并可用于判断系统是否趋于平衡、发散或呈现振荡行为。
===相图的定性分析===
在二维线性动力系统中，研究重点不仅在于求出解析解<math>\mathbf{x}(t)</math>，也在于通过相平面中的整体轨迹分布判断系统的定性行为。通过观察相轨迹的几何形状、流动方向以及其与平衡点的关系，可以在不显式求解微分方程的情况下，分析系统的稳定性、振荡性和长期演化趋势。例如，无阻尼弹簧振子在相平面中的轨迹通常表现为围绕平衡点的闭合曲线，反映了系统能量守恒和周期振荡的特征。
===平衡点附近的相平面结构===
在二维线性动力系统中，满足<math>A\mathbf{x}=\mathbf{0}</math>的状态称为系统的'''平衡点'''或'''不动点'''。对于连续时间系统<math>\dot{\mathbf{x}}=A\mathbf{x}</math>，若系统矩阵 A 非奇异，则原点 (0,0) 是唯一平衡点；若 A 奇异，则可能存在一条直线或更高维子空间上的平衡点。平衡点附近相轨迹的走向决定系统的局部稳定性，例如轨迹趋向平衡点表示稳定，远离平衡点表示不稳定，而围绕平衡点旋转则可能对应振荡行为。

===基于迹—行列式平面的分类===
对于二维线性系统，局部几何结构可以由矩阵<math>\displaystyle{ A }</math>的特征值完全刻画。设矩阵<math>\displaystyle{ A }</math>的两个特征值为<math>\displaystyle{ \lambda_1,\lambda_2 }</math>，其迹（Trace）与行列式（Determinant）分别为：

<math>\displaystyle{ T=\lambda_1+\lambda_2,\quad \Delta=\lambda_1\lambda_2 }</math>

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

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

因此，二维线性系统不动点的类型可统一由<math>\displaystyle{ T }</math>、<math>\displaystyle{ \Delta }</math>以及判别式<math>\displaystyle{ T^2-4\Delta }</math>来决定。
[[文件:二维线性系统不动点分类.png|缩略图|二维线性系统不动点的迹-行列式（Trace-Determinant）分类平面图]]

如右图所示，以行列式 [math]\Delta[/math] 为横轴，迹 [math]T[/math] 为纵轴，可以将系统所有的定性行为在二维平面上进行直观的几何划分：

* '''鞍点（Saddle Point）：''' 当 [math]\Delta < 0[/math]（即左半平面）时，特征值为一正一负的实数。此时系统在某个方向上吸引，在另一个方向上排斥，整体表现为不稳定。
* '''抛物线分界线：''' 曲线 [math]T^2 - 4\Delta = 0[/math] 是实数根与复数根的临界线。
** 在抛物线'''外部'''（[math]\Delta > 0[/math] 且 [math]T^2 - 4\Delta > 0[/math]），特征值为同号实数，对应于'''结点（Node）'''。
** 在抛物线'''内部'''（[math]\Delta > 0[/math] 且 [math]T^2 - 4\Delta < 0[/math]），特征值为共轭复数，对应于伴随旋转行为的'''焦点（Focus）'''。
* '''稳定与不稳定的分界：''' 对于右半平面（[math]\Delta > 0[/math]），横轴 [math]T = 0[/math] 是判断稳定性的关键。
** 当 [math]T < 0[/math]（下半平面），系统是'''渐近稳定'''的（稳定结点或稳定焦点），能量逐渐耗散，轨迹最终汇聚于原点。
** 当 [math]T > 0[/math]（上半平面），系统是'''不稳定'''的（不稳定结点或不稳定焦点），轨迹向外发散。
** 当 [math]T = 0[/math]（正半横轴），特征值为纯虚数，系统能量守恒，形成围绕原点的闭合周期轨道，即'''中心（Center）'''。下文提到的“相爱相杀”模型正属于此类。
* '''退化情况：''' 抛物线本身对应着特征值相等的退化结点或星点；而纵轴（[math]\Delta = 0[/math]）对应着存在零特征值的非孤立不动点（即存在一条由平衡点组成的直线）。

==经典实例：二维谐振子系统<ref name=":4" />==
为了更直观地理解线性系统的降维与重构，可以以经典的无阻尼弹簧振子为例<ref>J. R. Taylor, Classical Mechanics. Sausalito, CA, USA: University Science Books, 2005.</ref>。设一个质量为<math>\displaystyle{ m }</math>的物体连接在劲度系数为<math>\displaystyle{ k }</math>的弹簧上，其偏离平衡位置的位移记为<math>\displaystyle{ x }</math>。根据牛顿第二定律与胡克定律，该系统满足以下二阶常微分方程：

<math>\displaystyle{ m\ddot{x}+kx=0 }</math>

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

<math>\displaystyle{ \dot{v}=-\frac{k}{m}x }</math>

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

<math>\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{ x }</math>这一物理量是不够的，还必须掌握系统当下的速度<math>\displaystyle{ v }</math>。因此，描述该系统的自然状态空间不再是一维的直线，而是由位置与速度构成的二维'''相平面（Phase plane），即<math>\displaystyle{ (x,v) }</math>空间。'''
==跨学科案例：罗密欧与朱丽叶的恋爱动力学<ref name=":4" />==

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

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

<math>\displaystyle{ \dot{R}=aR+bJ }</math>

<math>\displaystyle{ \dot{J}=cR+dJ }</math>

写成矩阵形式即：

<math>\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{ a }</math>和<math>\displaystyle{ d }</math>描述个体对自身当前情感状态的内在反应；参数<math>\displaystyle{ b }</math>和<math>\displaystyle{ c }</math>描述双方对彼此情感变化的外部响应方式。不同参数符号和大小的组合，对应于不同的互动模式和长期演化结果。这一模型说明，即使系统变量来自社会或心理意义上的抽象量，其数学分析仍然可以采用与物理系统相同的线性动力学框架。
===特例分析：相爱相杀的循环===
作为一个典型情形，考虑如下系统：

<math>\displaystyle{ \dot{R}=bJ,\quad \dot{J}=cR }</math>

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

<math>\displaystyle{ A=\begin{pmatrix} 0 & b\\ c & 0 \end{pmatrix} }</math>

计算其迹与行列式可得：

<math>\displaystyle{ T=0,\quad \Delta=-bc>0 }</math>

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

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

线性系统

2026-06-25T20:44:05Z

张新宇：

线性系统

2026-05-15T10:18:08Z

张新宇：添加了引用的文献

在数学与科学中，'''线性系统（Linear System）'''是指系统的输出变化与输入变化成比例，并且系统内部各要素之间满足'''叠加原理（Superposition Principle）'''<ref name=":0">Alan V. Oppenheim, Alan S. Willsky, and S. Hamid Nawab, ''Signals and Systems'', 2nd ed., Prentice Hall, 1997.</ref><ref name=":1">S. Axler, Linear Algebra Done Right, 3rd ed. Cham, Switzerland: Springer, 2015.</ref>的系统。叠加原理包括两个基本性质：'''齐次性'''和'''可加性'''。也就是说，如果某一输入产生某一响应，那么该输入按比例放大或缩小后，系统响应也会按相同比例变化；如果多个输入同时作用于系统，则系统的总响应等于各个输入单独作用时响应的代数和。因此，线性系统的整体行为可以由局部行为组合得到，通常比非线性系统更加直观、规则且易于分析和控制。

从数学形式上看，线性系统通常可以表示为一组'''线性代数方程、线性微分方程或线性差分方程'''。在这些方程中，未知变量、未知函数及其导数只以一次形式出现，不包含未知量之间的乘积、高次幂、指数函数、三角函数等非线性关系。换句话说，线性系统的方程可以写成未知变量或未知函数的线性组合。例如，若两个输入分别产生两个输出，则这两个输入的线性组合也会产生对应输出的线性组合，这正是线性系统区别于非线性系统的核心特征<ref name=":0" />。

需要注意的是，对于线性微分方程而言，判断其是否线性的关键在于'''未知函数及其导数是否线性出现'''。即使方程中的自变量、已知系数或外部输入随时间、空间等因素发生非线性变化，只要未知函数本身及其导数之间不存在乘积、高次项或其他非线性形式，该方程仍然可以称为线性方程。例如，一个系统的参数可以随时间变化，但只要状态变量之间不发生非线性耦合，该系统仍属于线性系统。

由于线性系统具有结构清晰、数学性质良好、求解方法成熟等特点，它们被广泛应用于工程、物理、控制理论、通信系统和信号处理等领域。许多复杂的实际系统虽然本质上可能是非线性的，但在一定工作范围或小扰动条件下，可以通过线性化方法近似为线性系统，从而便于研究其稳定性、响应特性和控制规律。因此，线性系统不仅是一类重要的数学模型，也是理解、分析和设计复杂系统的基础工具<ref name=":0" /><ref name=":2">H. K. Khalil, Nonlinear Systems, 3rd ed. Upper Saddle River, NJ, USA: Prentice Hall, 2002.</ref>。
== 定义 ==
在数学中，若一个映射或算子 f(x) 满足以下两个性质，则称其为'''线性映射'''或'''线性算子'''<ref name=":1" />：

[math]\displaystyle{ f(x+y)=f(x)+f(y) }[/math]

[math]\displaystyle{ f(\alpha x)=\alpha f(x) }[/math]

其中，第一个性质称为'''可加性（Additivity）'''，表示两个输入共同作用时的输出等于它们分别作用时输出的和；第二个性质称为'''齐次性（Homogeneity）'''，表示输入按比例放大或缩小时，输出也按相同比例变化。二者合在一起通常称为'''叠加原理（Superposition Principle）<ref name=":0" /><ref name=":1" />'''，即

[math]\displaystyle{ f(\alpha x+\beta y)=\alpha f(x)+\beta f(y) }[/math]

因此，线性系统的本质在于：系统的整体响应可以由各个输入或状态变量的局部响应通过线性组合得到。

一般地，若一个系统可以写成

[math]\displaystyle{ f(x)=C }[/math]

其中 x 表示待求的未知对象，C 为已知项，并且 f(x) 是关于 x 的线性映射，则该系统称为'''线性系统（Linear System）'''；否则称为'''非线性系统（Nonlinear System）'''。当 C=0 时，该系统称为'''齐次线性系统'''；当 C≠0 时，则称为'''非齐次线性系统<ref name=":1" />'''。

这里的 x 可以是非常一般的数学对象，例如数、向量、函数，或满足一定边界条件的函数。因此，f(x) 也可以表示不同类型的线性运算，例如矩阵变换、积分算子、微分算子或差分算子。如果 f(x) 中包含对未知函数 x 的微分运算，则对应的系统就是'''线性微分方程系统'''<ref>W. E. Boyce, R. C. DiPrima, and D. B. Meade, Elementary Differential Equations and Boundary Value Problems, 11th ed. Hoboken, NJ, USA: Wiley, 2017.</ref>。
== 动力系统中的线性系统 ==
在动力学语境下，线性系统通常是研究多维连续动力系统的基本起点。一般地，n 维线性常系数微分方程组可以写成如下矩阵形式<ref name=":3">W. E. Boyce, R. C. DiPrima, and D. B. Meade, Elementary Differential Equations and Boundary Value Problems, 11th ed. Hoboken, NJ, USA: Wiley, 2017.</ref>：

<nowiki>[math]\displaystyle{ \frac{d\mathbf{x}}{dt}=A\mathbf{x} }[/math]</nowiki>

其中，[math]\displaystyle{ \mathbf{x}\in\mathbb{R}^n }[/math]表示系统的状态向量，A 为常数矩阵，用于描述系统内部各状态变量之间的线性耦合关系。该方程表明，系统状态的变化率由当前状态经过矩阵 A 的线性变换决定。

线性系统的一个核心优势在于，它们往往存在解析解。例如，对于上述常系数线性系统，其解可以形式化地写为

[math]\displaystyle{ \mathbf{x}(t)=e^{At}\mathbf{x}(0) }[/math]

其中，x(0) 为初始状态，[math]\displaystyle{ e^{At} }[/math] 为矩阵指数。通过分析矩阵 A 的特征值和特征向量，可以进一步研究系统的稳定性、振荡行为、相轨迹结构以及长期演化趋势<ref name=":3" />。
=== 向量场与相平面 ===
==== 向量场：方程的几何映射 ====
<nowiki>从前述谐振子的例子可以看出，当线性系统被抽象为 [math]\displaystyle{ \dot{\mathbf{x}}=A\mathbf{x} }[/math] 的形式后，它便超越了纯代数意义下的矩阵微分方程，映射为几何空间中的一个动态结构。对于二维空间中的任意状态点 [math]\displaystyle{ (x,v) }[/math]，矩阵乘法 [math]\displaystyle{ A\mathbf{x} }[/math] 都会赋予该点一个确定的变化速度与方向向量 [math]\displaystyle{ (\dot{x},\dot{v}) }[/math]。这些密布于空间各点、指示着系统瞬时演化趋势的箭头，共同构成了该系统的</nowiki><nowiki>'''</nowiki>向量场<nowiki>'''</nowiki>（Vector Field）。向量场直观地刻画了系统在全局各个状态下的“流速”与“流向”，是理解二维线性系统整体行为的基础<ref name=":4">S. H. Strogatz, 非线性动力学与混沌, 孙梅, 译. 北京, 中国: 机械工业出版社, 2016.</ref>。
==== 相平面与相轨迹 ====
由系统的全部状态变量所张成的几何空间，被称为系统的<nowiki>'''</nowiki>相空间<nowiki>'''</nowiki>（Phase Space）；在二维情况下，特称为<nowiki>'''</nowiki>相平面<nowiki>'''</nowiki>（Phase Plane）。当给定系统一个初始状态 [math]\displaystyle{ (x_0,v_0) }[/math] 时，随着时间 [math]\displaystyle{ t }[/math] 的推移，状态点会顺着向量场的指引在相平面上运动，其描绘出的连续曲线即为系统的<nowiki>'''</nowiki>相轨迹<nowiki>'''</nowiki>（Trajectory）<ref name=":4" />。
==== 定性分析 ====
在复杂系统的动力学研究中，对二维线性系统的研究重点往往不再局限于求解微分方程的显式时间解析解[math]\displaystyle{ \mathbf{x}(t) }[/math]，更重要的是俯瞰这些轨迹在相平面中的整体拓扑分布。通过考察轨迹的几何形状、流动趋势以及它们与特殊点（不动点）的相对关系，可以在不强求解析解的情况下，直接把握系统的宏观定性行为。

对于无阻尼弹簧振子而言，对应到相平面上，轨迹不会相交，而是表现为一系列围绕原点（平衡点）的闭合椭圆。这种首尾相连的封闭轨迹，在几何上完美诠释了该保守系统的能量守恒特性与周期性机械振荡<ref name=":4" />。
=== 不动点与稳定性 ===
<nowiki>在相平面的研究中，除了把握整体轨迹的几何拓扑，另一个核心问题是系统在平衡状态附近的局部演化行为。若系统在某一点处满足 [math]\displaystyle{ \dot{\mathbf{x}}=A\mathbf{x}=\mathbf{0} }[/math]，则称该点为系统的</nowiki><nowiki>'''</nowiki>不动点<nowiki>'''</nowiki>（Fixed point）或<nowiki>'''</nowiki>平衡点<nowiki>'''</nowiki>（Equilibrium point）<ref name=":4" />。对于非奇异的二维线性系统而言，原点 [math]\displaystyle{ (0,0) }[/math] 通常是最基本的不动点。
==== 稳定性的层级分类 ====
在线性系统中，常用的稳定性概念包括：
* 吸引的 (Attracting)：若从不动点附近出发的轨迹在 [math]\displaystyle{ t\to\infty }[/math] 时都收敛到该点。

* 李雅普诺夫稳定 (Lyapunov stable)：若任意足够靠近不动点的初始条件产生的轨迹，在未来始终保持在给定的较小邻域内。

* 渐近稳定 (Asymptotically stable)：若一个不动点同时满足上述两个条件。

* 中性稳定 (Neutrally stable)：轨迹始终保持在平衡点附近但不向其收敛（如无阻尼谐振子的闭合椭圆）。

* 不稳定 (Unstable)：无论扰动多么微小，轨迹都可能逐渐远离平衡点。
==== 基于特征值、迹与行列式的分类 ====
对于二维线性系统，局部几何结构可以由矩阵 [math]\displaystyle{ A }[/math] 的特征值完全刻画。设矩阵 [math]\displaystyle{ A }[/math] 的两个特征值为 [math]\displaystyle{ \lambda_1,\lambda_2 }[/math]，其迹（Trace）与行列式（Determinant）分别为：

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

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

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

因此，二维线性系统不动点的类型可统一由 [math]\displaystyle{ T }[/math]、[math]\displaystyle{ \Delta }[/math] 以及判别式 [math]\displaystyle{ T^2-4\Delta }[/math] 来决定<ref name=":4" />。[[文件:二维线性系统不动点分类.png|缩略图|链接=https://wiki.swarma.org/index.php/%E6%96%87%E4%BB%B6:%E4%BA%8C%E7%BB%B4%E7%BA%BF%E6%80%A7%E7%B3%BB%E7%BB%9F%E4%B8%8D%E5%8A%A8%E7%82%B9%E5%88%86%E7%B1%BB.png]]
== 在线性化和非线性系统分析中的作用 ==
尽管现实世界中的许多系统具有非线性特征，但线性系统在非线性系统分析中仍然具有基础性作用。一个重要原因是，许多非线性动力系统在平衡点或不动点附近的局部行为，可以通过相应的线性系统近似描述。

设一个连续时间非线性动力系统<ref name=":4" />为

<nowiki>[math]\displaystyle{ \frac{d\mathbf{x}}{dt}=\mathbf{F}(\mathbf{x}) }[/math]</nowiki>

其中，[math]\displaystyle{ \mathbf{x}\in\mathbb{R}^n }[/math] 为状态向量，[math]\displaystyle{ \mathbf{F}(\mathbf{x}) }[/math] 为非线性向量场。若 [math]\displaystyle{ \mathbf{x}^\ast }[/math] 是系统的不动点，即

[math]\displaystyle{ \mathbf{F}(\mathbf{x}^\ast)=\mathbf{0} }[/math]

则在该点附近可令

[math]\displaystyle{ \boldsymbol{\xi}=\mathbf{x}-\mathbf{x}^\ast }[/math]

通过泰勒展开，系统可近似为

<nowiki>[math]\displaystyle{ \frac{d\boldsymbol{\xi}}{dt}=J(\mathbf{x}^\ast)\boldsymbol{\xi} }[/math]</nowiki>

其中

<nowiki>[math]\displaystyle{ J(\mathbf{x}^\ast)=\left.\frac{\partial \mathbf{F}}{\partial \mathbf{x}}\right|_{\mathbf{x}=\mathbf{x}^\ast} }[/math]</nowiki>

称为该非线性系统在不动点处的<nowiki>'''雅可比矩阵'''</nowiki>（Jacobian matrix）。因此，非线性系统在不动点附近的局部行为，通常可以由一个线性常系数系统来近似刻画。

这种方法称为'''线性化'''。通过分析雅可比矩阵的特征值，可以判断不动点附近的局部稳定性。例如，若所有特征值的实部均小于零，则该不动点通常是局部稳定的；若存在特征值的实部大于零，则该不动点通常是不稳定的。若特征值为复数，则系统可能表现出振荡、螺旋收敛或螺旋发散等相空间行为。

因此，线性系统不仅本身具有独立的研究价值，也构成了研究非线性动力系统的重要基础。特别是在二维及更高维相空间中，线性系统可以用于刻画稳定节点、鞍点、焦点、中心等典型局部结构。这些结构为理解更复杂的非线性现象，如分岔、[https://wiki.swarma.org/index.php/%E6%9E%81%E9%99%90%E7%8E%AF 极限环]和混沌等，提供了基本分析框架<ref name=":2" /><ref name=":4" />。
== 经典实例：二维谐振子系统 ==
为了更直观地理解线性系统的降维与重构，可以以经典的无阻尼弹簧振子为例<ref>J. R. Taylor, Classical Mechanics. Sausalito, CA, USA: University Science Books, 2005.</ref>。设一个质量为 [math]\displaystyle{ m }[/math] 的物体连接在劲度系数为 [math]\displaystyle{ k }[/math] 的弹簧上，其偏离平衡位置的位移记为 [math]\displaystyle{ x }[/math]。根据牛顿第二定律与胡克定律，该系统满足以下二阶常微分方程：

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

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

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

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

[math]\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{ x }[/math] 这一物理量是不够的，还必须掌握系统当下的速度 [math]\displaystyle{ v }[/math]。因此，描述该系统的自然状态空间不再是一维的直线，而是由位置与速度构成的二维<nowiki>'''</nowiki>相平面<nowiki>'''</nowiki>（Phase plane），即 [math]\displaystyle{ (x,v) }[/math] 空间<ref name=":4" />。

== 跨学科案例：罗密欧与朱丽叶的恋爱动力学<ref name=":4" /> ==

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

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

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

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

写成矩阵形式即：

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

=== 特例分析：相爱相杀的循环 ===
作为一个典型情形，考虑如下系统：

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

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

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

计算其迹与行列式可得：

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

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

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

== 参考文献 ==
<references />

线性系统

2026-05-15T09:31:19Z

张新宇：加了一条参考文献，做测试

在数学与科学中，'''线性系统（Linear System）'''是指系统的输出变化与输入变化成比例，并且系统内部各要素之间满足'''叠加原理（Superposition Principle）'''<ref>Alan V. Oppenheim, Alan S. Willsky, and S. Hamid Nawab, ''Signals and Systems'', 2nd ed., Prentice Hall, 1997.</ref>的系统。叠加原理包括两个基本性质：'''齐次性'''和'''可加性'''。也就是说，如果某一输入产生某一响应，那么该输入按比例放大或缩小后，系统响应也会按相同比例变化；如果多个输入同时作用于系统，则系统的总响应等于各个输入单独作用时响应的代数和。因此，线性系统的整体行为可以由局部行为组合得到，通常比非线性系统更加直观、规则且易于分析和控制。

从数学形式上看，线性系统通常可以表示为一组'''线性代数方程、线性微分方程或线性差分方程'''。在这些方程中，未知变量、未知函数及其导数只以一次形式出现，不包含未知量之间的乘积、高次幂、指数函数、三角函数等非线性关系。换句话说，线性系统的方程可以写成未知变量或未知函数的线性组合。例如，若两个输入分别产生两个输出，则这两个输入的线性组合也会产生对应输出的线性组合，这正是线性系统区别于非线性系统的核心特征。

需要注意的是，对于线性微分方程而言，判断其是否线性的关键在于'''未知函数及其导数是否线性出现'''。即使方程中的自变量、已知系数或外部输入随时间、空间等因素发生非线性变化，只要未知函数本身及其导数之间不存在乘积、高次项或其他非线性形式，该方程仍然可以称为线性方程。例如，一个系统的参数可以随时间变化，但只要状态变量之间不发生非线性耦合，该系统仍属于线性系统。

由于线性系统具有结构清晰、数学性质良好、求解方法成熟等特点，它们被广泛应用于工程、物理、控制理论、通信系统和信号处理等领域。许多复杂的实际系统虽然本质上可能是非线性的，但在一定工作范围或小扰动条件下，可以通过线性化方法近似为线性系统，从而便于研究其稳定性、响应特性和控制规律。因此，线性系统不仅是一类重要的数学模型，也是理解、分析和设计复杂系统的基础工具。
== 定义 ==
在数学中，若一个映射或算子 f(x) 满足以下两个性质，则称其为'''线性映射'''或'''线性算子'''：

[math]\displaystyle{ f(x+y)=f(x)+f(y) }[/math]

[math]\displaystyle{ f(\alpha x)=\alpha f(x) }[/math]

其中，第一个性质称为'''可加性（Additivity）'''，表示两个输入共同作用时的输出等于它们分别作用时输出的和；第二个性质称为'''齐次性（Homogeneity）'''，表示输入按比例放大或缩小时，输出也按相同比例变化。二者合在一起通常称为'''叠加原理（Superposition Principle）'''，即

[math]\displaystyle{ f(\alpha x+\beta y)=\alpha f(x)+\beta f(y) }[/math]

因此，线性系统的本质在于：系统的整体响应可以由各个输入或状态变量的局部响应通过线性组合得到。

一般地，若一个系统可以写成

[math]\displaystyle{ f(x)=C }[/math]

其中 x 表示待求的未知对象，C 为已知项，并且 f(x) 是关于 x 的线性映射，则该系统称为'''线性系统（Linear System）'''；否则称为'''非线性系统（Nonlinear System）'''。当 C=0 时，该系统称为'''齐次线性系统'''；当 C≠0 时，则称为'''非齐次线性系统'''。

这里的 x 可以是非常一般的数学对象，例如数、向量、函数，或满足一定边界条件的函数。因此，f(x) 也可以表示不同类型的线性运算，例如矩阵变换、积分算子、微分算子或差分算子。如果 f(x) 中包含对未知函数 x 的微分运算，则对应的系统就是'''线性微分方程系统'''。
== 动力系统中的线性系统 ==
在动力学语境下，线性系统通常是研究多维连续动力系统的基本起点。一般地，n 维线性常系数微分方程组可以写成如下矩阵形式：

<nowiki>[math]\displaystyle{ \frac{d\mathbf{x}}{dt}=A\mathbf{x} }[/math]</nowiki>

其中，[math]\displaystyle{ \mathbf{x}\in\mathbb{R}^n }[/math]表示系统的状态向量，A 为常数矩阵，用于描述系统内部各状态变量之间的线性耦合关系。该方程表明，系统状态的变化率由当前状态经过矩阵 A 的线性变换决定。

线性系统的一个核心优势在于，它们往往存在解析解。例如，对于上述常系数线性系统，其解可以形式化地写为

[math]\displaystyle{ \mathbf{x}(t)=e^{At}\mathbf{x}(0) }[/math]

其中，x(0) 为初始状态，[math]\displaystyle{ e^{At} }[/math] 为矩阵指数。通过分析矩阵 A 的特征值和特征向量，可以进一步研究系统的稳定性、振荡行为、相轨迹结构以及长期演化趋势。
=== 向量场与相平面 ===
==== 向量场：方程的几何映射 ====
从前述谐振子的例子可以看出，当线性系统被抽象为 [math]\displaystyle{ \dot{\mathbf{x}}=A\mathbf{x} }[/math] 的形式后，它便超越了纯代数意义下的矩阵微分方程，映射为几何空间中的一个动态结构。对于二维空间中的任意状态点 [math]\displaystyle{ (x,v) }[/math]，矩阵乘法 [math]\displaystyle{ A\mathbf{x} }[/math] 都会赋予该点一个确定的变化速度与方向向量 [math]\displaystyle{ (\dot{x},\dot{v}) }[/math]。这些密布于空间各点、指示着系统瞬时演化趋势的箭头，共同构成了该系统的<nowiki>'''</nowiki>向量场<nowiki>'''</nowiki>（Vector Field）。向量场直观地刻画了系统在全局各个状态下的“流速”与“流向”，是理解二维线性系统整体行为的基础。
==== 相平面与相轨迹 ====
由系统的全部状态变量所张成的几何空间，被称为系统的<nowiki>'''</nowiki>相空间<nowiki>'''</nowiki>（Phase Space）；在二维情况下，特称为<nowiki>'''</nowiki>相平面<nowiki>'''</nowiki>（Phase Plane）。当给定系统一个初始状态 [math]\displaystyle{ (x_0,v_0) }[/math] 时，随着时间 [math]\displaystyle{ t }[/math] 的推移，状态点会顺着向量场的指引在相平面上运动，其描绘出的连续曲线即为系统的<nowiki>'''</nowiki>相轨迹<nowiki>'''</nowiki>（Trajectory）。
==== 定性分析 ====
在复杂系统的动力学研究中，对二维线性系统的研究重点往往不再局限于求解微分方程的显式时间解析解[math]\displaystyle{ \mathbf{x}(t) }[/math]，更重要的是俯瞰这些轨迹在相平面中的整体拓扑分布。通过考察轨迹的几何形状、流动趋势以及它们与特殊点（不动点）的相对关系，可以在不强求解析解的情况下，直接把握系统的宏观定性行为。

对于无阻尼弹簧振子而言，对应到相平面上，轨迹不会相交，而是表现为一系列围绕原点（平衡点）的闭合椭圆。这种首尾相连的封闭轨迹，在几何上完美诠释了该保守系统的能量守恒特性与周期性机械振荡。
=== 不动点与稳定性 ===
在相平面的研究中，除了把握整体轨迹的几何拓扑，另一个核心问题是系统在平衡状态附近的局部演化行为。若系统在某一点处满足 [math]\displaystyle{ \dot{\mathbf{x}}=A\mathbf{x}=\mathbf{0} }[/math]，则称该点为系统的<nowiki>'''</nowiki>不动点<nowiki>'''</nowiki>（Fixed point）或<nowiki>'''</nowiki>平衡点<nowiki>'''</nowiki>（Equilibrium point）。对于非奇异的二维线性系统而言，原点 [math]\displaystyle{ (0,0) }[/math] 通常是最基本的不动点。
==== 稳定性的层级分类 ====
在线性系统中，常用的稳定性概念包括：
* 吸引的 (Attracting)：若从不动点附近出发的轨迹在 [math]\displaystyle{ t\to\infty }[/math] 时都收敛到该点。

* 李雅普诺夫稳定 (Lyapunov stable)：若任意足够靠近不动点的初始条件产生的轨迹，在未来始终保持在给定的较小邻域内。

* 渐近稳定 (Asymptotically stable)：若一个不动点同时满足上述两个条件。

* 中性稳定 (Neutrally stable)：轨迹始终保持在平衡点附近但不向其收敛（如无阻尼谐振子的闭合椭圆）。

* 不稳定 (Unstable)：无论扰动多么微小，轨迹都可能逐渐远离平衡点。
==== 基于特征值、迹与行列式的分类 ====
对于二维线性系统，局部几何结构可以由矩阵 [math]\displaystyle{ A }[/math] 的特征值完全刻画。设矩阵 [math]\displaystyle{ A }[/math] 的两个特征值为 [math]\displaystyle{ \lambda_1,\lambda_2 }[/math]，其迹（Trace）与行列式（Determinant）分别为：

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

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

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

因此，二维线性系统不动点的类型可统一由 [math]\displaystyle{ T }[/math]、[math]\displaystyle{ \Delta }[/math] 以及判别式 [math]\displaystyle{ T^2-4\Delta }[/math] 来决定。[[文件:二维线性系统不动点分类.png|缩略图|链接=https://wiki.swarma.org/index.php/%E6%96%87%E4%BB%B6:%E4%BA%8C%E7%BB%B4%E7%BA%BF%E6%80%A7%E7%B3%BB%E7%BB%9F%E4%B8%8D%E5%8A%A8%E7%82%B9%E5%88%86%E7%B1%BB.png]]
== 在线性化和非线性系统分析中的作用 ==
尽管现实世界中的许多系统具有非线性特征，但线性系统在非线性系统分析中仍然具有基础性作用。一个重要原因是，许多非线性动力系统在平衡点或不动点附近的局部行为，可以通过相应的线性系统近似描述。

设一个连续时间非线性动力系统为

<nowiki>[math]\displaystyle{ \frac{d\mathbf{x}}{dt}=\mathbf{F}(\mathbf{x}) }[/math]</nowiki>

其中，[math]\displaystyle{ \mathbf{x}\in\mathbb{R}^n }[/math] 为状态向量，[math]\displaystyle{ \mathbf{F}(\mathbf{x}) }[/math] 为非线性向量场。若 [math]\displaystyle{ \mathbf{x}^\ast }[/math] 是系统的不动点，即

[math]\displaystyle{ \mathbf{F}(\mathbf{x}^\ast)=\mathbf{0} }[/math]

则在该点附近可令

[math]\displaystyle{ \boldsymbol{\xi}=\mathbf{x}-\mathbf{x}^\ast }[/math]

通过泰勒展开，系统可近似为

<nowiki>[math]\displaystyle{ \frac{d\boldsymbol{\xi}}{dt}=J(\mathbf{x}^\ast)\boldsymbol{\xi} }[/math]</nowiki>

其中

<nowiki>[math]\displaystyle{ J(\mathbf{x}^\ast)=\left.\frac{\partial \mathbf{F}}{\partial \mathbf{x}}\right|_{\mathbf{x}=\mathbf{x}^\ast} }[/math]</nowiki>

称为该非线性系统在不动点处的<nowiki>'''雅可比矩阵'''</nowiki>（Jacobian matrix）。因此，非线性系统在不动点附近的局部行为，通常可以由一个线性常系数系统来近似刻画。

这种方法称为'''线性化'''。通过分析雅可比矩阵的特征值，可以判断不动点附近的局部稳定性。例如，若所有特征值的实部均小于零，则该不动点通常是局部稳定的；若存在特征值的实部大于零，则该不动点通常是不稳定的。若特征值为复数，则系统可能表现出振荡、螺旋收敛或螺旋发散等相空间行为。

因此，线性系统不仅本身具有独立的研究价值，也构成了研究非线性动力系统的重要基础。特别是在二维及更高维相空间中，线性系统可以用于刻画稳定节点、鞍点、焦点、中心等典型局部结构。这些结构为理解更复杂的非线性现象，如分岔、极限环和混沌等，提供了基本分析框架。
== 经典实例：二维谐振子系统 ==
为了更直观地理解线性系统的降维与重构，可以以经典的无阻尼弹簧振子为例。设一个质量为 [math]\displaystyle{ m }[/math] 的物体连接在劲度系数为 [math]\displaystyle{ k }[/math] 的弹簧上，其偏离平衡位置的位移记为 [math]\displaystyle{ x }[/math]。根据牛顿第二定律与胡克定律，该系统满足以下二阶常微分方程：

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

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

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

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

[math]\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{ x }[/math] 这一物理量是不够的，还必须掌握系统当下的速度 [math]\displaystyle{ v }[/math]。因此，描述该系统的自然状态空间不再是一维的直线，而是由位置与速度构成的二维<nowiki>'''</nowiki>相平面<nowiki>'''</nowiki>（Phase plane），即 [math]\displaystyle{ (x,v) }[/math] 空间。

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

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

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

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

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

写成矩阵形式即：

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

=== 特例分析：相爱相杀的循环 ===
作为一个典型情形，考虑如下系统：

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

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

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

计算其迹与行列式可得：

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

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

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

== 参考文献 ==
[1] Alan V. Oppenheim, Alan S. Willsky, and S. Hamid Nawab, ''Signals and Systems'', 2nd ed., Prentice Hall, 1997.

线性系统

2026-04-24T03:53:55Z

张新宇：公式从latex转为了wiki风格

在数学与科学中，'''线性系统（Linear System）'''是指系统的输出变化与输入变化成比例，并且系统内部各要素之间满足'''叠加原理（Superposition Principle）'''的系统。叠加原理包括两个基本性质：'''齐次性'''和'''可加性'''。也就是说，如果某一输入产生某一响应，那么该输入按比例放大或缩小后，系统响应也会按相同比例变化；如果多个输入同时作用于系统，则系统的总响应等于各个输入单独作用时响应的代数和。因此，线性系统的整体行为可以由局部行为组合得到，通常比非线性系统更加直观、规则且易于分析和控制。

从数学形式上看，线性系统通常可以表示为一组'''线性代数方程、线性微分方程或线性差分方程'''。在这些方程中，未知变量、未知函数及其导数只以一次形式出现，不包含未知量之间的乘积、高次幂、指数函数、三角函数等非线性关系。换句话说，线性系统的方程可以写成未知变量或未知函数的线性组合。例如，若两个输入分别产生两个输出，则这两个输入的线性组合也会产生对应输出的线性组合，这正是线性系统区别于非线性系统的核心特征。

需要注意的是，对于线性微分方程而言，判断其是否线性的关键在于'''未知函数及其导数是否线性出现'''。即使方程中的自变量、已知系数或外部输入随时间、空间等因素发生非线性变化，只要未知函数本身及其导数之间不存在乘积、高次项或其他非线性形式，该方程仍然可以称为线性方程。例如，一个系统的参数可以随时间变化，但只要状态变量之间不发生非线性耦合，该系统仍属于线性系统。

由于线性系统具有结构清晰、数学性质良好、求解方法成熟等特点，它们被广泛应用于工程、物理、控制理论、通信系统和信号处理等领域。许多复杂的实际系统虽然本质上可能是非线性的，但在一定工作范围或小扰动条件下，可以通过线性化方法近似为线性系统，从而便于研究其稳定性、响应特性和控制规律。因此，线性系统不仅是一类重要的数学模型，也是理解、分析和设计复杂系统的基础工具。
== 定义 ==
在数学中，若一个映射或算子 f(x) 满足以下两个性质，则称其为'''线性映射'''或'''线性算子'''：

[math]\displaystyle{ f(x+y)=f(x)+f(y) }[/math]

[math]\displaystyle{ f(\alpha x)=\alpha f(x) }[/math]

其中，第一个性质称为'''可加性（Additivity）'''，表示两个输入共同作用时的输出等于它们分别作用时输出的和；第二个性质称为'''齐次性（Homogeneity）'''，表示输入按比例放大或缩小时，输出也按相同比例变化。二者合在一起通常称为'''叠加原理（Superposition Principle）'''，即

[math]\displaystyle{ f(\alpha x+\beta y)=\alpha f(x)+\beta f(y) }[/math]

因此，线性系统的本质在于：系统的整体响应可以由各个输入或状态变量的局部响应通过线性组合得到。

一般地，若一个系统可以写成

[math]\displaystyle{ f(x)=C }[/math]

其中 x 表示待求的未知对象，C 为已知项，并且 f(x) 是关于 x 的线性映射，则该系统称为'''线性系统（Linear System）'''；否则称为'''非线性系统（Nonlinear System）'''。当 C=0 时，该系统称为'''齐次线性系统'''；当 C≠0 时，则称为'''非齐次线性系统'''。

这里的 x 可以是非常一般的数学对象，例如数、向量、函数，或满足一定边界条件的函数。因此，f(x) 也可以表示不同类型的线性运算，例如矩阵变换、积分算子、微分算子或差分算子。如果 f(x) 中包含对未知函数 x 的微分运算，则对应的系统就是'''线性微分方程系统'''。
== 动力系统中的线性系统 ==
在动力学语境下，线性系统通常是研究多维连续动力系统的基本起点。一般地，n 维线性常系数微分方程组可以写成如下矩阵形式：

<nowiki>[math]\displaystyle{ \frac{d\mathbf{x}}{dt}=A\mathbf{x} }[/math]</nowiki>

其中，[math]\displaystyle{ \mathbf{x}\in\mathbb{R}^n }[/math]表示系统的状态向量，A 为常数矩阵，用于描述系统内部各状态变量之间的线性耦合关系。该方程表明，系统状态的变化率由当前状态经过矩阵 A 的线性变换决定。

线性系统的一个核心优势在于，它们往往存在解析解。例如，对于上述常系数线性系统，其解可以形式化地写为

[math]\displaystyle{ \mathbf{x}(t)=e^{At}\mathbf{x}(0) }[/math]

其中，x(0) 为初始状态，[math]\displaystyle{ e^{At} }[/math] 为矩阵指数。通过分析矩阵 A 的特征值和特征向量，可以进一步研究系统的稳定性、振荡行为、相轨迹结构以及长期演化趋势。
=== 向量场与相平面 ===
==== 向量场：方程的几何映射 ====
从前述谐振子的例子可以看出，当线性系统被抽象为 [math]\displaystyle{ \dot{\mathbf{x}}=A\mathbf{x} }[/math] 的形式后，它便超越了纯代数意义下的矩阵微分方程，映射为几何空间中的一个动态结构。对于二维空间中的任意状态点 [math]\displaystyle{ (x,v) }[/math]，矩阵乘法 [math]\displaystyle{ A\mathbf{x} }[/math] 都会赋予该点一个确定的变化速度与方向向量 [math]\displaystyle{ (\dot{x},\dot{v}) }[/math]。这些密布于空间各点、指示着系统瞬时演化趋势的箭头，共同构成了该系统的<nowiki>'''</nowiki>向量场<nowiki>'''</nowiki>（Vector Field）。向量场直观地刻画了系统在全局各个状态下的“流速”与“流向”，是理解二维线性系统整体行为的基础。
==== 相平面与相轨迹 ====
由系统的全部状态变量所张成的几何空间，被称为系统的<nowiki>'''</nowiki>相空间<nowiki>'''</nowiki>（Phase Space）；在二维情况下，特称为<nowiki>'''</nowiki>相平面<nowiki>'''</nowiki>（Phase Plane）。当给定系统一个初始状态 [math]\displaystyle{ (x_0,v_0) }[/math] 时，随着时间 [math]\displaystyle{ t }[/math] 的推移，状态点会顺着向量场的指引在相平面上运动，其描绘出的连续曲线即为系统的<nowiki>'''</nowiki>相轨迹<nowiki>'''</nowiki>（Trajectory）。
==== 定性分析 ====
在复杂系统的动力学研究中，对二维线性系统的研究重点往往不再局限于求解微分方程的显式时间解析解[math]\displaystyle{ \mathbf{x}(t) }[/math]，更重要的是俯瞰这些轨迹在相平面中的整体拓扑分布。通过考察轨迹的几何形状、流动趋势以及它们与特殊点（不动点）的相对关系，可以在不强求解析解的情况下，直接把握系统的宏观定性行为。

对于无阻尼弹簧振子而言，对应到相平面上，轨迹不会相交，而是表现为一系列围绕原点（平衡点）的闭合椭圆。这种首尾相连的封闭轨迹，在几何上完美诠释了该保守系统的能量守恒特性与周期性机械振荡。
=== 不动点与稳定性 ===
在相平面的研究中，除了把握整体轨迹的几何拓扑，另一个核心问题是系统在平衡状态附近的局部演化行为。若系统在某一点处满足 [math]\displaystyle{ \dot{\mathbf{x}}=A\mathbf{x}=\mathbf{0} }[/math]，则称该点为系统的<nowiki>'''</nowiki>不动点<nowiki>'''</nowiki>（Fixed point）或<nowiki>'''</nowiki>平衡点<nowiki>'''</nowiki>（Equilibrium point）。对于非奇异的二维线性系统而言，原点 [math]\displaystyle{ (0,0) }[/math] 通常是最基本的不动点。
==== 稳定性的层级分类 ====
在线性系统中，常用的稳定性概念包括：
* 吸引的 (Attracting)：若从不动点附近出发的轨迹在 [math]\displaystyle{ t\to\infty }[/math] 时都收敛到该点。

* 李雅普诺夫稳定 (Lyapunov stable)：若任意足够靠近不动点的初始条件产生的轨迹，在未来始终保持在给定的较小邻域内。

* 渐近稳定 (Asymptotically stable)：若一个不动点同时满足上述两个条件。

* 中性稳定 (Neutrally stable)：轨迹始终保持在平衡点附近但不向其收敛（如无阻尼谐振子的闭合椭圆）。

* 不稳定 (Unstable)：无论扰动多么微小，轨迹都可能逐渐远离平衡点。
==== 基于特征值、迹与行列式的分类 ====
对于二维线性系统，局部几何结构可以由矩阵 [math]\displaystyle{ A }[/math] 的特征值完全刻画。设矩阵 [math]\displaystyle{ A }[/math] 的两个特征值为 [math]\displaystyle{ \lambda_1,\lambda_2 }[/math]，其迹（Trace）与行列式（Determinant）分别为：

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

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

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

因此，二维线性系统不动点的类型可统一由 [math]\displaystyle{ T }[/math]、[math]\displaystyle{ \Delta }[/math] 以及判别式 [math]\displaystyle{ T^2-4\Delta }[/math] 来决定。[[文件:二维线性系统不动点分类.png|缩略图|链接=https://wiki.swarma.org/index.php/%E6%96%87%E4%BB%B6:%E4%BA%8C%E7%BB%B4%E7%BA%BF%E6%80%A7%E7%B3%BB%E7%BB%9F%E4%B8%8D%E5%8A%A8%E7%82%B9%E5%88%86%E7%B1%BB.png]]
== 在线性化和非线性系统分析中的作用 ==
尽管现实世界中的许多系统具有非线性特征，但线性系统在非线性系统分析中仍然具有基础性作用。一个重要原因是，许多非线性动力系统在平衡点或不动点附近的局部行为，可以通过相应的线性系统近似描述。

设一个连续时间非线性动力系统为

<nowiki>[math]\displaystyle{ \frac{d\mathbf{x}}{dt}=\mathbf{F}(\mathbf{x}) }[/math]</nowiki>

其中，[math]\displaystyle{ \mathbf{x}\in\mathbb{R}^n }[/math] 为状态向量，[math]\displaystyle{ \mathbf{F}(\mathbf{x}) }[/math] 为非线性向量场。若 [math]\displaystyle{ \mathbf{x}^\ast }[/math] 是系统的不动点，即

[math]\displaystyle{ \mathbf{F}(\mathbf{x}^\ast)=\mathbf{0} }[/math]

则在该点附近可令

[math]\displaystyle{ \boldsymbol{\xi}=\mathbf{x}-\mathbf{x}^\ast }[/math]

通过泰勒展开，系统可近似为

<nowiki>[math]\displaystyle{ \frac{d\boldsymbol{\xi}}{dt}=J(\mathbf{x}^\ast)\boldsymbol{\xi} }[/math]</nowiki>

其中

<nowiki>[math]\displaystyle{ J(\mathbf{x}^\ast)=\left.\frac{\partial \mathbf{F}}{\partial \mathbf{x}}\right|_{\mathbf{x}=\mathbf{x}^\ast} }[/math]</nowiki>

称为该非线性系统在不动点处的<nowiki>'''雅可比矩阵'''</nowiki>（Jacobian matrix）。因此，非线性系统在不动点附近的局部行为，通常可以由一个线性常系数系统来近似刻画。

这种方法称为'''线性化'''。通过分析雅可比矩阵的特征值，可以判断不动点附近的局部稳定性。例如，若所有特征值的实部均小于零，则该不动点通常是局部稳定的；若存在特征值的实部大于零，则该不动点通常是不稳定的。若特征值为复数，则系统可能表现出振荡、螺旋收敛或螺旋发散等相空间行为。

因此，线性系统不仅本身具有独立的研究价值，也构成了研究非线性动力系统的重要基础。特别是在二维及更高维相空间中，线性系统可以用于刻画稳定节点、鞍点、焦点、中心等典型局部结构。这些结构为理解更复杂的非线性现象，如分岔、极限环和混沌等，提供了基本分析框架。
== 经典实例：二维谐振子系统 ==
为了更直观地理解线性系统的降维与重构，可以以经典的无阻尼弹簧振子为例。设一个质量为 [math]\displaystyle{ m }[/math] 的物体连接在劲度系数为 [math]\displaystyle{ k }[/math] 的弹簧上，其偏离平衡位置的位移记为 [math]\displaystyle{ x }[/math]。根据牛顿第二定律与胡克定律，该系统满足以下二阶常微分方程：

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

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

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

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

[math]\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{ x }[/math] 这一物理量是不够的，还必须掌握系统当下的速度 [math]\displaystyle{ v }[/math]。因此，描述该系统的自然状态空间不再是一维的直线，而是由位置与速度构成的二维<nowiki>'''</nowiki>相平面<nowiki>'''</nowiki>（Phase plane），即 [math]\displaystyle{ (x,v) }[/math] 空间。

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

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

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

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

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

写成矩阵形式即：

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

=== 特例分析：相爱相杀的循环 ===
作为一个典型情形，考虑如下系统：

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

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

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

计算其迹与行列式可得：

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

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

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

线性系统

2026-04-24T03:34:59Z

张新宇：未完成的版本2-线性系统，只是基本上的框架完成，还在修改阶段

在数学与科学中，'''线性系统（Linear System）'''是指系统的输出变化与输入变化成比例，并且系统内部各要素之间满足'''叠加原理（Superposition Principle）'''的系统。叠加原理包括两个基本性质：'''齐次性'''和'''可加性'''。也就是说，如果某一输入产生某一响应，那么该输入按比例放大或缩小后，系统响应也会按相同比例变化；如果多个输入同时作用于系统，则系统的总响应等于各个输入单独作用时响应的代数和。因此，线性系统的整体行为可以由局部行为组合得到，通常比非线性系统更加直观、规则且易于分析和控制。

从数学形式上看，线性系统通常可以表示为一组'''线性代数方程、线性微分方程或线性差分方程'''。在这些方程中，未知变量、未知函数及其导数只以一次形式出现，不包含未知量之间的乘积、高次幂、指数函数、三角函数等非线性关系。换句话说，线性系统的方程可以写成未知变量或未知函数的线性组合。例如，若两个输入分别产生两个输出，则这两个输入的线性组合也会产生对应输出的线性组合，这正是线性系统区别于非线性系统的核心特征。

需要注意的是，对于线性微分方程而言，判断其是否线性的关键在于'''未知函数及其导数是否线性出现'''。即使方程中的自变量、已知系数或外部输入随时间、空间等因素发生非线性变化，只要未知函数本身及其导数之间不存在乘积、高次项或其他非线性形式，该方程仍然可以称为线性方程。例如，一个系统的参数可以随时间变化，但只要状态变量之间不发生非线性耦合，该系统仍属于线性系统。

由于线性系统具有结构清晰、数学性质良好、求解方法成熟等特点，它们被广泛应用于工程、物理、控制理论、通信系统和信号处理等领域。许多复杂的实际系统虽然本质上可能是非线性的，但在一定工作范围或小扰动条件下，可以通过线性化方法近似为线性系统，从而便于研究其稳定性、响应特性和控制规律。因此，线性系统不仅是一类重要的数学模型，也是理解、分析和设计复杂系统的基础工具。

== 定义 ==
在数学中，若一个映射或算子 f(x) 满足以下两个性质，则称其为'''线性映射'''或'''线性算子'''：

f(x+y)=f(x)+f(y),

f(αx)=αf(x).

其中，第一个性质称为'''可加性（Additivity）'''，表示两个输入共同作用时的输出等于它们分别作用时输出的和；第二个性质称为'''齐次性（Homogeneity）'''，表示输入按比例放大或缩小时，输出也按相同比例变化。二者合在一起通常称为'''叠加原理（Superposition Principle）'''，即

f(αx+βy)=αf(x)+βf(y).

因此，线性系统的本质在于：系统的整体响应可以由各个输入或状态变量的局部响应通过线性组合得到。

一般地，若一个系统可以写成

f(x)=C,

其中 x 表示待求的未知对象，C 为已知项，并且 f(x) 是关于 x 的线性映射，则该系统称为'''线性系统（Linear System）'''；否则称为'''非线性系统（Nonlinear System）'''。当 C=0 时，该系统称为'''齐次线性系统'''；当 C≠0 时，则称为'''非齐次线性系统'''。

这里的 x 可以是非常一般的数学对象，例如数、向量、函数，或满足一定边界条件的函数。因此，f(x) 也可以表示不同类型的线性运算，例如矩阵变换、积分算子、微分算子或差分算子。如果 f(x) 中包含对未知函数 x 的微分运算，则对应的系统就是'''线性微分方程系统'''。

== 动力系统中的线性系统 ==
在动力学语境下，线性系统通常是研究多维连续动力系统的基本起点。一般地，n 维线性常系数微分方程组可以写成如下矩阵形式：

dx/dt = Ax,

其中，x∈Rn 表示系统的状态向量，A 为常数矩阵，用于描述系统内部各状态变量之间的线性耦合关系。该方程表明，系统状态的变化率由当前状态经过矩阵 A 的线性变换决定。

线性系统的一个核心优势在于，它们往往存在解析解。例如，对于上述常系数线性系统，其解可以形式化地写为

x(t)=e^{At}x(0),

其中，x(0) 为初始状态，e^{At} 为矩阵指数。通过分析矩阵 A 的特征值和特征向量，可以进一步研究系统的稳定性、振荡行为、相轨迹结构以及长期演化趋势。

=== 向量场与相平面 ===

==== 向量场：方程的几何映射 ====
<nowiki>从前述谐振子的例子可以看出，当线性系统被抽象为 $\dot{\mathbf{x}} = A\mathbf{x}$ 的形式后，它便超越了纯代数意义下的矩阵微分方程，映射为几何空间中的一个动态结构。对于二维空间中的任意状态点 $(x, v)$，矩阵乘法 $A\mathbf{x}$ 都会赋予该点一个确定的变化速度与方向向量 $(\dot{x}, \dot{v})$。这些密布于空间各点、指示着系统瞬时演化趋势的箭头，共同构成了该系统的\textbf{向量场}（Vector Field）。向量场直观地刻画了系统在全局各个状态下的“流速”与“流向”，是理解二维线性系统整体行为的基础。</nowiki>

==== 相平面与相轨迹 ====
由系统的全部状态变量所张成的几何空间，被称为系统的\textbf{相空间}（Phase Space）；在二维情况下，特称为\textbf{相平面}（Phase Plane）。当给定系统一个初始状态 $(x_0, v_0)$ 时，随着时间 $t$ 的推移，状态点会顺着向量场的指引在相平面上游走，其描绘出的连续曲线即为系统的\textbf{相轨迹}（Trajectory）。

==== 定性分析 ====
在复杂系统的动力学研究中，对二维线性系统的研究重点往往不再局限于求解微分方程的显式时间解析解 $\mathbf{x}(t)$，更重要的是俯瞰这些轨迹在相平面中的整体拓扑分布。通过考察轨迹的几何形状、流动趋势以及它们与特殊点（不动点）的相对关系，可以在不强求解析解的情况下，直接把握系统的宏观定性行为。

对于无阻尼弹簧振子而言，对应到相平面上，轨迹不会相交，而是表现为一系列围绕原点（平衡点）的闭合椭圆。这种首尾相连的封闭轨迹，在几何上完美诠释了该保守系统的能量守恒特性与周期性机械振荡。

=== 不动点与稳定性 ===
<nowiki>在相平面的研究中，除了把握整体轨迹的几何拓扑，另一个核心问题是系统在平衡状态附近的局部演化行为。若系统在某一点处满足 $\dot{\mathbf{x}} = A\mathbf{x} = \mathbf{0}$，则称该点为系统的\textbf{不动点}（Fixed point）或\textbf{平衡点}（Equilibrium point）。对于非奇异的二维线性系统而言，原点 $(0,0)$ 通常是最基本的不动点。</nowiki>

==== 稳定性的层级分类 ====
在线性系统中，常用的稳定性概念包括：

* 吸引的 (Attracting)：若从不动点附近出发的轨迹在 $t \to \infty$ 时都收敛到该点。

* 李雅普诺夫稳定 (Lyapunov stable)：若任意足够靠近不动点的初始条件产生的轨迹，在未来始终保持在给定的较小邻域内。

* 渐近稳定 (Asymptotically stable)：若一个不动点同时满足上述两个条件。

* 中性稳定 (Neutrally stable)：轨迹始终保持在平衡点附近但不向其收敛（如无阻尼谐振子的闭合椭圆）。

* 不稳定 (Unstable)：无论扰动多么微小，轨迹都可能逐渐远离平衡点。

==== 基于特征值、迹与行列式的分类 ====
对于二维线性系统，局部几何结构可以由矩阵 $A$ 的特征值完全刻画。设矩阵 $A$ 的两个特征值为 $\lambda_1, \lambda_2$，其迹（Trace）与行列式（Determinant）分别为：

$$ T = \lambda_1 + \lambda_2, \quad \Delta = \lambda_1 \lambda_2 $$

根据特征方程 $\lambda^2 - T\lambda + \Delta = 0$，由求根公式可得：

<nowiki>$$ \lambda_{1,2} = \frac{T \pm \sqrt{T^2 - 4\Delta}}{2} $$</nowiki>

因此，二维线性系统不动点的类型可统一由 $T, \Delta$ 以及判别式 $T^2 - 4\Delta$ 来决定。
[[文件:二维线性系统不动点分类.png|缩略图]]

== 在线性化和非线性系统分析中的作用 ==
尽管现实世界中的许多系统具有非线性特征，但线性系统在非线性系统分析中仍然具有基础性作用。一个重要原因是，许多非线性动力系统在平衡点或不动点附近的局部行为，可以通过相应的线性系统近似描述。

设一个连续时间非线性动力系统为

dx/dt=F(x),

其中 x∈Rn 为状态向量，F(x) 为非线性向量场。若 x∗ 是系统的不动点，即

F(x∗)=0,

则在该点附近可令

ξ=x−x∗.

通过泰勒展开，系统可近似为

dξ/dt=J(x∗)ξ,

其中

J(x∗)=∂F/∂x|x=x∗

称为该非线性系统在不动点处的'''雅可比矩阵'''。因此，非线性系统在不动点附近的局部行为，通常可以由一个线性常系数系统来近似刻画。

这种方法称为'''线性化'''。通过分析雅可比矩阵的特征值，可以判断不动点附近的局部稳定性。例如，若所有特征值的实部均小于零，则该不动点通常是局部稳定的；若存在特征值的实部大于零，则该不动点通常是不稳定的。若特征值为复数，则系统可能表现出振荡、螺旋收敛或螺旋发散等相空间行为。

因此，线性系统不仅本身具有独立的研究价值，也构成了研究非线性动力系统的重要基础。特别是在二维及更高维相空间中，线性系统可以用于刻画稳定节点、鞍点、焦点、中心等典型局部结构。这些结构为理解更复杂的非线性现象，如分岔、极限环和混沌等，提供了基本分析框架。

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

$$ m\ddot{x} + kx = 0 $$

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

$$ \dot{v} = -\frac{k}{m}x $$

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

$$ \begin{pmatrix} \dot{x} \\ \dot{v} \end{pmatrix} = \begin{pmatrix} 0 & 1 \\ -\frac{k}{m} & 0 \end{pmatrix} \begin{pmatrix} x \\ v \end{pmatrix} $$

从该演化方程可以看出，要完全确定系统未来的状态，仅仅知道位置 $x$ 这个物理量是不够的，还必须掌握系统当下的速度 $v$。因此，描述该系统的自然状态空间不再是一维的直线，而是由位置与速度构成的二维\textbf{相平面}（Phase plane），即 $(x, v)$ 空间。

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

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

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

$$ \dot{R} = aR + bJ $$

$$ \dot{J} = cR + dJ $$

写成矩阵形式即：

$$ \begin{pmatrix} \dot{R} \\ \dot{J} \end{pmatrix} = \begin{pmatrix} a & b \\ c & d \end{pmatrix} \begin{pmatrix} R \\ J \end{pmatrix} $$

其中，参数 $a, d$ 描述个体对自身当前情感状态的内在反应；参数 $b, c$ 描述双方对彼此情感变化的外部响应方式。不同参数符号和大小的组合，对应于截然不同的“恋爱风格”和长期演化结局。这一模型深刻说明：即使系统的变量来自社会或心理意义上的抽象量，其数学分析仍然遵循与物理系统完全相同的线性动力学框架。

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

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

$$ \dot{R} = bJ, \quad \dot{J} = cR $$

其中，$b > 0, c < 0$。这表示罗密欧主要依据朱丽叶的态度来调整自己（正反馈）：朱丽叶越爱他，他越热情；而朱丽叶则对罗密欧的接近表现出回避倾向（负反馈）：罗密欧越热情，她越想后退。此时系统矩阵为：

$$ A = \begin{pmatrix} 0 & b \\ c & 0 \end{pmatrix} $$

计算其迹与行列式：

$$ T = 0, \quad \Delta = -bc > 0 $$

根据前文的分类判别式，该系统的特征值为一对纯虚数，原点对应于一个\textbf{中心（Center）}。在相平面上，轨迹表现为围绕原点的闭合曲线。这意味着两人的情感状态既不会收敛到平淡的静止平衡，也不会演变为无限发散的狂热，而是呈现出周期性的“爱恨循环”。

这一案例完美地展示了动力学中的“同构性”：二维线性系统的不动点分类不仅能够解释弹簧振子等经典物理系统中的能量守恒与周期振动，也同样适用于刻画社会学中抽象的反作用力。在数学上，这两类问题都归结为对矩阵特征值、迹与行列式的分析；而在几何上，它们都可以通过相平面中的轨迹拓扑得到最直观的解答。

文件:二维线性系统不动点分类.png

2026-04-24T03:28:56Z

张新宇：

二维线性系统不动点分类的 $T-\Delta$ (迹-行列式) 参数平面图

线性系统

2026-04-24T03:22:31Z

张新宇：未完成的版本1-线性系统