通过 Householder 的 QR:计算复杂度较低的变体?

计算科学 线性代数 算法 线性求解器
2021-12-26 23:34:59

我是一个概率论者,需要对一个相当大的线性最小二乘问题进行一些计算,所以我正在尝试尽可能优化计算。

在使用 Householder 反射计算超定系统的矩阵表示的 QR 分解时,常用算法建议您计算一些和超平面,遵循几何动机。AAx=buuuUU

这产生了 Householder 反射器H=I2u,uuu

现在,为了消除 A 列的子列从左侧乘以AAH

  1. 的计算需要次乘法。uum2
  2. 的直接计算需要次乘法。HAm2n

乘积等价于(或扩展)为,它需要乘法——通过关联性。AHA2u,uu(uA)2mn

我可能遗漏了一些相当明显的东西。不过,有没有办法避免每次迭代都进行这种计算,并从次触发器中降低复杂度?2mn

0个回答
没有发现任何回复~