“因此从任意初值收敛到 $x $ 当且仅当迭代矩阵的谱半径满足 $\rho(G)<1$。谱半径不等于任意指定的矩阵范数;若某个相容范数有 $ G <1$,它是方便的充分条件,而不是必要条件。”
形式陈述 ​
矩阵空间上的范数首先满足正定性、绝对齐次性与三角不等式。若还对可乘矩阵满足
则称它次乘;若给定输入、输出向量范数后满足
它正是抽象算子范数在有限维坐标中的形式,表示单位输入可能遭受的最大放大。由同一向量范数在定义域和值域诱导的矩阵范数自动相容且次乘。
对
以及
前两式分别是最大绝对列和与最大绝对行和,可以直接计算;
也次乘,并满足
直觉 ​
向量范数规定输入和输出各自怎样量尺寸,诱导矩阵范数就在这套尺子下寻找最容易被放大的方向。换一种向量范数,同一矩阵的“最大增益”也会改变;因此没有脱离尺度选择的唯一矩阵大小。
Frobenius 范数采用另一幅图像:它把矩阵的全部元素,等价地全部奇异值,当作一个长向量计算总能量。它容易求、对酉变换不变,却不直接回答某个单位输入最多会被放大多少。把两类范数都统称为 operator norm,会抹掉这个区别。
例子与边界 ​
取
最大列和与最大行和分别给出
而
的最大特征值为
对方阵,谱半径
的谱半径恒为
推论与应用 ​
矩阵范数为线性系统条件数、残差到前向误差的转换、分解的后向误差与迭代收敛界提供共同尺度。次乘性允许把一串算子作用的误差逐层控制,而诱导范数的相容性保证矩阵界确实能约束向量误差。
奇异值分解解释
参考资料
- Lloyd N. Trefethen and David Bau III, Numerical Linear Algebra, SIAM, 1997, Lectures 3 and 4.
- Nicholas J. Higham, Accuracy and Stability of Numerical Algorithms, 2nd ed., SIAM, 2002, Ch. 6.
- LAPACK Users’ Guide, Basic Linear Algebra Subprograms.