对角矩阵的压缩存储
什么是对角矩阵
-
定义
若一个n阶方阵A满足其所有非零元素都集中在以主对角为中心的带状区域中,则称其为n阶对角矩阵(diagonal matrix)。 -
图片解释
矩阵的压缩
1,当带宽b=1时
只有一条带子,像上面那个图一样,第一行(列)最后一行(列)都只有两个元素,
按照行存储 也就是一行一行的存储
a[i][j]=2(i-1)+j
2,当b不等于1.且b小于n/2
-
当1<=i<=b时,前几行个数是递增的b+1到2*b+1
1/2*(i-1)*(i+2*b)+j
-
当b+1<=i<=n-b ,,每行个数相同均为2*b+1
k=2*h*i+j-1/2*h*(h+3)
-
当n-h+1<=i<=n,个数递减从2*b+1到b+1
k=1/2*(n-i+1)(i-n-2*h)+j+h*(2*n-h)