Диаграмма извлечения блока из матрицы

Question

У меня есть матрица njxnj, состоящая из nxn-матриц. Я хочу извлечь диагональные j блоков из nxn-матриц. то есть я хочу, чтобы извлечь диагональ (при п = 2, J = 4):

Что бы быть наиболее эффективный способ сделать это?

Answer 1

Чтобы проиндексировать элементы, вы можете использовать blkdiag, чтобы создать соответствующую маску.

%your parameters 
n=2 
j=4 
%some example matrix 
M=magic(n*j); 
%create the input for blkdiag, j matrices of size n 
h=repmat({true(n)},j,1) 
%use blkdiag to select the elements 
M(logical(blkdiag(h{:})))