r - 做一个 block 矩阵0s

Question

我有一个大矩阵 (2160x2160)，首先我想定义一个分块矩阵(比如 20x20)并用 0 填充对角分块矩阵。我正在为如何实现这一目标而苦苦挣扎。

例如:

From the matrix

         [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9]
  [1,]    1      2     3     4     5     6     7     8     9
  [2,]    11     12    13    14    15    16    17    18    19
  [3,]    21     22    23    24    25    26    27    28    29
  [4,]    31     32    33    34    35    36    37    38    39
  [5,]    41     42    43    44    45    46    47    48    49
  [6,]    51     52    53    54    55    56    57    58    59
  [7,]    61     62    63    64    65    66    67    68    69
  [8,]    71     72    73    74    75    76    77    78    79
  [9,]    81     82    83    84    85    86    87    88    89

我首先想将 block 矩阵定义为 3x3，然后用零填充原始矩阵的对角线 block ，我最终预期的示例结果是:

         [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9]
  [1,]    0      0     0     4     5     6     7     8     9
  [2,]    0      0     0     14    15    16    17    18    19
  [3,]    0      0     0     24    25    26    27    28    29
  [4,]    31     32    33     0     0     0    37    38    39
  [5,]    41     42    43     0     0     0    47    48    49
  [6,]    51     52    53     0     0     0    57    58    59
  [7,]    61     62    63    64    65    66     0     0     0
  [8,]    71     72    73    74    75    76     0     0     0
  [9,]    81     82    83    84    85    86     0     0     0

我已经能够通过 blockmatrix 库的函数 as.blockmatrix 来定义 block 矩阵。我想我应该打破子矩阵中的 block ，用 0 填充我需要的 block ，然后重新构建我的原始矩阵，改变对角线，但我很难做到这一点。也欢迎任何其他实现此目的的想法(在 blockmatrix 库的潜在解决方案之外)

我还查看了 Matrix 包中的 bdiag 函数以及 magic 下的 adiag 函数，但没有任何想法。

预先感谢您的帮助!

Answer 1

设 M 为您的矩阵，零 block 的大小为 n x n。然后

M[(col(M) - 1) %/% n - (row(M) - 1) %/% n == 0] <- 0

给出结果，其中%/%表示整除。使用 (col(M) - 1) %/% n - (row(M) - 1) %/% n 给出一个 block 矩阵，其中对角线上的 block 仅包含零。