qr-decomposition相关内容

C ++中的薄QR分解

是否存在易于使用的c ++库,用于矩形矩阵的“薄” QR分解? Eigen似乎仅支持完整的Q矩阵。我可以取一个完整的Q并丢弃一些列,但是不计算它们开头会更有效吗? 解决方案 Newmat 完全满足您的要求。 要将A分解为QR,您可以执行以下操作: 矩阵Q = A; UpperTriangularMatrix R; QRZ(Q,R) 如果A为3x5矩 ..
发布时间:2020-10-27 23:36:06 C/C++开发

R和Matlab中的qr函数

我有一个关于将matlab函数转换为R的问题,我希望有人可以提供帮助. 在matlab和R中使用的标准QR分解称为qr().据我了解,在两种语言中执行qr分解的标准方法是: Matlab: [Q,R] = qr(A) 满足QR = A R: z ..
发布时间:2020-05-06 14:39:27 其他开发

从QR分解中获得帽子矩阵以进行加权最小二乘回归

我正在尝试扩展程序包McSptial的lwr()函数,该函数适合将weigthed回归作为非参数估计.在lwr()函数的核心中,它使用solve()而不是QR分解来反转矩阵,从而导致数值不稳定.我想更改它,但无法弄清楚以后如何从QR分解中获得帽子矩阵(或其他导数). 有数据: set.seed(0); xmat ..
发布时间:2020-04-30 12:22:55 其他开发

如何从QR分解输出中获得Q?

来自LAPACK的 DGEQRF和SGEQRF 返回QR的Q部分以打包格式分解.拆箱似乎需要O(k^3)步骤(k个低排名产品),而且似乎不是很简单.另外,对我来说,k顺序乘法的数值稳定性还不清楚. LAPACK是否包含用于解包Q的子例程,如果没有,我应该怎么做? 解决方案 是的,LAPACK确实提供了一个从基本反射器(即DGEQRF返回的数据部分)中检索Q的例程,称为* DORGQ ..
发布时间:2020-04-30 12:04:07 其他开发