spdiags函数是提取和创建稀疏带和对角矩阵的函数. 用法:[B,D ] = spdiags(A),从矩阵A中取出所有非零对角元素,并保存在矩阵B中,向量D表示非零元素的对角线位置.B = spdiags(A,D ),从矩阵A中取出由D指定的对角线元素,并保存在矩阵B中.A = spdiags(B,D ,m, n),产生一个m×n稀疏矩阵A,其元素是B中的列元素放在由D指定的对角线位置上.例如1:A = [11 0 13 00 22 0 240 0 33 041 0 0 440 52 0 00 0 63 00 0 0 74]这里 m = 7, n = 4, p = 3.[B,d] = spdiags(A)B = 41 11 052 22 063 33 1374 44 24d = -302例如2:n=3; e = ones(n,1);A = spdiags([e -2*e e], -1:1, n, n)A = (1,1) -2(2,1) 1(1,2) 1(2,2) -2(3,2) 1(2,3) 1(3,3) -2例如3:spdiags(B,d,7,4)ans = (1,1) 11(4,1) 41(2,2) 22(5,2) 52(1,3) 13(3,3) 33(6,3) 63(2,4) 24(4,4) 44(7,4) 74.用法:[B,D ] = spdiags(A),从矩阵A中取出所有非零对角元素,并保存在矩阵B中,向量D表示非零元素的对角线位置.B = spdiags(A,D ),从矩阵A中取出由D指定的对角线元素,并保存在矩阵B中.A = spdiags(B,D ,m, n),产生一个m×n稀疏矩阵A,其元素是B中的列元素放在由D指定的对角线位置上.例如1:A = [11 0 13 00 22 0 240 0 33 041 0 0 440 52 0 00 0 63 00 0 0 74]这里 m = 7, n = 4, p = 3.[B,d] = spdiags(A)B = 41 11 052 22 063 33 1374 44 24d = -302例如2:n=3; e = ones(n,1);A = spdiags([e -2*e e], -1:1, n, n)A = (1,1) -2(2,1) 1(1,2) 1(2,2) -2(3,2) 1(2,3) 1(3,3) -2例如3:spdiags(B,d,7,4)ans = (1,1) 11(4,1) 41(2,2) 22(5,2) 52(1,3) 13(3,3) 33(6,3) 63(2,4) 24(4,4) 44(7,4) 74