矩阵合并

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

3*3矩阵A
2*2矩阵B
求矩阵diag（A, B）的代码？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Dia 矩阵

c(diag(A),diag(B))?

ltx5151 发表于 2013-2-26 10:01
c(diag(A),diag(B))?

感谢回复。但是不行，出来的是一行的向量。

diag(c(diag(A),diag(B)))

qoiqpwqr 发表于 2013-2-26 11:15
diag(c(diag(A),diag(B)))

感谢回复。可是还是不行。因为矩阵A,B不一定是对角阵，所以diag(c(diag(A),diag(B)))会丢掉A,B的某些元素。

子衿1219 发表于 2013-2-26 11:28
感谢回复。可是还是不行。因为矩阵A,B不一定是对角阵，所以diag(c(diag(A),diag(B)))会丢掉A,B的某些元素 ...

你想要的结果是什么？给个例子

你的意思是创建一个分块矩阵C，矩阵A、B位居这个分块矩阵的对角位置。如此的话，试试下面这个（比较笨）：

AA=rbind(A,matrix(0,2,3);
BB=rbind(matrix(0,3,2),B);
CC=cbind(AA,BB)

flamingstone 发表于 2013-2-26 15:12
你的意思是创建一个分块矩阵C，矩阵A、B位居这个分块矩阵的对角位置。如此的话，试试下面这个（比较笨）：
...

嗯。谢谢了。我现在就是用的这种“笨”方法。谢啦。

这样也行

> C <- matrix(0, 5, 5)
> C[1:3, 1:3] <- A
> C[4:5, 4:5] <- B

rbind(cbind(a,matrix(0,nrow(a),ncol(b))),  cbind(matrix(0,nrow(b),ncol(a)),b))