Scipy中的稀疏矩阵--介绍

稀疏矩阵介绍

1 介绍

（稠密）矩阵是：

• 数学对象
• 存储二维数组的数据结构

 重要特质：

• 所有元素内存一次性分配，在Numpy ndarray中通常是一个连续的块
• 快速访问单个元素

1.1 为什么使用稀疏矩阵？

• 内存增长速度为n^2
• 例子（双精度矩阵）：

>>> import numpy as np
>>> import matplotlib.pyplot as plt
>>> x = np.linspace(0, 1e6, 10)
>>> plt.plot(x, 8.0 * (x**2) / 1e6, lw=5)
[<matplotlib.lines.Line2D object at ...>]
>>> plt.xlabel(’size n’)
<matplotlib.text.Text object at ...>
>>> plt.ylabel(’memory [MB]’)
<matplotlib.text.Text object at ...>

1.2 稀疏矩阵VS稀疏矩阵存储方案

• 稀疏矩阵是一个大部分元素为空的矩阵
• 存储所有的零元素是非常浪费的->仅存储非零元
• 考虑压缩
• 优势：节约大量内存
• 弊端：依赖于实际的存储方案，有时候不可行

1.3 典型应用

• 求解偏微分方程（PDEs）：有限元方法、工程机械，电工，物理等
• 图理论
• 等等

1.4 先决条件

最新版本的：

• numpy
• scipy
• matplotlib（可选）
• ipython（功能强大）

1.5 稀疏结构可视化

• matplotlib中的spy()
• 图例：