微分方程 · 相圖

相圖

相圖將微分方程組的每個解畫成變數平面上的曲線。對線性系統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改變後,它顯示的是新系統的解,而非舊解的歷史紀錄。

操作方式

  • 在相平面空白處按一下即可新增起點。靠近直線解時,起點會吸附到直線上。

相關內容

延伸閱讀: 維基百科: Phase portrait (英文).