我有一个 dgTMatrix 类型的稀疏矩阵，它是用 R 中的矩阵dgTMatrix构建的，我正在尝试使用 QR 分解求解Ax=b中的向量x ，但它无法正常工作。 例如下面我正在求解随机A, b ，您可以看到该方法适用于Ax-b = 0的范数 我A和b是这样的： 当我对A进行 QR 分解时，我得到以 ...
在 R 中的稀疏矩阵上应用 Matrix::qr() 时，output 与 base::qr 有很大不同。 有V、beta、p、R、q但没有rank和pivot。 下面是一个小示例代码。 我想检测 A 稀疏矩阵的线性相关列，这需要 pivot 和秩。 我应该如何获得这些信息？ 还有另一个相关但未回 ...
我需要在某些特征相关的数据集上应用 Python 线性回归。 Scikit 的线性回归使用奇异值分解来最小化平方误差。 但是，这种方法没有考虑变量之间的相关性。 R 使用 QR 分解，这是一种使用“lm”方法对协变量进行去相关的技术。 相关变量的系数作为 NaN 值返回，这对于我想要做的事情很 ...
我有大量矩阵，我需要对其执行 QR 分解并存储生成的 Q 矩阵（归一化以使 R 矩阵在其对角线上具有正数）。 除了使用qr()函数之外，还有其他方法吗？ 这是工作示例： 所以计算 62,500 QR 分解需要 3.07 秒。 我想知道是否有更快的东西？ ...
我在处理排名不足的情况时非常依赖qr()函数，但最近遇到了一些无法正常工作的示例。 考虑下面的矩阵badX ： badX <- structure(c(-1.641906809157e-10, 0, 0, 0, 0, -0.5, 0, 0, -1.10482935525559e-16, ...
我用拉帕克。 我正在尝试在 C 中对复杂数据进行 QR 分解。 为此，我编写了 function （基于 Haatschii 代码How to get the Q from the QR factorization output? ）： 值得注意的是，Alireza 的作者也遇到了 functio ...
我正在尝试使用 QR 方法找到矩阵 A 的特征向量。 我找到了对应于最大特征值的特征值和特征向量。 如何在不使用 numpy.linalg.eig 的情况下找到特征向量的 rest？ 接下来我检查条件： 检查条件： ...
随着 Vivado HLS 安装，提供了一个线性代数库。 在这个库中有这个函数 QRF_ALT，它应该是 Qrf-basic 的高吞吐量版本：使用给定旋转的 QR 分解。 我的问题是它如何实现高吞吐量？ 从代码附带的简短评论中，我的猜测是它使用了批量矩阵行的并行处理，但是它是如何调度它的，我无法理 ...
数组中的所有标准化元素都设置为 0。 我正在尝试为阵列创建反射器。 当我通过 Python 中的第一个元素(x[1:] /= x[0], x[0] = 1)对数组中的所有元素进行归一化时，我只会在第一个元素中得到全零和 1。 结果应如下所示(u = [1, 4/9, 1/9, 3/9, 1/9]) ...
我正在尝试通过家庭反射器实现 QR 分解。 在一个非常简单的数组上尝试这个时，我得到了奇怪的数字。 任何人都可以告诉我，为什么在 function 定义的最后一行在 vec 和 vec.T 之间使用 @ vs * 运算符可以获得主要奖励积分。 截至今天早上，这已经难倒了两位数学/计算机科学博士。 ...
我需要使矩阵A所有其他列正交于其列j 。 我使用以下算法： 这个想法来自用Gram-Schmidt过程进行QR分解 。 但是由于Gram-Schmidt过程，该代码并未经过优化且不稳定。 Numpy是否提供任何方法来计算这些向量的正交投影？ 带户主矩阵 我 ...
我试图用 QR 分解和 linalg.solve 解决一个超定系统，但我得到的错误是 LinAlgError: 数组的最后 2 个维度必须是 square 。 当 R 数组不是正方形时会发生这种情况，对吗？ 代码看起来像这样 有没有办法以更有效的方式为任意 A 维度编写此代码？ 如果没有，我 ...
我目前正在尝试为矩形矩阵实现基于 Householder 的 QR 分解，如http://eprints.ma.man.ac.uk/1192/1/qrupdating_12nov08.pdf （第 3、4、5 页）中所述。 显然我弄错了一些伪代码，因为 (1) 我的结果不同于numpy.qr.li ...
我对QR分解方法有疑问。 我使用dgeqrf子例程进行分解，但编译器中没有错误，但此后出现了问题。 我还没发现错误在哪里。 另一个问题是，如果A矩阵为零，则A = Q * R =>可以分解为零或失去秩。 ...
我有以下算法实现 我的问题是在每次迭代（可能是每次重新启动）时获得一个向量（可能是矩阵） vetnorm_r ，其中包含残差的所有范数（作为输出）。 我不知道如何建立这个向量或矩阵。 ％输入A REAL非对称正定矩阵％x REAL初始猜测向量％b REAL右侧向量％M REAL ...
我正在尝试求解A * x = b，其中A具有复杂的值并且很密集。 我使用了cuSolverDN库中的cusolverDnCgeqrf（）方法对一组线性方程组进行QR分解。 但是，我想多次执行此操作以加快处理速度。 是否有此方法的“批处理”版本？ 还是我可以使用另一个CUDA库？ ...
我目前正在R中构建一个应用程序，以计算QR矩阵分解，QR非负矩阵分解和计算ICA。 目前，我正在完成第一个任务。 我收到以下错误： source（“ trial.R”）source（“ trial.R”）中的错误：trial.R：153：0：输入151的意外结束：152：^ ...
我有三个关于张量因子分解的问题。 张量因子分解（分解）的情况（或应用）是什么？ 这有可能成为未来的主流技术？ 你如何使用它？ ...
避免阵列分配有利于提高性能。 但是，我还没有理解什么是可以执行矩阵A的QR分解的最有效方式。 （注意：需要Q和R矩阵） 简单地使用Q, R = qr(A)可能不是最好的想法，因为它分配Q和R两者，其中两者都可以重新分配。 函数qrfact允许以打包格式存储分解 。 但是，我之后 ...
我正在手动尝试计算回归系数，而不是对数据使用任何默认值http://people.sc.fsu.edu/~jburkardt/datasets/regression/x31.txt 这是我的代码，可以正确产生满足A = QR的Q＆R。 但是我找不到系数作为问与答产生问题的维度。 有什么 ...