该项目利用Householder反射矩阵实现对任意维度的实数矩阵进行QR分解。其核心原理是通过构造一系列对称且正交的Householder矩阵,依次消去目标矩阵每列在对角线以下的元素,最终将其逐步化简为上三角矩阵R。与此同时,所有步骤中产生的反射矩阵的累积乘积构成了分解得到的正交矩阵Q。
这种方法在数值计算领域具有极高的稳定性,相比于Gram-Schmidt正交化方法,Householder QR分解能更有效地抑制浮点运算中的舍入误差积累,特别是在处理病态矩阵时表现出色。程序内部实现了自动化的单位反射向量