第2章MATLAB数值计算.ppt

2.4.1 稀疏矩阵的建立 1．使用sparse函数产生稀疏矩阵 sparse函数用于创建稀疏矩阵，或将1个全元素矩阵直接转换成为稀疏矩阵。 语法： sparse(i,j,s,m,n) %直接创建稀疏矩阵 sparse(p) %由全元素矩阵P转换为稀疏矩阵 说明：i、j是非0元素的行、列下标；s是非0元素所形成的向量；m、n是s的行、列维数，可省略；i、j、s都是长度相同的向量，生成矩阵的元素s(k)下标分别是i(k)和j(k)；P为全元素矩阵。 【例2.21】 产生稀疏矩阵。 a=eye(3); a(4,:)=[ ?5?2 ?3] a = 1 0 0 0 1 0 0 0 1 ?5 ?2 ?3 b=sparse(a) %创建稀疏矩阵 b = (1,1) 1 (4,1) ?5 (2,2) 1 (4,2) ?2 (3,3) 1 (4,3) ?3 c=sparse([1 4 2 4 3 4],[1 1 2 2 3 3],[1 ?5 1 ?2 1 ?3]); %创建与b相同的稀疏矩阵 程序说明：sparse的前两个参数向量分别表示稀疏矩阵元素的行和列下标，第一个元素的下标是(1,1)，第三个参数向量是稀疏矩阵元素，第一个元素1。 与sparse函数相反，full函数可将稀疏矩阵转变为全元素矩阵。 语法： full(p) %将稀疏矩阵P转变为全元素矩阵 2．用spdiags函数创建稀疏矩阵 spdiags函数用对角线元素构建1个稀疏矩阵。 语法： spdiags(D,k,m,n) 说明：矩阵D的每一列代表矩阵的对角线向量；k代表对角线的位置（0代表主对角线，?1代表向下位移一单位的次对角线，1代表向上位移一单位的次对角线，依此类推）；m、n分别代表矩阵的行、列维数。 【例2.21续】 用spdiags函数创建稀疏矩阵。 D=[3 2 9;2 4 9;1 1 4] D = 3 2 9 2 4 9 1 1 4 d=[0 1 2]; s=spdiags(D,d,4,3) %构成4行3列的稀疏矩阵 s = (1,1) 3 (1,2) 4 (2,2) 2 (1,3) 4 (2,3) 1 (3,3) 1 full(s) %转换成全元素矩阵 ans = 3 4 4 0 2 1 0 0 1 0 0 0 可以使用“whos”命令比较【例2.21】矩阵A和B所占用的内存大小。 whos Name Size Byte Class a 4?3 96 double array b 4?3 88 sparse array Grand total is 18 elements using 184 Byte 可以看出稀疏矩阵B占用的内存为88字节，比全元素矩阵A占用的96少。如果稀疏矩阵B含非0元素更少则占用内存字节数更少。 对于1个只包含实数m?n的稀疏矩阵，含有nnz个非零元素，MATLAB使用3个内部数组储存此稀疏矩阵的信息。 （1）第1个数组：以double方式储存nnz个非零元素，使用的空间为8?nnz（Byte）。 （2）第2个数组：以整数方式储存nnz个非零元素，每个元素的行下标使用的空间为4?nnz（Btye）。 （3）第3个数组：以整数方式储存n个列的每个列的起始指针，使用的空间为4?n（Byte）。 （4）如果是复数稀疏矩阵，则需要第4个数组，以double方式储存nnz个非零元素的虚数部分。 MATLAB提供了如下几个返回稀疏矩阵元素个数的函数。 ① nnz：可返回稀疏矩阵的非零元素个数。 ② nonzeros：返回1个包含所有非零元素的列向量。 ③ nzmax：返回最大的非零元素个数，当nnznzmax时，MATLAB会动态调整以便增加内存给nzmax，用于储存新增的非零元素。 【例2.21续】 查看稀疏矩阵的非零元素。 nnz(b) %得出非零元素个数 ans = 6 nonz