微分方程 · 相图

相图

相图把微分方程组的每个解画成变量平面上的一条曲线。对线性方程组x′ = Ax,A的迹和行列式给出它的特征值;除非两个特征值相等,特征值就决定了这是哪一类图。

相图与本页画面

相图把微分方程组的解画成变量平面上的曲线,每个起点引出一条。本页画的是线性方程组x′=Axx' = Ax,AA是2×22 \times 2矩阵。在迹–行列式地图上移动点(τ,Δ)(\tau, \Delta),画面会经过鞍点、结点、螺旋和中心;旁边的特征值平面画出每幅图背后的两个特征值。

方程组、它的迹和行列式
(xy)′=(abcd)(xy),τ=a+d,Δ=ad−bc\begin{pmatrix} x \\ y \end{pmatrix}' = \begin{pmatrix} a & b \\ c & d \end{pmatrix} \begin{pmatrix} x \\ y \end{pmatrix}, \qquad \tau = a + d, \quad \Delta = ad - bc

每条曲线都顺着方向场走。方向场的箭头一样长,只表示方向;快慢要看等时点,运动慢的地方点就密。有些解沿着过原点的直线走。如果Av=λvAv = \lambda v,那么x(t)=eλtvx(t) = e^{\lambda t} v就是一个解,因为它的导数λeλtv\lambda e^{\lambda t} v正好等于AA作用在它上面。所以这些直线就是AA的特征方向,λ\lambda是沿这条线的增长率:λ>0\lambda > 0时往外走,λ<0\lambda < 0时往里走,λ=0\lambda = 0时不动。

特征值是λ2−τλ+Δ=0\lambda^2 - \tau\lambda + \Delta = 0的根,所以τ\tau是两根之和,Δ\Delta是两根之积,τ2−4Δ\tau^2 - 4\Delta是判别式。地图上每条分界线都对应特征值平面上的一件事。Δ=0\Delta = 0以下,两个特征值异号,画面是鞍点;在这条线上有一个特征值是0,一整条直线上的点都不动。在抛物线τ2=4Δ\tau^2 = 4\Delta上两个特征值相遇,抛物线以上变成一对共轭复数。Δ>0\Delta > 0时,τ\tau的正负区分衰减和增长;在τ=0\tau = 0上,这对复数落在虚轴上。所以τ\tau和Δ\Delta决定了特征值,也决定了相图的类型;只有特征值重合时,还要看AA是不是单位矩阵的倍数。

特征值是复数α±βi\alpha \pm \beta i(β>0\beta > 0)时,每个非零解绕原点转一圈用时2π/β2\pi/\beta,这段时间里整体缩放e2πα/βe^{2\pi\alpha/\beta}倍。所以相隔一圈的等时点落在从原点出发的同一条射线上,每个都按同一个倍数缩放。α=0\alpha = 0时倍数是1,每条非零轨道都闭合成椭圆,这就是中心。取Δ=1\Delta = 1、A=(01−1τ)A = \begin{pmatrix} 0 & 1 \\ -1 & \tau \end{pmatrix},方程组就是阻尼振子x′′+2ζx′+x=0x'' + 2\zeta x' + x = 0,其中ζ=−τ/2\zeta = -\tau/2。

穿过抛物线时,两条直线解像剪刀一样合成一条,随后消失,画面变成螺旋。除了抛物线的顶点,抛物线上一般的矩阵只有这一条直线解,其他解进出原点时都与它相切。在(τ,Δ)=(−2,1)(\tau, \Delta) = (-2, 1)处,矩阵−I-I的特征值相同,每个方向却都是直线解。换基A↦PAP−1A \mapsto PAP^{-1}会改变画面的形状,但保持τ\tau、Δ\Delta、特征值和类型。单位矩阵的倍数换基后不变。

值得注意

  • 结点的两条直线解一般并不垂直。不在这两条直线上的解,在原点附近沿慢的那条进出,走远后转向快的那条。
  • 蓝色曲线是过起点的整条解。黄色路径是放下起点以来走过的一段,所以AA改变后,它显示的是新的解,不是旧解的历史。

操作

  • 在相平面空白处按一下可以加一个起点。靠近直线解时,它会落到直线上。

相关概念

延伸阅读:维基百科:相图 (动态系统)。