假设我们有一个订单A
的矩阵k2*k
,其中k=k1+k2
为k1=k-k2
。矩阵A
的特征是第一k1
列为零,剩余k2*k2
为单位矩阵。如何在R中创建这样的自定义矩阵?
注意:
对于k1
和k2
的较小尺寸,这很容易。但我正在寻找一个自动命令来处理更大尺寸的k1
和k2
。
答案 0 :(得分:2)
k <- 20
k2 <- 15
k1 <- k - k2
diagonal <- diag(k2)
zeros <- matrix(0, nrow = k2, ncol = k1)
result <- cbind(zeros, diagonal)