正交矩阵
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
正交矩阵
正交矩阵是实数特殊化的酉矩阵,因此总是正规矩阵。
尽管我们在这里只考虑实数矩阵,这个定义可用于其元素来自任何域的矩阵。
正交矩阵毕竟是从内积自然引出的,对于复数的矩阵这导致了归一要求。
目录
定义 1
n阶实矩阵 A称为正交矩阵,如果:A×A′=E(E为单位矩阵,A'表示“矩阵A的转置矩阵”。
)若A为正交阵,则下列诸条件是等价的:
1) A 是正交矩阵
2) A×A′=E(E为单位矩阵)
3) A′是正交矩阵
4) A的各行是单位向量且两两正交
5) A的各列是单位向量且两两正交
6) (Ax,Ay)=(x,y) x,y∈R
正交矩阵通常用字母Q表示。
举例:A=[r11 r12 r13;r21 r22 r23;r31 r32 r33]
下面是一些小正交矩阵的例子和可能的解释。
恒等变换。
旋转16.26°。
针对x轴反射。
旋转反演(rotoinversion): 轴 (0,-3/5,4/5),角度90°。
置换坐标轴。
编辑本段
基本构造
低维度
最简单的正交矩阵是1×1 矩阵 [1] 和 [−1],它们可分别解释为恒等和实数线针对原点的反射。
如下形式的2×2 矩阵
它的正交性要求满足三个方程
矩阵性质
实数方块矩阵是正交的,当且仅当它的列形成了带有普通欧几里得点积的欧几里得空间R的正交规范基,它为真当且仅当它的行形成R的正交基。
假设带有正交(非正交规范)列的矩阵叫正交矩阵可能是诱人的,但是这种矩阵没有特殊价值而没有特殊名字;他们只是MM = D,D是对角矩阵。
任何正交矩阵的行列式是 +1 或−1。
这可从关于行列式的如下基本事实得出:
反过来不是真的;有 +1 行列式不保证正交性,即使带有正交列,可由下列反例证实。
对于置换矩阵,行列式是 +1 还是−1 匹配置换是偶还是奇的标志,行列式是行的交替函数。
比行列式限制更强的是正交矩阵总可以是在复数上可对角化来展示特征值的完全的集合,它们全都必须有(复数)绝对值 1。
群性质
正交矩阵的逆是正交的,两个正交矩阵的积是正交的。
事实上,所有
n×n正交矩阵的集合满足群的所有公理。
它是n(n−1)/2 维的紧致李群,叫做正交群并指示为O(n)。
行列式为 +1 的正交矩阵形成了路径连通的子群指标为 2 的O(n)
正规子群,叫做旋转的特殊正交群SO(n)。
商群O(n)/SO(n) 同构于O(1),带有依据行列式选择 [+1] 或 [−1] 的投影映射。
带有行列式−1 的正交矩阵不包括单位矩阵,所以不形成子群而只是陪集;它也是(分离的)连通的。
所以每个正交群被分为两个部分;因为投影映射分裂,O(n) 是SO(n) 与O(1)的半直积。
用实用术语说,一个相当的陈述是任何正交矩阵可以通过采用一个旋转矩阵并可能取负它的一列来生成,如我们在2×2 矩阵中看到的。
如果n是奇数,则半直积实际上是直积,任何正交矩阵可以通过采用一个旋转矩阵并可能取负它的所有列来生成。
现在考虑 (n+1)×(n+1) 右底元素等于 1 的正交矩阵。
最后一列(和最后一行)的余下元素必须是零,而任何两个这种矩阵的积有同样的形式。
余下的矩阵是n×n正交矩阵;因此O(n) 是O(n+1) (和所有更高维群)的子群。
因为 Householder 正交矩阵形式的基本反射可把任何正交矩阵简约
成这种约束形式,一系列的这种反射可以把任何正交矩阵变回单位矩阵;
因此正交群是反射群。
最后一列可以被固定为任何单位向量,并且每种选择给出不同的O(n) 在O(n+1) 中的复本;以这种方式O(n+1) 是在单位球S与纤维O(n) 上的丛。
类似的,SO(n) 是SO(n+1) 的子群;任何特定正交矩阵可以使用类似过程通过 Givens 平面旋转来生成。
丛结构持续: SO(n) ↪SO(n+1) → S。
一个单一旋转可以在最后一列的第一行生成一个零,而n−1 次旋转序列将置零n×n旋转矩阵的除了最后一列的最后一行的所有元素。
因为平面是固定的,每次旋转只有一个自由度,就是它的角度。
通过归纳,SO(n) 因此有
自由度,O(n) 也是。
置换矩阵简单一些;它们不形成李群,只是一个有限群,n! 次对称群Sn。
通过同类的讨论,Sn是Sn+1 的子群。
偶置换生成行列式 +1 的置换矩阵的子群,n!/2 次交错群。
规范形式
更广泛的说,任何正交矩阵的效果分离到在正交二维空间上的独立动作。
就是说,如果Q是狭义正交的,则你可以找到(旋转)改变基的一个正交矩阵P,把Q带回到分块对角形式:
(n偶数), (n奇数)。
这里的矩阵R1,...,Rk是2×2 旋转矩阵,而余下的元素是零。
作为例外,一个旋转块可以是对角的,±I。
因此如果需要的话取负一列,并注意2×2 反射可对角化为 +1 和−1,任何正交矩阵可变为如下形式
, 矩阵R1,…,Rk给出位于复平面中单位圆上的特征值的共轭对;所以这个分解复合确定所有带有绝对值 1 的特征值。
如果n是奇数,至少有一个实数特征值 +1 或−1;对于3×3 旋转,关联着 +1 的特征向量是旋转轴。
编辑本段
数值线性代数
利益
数值分析自然的利用了正交矩阵的很多数值线性代数的性质。
例如,经常需要计算空间的正交基,或基的正交变更;二者都采用了正交矩阵的形式。
有行列式±1 和所有模为 1 的特征值是对数值稳定性非常有利的。
一个蕴涵是条件数为 1 (这是极小的),所以在乘以正交矩阵的时候错误不放大。
很多算法为此使用正交矩阵如 Householder反射和 Givens旋转。
有帮助的不只是正交矩阵是可逆的,还有它的逆矩阵本质上是免花费的,只需要对换索引(下标)。
置换是很多算法成功的根本,包括有局部定支点(partial pivoting)的运算繁重的高斯消去法(这里的置换用来定支点)。
但是它们很少明显作为矩阵出现;它们的特殊形式允许更有限的表示,比如n个索引的列表。
同样的,使用 Householder 和 Givens 矩阵的算法典型的使用特殊方法的乘法和存储。
例如,Givens 旋转只影响它所乘的矩阵的两行,替代完全的n次的矩阵乘法为更有效的n次运算。
在使用这些反射和旋转向矩阵介入零的时候,腾出的空间足够存储充足的数据来重生成这个变换。
分解
一些重要的矩阵分解(Golub & Van Loan, 1996)涉及到了正交矩阵,包括:
QR分解M = QR, Q正交,R上三角。
奇异值分解M = UΣV, U和V 正交,Σ非负对角。
谱分解S = QΛQ, S对称,Q正交,Λ对角。
极分解M = QS, Q正交,S对称非负确定。