正交矩阵

╰半橙微兮° 2022-03-17 09:23 240阅读 0赞

分享一下我老师大神的人工智能教程。零基础！通俗易懂！风趣幽默！还带黄段子！希望你也加入到我们人工智能的队伍中来！[https://blog.csdn.net/jiangjunshow][https_blog.csdn.net_jiangjunshow]

正交[矩阵][Link 1]是实数特殊化的酉矩阵，因此总是[正规矩阵][Link 2]。尽管我们在这里只考虑实数矩阵，这个定义可用于其元素来自任何域的矩阵。正交矩阵毕竟是从内积自然引出的，对于复数的矩阵这导致了归一要求。

## 定义 ##

定义 1

如果：AA'=E（E为单位矩阵，A'表示“矩阵A的转置矩阵”。）或A′A=E，则n阶实矩阵 A称为正交矩阵， 若A为正交阵，则满足以下条件:

1) A 是正交矩阵

2) AA′=E（E为单位矩阵）（\#add它的转置矩阵是它的逆矩阵，这是很重要的）

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\]

则有：

r11^2+r12^2+r13^2=r21^2+r22^2+r23^2=r31^2+r32^2+r33^2=1

r11\*r12+r21\*r22+r31\*r32=0等性质

正交方阵是欧氏空间中标准正交基到标准正交基的过渡矩阵。

在 [矩阵论][Link 3]中， [实数][Link 4] **正交矩阵**是 [方块矩阵][Link 5] Q，它的 [转置矩阵][Link 6]是它的 [逆矩阵][Link 7]:

，如果正交矩阵的 [行列式][Link 8]为 +1，则我们称之为 **特殊正交矩阵**:

## 概述 ##

要看出与内积的联系，考虑在 *n* 维实数 [内积空间][Link 9]中的关于正交基写出的向量 **v**。 **v** 的长度的平方是 **vv**。如果矩阵形式为 *Q* **v** 的线性变换保持了向量长度，则

所以有限维线性 [等距同构][Link 10]，比如 [旋转][Link 11]、 [反射][Link 12]和它们的组合，都产生正交矩阵。反过来也成立: 正交矩阵蕴涵了正交变换。但是， [线性代数][Link 13]包括了在既不是有限维的也不是同样维度的空间之间的 [正交变换][Link 14]，它们没有等价的正交矩阵。

有多种原由使正交矩阵对理论和实践是重要的。 *n*× *n* 正交矩阵形成了一个 [群][Link 15]，即指示为 *O*( *n*) 的 [正交群][Link 16]，它和它的子群广泛的用在数学和物理科学中。例如，分子的 [点群][Link 17]是 *O*(3) 的子群。因为浮点版本的正交矩阵有有利的性质，它们是字数值线性代数中很多算法比如 QR分解的关键，通过适当的规范化， [离散余弦变换][Link 18] (用于 [MP3][] 压缩)可用正交矩阵表示。

## 例子 ##

下面是一些小正交矩阵的例子和可能的解释。

恒等变换。 旋转 16.26°。 针对 *x* 轴反射。 旋转反演(rotoinversion): 轴 (0,-3/5,4/5)，角度90°。 置换坐标轴。

## 基本构造 ##

### 低维度 ###

最简单的正交矩阵是 1×1 矩阵 \[1\] 和 \[−1\]，它们可分别解释为恒等和实数线针对原点的反射。

如下形式的 2×2 矩阵

它的正交性要求满足三个方程

在考虑第一个方程时，不丢失一般性而设 *p* = cos θ, *q* = sin θ；因此要么 *t* = − *q*, *u* = *p* 要么 *t* = *q*, *u* = − *p*。我们可以解释第一种情况为旋转 θ(θ = 0 是单位矩阵)，第二个解释为针对在角 θ/2 的直线的反射。

旋转 反射 在 45°的反射对换 *x* 和 *y*；它是 [置换矩阵][Link 19]，在每列和每行带有一个单一的 1(其他都是 0):

单位矩阵也是置换矩阵。

反射是它自己的逆，这蕴涵了反射矩阵是 [对称][Link 20]的(等于它的转置矩阵)也是正交的。两个旋转矩阵的积是一个旋转矩阵，两个反射矩阵的积也是旋转矩阵。

### 更高维度 ###

不管维度，总是可能把正交矩阵按纯旋转与否来分类，但是对于 3×3 矩阵和更高维度矩阵要比反射复杂多了。例如，

和 表示通过原点的 [反演][Link 21]和关于 *z* 轴的旋转反演(逆时针旋转90°后针对 *x*\- *y*平面反射，或逆时针旋转 270°后对原点反演)。

旋转也变得更加复杂；它们不再由一个角来刻画，并可能影响多于一个平面子空间。尽管经常以一个轴和角来描述 3×3 旋转矩阵，在这个维度旋转轴的存在是偶然的性质而不适用于其他维度。

但是，我们有了一般适用的基本建造板块如置换、反射、和旋转。

### 基本变换 ###

最基本的置换是换位(transposition)，通过交换单位矩阵的两行得到。任何 *n*× *n* 置换矩阵都可以构造为最多 *n*−1 次换位的积。 构造自非零向量 **v** 的 Householder反射为

这里的分子是对称矩阵，而分母是 **v** 的平方量的一个数。这是在垂直于 **v** 的超平面上的反射(取负平行于 **v** 任何向量分量)。如果 **v** 是单位向量，则 *Q* = *I*−2 **vv** 就足够了。Householder 反射典型的用于同时置零一列的较低部分。任何 *n*× *n* 正交矩阵都可以构造为最多 *n* 次这种反射的积。

Givens旋转作用于由两个坐标轴所生成的二维(平面)子空间上，按选定角度旋转。它典型的用来置零一个单一的次对角线元素(subdiagonal entry)。任何 *n*× *n* 的旋转矩阵都可以构造为最多 *n*( *n*−1)/2 次这种旋转的积。在 3x3 矩阵的情况下，三个这种旋转就足够了；并且通过固定这个序列，我们可以用经常叫做 [欧拉角][Link 22]的三个角来(尽管不唯一)描述所有 3×3 旋转矩阵。

雅可比旋转有同 Givens 旋转一样的形式，但是被用做 [相似变换][Link 23]，选择来置零 2×2 子矩阵的两个远离对角元素(off-diagonal entry)。

## 性质 ##

### 矩阵性质 ###

实数方块矩阵是正交的，当且仅当它的列形成了带有普通欧几里得 [点积][Link 24]的 [欧几里得空间][Link 25] **R** 的正交规范基，它为真当且仅当它的行形成 **R** 的正交基。假设带有正交(非正交规范)列的矩阵叫正交矩阵可能是诱人的，但是这种矩阵没有特殊价值而没有特殊名字；他们只是 *MM* = *D*， *D* 是 [对角矩阵][Link 26]。

任何正交矩阵的行列式是 +1 或 −1。这可从关于行列式的如下基本事实得出:

反过来不是真的；有 +1 行列式不保证正交性，即使带有正交列，可由下列反例证实。

对于置换矩阵，行列式是 +1 还是 −1 匹配置换是偶还是奇的 [标志][Link 27]，行列式是行的交替函数。

比行列式限制更强的是正交矩阵总可以是在 [复数][Link 28]上可对角化来展示 [特征值][Link 29]的完全的集合，它们全都必须有(复数) [绝对值][Link 30] 1。

### 群性质 ###

正交矩阵的逆是正交的，两个正交矩阵的积是正交的。事实上，所有 *n*× *n* 正交矩阵的集合满足群的所有公理。它是 *n*( *n*−1)/2 维的 [紧致][Link 31] [李群][Link 32]，叫做正交群并指示为 *O*( *n*)。

行列式为 +1 的正交矩阵形成了路径连通的子群指标为 2 的 *O*( *n*) [正规子群][Link 33]，叫做旋转的特殊正交群 *SO*( *n*)。 [商群][Link 34] *O*( *n*)/ *SO*( *n*) 同构于 *O*(1)，带有依据行列式选择 \[+1\] 或 \[−1\] 的投影映射。带有行列式 −1 的正交矩阵不包括单位矩阵，所以不形成子群而只是 [陪集][Link 35]；它也是(分离的)连通的。所以每个正交群被分为两个部分；因为投影映射 [分裂][Link 36]， *O*( *n*) 是 *SO*( *n*) 与 *O*(1)的 [半直积][Link 37]。用实用术语说，一个相当的陈述是任何正交矩阵可以通过采用一个旋转矩阵并可能取负它的一列来生成，如我们在 2×2 矩阵中看到的。如果 *n* 是奇数，则半直积实际上是 [直积][Link 38]，任何正交矩阵可以通过采用一个旋转矩阵并可能取负它的所有列来生成。

现在考虑 ( *n*\+1)×( *n*\+1) 右底元素等于 1 的正交矩阵。最后一列(和最后一行)的余下元素必须是零，而任何两个这种矩阵的积有同样的形式。余下的矩阵是 *n*× *n* 正交矩阵；因此 *O*( *n*) 是 *O*( *n*\+1) (和所有更高维群)的子群。

因为 Householder 正交矩阵形式的基本反射可把任何正交矩阵简约成这种约束形式，一系列的这种反射可以把任何正交矩阵变回单位矩阵；因此正交群是反射群。最后一列可以被固定为任何单位向量，并且每种选择给出不同的 *O*( *n*) 在 *O*( *n*\+1) 中的复本；以这种方式 *O*( *n*\+1) 是在单位球 *S* 与纤维 *O*( *n*) 上的 [丛][Link 39]。

类似的， *SO*( *n*) 是 *SO*( *n*\+1) 的子群；任何特定正交矩阵可以使用类似过程通过 Givens 平面旋转来生成。丛结构持续: *SO*( *n*) ↪ *SO*( *n*\+1) → *S*。一个单一旋转可以在最后一列的第一行生成一个零，而 *n*−1 次旋转序列将置零 *n*× *n* 旋转矩阵的除了最后一列的最后一行的所有元素。因为平面是固定的，每次旋转只有一个自由度，就是它的角度。通过归纳， *SO*( *n*) 因此有

自由度， *O*( *n*) 也是。

置换矩阵简单一些；它们不形成李群，只是一个有限群， *n*! 次 [对称群][Link 40] *Sn*。通过同类的讨论， *Sn* 是 *Sn*\+1 的子群。偶置换生成行列式 +1 的置换矩阵的子群， *n*!/2 次交错群。

### 规范形式 ###

更广泛的说，任何正交矩阵的效果分离到在正交二维空间上的独立动作。就是说，如果 *Q* 是狭义正交的，则你可以找到(旋转)改变基的一个正交矩阵 *P*，把 *Q* 带回到分块对角形式:

( *n* 偶数)， ( *n* 奇数)。 这里的矩阵 *R*1,..., *Rk* 是 2×2 旋转矩阵，而余下的元素是零。作为例外，一个旋转块可以是对角的， ± *I*。因此如果需要的话取负一列，并注意 2×2 反射可对角化为 +1 和 −1，任何正交矩阵可变为如下形式

, 矩阵 *R*1,…, *Rk* 给出位于 [复平面][Link 41]中单位圆上的特征值的共轭对；所以这个分解复合确定所有带有绝对值 1 的特征值。如果 *n* 是奇数，至少有一个实数特征值 +1 或 −1；对于 3×3 旋转，关联着 +1 的特征向量是旋转轴。

## 数值线性代数 ##

### 利益 ###

[数值分析][Link 42]自然的利用了正交矩阵的很多数值线性代数的性质。例如，经常需要计算空间的正交基，或基的正交变更；二者都采用了正交矩阵的形式。有行列式 ±1 和所有模为 1 的特征值是对数值稳定性非常有利的。 一个蕴涵是 [条件数][Link 43]为 1 (这是极小的)，所以在乘以正交矩阵的时候错误不放大。很多算法为此使用正交矩阵如 Householder反射和 Givens旋转。有帮助的不只是正交矩阵是可逆的，还有它的逆矩阵本质上是免花费的，只需要对换索引(下标)。

置换是很多算法成功的根本，包括有局部定支点(partial pivoting)的运算繁重的 [高斯消去法][Link 44](这里的置换用来定支点)。但是它们很少明显作为矩阵出现；它们的特殊形式允许更有限的表示，比如 *n* 个索引的列表。

同样的，使用 Householder 和 Givens 矩阵的算法典型的使用特殊方法的乘法和存储。例如，Givens 旋转只影响它所乘的矩阵的两行，替代完全的 *n* 次的 [矩阵乘法][Link 45]为更有效的 *n* 次运算。在使用这些反射和旋转向矩阵介入零的时候，腾出的空间足够存储充足的数据来重生成这个变换。

### 分解 ###

一些重要的 [矩阵分解][Link 46](Golub & Van Loan, 1996)涉及到了正交矩阵，包括:

QR分解 *M* = *QR*, *Q* 正交， *R* 上三角。 [奇异值分解][Link 47] *M* = *UΣV*, *U* 和 *V* 正交， *Σ* 非负对角。 谱分解 *S* = *QΛQ*, *S* 对称， *Q* 正交， *Λ* 对角。 极分解 *M* = *QS*, *Q* 正交， *S* 对称非负确定。

[https_blog.csdn.net_jiangjunshow]: https://blog.csdn.net/jiangjunshow/article/details/77338485
[Link 1]: http://baike.baidu.com/view/10337.htm
[Link 2]: http://baike.baidu.com/view/3122973.htm
[Link 3]: http://baike.baidu.com/view/2092921.htm
[Link 4]: http://baike.baidu.com/view/14749.htm
[Link 5]: http://baike.baidu.com/view/2136856.htm
[Link 6]: http://baike.baidu.com/view/1659988.htm
[Link 7]: http://baike.baidu.com/view/597891.htm
[Link 8]: http://baike.baidu.com/view/111348.htm
[Link 9]: http://baike.baidu.com/view/454412.htm
[Link 10]: http://baike.baidu.com/view/972495.htm
[Link 11]: http://baike.baidu.com/view/131763.htm
[Link 12]: http://baike.baidu.com/view/77128.htm
[Link 13]: http://baike.baidu.com/view/32243.htm
[Link 14]: http://baike.baidu.com/view/689186.htm
[Link 15]: http://baike.baidu.com/view/48541.htm
[Link 16]: http://baike.baidu.com/view/2627450.htm
[Link 17]: http://baike.baidu.com/view/2430960.htm
[Link 18]: http://baike.baidu.com/view/391624.htm
[MP3]: http://baike.baidu.com/view/1310.htm
[Link 19]: http://baike.baidu.com/view/3150602.htm
[Link 20]: http://baike.baidu.com/view/90573.htm
[Link 21]: http://baike.baidu.com/view/750178.htm
[Link 22]: http://baike.baidu.com/view/90087.htm
[Link 23]: http://baike.baidu.com/view/689137.htm
[Link 24]: http://baike.baidu.com/view/2744555.htm
[Link 25]: http://baike.baidu.com/view/42815.htm
[Link 26]: http://baike.baidu.com/view/1445899.htm
[Link 27]: http://baike.baidu.com/view/28475.htm
[Link 28]: http://baike.baidu.com/view/10078.htm
[Link 29]: http://baike.baidu.com/view/689250.htm
[Link 30]: http://baike.baidu.com/view/220956.htm
[Link 31]: http://baike.baidu.com/view/686985.htm
[Link 32]: http://baike.baidu.com/view/53493.htm
[Link 33]: http://baike.baidu.com/view/1004664.htm
[Link 34]: http://baike.baidu.com/view/1962683.htm
[Link 35]: http://baike.baidu.com/view/1004659.htm
[Link 36]: http://baike.baidu.com/view/552808.htm
[Link 37]: http://baike.baidu.com/view/2988681.htm
[Link 38]: http://baike.baidu.com/view/360052.htm
[Link 39]: http://baike.baidu.com/view/241580.htm
[Link 40]: http://baike.baidu.com/view/1536658.htm
[Link 41]: http://baike.baidu.com/view/795684.htm
[Link 42]: http://baike.baidu.com/view/295760.htm
[Link 43]: http://baike.baidu.com/view/1014731.htm
[Link 44]: http://baike.baidu.com/view/1174451.htm
[Link 45]: http://baike.baidu.com/view/2455255.htm
[Link 46]: http://baike.baidu.com/view/2362946.htm
[Link 47]: http://baike.baidu.com/view/3068725.htm