豪斯霍尔德变换

豪斯霍尔德变换Householder transformation)或译“豪斯霍德转换”[1],又称初等反射Elementary reflection),最初由A.C Aitken在1932年提出[2]阿尔斯通·斯科特·豪斯霍尔德在1958年指出了这一变换在数值线性代数上的意义[3]。这一变换将一个向量变换为由一个超平面反射的镜像,是一种线性变换。其变换矩阵被称作豪斯霍尔德矩阵,在一般内积空间中的类比被称作豪斯霍尔德算子英语Householder operator。超平面的法向量被称作豪斯霍尔德向量。

定义

 
豪斯霍尔德变换示意图:向量x在豪斯霍尔德向量v的超平面 上的镜像是HxH是豪斯霍尔德矩阵。

如果   给出为单位向量 单位矩阵,则描述上述线性变换的是 豪斯霍尔德矩阵  表示向量  共轭转置

 

性质

豪斯霍尔德矩阵 有如下性质:

  • 它是对称矩阵,即  
  • 它是正交矩阵,即  
  • 它是埃尔米特矩阵,即  
  • 它是对合的,即  

进一步的,  实际上按上面描述的那样反射了   (用它的位置向量   来识别),因为

 

这里的   表示内积。注意   等于从 X 到超平面的距离。

应用

豪斯霍尔德变换可以将向量的某些元素置零,同时保持该向量的范数不变。例如,将非零列向量 变换为单位基向量 乘以一个常数的豪斯霍尔德矩阵为

 

其中豪斯霍尔德向量 满足:

 

Dubrulle 在2000年给出了将豪斯霍尔德变换应用于生成一个一般的稀疏向量的一个数值稳定的算法[4]

对一个矩阵的各个列向量逐一进行相应的豪斯霍尔德变换,可以将这个矩阵变换为上海森伯格矩阵上三角矩阵等形式[5]。后者就是QR分解的豪斯霍尔德算法。

参考文献

  1. ^ 胡家彰. MIMO通讯系统之低复杂度天线选择页面存档备份,存于互联网档案馆
  2. ^ H.W. Turnbull, A.C. Aitken, An Introduction to the Theory of Canonical Matrices, Blackie, London: Glasgrow, 1932
  3. ^ Alston S. Householder, Unitary Triangularization of a Nonsymmetric Matrix, Journal ACM, 5 (4), 1958, 339-342. DOI:10.1145/320941.320947
  4. ^ A.A. Dubrulle, Householder Transformations Revisited, SIAM Journal on Matrix Analysis and Applications, 2001
  5. ^ David D. Morrison, Remarks on the Unitary Triangularization of a Nonsymmetric Matrix, Journal ACM, 7 (2), 1960, 185-186. DOI:10.1145/321021.321030

参见

外部链接