2011-07-04から1日間の記事一覧

ブロック対角化に関するメモ

R

RjpWikiの記事(こちら) L <- list(matrix(seq(4),4,4),matrix(seq(9),3,3),matrix(seq(25),5,5)) library(Matrix) bdiag(L) M<- as.matrix(bdiag(L)) N<-length(M[,1]) # 対角状態を崩す M2<-M[sample(1:N),sample(1:N)] M3<-sign(M2) par(ask=TRUE) heatmap…