伪谱：看见脆弱性，最终的收获

当特征值误导你时

鲍尔-法伊克已警告我们：非正规矩阵的特征值在微小扰动下可能剧烈移动。因此 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. 用盖尔什戈林定位特征值，用外尔或鲍尔-法伊克界定其移动，用伪谱可视化脆弱性。