访问稀疏矩阵

以下多条命令可以提供有关稀疏矩阵的非零元素的概要信息：

nnz 返回稀疏矩阵中的非零元素数。

nonzeros 返回包含稀疏矩阵的所有非零元素的列向量。

nzmax 返回为稀疏矩阵的非零项分配的存储空间量。

要尝试上述中的一些命令，请加载提供的稀疏矩阵 west0479，该矩阵是 Harwell-Boeing 集合之一。

该矩阵为八个阶段的化工精馏塔建模。

尝试以下命令。

注意

使用 Ctrl+C 随时停止列出 nonzeros。

请注意，最初在默认情况下，nnz 与 nzmax 的值相同。即，非零元素数等于为非零元素分配的存储位置数。但是，如果将其他的数组元素置零，MATLAB® 不会动态释放内存。将某些矩阵元素的值更改为零时会更改 nnz 的值，但不会更改 nzmax 的值。

但是，可以根据需要将尽可能多的非零元素添加到矩阵中。不受 nzmax 原始值的限制。

对于任何矩阵，无论是满矩阵还是稀疏矩阵，find 函数都会返回非零元素的索引和值。其语法是

find 返回向量 i 中的非零值的行索引、向量 j 中的列索引以及向量 s 中的自身非零值。下面的示例使用 find 查找稀疏矩阵中的非零索引和值。sparse 函数同时使用 find 输出和矩阵大小重新创建矩阵。

由于稀疏矩阵是以压缩稀疏列格式存储的，因此为稀疏矩阵进行索引的相关成本与为满矩阵进行索引的相关成本不同。在只需更改稀疏矩阵中的若干元素时，这类成本可忽略不计，因此，在这类情况下，通常使用常规数组索引来重新分配值：

如果线性索引超过 2^48-1（即当前矩阵中允许的元素数上限），使用线性索引在大型稀疏矩阵中访问或指定元素将失败。

要访问其线性索引大于 intmax 的元素，请使用数组索引：

尽管在稀疏矩阵中进行索引以更改单个元素的成本可忽略不计，但该成本在循环环境下会增加，而且在大型矩阵中该操作可能会使执行速度变得很慢。因此，在需要更改大量稀疏矩阵元素的情况下，最好使用向量化方法而不要使用循环方法来执行该操作。例如，考虑稀疏单位矩阵：

如果为稀疏矩阵预分配内存，然后以类似的逐个元素的方式填充，会使对稀疏数组进行索引产生大量开销：

构建索引和值向量则无需为稀疏数组进行索引，因此这种方法的速度快得多：

因此，最好使用构造函数（例如 sparse 或 spdiags 函数）一次构造所有稀疏矩阵。

例如，假定需要稀疏形式的坐标矩阵 C：

C   =  (4000−10400−10040−101010141−14)

使用 sparse 函数，以及行下标、列下标和值组成的三联对组，直接构造该五列矩阵：

以图的形式查看非零元素在稀疏矩阵内的分布通常很有用。MATLAB spy 函数生成稀疏结构的模板视图，其中图表中的每点代表一个非零数组元素的位置。

例如：

加载提供的稀疏矩阵 west0479，该矩阵是 Harwell-Boeing 集合之一。

查看稀疏结构体。

sparse