6.3 複特徵值與實應用#
前言：看見看不見的振盪#
對矩陣 \(\begin{bmatrix} 0 & -1 \\ 1 & 0 \end{bmatrix}\) （旋轉 90°）解特徵方程： \(\lambda^2 + 1 = 0\) ， \(\lambda = \pm i\) ——純虛數！實矩陣竟有複特徵值？是的，而且這不是麻煩，是禮物：複特徵值攜帶振盪資訊。實部控制增長/衰減，虛部控制振盪頻率。彈簧的簡諧運動、交流電路、螺旋星系、量子相位，都是複特徵值的實體化身。本節建立複特徵值理論，並展示它如何在純實數世界中工作。
歷史前因後果#
虛數的合法化之路#
- 1545 年，卡爾達諾（Gerolamo Cardano, 1501–1576）在《大術》（Ars Magna）解三次方程時首次被迫操作 \(\sqrt{-15}\) ，稱之為「心理折磨」（mental tortures）。1572 年邦貝利（Rafael Bombelli, 1526–1572）建立虛數的運算規則，讓三次方程的「 casus irreducibilis」（不可約情形）得以計算。
- 1748 年，歐拉（Leonhard Euler, 1707–1783）發表 \(e^{i\theta} = \cos\theta + i\sin\theta\) ——複指數與振盪的橋梁。十九世紀柯西與高斯完成複分析基礎，虛數終於「合法」。
- 1858 年，凱萊的矩陣論允許複特徵值；柯西（1829）與厄米特（Charles Hermite, 1855）證明：對稱（自伴）矩陣特徵值必為實數——複特徵值只出現在非對稱矩陣。
- 1892 年，Poincaré 在天體力學中系統使用「特徵指數」（複數）分類軌道穩定性；Liapunov（1892《運動穩定性通論》）以特徵值實部建立穩定性理論——現代控制理論的基礎。
-
量子力學（1925–26）：狀態空間是複希爾伯特空間，相位 \(e^{i\theta}\) 是干涉現象的載體——複數從「方便的技巧」變成自然的必需品。
-
解決了什麼問題：振盪系統（機械、電路、軌道）的統一描述；實係數多項式虛根成對出現的解釋。
- 引發了什麼新問題：複矩陣的「實化身」——實部虛部如何安排（實 Jordan 形）；複矩陣的正定性與酉對角化（第 8 章前置）。
複特徵值的基本理論#
為什麼實矩陣會有複特徵值？#
特徵多項式是實係數多項式，但實係數多項式的根可以是複數。共軛成對定理：若 \(\lambda = a + bi\) 是 \(A\) 的特徵值（ \(A\) 實），則 \(\bar\lambda = a - bi\) 也是，且特徵向量也共軛： \(\mathbf{v}\) 對應 \(\lambda\) ，則 \(\bar{\mathbf{v}}\) 對應 \(\bar\lambda\) 。
證明： \(A\mathbf{v} = \lambda\mathbf{v}\) ，兩邊取共軛（ \(A\) 實， \(\bar A = A\) ）： \(A\bar{\mathbf{v}} = \bar\lambda\bar{\mathbf{v}}\) 。∎
幾何重數也成對，因此複特徵值不影響可對角化判準（對 \(n\) 計數時成對出現）。
複特徵值的幾何：旋轉–縮放#
把 \(\lambda = a + bi\) 寫成極式：
\(A^k\) 的作用中， \(\lambda^k = r^ke^{ik\theta}\) ：半徑按 \(r^k\) 縮放、角度按 \(k\theta\) 旋轉。
| 特徵值 | 長期行為 |
|---|---|
| $ | \lambda |
| $ | \lambda |
| $ | \lambda |
| \(\lambda\) 實 | 沿直線縮放（6.1–6.2 情形） |
| \(\lambda\) 純虛 | 純振盪、不衰減（簡諧） |
實部與虛部的分工（微分方程）#
\(\mathbf{x}' = A\mathbf{x}\) ， \(\lambda = \alpha \pm i\beta\) 的解：
- \(\alpha = \mathrm{Re}\,\lambda\) ：增長/衰減率（ \(\alpha < 0\) 阻尼、 \(\alpha = 0\) 無阻尼、 \(\alpha > 0\) 自激）；
- \(\beta = \mathrm{Im}\,\lambda\) ：振盪角頻率（週期 \(2\pi/\beta\) ）。
穩定性判準（Liapunov）：連續系統穩定 ⇔ 所有 \(\mathrm{Re}(\lambda) < 0\) ；離散系統穩定 ⇔ 所有 \(|\lambda| < 1\) 。
二維實空間中的複特徵對#
\(2\times2\) 實矩陣 \(A\) 有複特徵值 \(\alpha \pm i\beta\) （ \(\beta \ne 0\) ）時，取 \(\mathbf{v} = \mathbf{p} + i\mathbf{q}\) （ \(\mathbf{p}, \mathbf{q}\) 實），則實解為：
幾何：軌跡是螺旋線——在由 \(\mathbf{p}, \mathbf{q}\) 張成的平面上旋轉（角速 \(\beta\) ），同時按 \(e^{\alpha t}\) 膨脹/收縮。
例題#
例題 6.3-1（純旋轉）： \(R = \begin{bmatrix} \cos\theta & -\sin\theta \\ \sin\theta & \cos\theta \end{bmatrix}\) 。
解：特徵方程 \((\cos\theta - \lambda)^2 + \sin^2\theta = 0\) ， \(\lambda = \cos\theta \pm i\sin\theta = e^{\pm i\theta}\) 。 \(|\lambda| = 1\) ——純旋轉，無增長。特徵向量（複）：對 \(\lambda = e^{i\theta}\) ， \((R - \lambda I)\mathbf{v} = \mathbf{0}\) 給 \(\mathbf{v} = (1, -i)\) （驗證： \(R(1, -i) = (\cos\theta + i\sin\theta,\ \sin\theta - i\cos\theta) = e^{i\theta}(1, -i)\) ✓）。
解讀：旋轉沒有實的「不動方向」（除非 \(\theta = 0, \pi\) ）——這解釋了為什麼直觀找不到不動箭：不動方向「住在複數裡」。 \(\theta = 90°\) 時 \(\lambda = \pm i\) ， \(R^4 = I\) （轉四圈回原點）： \(i^4 = 1\) ✓——複數運算完美預測矩陣行為。
例題 6.3-2（阻尼振盪）：解 \(\mathbf{x}' = A\mathbf{x}\) ， \(A = \begin{bmatrix} 0 & 1 \\ -5 & -2 \end{bmatrix}\) 。
解：特徵方程 \(\lambda^2 + 2\lambda + 5 = 0\) ， \(\lambda = \frac{-2 \pm \sqrt{4 - 20}}{2} = -1 \pm 2i\) 。
對 \(\lambda = -1 + 2i\) ： \((A - \lambda I)\mathbf{v} = \begin{bmatrix} 1-2i & 1 \\ -5 & 1-2i \end{bmatrix}\mathbf{v} = \mathbf{0}\) 。取 \(\mathbf{v} = (1, -1+2i)\) ，即 \(\mathbf{p} = (1, -1)\) 、 \(\mathbf{q} = (0, 2)\) 。
實解（ \(\alpha = -1, \beta = 2\) ）：
判讀： \(\alpha = -1 < 0\) 衰減； \(\beta = 2\) 振盪角頻率（週期 \(\pi\) ）。阻尼振盪——像盪鞦韆有空氣阻力。軌跡是收縮螺旋。這是二階方程 \(x'' + 2x' + 5x = 0\) （欠阻尼： \(4 < 20\) ）的狀態空間版本。
例題 6.3-3（離散系統：螺旋外擴）： \(A = \begin{bmatrix} 0.5 & -0.8 \\ 0.8 & 0.5 \end{bmatrix}\) ，討論 \(A^k\) 的長期行為。
解： \(A = rR_\phi\) 形式（ \(0.5^2 + 0.8^2 = 0.89 < 1\) ）： \(r = \sqrt{0.89} \approx 0.9434\) ， \(\phi = \arctan(0.8/0.5) \approx 58°\) 。特徵值 \(\lambda = 0.5 \pm 0.8i\) ， \(|\lambda| = 0.9434 < 1\) 。
行為：每次迭代縮到 \(94.34\%\) 、旋轉 \(58°\) 。長期收縮螺旋趨於零——穩定。若把 \(0.5\) 換成 \(1.05\) ： \(r = \sqrt{1.1025 + 0.64} > 1\) ，發散螺旋——不穩定。判準一目了然： \(|\lambda|\) 與 1 比大小。
例題 6.3-4（PageRank 之外的網路：人口振盪）： LESLIE 矩陣（年齡結構人口模型，Leslie 1945）：
特徵方程： \(-\lambda^3 + 2\lambda + 0.375 = 0\) ……（展開： \(\det(L - \lambda I) = -\lambda^3 + 4\cdot0.5\cdot\lambda + 3\cdot0.5\cdot0.25 = -\lambda^3 + 2\lambda + 0.375\) ）。解 \(\lambda^3 - 2\lambda - 0.375 = 0\) ：試 \(\lambda \approx 1.5\) ： \(3.375 - 3 - 0.375 = 0\) ✓。故 \(\lambda_1 = 1.5\) ，其餘兩根由 \(\lambda^2 + 1.5\lambda + 0.25 = 0\) ： \(\lambda = \frac{-1.5 \pm \sqrt{2.25 - 1}}{2} = \frac{-1.5 \pm \sqrt{1.25}}{2} \approx -0.192, -1.059\) 。
判讀：主特徵值 \(\lambda_1 = 1.5 > 1\) 且唯一最大模——人口每年成長 50%（Perron–Frobenius 定理保證其他模式相對消失）。年齡分布收斂到 \(\lambda_1\) 的特徵向量比例。若主特徵值是複數，人口將出現世代振盪——生育高峰的循環現象在數學上就是複主特徵值。
什麼矩陣永遠沒有複特徵值？#
| 矩陣類型 | 特徵值性質 | 原因 |
|---|---|---|
| 對稱（ \(A = A^T\) ） | 全實 | 譜定理（第 8 章），柯西 1829 |
| 上三角 | 對角線元素 | 特徵多項式 = \(\prod(t - a_{ii})\) |
| 旋轉（正交， \(\det = 1\) ） | \(1, e^{\pm i\theta}\) | 保長度的推論 |
| 正矩陣 | 有最大正實特徵值 | Perron–Frobenius（1907/1912） |
| 隨機矩陣（列和 1） | \(\lambda = 1\) 必為特徵值 | \(\mathbf{1}\) 是 \(A^T\) 的特徵向量 |
柯西的直覺版：對稱矩陣對二次型 \(\mathbf{x}^TA\mathbf{x}\) 恆為實值；若 \(A\mathbf{v} = \lambda\mathbf{v}\) （複 \(\mathbf{v}\) ），取共軛內積可得 \(\lambda\|\mathbf{v}\|^2 = \bar{\mathbf{v}}^TA\mathbf{v}\) 為實數，故 \(\lambda \in \mathbb{R}\) （完整證明見 8.1 節）。
與其他領域的關聯#
- 微分方程：二階方程 \(x'' + c x' + kx = 0\) 的三種阻尼（過阻尼、臨界、欠阻尼）對應特徵根實、重根、複——特徵值是「振動物理學」的通用語言。
- 電路學：RLC 電路的自然響應由複特徵值決定；阻抗 \(Z = R + i\omega L\) 是複數的物理化身。
- 控制工程：極點（特徵值）在複平面的位置決定響應品質；控制設計 = 移動極點。Nyquist 圖、根軌跡法都是複特徵值的圖像化。
- 量子力學：演化 \(e^{-i\hat{H}t/\hbar}\) ——虛部在指數上，能量實特徵值給出相位振盪 \(e^{-iEt/\hbar}\) ；干涉、隧穿都靠複相位。
- 訊號處理：DFT 的旋轉因子 \(e^{2\pi i k/N}\) ；Z-變換的極點（ \(|\lambda|\) 與 1 的關係）決定濾波器穩定性。
- 混沌與分形：Mandelbrot 集合的定義 \(z \to z^2 + c\) 的逃逸判準，本質是「複動力系統的特徵值模」迭代。
本節小結#
- 實矩陣的複特徵值共軛成對（ \(\lambda, \bar\lambda\) ），特徵向量亦然。
- \(\lambda = re^{i\theta}\) ：模 \(r\) 控制增長/衰減、幅角 \(\theta\) 控制旋轉；微分方程中 \(\mathrm{Re}\,\lambda\) 是阻尼、 \(\mathrm{Im}\,\lambda\) 是頻率。
- 複特徵對在實空間的化身是螺旋軌跡：旋轉平面（由 \(\mathbf{p}, \mathbf{q}\) 張成）+ 徑向縮放。
- 穩定性判準：連續 \(\mathrm{Re}\,\lambda < 0\) 、離散 \(|\lambda| < 1\) ；主特徵值決定長期宿命。
- 對稱矩陣特徵值必實（柯西 1829）——第 8 章譜定理的預告。
預告：既然對稱矩陣特徵值全實、且行為如此良好，它值專章對待。第 7 章先建立內積空間的幾何工具（長度、角度、正交），第 8 章用它们征服對稱矩陣。