“稠密数值计算通常用Householder 反射。第 $k$ 步对 $A {k:m,k}$ 构造反射,左乘后把对角线下方元素清零;一串反射向量隐式表示 $Q$,尾部更新最终产生 $R$。实 $…”
形式陈述
在标准内积空间
给定
其中
为什么正好能消去尾部?令
对任意向量
以
直觉
Gram–Schmidt 一次去掉一个已有方向的分量,Householder 反射则选择一面镜子,让整根向量在一次变换后与坐标轴重合。镜子保持长度和夹角,所以在连续施加反射时,不会主动放大
向量
例子与边界
对
这里
若
还需区分代数公式和安全实现。输入为零向量时已经无需消去,应采用恒等变换;不能代入
反射不是任意旋转。实 Householder 矩阵的行列式为
推论与应用
QR 分解逐列使用 Householder 反射,把每列对角线下方的整段向量清零。相同原语还用于 Hessenberg 化、双对角化与许多特征值、奇异值算法,因为酉变换既保留范数,又适合用紧凑向量实现。
Gram–Schmidt仍是理解正交化与张成空间保持的直接方法;对于稠密数值 QR,Householder 通常提供更可靠的正交性。两种方法解决同一精确结构目标,却拥有不同计算图和舍入路径。
同一个实内积反射公式也能约束离散几何。有限结晶根系要求一组非零法向在每个根反射下保持不变,并要求反射系数
参考资料
- 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.