矩陣指數計算機

矩陣指數計算機 用於計算 \(e^{At}\)，即齊次線性系統 \(x'(t)=Ax(t)\) 的狀態轉移矩陣。它專為線性代數、控制理論、微分方程、馬可夫鏈生成器以及任何由常數矩陣驅動連續時間演化的模型而設計。

矩陣指數的意義

對於純量 \(a\)，指數 \(e^{at}\) 求解 \(x'=ax\)。對於方陣 \(A\)，將數值的冪替換為矩陣的冪後，同樣的概念依然適用：

結果並非通過對 \(A\) 的每個條目求指數得到。冪運算 \(A^2,A^3,\ldots\) 中的矩陣乘法捕捉了變量之間的耦合，這正是線性 ODE 系統所需要的。

求解線性 ODE 系統

若 \(A\) 為常數且 \(x(0)=x_0\)，則初值問題的解為：

這就是為什麼 \(e^{At}\) 常被稱為 狀態轉移矩陣 或 基本矩陣解。每一列顯示了標準基底狀態在時間 \(t\) 後移動到的位置。

如何使用矩陣指數計算機

1. 輸入矩陣 A： 每行輸入一行，條目之間使用空格或逗號分隔。
2. 選擇時間 t： 正值用於正向演化，負值用於反向演化。
3. 求解 ODE 時添加 x(0)： 向量的條目數必須與矩陣維度相同。
4. 計算並檢查： 讀取 \(e^{At}\)、可選的 \(x(t)\)、跡恆等式，以及當 A 為 2×2 時的 2D 動畫。

數值方法

本計算機使用 縮放與平方（scaling and squaring） 配合 13 階 Padé 近似法。具體而言，它首先將 \(At\) 縮放為較小的矩陣，計算有理近似值，然後重複平方結果以回到原始時間尺度。這比簡單地截斷泰勒級數更穩定。

重要恆等式：體積縮放

矩陣指數的行列式有一個簡潔的跡公式：

對於 2D 系統，這描述了流動下的面積縮放；對於 3D 系統，它描述了體積縮放。負跡傾向於收縮體積，而正跡則擴張體積。

何時使用此工具

常見問題

計算機能處理不可對角化的矩陣嗎？

可以。Padé 方法直接計算 \(e^{At}\)，因此不需要對角化。只要數值保持在穩定限制內，Jordan 塊和重根特徵值都是有效的輸入。

為什麼 ||At|| 有限制？

極大的 \(\|At\|_1\) 值可能導致巨大的指數條目或浮點溢出。計算機保持一個保守的界限，以便用戶獲得可靠且瀏覽器友好的結果，而不是誤導性的無窮大。

這會產生符號公式嗎？

此工具專注於數值矩陣指數和 ODE 狀態值。對於符號閉合形式、對角化和 Jordan 標準型工作流程，請使用專用的特徵值或 Jordan 標準型計算機。