偽譜：看見脆弱性，最終的收穫

當特徵值誤導你時

鮑爾-法伊克已警告我們：非正規矩陣的特徵值在微小擾動下可能劇烈移動。因此 rho(A) < 1 不再保證 ||A^k|| 對中等的 k 保持小——在譜半徑公式承諾的漸近衰減開始前，瞬態增長可能極大。僅憑譜是錯誤的診斷工具。

修正方法是用偽譜代替譜（即 A - z*I 恰好奇異的集合）：A - z*I 幾乎奇異的 z 的集合。對 epsilon > 0，epsilon-偽譜是 { z : ||(A - z*I)^-1|| > 1/epsilon }，等價地 { z : sigma_min(A - z*I) < epsilon }。

把豫解式範數當作地圖

量 ||(A - z*I)^-1|| 是豫解式範數。它恰在特徵值處飆升至無窮。對正規矩陣，飆升很尖銳，epsilon-偽譜只是每個特徵值周圍的 epsilon-圓盤。對非正規矩陣，豫解式在遠離譜處仍很大，於是偽譜向外膨脹——一幅潛在不穩定的圖像。

Probe the resolvent on a grid of complex points z:

  for each z in grid:
      M = A - z*I
      s = sigma_min(M)          # smallest singular value (SVD)
      resolvent_norm = 1 / s    # = ||(A - zI)^-1||_2
      z is in eps-pseudospectrum  <=>  s < eps

Normal A:  s small ONLY near true eigenvalues  -> tight discs
Non-normal A:  s small over a WIDE region       -> bulging set

Key identity:  ||(A - zI)^-1||_2 = 1 / sigma_min(A - zI)

收穫：你能信賴的秩與穩定性

同樣的 sigma_min 思想修正了另一個脆弱概念：秩。精確秩計數非零奇異值，但在浮點中沒有東西恰好為零。數值秩計數高於某容差的奇異值，因此它是該容差內最近矩陣的秩——用擾動思維使秩變得穩健。

1. 用合適的範數度量大小，並利用等價性自由切換。
2. 用次可乘性與諾伊曼級數來擾動逆矩陣並鏈式傳遞誤差。
3. 在動手計算前，讀取 kappa(A) 以預測相對誤差會被放大得多嚴重。
4. 用蓋爾什戈林定位特徵值，用外爾或鮑爾-法伊克界定其移動，用偽譜可視化脆弱性。