20得票3回答
改善一个状况不佳的矩阵

我有一个矩阵条件很差,其 rcond() 接近于零,因此,该矩阵的逆矩阵不正确。我尝试使用 pinv(),但仍无法解决问题。这是我计算矩阵逆的方法:X = (A)\(b); 我寻找了解决这个问题的方法,并在这个链接(最后一个解决方案)中找到了改进矩阵的方法。 那里的解决方案建议使用以下内容:A...

20得票1回答
使用Python计算非常大稀疏矩阵的特征值

我有一个非常大的稀疏矩阵,它表示马尔可夫链中的转移矩阵,即矩阵的每一行之和等于一,并且我想找到第一个特征值及其相应的向量,使其小于一。我知道特征值在区间[-1,1]内有界,并且它们都是实数(非复数)。 我试图使用Python的scipy.sparse.eigs函数来计算这些值,然而,函数的一个...

20得票1回答
如何在Eigen中将稀疏矩阵转换为密集矩阵

有没有一种简单快速的方法将稀疏矩阵转换为双精度密集矩阵? 因为我的 SparseMatrix 经过一些矩阵乘积之后不再是稀疏的,而变成了密集的。 我还有一个问题:Eigen库具有极佳的性能,这是怎么实现的呢?我不明白为什么,因为只有头文件,没有编译源码。

19得票1回答
Scipy稀疏矩阵CSR转换为TensorFlow SparseTensor - Mini-Batch梯度下降

我有一个来自SVM-Light格式的稀疏TF-IDF特征矩阵创建的Scipy稀疏CSR矩阵。 特征数量非常庞大且稀疏,因此必须使用SparseTensor,否则速度太慢。例如,特征数为5,样本文件可能如下所示:0 4:1 1 1:3 3:4 0 5:1 0 2:1 解析后,训练集如下:trai...

19得票5回答
在非常大的稀疏矩阵上应用PCA

我正在使用R进行文本分类任务,并获得了一个大小为22490 x 120,000的文档-词条矩阵(只有400万个非零条目,不到1%的条目)。现在我想通过利用PCA(主成分分析)来降低维度。不幸的是,R无法处理这个巨大的矩阵,因此我把这个稀疏矩阵存储在“Matrix Market Format”文...

18得票14回答
如何高效存储具有高度冗余值的矩阵

我有一个非常大的矩阵(100M行x 100M列),其中有很多相邻重复的值。例如:8 8 8 8 8 8 8 8 8 8 8 8 8 8 4 8 8 1 1 1 1 1 8 8 8 8 8 4 8 8 1 1 1 1 1 8 8 8 8 8 4 8 8 1 1 1 1 1 8 8 8 8 8 4...

18得票2回答
另一个生命游戏问题(无限网格)?

我一直在玩康威生命游戏,最近发现了一些非常快的实现方式,比如Hashlife和Golly。(在此下载Golly - http://golly.sourceforge.net/) 有一件事情让我困惑不解,那就是编程人员是如何实现无限的网格的呢?我们不能保留任何无限的数组,如果你运行Golly并...

18得票2回答
稀疏矩阵的内存高效创建

我有一个由6000个独特字符串组成的列表,其中包含50000个字符串向量的各种组合。 目标:我想将它们转换为“相对频率”(table(x)/length(x)),并将它们存储在稀疏矩阵中。低内存消耗比速度更重要。目前内存是瓶颈。(尽管源数据大约有50 mb,目标格式的数据大小仅为10mb--...

18得票4回答
在glmnet中出现的R错误:外部函数调用中出现NA/NaN/Inf。

我正在尝试使用glmnet创建模型(目前使用cv找到lambda值),但是出现了错误NA/NaN/Inf in foreign function call (arg 5)。我认为这与我的数据集中的NA值有关,因为当我删除所有带有NA的数据点时,命令成功运行。 我原本以为glmnet可以处理N...

18得票2回答
将Eigen稠密矩阵转换为稀疏矩阵

如何将 Eigen::Matrix<double,Dynamic,Dynamic> 转换为 Eigen::SparseMatrix<double>?我正在寻找一种更好的方法,而不是通过迭代密集矩阵来实现。