1. 程式人生 > >線性代數之——微分方程和 exp(At)

線性代數之——微分方程和 exp(At)

本節的核心是將常係數微分方程轉化為線性代數問題。

d u d t =

λ u u ( t )
= C e λ t \frac{du}{dt}=\lambda u \quad 的解為 \quad u(t) = Ce^{\lambda t}

代入 t = 0 t=0 ,可得 u ( 0 ) = C u(0) = C ,因此有 u ( t ) = u ( 0 ) e λ t u(t) = u(0)e^{\lambda t} 。這是隻有一個變數的情況,線上性代數裡,我們擴充套件到 n n 個方程的情況。

d u d t = A u u ( 0 ) t = 0 \frac{d\boldsymbol u}{dt}=A \boldsymbol u \quad 初始條件為向量 \quad \boldsymbol u(0)_{t=0}

注意,這裡 A A 是常矩陣,不隨時間而改變。而且這些方程是線性的,如果 u ( t ) \boldsymbol u(t) v ( t ) \boldsymbol v(t) 都是方程組的解,那麼它們的線性組合 C u ( t ) + D v ( t ) C\boldsymbol u(t)+D\boldsymbol v(t) 也是解,我們需要 n n 個這樣的常數來匹配方程組的初始條件。

1. d u d t = A u \frac{d\boldsymbol u}{dt}=A \boldsymbol u 的解

其中一個解是 e λ t x e^{\lambda t} \boldsymbol x λ \lambda 是矩陣 A A 的特徵值,而 x \boldsymbol x 是特徵向量。將這個解代入原方程,利用 A x = λ x A\boldsymbol x=\lambda \boldsymbol x 可得

d u d t = λ e λ t x = A e λ t x = A u \frac{d\boldsymbol u}{dt} = \lambda e^{\lambda t} \boldsymbol x = A e^{\lambda t} \boldsymbol x=A \boldsymbol u

這個解的所有部分都有 e λ t e^{\lambda t} ,當 λ > 0 \lambda>0 時,解會增長;當 λ < 0 \lambda<0 時,解會衰減。而當 λ \lambda 為虛數時,則它的實部決定解是增長還是衰減。

  • 例 1

求解 d u d t = A u = [ 0 1 1 0 ] u u 0 = [ 4 2 ] \frac{d\boldsymbol u}{dt}=A \boldsymbol u = \begin{bmatrix}0&1 \\ 1&0\end{bmatrix}\boldsymbol u,\boldsymbol u_0 = \begin{bmatrix}4 \\ 2\end{bmatrix}

矩陣 A A 的特徵值為 1 和 -1,特徵向量為 (1, 1) 和 (1, -1),因此兩個純指數解為:

u 1 ( t ) = e λ 1 t x 1 = e t [ 1 1 ] \boldsymbol u_1(t) = e^{\lambda_1 t} \boldsymbol x_1 = e^t\begin{bmatrix}1 \\ 1\end{bmatrix}

u 2 ( t ) = e λ 2 t x 2 = e t [ 1 1 ] \boldsymbol u_2(t) = e^{\lambda_2 t} \boldsymbol x_2 = e^{-t}\begin{bmatrix}1 \\ -1\end{bmatrix}

這些 u \boldsymbol u 依然是矩陣的特徵向量,它們滿足 A u 1 = u 1 A\boldsymbol u_1 = \boldsymbol u_1 A u 2 = u 2 A\boldsymbol u_2 = -\boldsymbol u_2 ,只不過是係數隨著 t t 改變罷了。方程組的全解為這些特解的線性組合。

利用初始條件我們可以確定出係數 C C D D

因此,我們可以通過以下三個步驟來求解 d u d t = A u \frac{d\boldsymbol u}{dt}=A \boldsymbol u

  • u 0 \boldsymbol u_0 寫成特徵向量的線性組合, u 0 = c 1 x 1 +