Как распределить записи матрицы в большую с использованием продуктов матрицы

Учитывая матрицу A, где A есть

x_11  x_12  x_13 
x_21  x_22  x_23
x_31  x_32  x_33
x_41  x_42  x_43

Как я могу эффективно создать эту вторую матрицу, используя продукты матрицы?

x_11    0   0   x_21    0   0   x_31    0   0   x_41    0   0
0   x_12    0   0   x_22    0   0   x_32    0   0   x_42    0
0   0   x_13    0   0   x_23    0   0   x_33    0   0   x_43

2 ответа

Решение

Скажите, что ваша оригинальная матрица m:

m <- matrix(1:12, nrow=4)

Теперь немного жонглирования row а также col:

m2 <- matrix(0, nrow=3, ncol=length(m))
m2[((col(m2) - 1) %% 3 + 1) == row(m2)] <- m
m2

Вот один из способов:

m <- as.matrix(read.table(text='11  12  13 
21  22  23
31  32  33
41  42  43'))

do.call(cbind, lapply(split(m, seq(nrow(m))), diag))
#      [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12]
# [1,]   11    0    0   21    0    0   31    0    0    41     0     0
# [2,]    0   12    0    0   22    0    0   32    0     0    42     0
# [3,]    0    0   13    0    0   23    0    0   33     0     0    43
Другие вопросы по тегам