matlab

1．带状（对角）稀疏矩阵

函数 spdiags

格式 [B,d] = spdiags(A) %从矩阵A中提取所有非零对角元素，这些元素保存在矩阵B中，向量d表示非零元素的对角线位置。

B = spdiags(A,d)     %从A中提取由d指定的对角线元素，并存放在B中。

A = spdiags(B,d,A) %用B中的列替换A中由d指定的对角线元素，输出稀疏矩阵。

A = spdiags(B,d,m,n) %产生一个m×n稀疏矩阵A，其元素是B中的列元素放

在由d指定的对角线位置上。

例1-100

>>A = [11   0   13    0

       0   22    0   24

         0   0   33    0

      41    0    0   44

     0   52    0    0

        0    0   63    0

        0   0    0   74];

>>[B,d] = spdiags(A)

B =

    41    11     0

    52    22    0

    63    33    13

    74    44    24

d =

    -3      %表示B的第1列元素在A中主对角线下方第3条对角线上

     0      %表示B的第2列在A的主对角线上

     2      %表示B的第3列在A的主对角线上方第2条对角线上

例1-101

>> B=[1 2 3 4

      5 6 7 8

      9 10 11 12

      13 14 15 16];

>> d=[-2 0 1 3];

>> A=spdiags(B,d,4,4);

>> full(A)

ans =

     2     7     0    16

     0     6    11     0

     1     0    10    15

     0     5     0    14