对角化大型矩阵的工具

你试过SVD吗, svds 在matlab中实现稀疏矩阵。

你考虑过 following property : exp(A*t)=L^(-1){(sI-A)^(-1)} 其中L^(-1)是拉普拉斯逆变换?-前提是可以反转(sI-A)

如果您可以访问64位机器,并且使用64位支持编译的八度音阶,那么您可能可以绕过这个问题。

另外,我不知道您在什么平台上运行所有这些,但是在基于UNIX的系统中,您可以使用 ulimit

例如,您可以运行

ulimit -u unlimited

ulimit -s [stacksize]

增加堆栈大小限制。

或者,如果你的矩阵是这样构造的:A=[B zeros;zeros C],那么你可以把B和C分别对角化,然后把它们放在一个矩阵中。我猜你也可以为eig做类似的事情。

igraph package 和功能 arpack 哪个是接口 ARPACK library