如何结合两个矩阵在R中创建块对角矩阵？
块对角矩阵表示在最后一个元素末尾将一个矩阵添加到另一个矩阵。例如，如果我们有一个具有九个值的矩阵，另一个矩阵也具有九个值，则第二个矩阵将被添加到第一个矩阵，并且第一个矩阵下方的元素将为零，而第二个矩阵上方的元素也将为零。
示例
M1<-matrix(1:9,ncol=3) M1
输出结果
[,1] [,2] [,3] [1,] 1 4 7 [2,] 2 5 8 [3,] 3 6 9
示例
M2<-matrix(10:18,ncol=3) M2
输出结果
[,1] [,2] [,3] [1,] 10 13 16 [2,] 11 14 17 [3,] 12 15 18
示例
rbind(cbind(M1,matrix(0,nrow=nrow(M1),ncol=ncol(M2))),cbind(matrix(0,nrow=nrow(M2),ncol=ncol(M1)),M2))
输出结果
[,1] [,2] [,3] [,4] [,5] [,6] [1,] 1 4 7 0 0 0 [2,] 2 5 8 0 0 0 [3,] 3 6 9 0 0 0 [4,] 0 0 0 10 13 16 [5,] 0 0 0 11 14 17 [6,] 0 0 0 12 15 18
示例
M3<-matrix(1:25,nrow=5) M3
输出结果
[,1] [,2] [,3] [,4] [,5] [1,] 1 6 11 16 21 [2,] 2 7 12 17 22 [3,] 3 8 13 18 23 [4,] 4 9 14 19 24 [5,] 5 10 15 20 25
示例
M4<-matrix(26:50,ncol=5) M4
输出结果
[,1] [,2] [,3] [,4] [,5] [1,] 26 31 36 41 46 [2,] 27 32 37 42 47 [3,] 28 33 38 43 48 [4,] 29 34 39 44 49 [5,] 30 35 40 45 50
示例
rbind(cbind(M3,matrix(0,nrow=nrow(M4),ncol=ncol(M4))),cbind(matrix(0,nrow=nrow(M4),ncol=ncol(M3)),M4))
输出结果
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 1 6 11 16 21 0 0 0 0 0 [2,] 2 7 12 17 22 0 0 0 0 0 [3,] 3 8 13 18 23 0 0 0 0 0 [4,] 4 9 14 19 24 0 0 0 0 0 [5,] 5 10 15 20 25 0 0 0 0 0 [6,] 0 0 0 0 0 26 31 36 41 46 [7,] 0 0 0 0 0 27 32 37 42 47 [8,] 0 0 0 0 0 28 33 38 43 48 [9,] 0 0 0 0 0 29 34 39 44 49 [10,] 0 0 0 0 0 30 35 40 45 50
示例
M5<-matrix(sample(1:5,25,replace=TRUE),nrow=5)M5
输出结果
[,1] [,2] [,3] [,4] [,5] [1,] 2 2 3 1 2 [2,] 2 4 1 5 4 [3,] 4 1 5 4 2 [4,] 2 4 4 2 5 [5,] 5 1 2 5 3
示例
M6<-matrix(sample(6:10,25,replace=TRUE),ncol=5) M6
输出结果
[,1] [,2] [,3] [,4] [,5] [1,] 10 9 9 6 6 [2,] 8 7 10 9 8 [3,] 8 9 9 8 9 [4,] 7 8 6 8 7 [5,] 8 9 6 8 7
示例
rbind(cbind(M5,matrix(0,nrow=nrow(M6),ncol=ncol(M6))),cbind(matrix(0,nrow=nrow(M6),ncol=ncol(M5)),M6))
输出结果
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 3 1 3 1 3 0 0 0 0 0 [2,] 4 2 3 1 5 0 0 0 0 0 [3,] 5 4 3 5 4 0 0 0 0 0 [4,] 1 3 5 1 5 0 0 0 0 0 [5,] 1 2 2 4 5 0 0 0 0 0 [6,] 0 0 0 0 0 8 8 7 6 10 [7,] 0 0 0 0 0 6 6 8 10 7 [8,] 0 0 0 0 0 9 6 7 10 9 [9,] 0 0 0 0 0 10 6 10 6 7 [10,] 0 0 0 0 0 10 10 9 7 6
示例
M7<-matrix(round(rnorm(36),0),nrow=6) M7
输出结果
[,1] [,2] [,3] [,4] [,5] [,6] [1,] 1 0 2 -1 1 0 [2,] -2 0 2 0 0 -1 [3,] -1 0 0 0 0 0 [4,] 0 0 1 1 1 0 [5,] 1 -1 0 -1 2 1 [6,] 0 1 0 0 2 -1
示例
M8<-matrix(round(runif(36,2,5),0),nrow=6) M8
输出结果
[,1] [,2] [,3] [,4] [,5] [,6] [1,] 3 4 3 3 3 4 [2,] 2 2 2 3 4 5 [3,] 3 2 3 4 4 3 [4,] 3 3 5 2 3 2 [5,] 3 5 4 4 5 2 [6,] 4 2 2 4 2 3
示例
rbind(cbind(M7,matrix(0,nrow=nrow(M8),ncol=ncol(M8))),cbind(matrix(0,nrow=nrow(M8),ncol=ncol(M7)),M8))
输出结果
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 0 0 0 1 -1 1 0 0 0 0 0 0 [2,] -1 1 2 1 -1 1 0 0 0 0 0 0 [3,] 0 0 1 0 1 1 0 0 0 0 0 0 [4,] -1 0 -2 0 1 -3 0 0 0 0 0 0 [5,] 0 1 1 0 1 0 0 0 0 0 0 0 [6,] 0 0 0 0 0 1 0 0 0 0 0 0 [7,] 0 0 0 0 0 0 2 4 2 2 3 4 [8,] 0 0 0 0 0 0 4 4 5 3 3 4 [9,] 0 0 0 0 0 0 3 5 3 4 5 3 [10,] 0 0 0 0 0 0 3 4 5 4 2 5 [11,] 0 0 0 0 0 0 4 3 3 3 4 2 [12,] 0 0 0 0 0 0 5 3 5 2 3 5