Matlab - 计算大稀疏（A'*A）矩阵的最大特征值

Question

我有一个大的 (400K*400K) 稀疏矩阵，我需要计算 A'*A 的最大特征值。

问题是由于内存问题，Matlab 甚至无法计算 A'。
我还尝试了[a,b,c] = find(A)，然后通过创建转置稀疏矩阵进行转置，但尽管find() 有效，但创建 sprase 却没有。

有没有很好的解决方案？它可以在 matlab 函数中，也可以在其他技术中计算这种乘法的最大特征值。

谢谢。

Answer 1

如果 A 是稀疏的，请参阅此 thread 和 this documentation 中的一些讨论（基本上是逐个进行）以了解转置它的方法等。

但是现在你需要计算B=A'*A。问题是，它仍然稀疏吗？假设是这样，使用链接中提到的先前技术应该没有问题。

那么在你获得B=A'*A之后，使用eigs

eigs(B,1)

求最大量级特征值。