查看“︁线性动态系统”︁的源代码

{{unreferenced|time=2017-01-09T12:05:48+00:00}}

'''线性动力系统'''是演化函数满足[[线性关系]]的[[动力系统]]。一般来说，动力系统没有[[解析解]]，而线性动力系统则可以精确求解，且具有丰富的数学性质。线性系统还可用于理解一般动力系统的定性行为，方法是计算系统[[平衡点 (数学)|平衡点]]，并将其近似为每个平衡点周围的线性系统。
==概述==

线性动力系统中，状态向量（<math>N</math>维向量，记作<math>\mathbf{x}</math>）的变化等于常数矩阵（记作<math>\mathbf{A}</math>）乘以
<math>\mathbf{x}</math>。这种变化有两种形式：[[流 (数学)|流]]，即<math>\mathbf{x}</math>随时间连续变化

:<math>
\frac{d}{dt} \mathbf{x}(t) = \mathbf{A} \mathbf{x}(t)
</math>

或映射，其中
<math>\mathbf{x}</math>随时间离散变化

:<math>
\mathbf{x}_{m+1} = \mathbf{A} \mathbf{x}_{m}
</math>

这些方程在以下意义上是线性的：如若
<math>\mathbf{x}(t)</math>、<math>\mathbf{y}(t)</math> 
是两个有效解，那么它们的任意[[线性组合]]也是有效解，如
<math>\mathbf{z}(t) \ \stackrel{\mathrm{def}}{=}\  \alpha \mathbf{x}(t) + \beta \mathbf{y}(t)</math> 
，其中<math>\alpha</math>、<math>\beta</math>
为任意[[标量]]。矩阵<math>\mathbf{A}</math>不必对称。

线性动力系统可精确求解，而大多数非线性动力系统不能。偶尔，非线性系统也可将变量变换为线性系统来精确求解。此外，（几乎）任何非线性系统的解都可通过其定点附近的等效线性系统很好地逼近。因此，理解线性系统及其解法是理解更复杂的非线性系统的关键一步。

==线性动力系统的解==

若初向量<math>\mathbf{x}_{0} \ \stackrel{\mathrm{def}}{=}\  \mathbf{x}(t=0)</math>
与矩阵<math>\mathbf{A}</math>的一个右特征向量<math>\mathbf{r}_{k}</math>对齐，那么动力系统就是

:<math>
\frac{d}{dt} \mathbf{x}(t) = 
\mathbf{A} \mathbf{r}_{k} = \lambda_{k} \mathbf{r}_{k}
</math>

其中<math>\lambda_{k}</math>是相应的特征值；方程的解为
:<math>
\mathbf{x}(t) = 
\mathbf{r}_{k} e^{\lambda_{k} t}
</math>
可以通过替换来确认。

若<math>\mathbf{A}</math>是[[可对角化矩阵]]，那么<math>N</math>维空间中的任何向量都可用矩阵<math>\mathbf{A}</math>的左右特征向量（记作<math>\mathbf{l}_{k}</math>）的[[线性组合]]表示。
:<math>
\mathbf{x}_{0} = 
\sum_{k=1}^{N} 
\left( \mathbf{l}_{k} \cdot \mathbf{x}_{0} \right)
\mathbf{r}_{k}
</math>

因此<math>\mathbf{x}(t)</math>的一般解是右特征向量各解的线性组合
:<math>
\mathbf{x}(t) = 
\sum_{k=1}^{n} 
\left( \mathbf{l}_{k} \cdot \mathbf{x}_{0} \right)
\mathbf{r}_{k} e^{\lambda_{k} t}
</math>

离散情形同理。

== 二维中的分类 ==
[[Image:LinDynSysTraceDet.jpg|right|thumb|400px|非线性系统的线性近似：根据[[雅各布矩阵]]的迹和行列式（平衡点附近系统的线性化）对二维定点进行分类。]]
[[特征多项式]]det('''A''' - &lambda;'''I''')的根就是'''A'''的特征值。根<math>\lambda_n</math>的符号与相互关系可确定动力系统
:<math>
\frac{d}{dt} \mathbf{x}(t) = \mathbf{A} \mathbf{x}(t).
</math>
的稳定性。对于2维系统，特征多项式形式为<math>\lambda^2-\tau\lambda+\Delta=0</math>，其中<math>\tau</math>是[[迹]]，<math>\Delta</math>是'''A'''的[[行列式]]。因此，这两个根的形式是
:<math>\lambda_1=\frac{\tau+\sqrt{\tau^2-4\Delta}}{2}</math>
:<math>\lambda_2=\frac{\tau-\sqrt{\tau^2-4\Delta}}{2}</math>,
且<math>\Delta=\lambda_1\lambda_2</math>，<math>\tau=\lambda_1+\lambda_2</math>。因此，如果<math>\Delta<0</math>，则特征值符号相反，定点为鞍点。若<math>\Delta>0</math>，则特征值符号相同。因此，若<math>\tau>0</math>则两个特征值都为正，该点不稳定；若<math>\tau<0</math>则两个特征值都为负，该点稳定。判别式将说明该点是节点型还是螺旋型（即特征值是实值还是复值）。

==相關條目==
*[[線性系統]]
*[[动力系统]]
*{{le|動力系統及微分方程主題列表|List of dynamical systems and differential equations topics}}

*[[矩阵微分方程]]
[[Category:动力系统]]