Skip to content

定义Definition

Householder 反射

Householder reflection · Householder transformation

用一个向量紧凑表示酉反射,把整段向量化到单一坐标方向,并作为稳定 QR 的基本原语。

形式陈述 ​

在标准内积空间 Fm 中,对非零向量 v,定义 Householder 反射

H=I−2vv∗v∗v.

矩阵 Pv=vv∗/(v∗v) 是到 span{v} 的正交投影,满足 Pv∗=Pv 与 Pv2=Pv。因此

H∗=H,H∗H=H2=I.

H 在 v 方向上乘以 −1,在 v⊥ 上保持不变;实数情形是关于超平面 v⊥ 的镜面反射,行列式为 −1。

给定 x≠0,令

phase(x1)={x1/|x1|,x1≠0,1,x1=0,α=−phase(x1)‖x‖2,

其中 e1=(1,0,…,0)T 是第一坐标向量。取 v=x−αe1,便有

Hx=αe1.

为什么正好能消去尾部?令 r=‖x‖2。上述选择给出 v∗x=r(r+|x1|) 与 v∗v=2r(r+|x1|),所以 2v(v∗x)/(v∗v)=v,从而 Hx=x−v=αe1。负相位选择还使 x1 与 −α 同向相加,避免在构造 v 时让两个接近的数相减。输入向量 x 后,算法只需输出 v、α 和必要的缩放,而不必形成稠密 H。

对任意向量 y,可按

Hy=y−2vv∗yv∗v

以 Θ(m) 成本应用反射;对尾部矩阵则是一次秩一更新。完成准则是 v∗v 有效、结果有限,并核对目标向量除首分量外已降到舍入尺度。显式构造 m×m 的 H 再做矩阵乘法会把线性存储和秩一更新膨胀成不必要的稠密工作。

直觉

Gram–Schmidt 一次去掉一个已有方向的分量,Householder 反射则选择一面镜子,让整根向量在一次变换后与坐标轴重合。镜子保持长度和夹角,所以在连续施加反射时,不会主动放大 2-范数误差;被“清零”的尾部也可以直接成为三角结构。

向量 v 不是额外的大矩阵,而是这面镜子的法向。保存一串反射向量,就等于紧凑保存一个复杂正交变换;只有读者确实需要 Q 的每个元素时,才值得把它们显式展开。

例子与边界

对 x=(3,4)T,取 α=−5 和 v=x−αe1=(8,4)T,得到

H=(−3/5−4/5−4/53/5),Hx=(−50).

这里 v∗v=80、v∗x=40,秩一公式直接给出 Hx=x−v=(−5,0)T,不必先写出矩阵。反射保持 ‖x‖2=5,却一次消去了第二分量。若后续只需把 H 作用到别的列,使用 v 的秩一公式即可。

若 x 已非常接近正的 e1 方向,却错误选择 α=+‖x‖2,则 v=x−αe1 的首分量由两个近数相减得到,可能丢失大部分有效位;相反相位使其相加并保持可表示尺度。复数情形不能只取实数符号,还必须使用 x1 的相位,否则目标首分量与消去公式不一致。

还需区分代数公式和安全实现。输入为零向量时已经无需消去,应采用恒等变换;不能代入 v=0 后计算 0/0。极大或极小分量也可能使直接求平方和溢出或下溢,因此范数与反射向量应使用缩放计算。精确的 H 保范数,并不意味着浮点生成和应用 H 毫无舍入误差。

反射不是任意旋转。实 Householder 矩阵的行列式为 −1,而纯旋转行列式为 1;两个反射的乘积才可能形成旋转。这个代数边界不妨碍它作为正交变换原语,却不能在几何描述中省略。

推论与应用

QR 分解逐列使用 Householder 反射,把每列对角线下方的整段向量清零。相同原语还用于 Hessenberg 化、双对角化与许多特征值、奇异值算法,因为酉变换既保留范数,又适合用紧凑向量实现。

Gram–Schmidt仍是理解正交化与张成空间保持的直接方法;对于稠密数值 QR,Householder 通常提供更可靠的正交性。两种方法解决同一精确结构目标,却拥有不同计算图和舍入路径。

同一个实内积反射公式也能约束离散几何。有限结晶根系要求一组非零法向在每个根反射下保持不变,并要求反射系数 2(α,β)/(α,α) 为整数。这里输入不仅是镜面,还保留法向长度;B3 与 C3 有同样的Householder反射,却有不同的根长数据。数值QR关注稳定地消去坐标,根系则用这份整性和有限闭包来认证反射词。

参考资料
  • Lloyd N. Trefethen and David Bau III, Numerical Linear Algebra, SIAM, 1997, Lecture 10.
  • Gene H. Golub and Charles F. Van Loan, Matrix Computations, 4th ed., Johns Hopkins University Press, 2013, §5.1.
  • Tobin A. Driscoll and Richard J. Braun, Fundamentals of Numerical Computation, Python online edition, §3.4.1 Householder reflections, accessed 2026.
  • LAPACK Users’ Guide, 3rd ed., SIAM, 1999, QR Factorization.
关系图谱15 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系