在r

Question

我想形成一个100x100的块对角矩阵，其对角线实体为5，（k = 1,2，（（2k-1）2k）和（2k（2k-1））实体与0.5相同。 ..50，所有其他实体均为零。

Answer 1

使用base R，您可以使用kronecker()创建块对角矩阵，即

M <- kronecker(diag(50),matrix(c(5,0.5,0.5,5),2,2))

Answer 2

> library(Matrix)
> M <- as.matrix(bdiag(lapply(seq_len(50), function(X) matrix(c(5, 0.5, 0.5, 5), 2))))
> dim(M)