Skip to content

060 复特征值 Complex-Valued Eigenvalues

这一节继续讨论 个方程组成的常系数线性齐次方正组 其中系数矩阵 是实数矩阵。如果解的形式是 ,那么 是系数方程 的特征值和特征向量。 的特征值 是特征方程 的根,相应的非零特征向量满足 如果 的系数是实数,多项式方程 的系数是实数,任意复数根会共轭出现。如果 的特征向量,其中 是实数,那么 也是 的特征向量。在正式讨论之前先看个例子。

例 1 求方程组 的实数基础解系。给出相图和一些解的曲线。

解:下图是方程组 的方向场。解的轨迹顺时针向原点旋转。

假定解的形式是 那么得到线性方程组 特征方程是 因此特征值是 ,代入 得到对应的特征向量是 从上面可以看出特征向量 也互为复共轭。因此方程组 的基础解析是 为了获得实数解,可以取 的实部和虚部。那么 因此方程组 的实数解是 它们的朗斯基是 不会为零,因为 是线性无关的,是方程组 的基础解系。

的图分别在下图中黑色实线和虚线表示。

因为 因此 分别会通过 。方程组 的其他解是这两个解的线性组合,其中一些如上图所示。当 时,每个解都沿着螺旋线趋于原点,并绕着原点旋转无数次,这是因为解 是一个衰减的指数和正弦、余弦函数的乘积。下图是 变化的图像如下(b)所示,每一条曲线代表一种随时间衰减的振荡。 变化的图像如(c)所示。

上面的(a)图是一个典型的特征值是共轭复数且实部为负数的 的方程组 的解的图像。原点称为螺旋点(spiral point),随着 的增加,所有轨迹都趋于原点,因此是渐进稳定的。如果复数的实部是正数,图像也和(a)类似,不过轨迹方向是远离原点的,因此轨迹会变得无穷大,此时原点是不稳定的。

如果实部为零,那么既不会趋于原点也不会远离原点,轨迹是周而复始的绕着原点的曲线。后续会有这样的例子。此时原点称为中心(center),也被称为是稳定的,但不是渐进稳定的。这三种情况下,运动方向可能是顺时针,也可能是逆时针,这依赖于系数矩阵

上面的图是计算机辅助作图,但是人工可以给出草图。假设特征值是 ,如果 ,向内旋转,,向外旋转,如果 ,那么是闭合曲线。为了确定是顺时针还是逆时针,只需要确定一个点的方向即可。比如 中我们选择 ,那么 。因此该点处切线方向 确定了, 方向上是正的,因此轨迹从第二象限到第一象限,那么是顺时针运动。

回到线性微分方程组 假定有一组共轭复数特征值 ,相应的特征向量是 也是共轭的。 满足 是实数矩阵,取上面方程的共轭得到 其中 的共轭,那么 。那么解是 它们彼此复共轭。和例 1 一样,可以取 的实部和虚部,得到方程 的实数解。

,其中 是实数,那么 分离实部和虚部得到 ,那么 是方程 的实数解。 是线性无关的。下面给出证明。

首先证明 是线性无关的,使用 来表示 代入 得到 那么 由于 线性无关,那么 因此 ,这就得到了 是线性无关的。

下面证明 是线性无关的,那么需要证明 。按照 合并同类项得到 因为 线性无关,那么 解这个方程组得到 ,因此 线性无关。

假定矩阵 有共轭特征值 是实数特征值且不相等,特征向量是 。那么方程组 的通解是 其中 给出。这里要强调 矩阵都是实数,因为这样复数特征值和特征向量才是共轭的。

对于 的实系数方程组,有如下三种情况:

  • 特征值是实数且符号相反, 是鞍点
  • 特征值是实数且符号相同,且不相等, 是节点
  • 特征值是复数且实部不为零, 是螺旋点

其他可能情况是上述两种情况之间的过渡状态,不过现实应用中不常遇到。比如零特征值对应鞍点和节点之间。纯虚特征值的情况介于渐进稳定螺旋点和不稳定螺旋点之间。实数且相等的情况介于节点和螺旋点之间。

例 2 方程组 包含参数 ,描述系统的解是如何依赖于 的,特别是找到 的分叉点,也就是轨迹的行为发生显著变化的点。

解:轨迹的行为完全由特征值确定。特征方程是 那么特征值是 根据上式可知当 是复共轭特征值,其他时候特征值是实数。

是两个分叉点,特征值从实数变为复数,或者反之。对于 ,特征值都是负数,轨迹趋于原点(渐进稳定点)。 时,特征值都是正数,轨迹远离原点(非稳定点)。 时,特征值是复数,轨迹是螺旋线。当 时,特征值实部是负数,螺旋线向内,当 时,特征值实部是正数,螺旋线向内。

第三个分叉点是 ,螺旋线的旋转方向发生了变化。当 时,轨迹是绕着原点的封闭曲线。在分叉点 时,特征值是实数并且相等。这种情况下原点还是节点,但是相图和 7.5 有所不同,在 7.8 小节讨论。

多弹簧系统

下面讨论 7.1 小节开始时由两个物体三个弹簧组成的系统。假定没有外力,那么 那么方程组是 这个方程组可以当做两个二阶方程的方程组求解,也可以变成四个一阶方程组成的方正组。令 ,那么 可以写作

例 3 假定 ,那么这些方程是 分析运行,画出典型轨迹图。

解:令 ,那么 写成矩阵形式是 两个物体相对平衡位置的位移, 是它们的速度。假定解 ,其中 是矩阵 的特征值, 是对应的特征向量。可以手算,也可以使用计算机辅助计算。得到 的特征方程是 因此特征值是 ,对应的特征向量是 是共轭解,可以找到其实部和虚部得到实数解。有 类似的 线性无关,因此是基础解系。方程 的通解是 其中 是任意常量。

这个方程组的相空间是四维的,每一解由一组特定的 确定,对应着四维空间的一条轨迹。 的解的周期是 ,因此轨迹是闭合曲线。不管 是从哪个点开始,在 时回到起点然后循环。这里无法给出四维空间的区间,而是给出在 平面或者 平面上的投影。

首先分析 的前两项,其频率是 1,周期是 。注意 ,这意味着两个物体始终同向运动,位移也同向,第二个物体的位移和速度都是第一个物体的三分之二。聚焦于 ,如果画出 相对于 的曲线,得到的是振幅为 3 和 2 的余弦曲线,如下图(a)所示。曲线在 平面的投影是半径为 3 的圆,如下图(b)所示,从 点开始顺时针运动,周期是 。(b)图也展示了曲线在 平面的投影,是半径为 2 的圆,从 开始顺时针运动,周期也是 。原点是中心点。分析 或者是 的线性组合可以得到类似的图像。

的后两项描述了频率为 1 周期为 的运动。类似的,。这说明两个物体运动反向,第二个物体的位移和速度都是第一个物体的三分之四。和之前类似,分析 ,画出 相对于 的曲线,如下图(a)所示。相位差是 的振幅是 的三分之四。(b)图展示了两个物体在各自相平面的轨迹投影。两个图像都是椭圆,内圈从 开始,外圈从 开始,顺时针运动,周期是 。原点是中心。分析 或者是 的线性组合可以得到类似的图像。

前面两段描述了两个物体组成的系统的基本模式(fundamental mode)。每一个都需要特殊的初始条件。比如为了得到频率为 1 的情况, 中的 应该为零。这意味着初始条件是 。类似的,当 为零时得到的是频率为 2 的模式,此时初始条件是

如果是更一般的初始条件,解是两种基本模式的组合。比如下图(a)是 相对 的变化。

下图(b)是轨迹在 上的投影。

后面这张图有一点误导性,轨迹的投影有相交的点,但是这在四维空间中是不会发生的,因为这违反了解的唯一性:不能有两个不同的解从同一个点开始。