Skip to content

正交矩阵与酉矩阵

Orthogonal matrix · Unitary matrix · 酉矩阵

在实或复内积空间中保持内积的方阵,其逆分别等于转置或共轭转置。

形式陈述

实方阵 QRn×n 称为正交矩阵,若

QTQ=In.

复方阵 UCn×n 称为酉矩阵,若

UU=In,

其中 U=UT 是共轭转置。因为矩阵是方阵,左逆同时是右逆,所以

Q1=QT,U1=U.

这些条件等价于列向量构成标准内积下的标准正交基;也等价于行向量构成标准正交基。对任意向量 x,y

Ux,Uy=xUUy=xy=x,y.

因而酉矩阵保持内积、Euclidean 范数、距离与夹角。实正交矩阵是酉矩阵在实数域上的对应物。

正交矩阵群记作 O(n)。由行列式恒等式

det(QTQ)=det(Q)2=1

可知 detQ=±1;其中行列式为 1 的子群 SO(n) 表示保定向的旋转,行列式为 1 的元素包含反射。酉矩阵的特征值都位于单位圆上,因为若 Uv=λvv0,则 v=Uv=|λ|v

直觉

正交或酉变换可以旋转、反射以及在复空间中改变相位,但不会伸长任何方向。一般矩阵会把单位球变成椭球;正交矩阵把单位球原样搬到自身。逆变换无需重新求解线性方程,只要转置或共轭转置,就能把这次刚性运动撤销。

“列正交”意味着不同输入坐标被送到互相垂直且长度为一的输出方向。方阵条件保证这些方向填满整个空间;少一列时仍可保持定义域中的长度,却不是整个陪域上的可逆刚性运动。

例子与计算

二维旋转矩阵

Qθ=(cosθsinθsinθcosθ)

满足 QθTQθ=IdetQθ=1。取 θ=π/2,向量 (2,1)T 被送到 (1,2)T,两者范数都为 5,而转置矩阵正好旋转 π/2 将它送回。

Householder 反射

H=I2vvTvTv,v0,

满足 HT=HH2=I,所以正交且 detH=1。它把 v 方向反向,保持 v 中的每个向量不变。

复数例子 U=diag(1,i) 满足 UU=I,但 UTU=diag(1,1)I。这说明复空间必须使用共轭转置;只做普通转置会丢失 Hermitian 内积中的共轭。

边界与失败情形

ARm×nm>n,条件 ATA=In 只说明其列正交,A 是从 Rn 到列空间的等距嵌入。此时 AAT 是投影到列空间的矩阵,通常不等于 Im,也不存在双侧逆。因此不能把长方矩阵称为正交矩阵而沿用 A1=AT

近似满足 Q^TQ^I 不等于定义上正交。浮点 QR 算法产生的矩阵会有正交性残差

IQ^TQ^,

其大小取决于算法与舍入误差。数值稳定性是实现性质,不包含在代数定义中;经典 Gram–Schmidt 与 Householder QR 即使目标分解相同,有限精度表现也不同。

保持行列式绝对值为一并不足以正交。例如 A=diag(2,1/2) 的行列式为一,却把坐标方向分别拉伸两倍和压缩一半,故 ATAI。保持某一个向量的范数也不足以推出保持所有内积。

推论与应用

正交相似变换 AQTAQ 保持特征多项式,并且在 Euclidean 范数下不放大向量误差。实对称矩阵可被正交对角化,复正规矩阵可被酉对角化;这里的正交/酉基使谱分量之间不发生几何畸变。

Givens 旋转一次消去选定的一个矩阵元素,Householder 反射一次处理整段向量,两者都是结构化正交变换。QR 分解通过它们把矩阵化为上三角形,同时尽量避免条件数因变换本身恶化。离散 Fourier 变换在适当归一化后则是酉矩阵。

参考资料
  • Roger A. Horn and Charles R. Johnson, Matrix Analysis, 2nd ed., Cambridge University Press, 2013,§2.1–2.5。
  • Gene H. Golub and Charles F. Van Loan, Matrix Computations, 4th ed., Johns Hopkins University Press, 2013,§5.1–5.2。
  • Sheldon Axler, Linear Algebra Done Right, 3rd ed., Springer, 2015,Ch. 7。