R有一个qr()
函数,它使用LINPACK或LAPACK执行QR分解(根据我的经验,后者快5% )。返回的主要对象是一个包含在上三角矩阵R中的矩阵"qr“(即R=qr[upper.tri(qr)]
)。到目前一切尚好。qr的下三角部分包含“紧凑形式”的Q。可以使用qr.Q()
从qr分解中提取Q。我想找出qr.Q()
的倒数。换句话说,我确实有Q和R,并想把它们放在一个"qr“对象中。R是微不足道的,但Q不是。我们的目标是应用qr.solve()
,这比大型系统上的solve()
快得多。
https://stackoverflow.com/questions/3031215
复制相似问题