矩阵范数与扰动

诺伊曼级数

诺伊曼级数是几何级数 1/(1 - r) = 1 + r + r^2 + ... 的矩阵版本。它把 I - A 的逆表示为幂的无穷和:(I - A)^-1 = I + A + A^2 + A^3 + ...。当级数收敛时,它既证明逆存在,又给出逼近它的构造性方法。

收敛有一个干净的充分条件:若在任一次乘范数下 ||A|| < 1,则部分和构成柯西序列(因为 ||A^k|| <= ||A||^k 几何式衰减),于是级数收敛,其和为 (I - A)^-1。精确条件是谱的:级数收敛当且仅当 rho(A) < 1,即谱半径小于一。

实用之处在于逆的扰动。若你已知 B^-1 而把 B 略微扰动为 B + E,写 B + E = B(I + B^-1 E) 并展开:(B + E)^-1 = (I + B^-1 E)^-1 B^-1 = (I - B^-1 E + ...) B^-1,当 ||B^-1 E|| < 1 时有效。这就是你估计线性系统的解在矩阵小变化下如何移动的方法。

它也直接界定求逆的条件。从级数得到 ||(I - A)^-1|| <= 1 / (1 - ||A||),它量化了 I - A 在其逆爆炸之前能漂移到多接近奇异——这是条件数估计与迭代求解器收敛背后的解析引擎。

(I - A)^-1 = I + A + A^2 + ..., ||(I - A)^-1|| <= 1/(1 - ||A||)

当 ||A|| < 1 时收敛;该界表明当 ||A|| 趋近 1 时逆爆炸(I - A 趋于奇异)。

范数条件 ||A|| < 1 是充分而非必要的;谱条件 rho(A) < 1 才是精确的。非正规的 A 可以有 rho(A) < 1 却 ||A|| > 1,于是即便简单的范数检验失败,级数仍收敛。

又称
operator geometric seriesNeumann expansion