MATLAB 5手册009.PDFVIP

下载 第9章 稀 疏 矩 阵 在许多问题中提到了含有大量 0元素的矩阵，这样的矩阵称为稀疏矩阵。比如求解普通或 者部分微分方程的数值解。为了节省存储空间和计算时间， M AT L A B考虑到矩阵的稀疏性， 在对它运算时有特殊的命令。 9.1 矩阵为什么稀疏 一个稀疏矩阵中有许多元素等于零，这便于矩阵的计算和保存。如果 M AT L A B把一个矩 阵当作稀疏矩阵，那么只需在 m ×3 的矩阵中存储m个非零项。第 1列是行下标，第 2 列是列下 标，第 3列是非零元素值，不必保存零元素。如果存储一个浮点数要 8个字节，存储每个下标 要4个字节，那么整个矩阵在内存中存储需要 1 6 ×m个字节。 ■ 例9 . 1 A = e y e ( 1 0 0 0 ) ; 得到一个 1 0 0 0 ×1 0 0 0的单位矩阵，存储它需要8 Mb 空间。如果使用命令： B = s p e y e ( 1 0 0 0 ) ; 用一个 1 0 0 0 ×3 的矩阵来代表，每行包含有一个行下标、列下标和元素本身。现在只需 1 6 K b的空间就可以存储 1 0 0 0 ×1 0 0 0的单位矩阵，它只需要满单位矩阵的 0 . 2 %存储空间。对于 许多的广义矩阵也可这样来作。 ■ 稀疏矩阵的计算速度更快，因为 M AT L A B 只对非零元素进行操作，这是稀疏矩阵的第二 个突出的优点。 ■ 例9 . 2 假设矩阵A 、B和例9 . 1 中的矩阵一样。计算2 *A 需要一百万次的浮点运算，而计算 2 *B 只 需要2 0 0 0次浮点运算。 ■ 因为M AT L A B 不能自动创建稀疏矩阵，所以要用特殊的命令来得到稀疏矩阵，在下一节 中将给出这些命令。前面章节中的算术和逻辑运算都适用于稀疏矩阵。 9.2 创建和转换稀疏矩阵 在M AT L A B 中，用命令s p a r s e来创建一个稀疏矩阵。 命令集8 7 创建稀疏矩阵 s p a r s e ( A ) 由非零元素和下标建立稀疏矩阵 A 。如果A 已是一个稀疏 矩阵，则返回A本身。 s p a r s e ( m , n ) 生成一个m ×n 的所有元素都是0 的稀疏矩阵。 1 2 4 M A TLAB 5 手册 下载 s p a r s e ( u , v , a ) 生成一个由长度相同的向量 u ，v和a定义的稀疏矩阵。其 中u 和v 是整数向量， a是一个实数或者复数向量。 (u , v ) i i 对应值 a ，如果a 中有零元素，则将这个元素排除在外。 i 稀疏矩阵的大小为m a x (u) ×m a x (v) 。 s p a r s e ( u , v , a , m , n ) 生成一个m ×n 的稀疏矩阵，(u , v )对应值a 。向量u ，v和a