正交矩阵
数学术语
如果AAT=E(E为单位矩阵,AT表示“矩阵A的转置矩阵”)或ATA=E,则n阶实矩阵A称为正交矩阵。正交矩阵是实数特殊化的酉矩阵,因此总是属于正规矩阵。尽管在这里只考虑实数矩阵,但这个定义可用于其元素来自任何域的矩阵。正交矩阵毕竟是从内积自然引出的,所以对于复数的矩阵则导致了归一要求。正交矩阵不一定是实矩阵。实正交矩阵(即该正交矩阵中所有元都是实数)可以看作是一种特殊的酉矩阵,但也存在一种复正交矩阵,这种复正交矩阵不是酉矩阵。
定义
如果:AAT=E(E为单位矩阵,AT表示“矩阵A的转置矩阵”。)或ATA=E,则n阶实矩阵A称为正交矩阵,若A为正交阵,则满足以下条件:
1. AT是正交矩阵
2. (E为单位矩阵)
3. AT的各行是单位向量且两两正交
4. AT的各列是单位向量且两两正交
5. (Ax,Ay)=(x,y)x,y∈R
6. |A|=1或-1
7.
8. 正交矩阵通常用字母Q表示。
9. 举例:
若A=[r11r12r13;r21r22r23;r31r32r33],则有:
定理
矩阵论中,实数正交矩阵是方块矩阵Q,它的转置矩阵是它的逆矩阵,如果正交矩阵的行列式为+1,则称之为特殊正交矩阵。
1.方阵A正交的充要条件是A的行(列)向量组是单位正交向量组
2.方阵A正交的充要条件是A的n个行(列)向量是n维向量空间的一组标准正交基
3.A是正交矩阵的充要条件是:A的行向量组两两正交且都是单位向量
4.A的列向量组也是正交单位向量组。
5.正交方阵是欧氏空间中标准正交基到标准正交基的过渡矩阵。
举例
下面是一些小正交矩阵的例子和可能的解释。
恒等变换
恒等变换就是把一个解析式变成与它恒等的另一个解析式.使用恒等变换往往是在碰到的问题比较繁杂、一时难以下手的时候,通过恒等变换把要解决的问题简化,由未知到已知,最终解决问题.所以,恒等变换的特点就是:将复杂的问题通过表达形式的变形转化成容易解决的简单问题。
X轴反射
反射变换(refIection)又称为镜像反射或镜像变换,类似于一个对象在一面镜子中的影子。二维平面上给定一条直线可以做关于直线的镜像反射;三维空间中,给定一个平面,可以做关于这个平面的镜像反射。对于矩阵变换如图1所示。
关于正交矩阵其他解释还有:旋转反演(rotoinversion):轴(0,-3/5,4/5),角度90°;置换坐标轴等。
基本构造
低维度构造
最简单的正交矩阵是1×1矩阵[1]和[−1],它们可分别解释为恒等和实数线针对原点的反射。
它的正交性要求满足三个方程,在考虑第一个方程时,不丢失一般性而设p=cosθ,q=sinθ;因此要么t=−q,u=p要么t=q,u=−p。可以解释第一种情况为旋转θ(θ=0是单位矩阵),第二个解释为针对在角θ/2的直线的反射。
旋转反射在45°的反射对换x和y;它是置换矩阵,在每列和每行带有一个单一的1(其他都是0):单位矩阵也是置换矩阵。
反射是它自己的逆,这蕴涵了反射矩阵是对称的(等于它的转置矩阵)也是正交的。两个旋转矩阵的积是一个旋转矩阵,两个反射矩阵的积也是旋转矩阵。
更高维度构造
如果不管维度,总是有可能把正交矩阵按纯旋转与否来分类的,但是对于3×3矩阵和更高维度矩阵要比反射复杂多了。例如,表示通过原点的反演和关于z轴的旋转反演(逆时针旋转90°后针对x-y平面反射,或逆时针旋转270°后对原点反演)。旋转也变得更加复杂;它们不再由一个角来刻画,并可能影响多于一个平面子空间。尽管经常以一个轴和角来描述3×3旋转矩阵,在这个维度旋转轴的存在是偶然的性质而不适用于其他维度。但是,一般适用的基本建造板块如置换、反射、和旋转可以满足这些情况。
基本变换
正交矩阵的最基本置换是换位(transposition),通过交换单位矩阵的两行得到。任何n×n置换矩阵都可以构造为最多n−1次换位的积。构造自非零向量v的Householder反射,这里的分子是对称矩阵,而分母是v的平方量的一个数。这是在垂直于v的超平面上的反射(取负平行于v任何向量分量)。如果v是单位向量,则Q=I−2vv就足够了。Householder反射典型的用于同时置零一列的较低部分。任何n×n正交矩阵都可以构造为最多n次这种反射的积。
Givens旋转作用于由两个坐标轴所生成的二维(平面)子空间上,按选定角度旋转。它典型的用来置零一个单一的次对角线元素(subdiagonalentry)。任何n×n的旋转矩阵都可以构造为最多n(n−1)/2次这种旋转的积。在3x3矩阵的情况下,三个这种旋转就足够了;并且通过固定这个序列,可以用经常叫作欧拉角的三个角来(尽管不唯一)描述所有3×3旋转矩阵。
雅可比旋转有同Givens旋转一样的形式,但是被用做相似变换,选择来置零2×2子矩阵的两个远离对角元素(off-diagonalentry) 。
基本特性
矩阵性质
实数方块矩阵是正交的,当且仅当它的列形成了带有普通欧几里得点积欧几里得空间R的正交规范基,它为真当且仅当它的行形成R的正交基。假设带有正交(非正交规范)列的矩阵叫正交矩阵可能是诱人的,但是这种矩阵没有特殊价值而没有特殊名字;他们只是MM=D,D是对角矩阵
1.逆也是正交阵;
2.积也是正交阵;
3.行列式的值为正1或负1。
任何正交矩阵的行列式是+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的特征向量是旋转轴
理论实践
有多种原由使正交矩阵对理论和实践是重要的。n×n正交矩阵形成了一个群,即指示为O(n)的正交群,它和它的子群广泛的用在数学和物理科学中。例如,分子的点群是O(3)的子群。因为浮点版本的正交矩阵有有利的性质,它们是字数值线性代数中很多算法比如QR分解的关键,通过适当的规范化,离散余弦变换(用于MP3压缩)可用正交矩阵表示 。
分析算法
利益
数值分析自然的利用了正交矩阵的很多数值线性代数的性质。例如,经常需要计算空间的正交基,或基的正交变更;二者都采用了正交矩阵的形式。有行列式±1和所有模为1的特征值是对数值稳定性非常有利的。一个蕴涵是条件数为1(这是极小的),所以在乘以正交矩阵的时候错误不放大。很多算法为此使用正交矩阵如Householder反射和Givens旋转。有帮助的不只是正交矩阵是可逆的,还有它的逆矩阵本质上是免花费的,只需要对换索引(下标)。
置换是很多算法成功的根本,包括有局部定支点(partialpivoting)的运算繁重的高斯消去法(这里的置换用来定支点)。但是它们很少明显作为矩阵出现;它们的特殊形式允许更有限的表示,比如n个索引的列表。
同样的,使用Householder和Givens矩阵的算法典型的使用特殊方法的乘法和存储。例如,Givens旋转只影响它所乘的矩阵的两行,替代完全的n次的矩阵乘法为更有效的n次运算。在使用这些反射和旋转向矩阵介入零的时候,腾出的空间足够存储充足的数据来重生成这个变换。
分解
一些重要的矩阵分解(Golub&VanLoan,1996)涉及到了正交矩阵,包括:
参考资料
最新修订时间:2024-10-11 21:04
目录
概述
定义
参考资料