稀疏矩阵的线性代数

对于稀疏矩阵来说，其线性代数操作可以使用 scipy.sparse.linalg 实现：

In [1]:

import scipy.sparse.linalg

矩阵操作

scipy.sparse.linalg.inv
- 稀疏矩阵求逆
scipy.sparse.linalg.expm
- 求稀疏矩阵的指数函数

矩阵范数

scipy.sparse.linalg.norm
- 稀疏矩阵求范数

线性方程组求解

提供了一系列求解方法： http://docs.scipy.org/doc/scipy/reference/sparse.linalg.html#solving-linear-problems

主要使用的是迭代方法求解。

特征值分解和奇异值分解

对于特别大的矩阵，原来的方法可能需要太大的内存，考虑使用这两个方法替代：

scipy.sparse.linalg.eigs
- 返回前 k 大的特征值和特征向量
scipy.sparse.linalg.svds
- 返回前 k 大的奇异值和奇异向量

所有的这些操作既可以在稀疏矩阵上使用，也可以在普通矩阵上使用。