相平面
相平面(phase plane)是在應用數學(特別是非線性系統)中,視覺化的展示特定微分方程特徵的方式。相平面是一個由二個狀態變數為座標軸組成的平面,例如說(x, y)或(q, p)等。相平面是多維度相空間在二維空間中的例子。
相平面法(phase plane method)是指用繪圖的方式,來確認微分方程的解中是否存在極限環。
微分方程的解可以形成函數族。用繪圖的方式,可以畫在二維的相平面上,類似二維的向量場。向量會表示某一點對應特定參數(例如時間)的導數,也就是(dx/dt, dy/dt),會繪製在對應的點上,以箭頭表示。若有夠多的點,就可以分析此區域內的系統行為,若有極限環,也可以識別出來。
整個場即可形成相圖,在流線上的特定路徑(一個永遠和向量相切的路徑)即為相路徑(phase path)。向量場上的相表示微分方程所說明的系統隨時間的演化。
相平面可以用來解析物理系統的行為,特別是振盪系統,例如獵食者-獵物模型(可參考洛特卡-沃爾泰拉方程)。這些模型中的相路徑可能是向內旋轉,慢慢趨近0,也可能是向外旋轉,慢慢趨近無限大,或是接近中性的平衡位置,此情形稱為centre,路徑可能是圓形、橢圓或是其他形狀。在判斷其系統是否穩定時很有用[1]。
另一個振盪系統的例子是一些多步的化學反應,其中有些會有化學平衡,不是完全反應。此情形下可以將反應物及生成物濃度(或質量)的變化利用微分方程來建模,可以對其化學動力學有更清楚的瞭解[2]。
線性系統的例子
[編輯]可以整理為矩陣方程式:
其中A是2 × 2的係數矩陣,而x = (x, y)是二個自變量組成的座標向量。
此系統可以解析求解,例如積分下式[3]:
用特徵值求解
[編輯]上述方程常見的解法是用右邊矩陣型式的係數,利用特徵值求解,利用以下的行列式求得:
以及以下的特徵向量
特徵值表示指數項的冪次,而特徵向量為其係數。若將解寫成代數型式,可以表示為幾個指數項配合對應係數的和。因為特徵向量的唯一性,每一個此方式得到的解會有待確定的係數c1, c2, ... cn.
通解為:
其中λ1和λ2是特徵值,而(k1, k2), (k3, k4)是基礎特徵向量。係數c1和c2和特徵向量的唯一性有關,只有在初始條件已知時才能求解。
上述行列式可以得到以下的特徵多項式:
是以下型式的一元二次方程:
其中;
("tr"表示矩陣的跡)以及
特徵值的顯式解可以由求得二次方程而得:
其中
特徵向量及節點
[編輯]特徵向量及節點會決定相路徑的輪廓,可以用圖像來描繪動態系統的解,如下所述:
要畫相平面時,會先畫對應二個特徵向量的直線(表示系統既不趨近直線,也不遠離直線的穩定條件)。之後相平面就會用有箭頭的實線代替有向量場上每點的箭頭。特徵值的正負號會影響的相平面的特點:
以上說明可以用微分方程解中指數解的行為來理解。
重覆的特徵值
[編輯]若是二個特徵值為相同的實數,會需要透過一個未知的向量以及第一個特徵向量來求解係數矩陣,產生第二個解。不過若系統對應,也可以用正交的特徵向量來產生第二個解。
有虛部的特徵值
[編輯]有虛部的複數特徵值,表示其解包括有正弦及餘弦函數(可以表示為冪次為複數的次數)。此情形比較簡單,只需要一個特徵值以及一個特徵向量就可以產生系統的完整解。
相關條目
[編輯]參考資料
[編輯]- ^ 1.0 1.1 1.2 1.3 D.W. Jordan; P. Smith. Non-Linear Ordinary Differential Equations: Introduction for Scientists and Engineers 4th. Oxford University Press. 2007. ISBN 978-0-19-920825-8.
- ^ K.T. Alligood; T.D. Sauer; J.A. Yorke. Chaos: An Introduction to Dynamical Systems. Springer. 1996. ISBN 978-0-38794-677-1.
- ^ W.E. Boyce; R.C. Diprima. Elementary Differential Equations and Boundary Value Problems 4th. John Wiley & Sons. 1986. ISBN 0-471-83824-1.