Skip to content

矩阵范数与诱导算子范数

Matrix norm · Induced matrix norm · Operator norm of a matrix

用诱导范数和常用可计算矩阵范数度量线性映射的放大能力,并区分算子范数、Frobenius 范数与谱半径。

形式陈述

矩阵空间上的范数首先满足正定性、绝对齐次性与三角不等式。若还对可乘矩阵满足

ABAB,

则称它次乘;若给定输入、输出向量范数后满足 AxYAxX,则称它与这组向量范数相容。选择向量范数后,由矩阵代表的线性映射自然得到诱导范数

AXY=supx0AxYxX=supxX=1AxY.

它正是抽象算子范数在有限维坐标中的形式,表示单位输入可能遭受的最大放大。由同一向量范数在定义域和值域诱导的矩阵范数自动相容且次乘。

A=(aij)Fm×n,最常用的三种诱导范数为

A1=max1jni=1m|aij|,A=max1imj=1n|aij|,

以及

A2=σmax(A).

前两式分别是最大绝对列和与最大绝对行和,可以直接计算;2-范数由最大奇异值给出。Frobenius 范数

AF=(i,j|aij|2)1/2=(iσi2)1/2

也次乘,并满足 Ax2AFx2,但在维数大于一时不是某个同空间向量范数的诱导矩阵范数:任何这样的诱导范数都有 I=1,而 InF=n

直觉

向量范数规定输入和输出各自怎样量尺寸,诱导矩阵范数就在这套尺子下寻找最容易被放大的方向。换一种向量范数,同一矩阵的“最大增益”也会改变;因此没有脱离尺度选择的唯一矩阵大小。

Frobenius 范数采用另一幅图像:它把矩阵的全部元素,等价地全部奇异值,当作一个长向量计算总能量。它容易求、对酉变换不变,却不直接回答某个单位输入最多会被放大多少。把两类范数都统称为 operator norm,会抹掉这个区别。

例子与边界

A=(1230).

最大列和与最大行和分别给出

A1=4,A=3,

AF=14。由于

AA=(10224)

的最大特征值为 7+13,所以 A2=7+133.2566。这些数都正确,只是在回答不同的尺度问题,不能省略下标后直接比较。

对方阵,谱半径 ρ(A)=maxi|λi| 不超过任意诱导范数,但通常不等于它。非正规矩阵

BM=(1M01)

的谱半径恒为 1,而 BM2|M| 增大。这说明特征值只描述不变方向上的渐近伸缩,矩阵范数还会看见非正交方向之间的瞬时放大。

推论与应用

矩阵范数为线性系统条件数、残差到前向误差的转换、分解的后向误差与迭代收敛界提供共同尺度。次乘性允许把一串算子作用的误差逐层控制,而诱导范数的相容性保证矩阵界确实能约束向量误差。

奇异值分解解释 2-范数和 Frobenius 范数的谱表达;选择 1-范数或无穷范数时,则可使用矩阵元素直接计算或估计。工程报告必须写明使用的范数,因为“条件数大”或“扰动小”在没有下标时并不是完整结论。

参考资料
  • 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.