Создать массив единиц в указанном порядке

Question

У меня есть матрица положение:

positionMatrix = [1 2 3; 1 3 2; 2 1 3]; 

Я хочу простую реализацию (не для цикла), что бы генерирует массивы следующим образом:

% there are 3 lines in positionMatrix, so it should generates 3 arrays of ones 
array 1 should be [1 0 0; 0 1 0; 0 0 1] %from positionMatrix 1 2 3 
array 2 should be [1 0 0; 0 0 1; 0 1 0] %from positionMatrix 1 3 2 
array 3 should be [0 1 0; 1 0 0; 0 0 1] %from positionMatrix 2 1 3 

positionMatrix мог - M x N (с M не равным N).

Answer 1

Генерация выходного сигнала в виде единого массива 3D

[M N] = size(positionMatrix); 
mx = max(positionMatrix(:)); % max column index 
out = zeros([N mx M]); 
out(sub2ind(size(out), ... 
     repmat(1:N, [M 1]),... 
     positionMatrix, ... 
     repmat((1:M)', [1 N]))) = 1; 
out(:,:,1) = 
1  0  0 
0  1  0 
0  0  1 
out(:,:,2) = 
1  0  0 
0  0  1 
0  1  0 
out(:,:,3) = 
0  1  0 
1  0  0 
0  0  1 

, если вы хотите, чтобы каждый выход матрицы в другую ячейку, которую можно использовать mat2cell

>> mat2cell(out, N, mx, ones(1,M)) 

Answer 2

Не уверен, что это именно то, что вам нужно, но вот что-то, что приближается к требованиям:

positionMatrix = [1 2 3; 1 3 2; 2 1 3]; 

myArray = [positionMatrix == 1, positionMatrix == 2, positionMatrix == 3] 

Это, конечно, предполагает, что ваш positionMatrix будет расти в количестве строк, но нет (много) в количестве столбцов.

Answer 3

Здесь я иду с accumarray снова. На самом деле, это один довольно интуитивное с accumarray, если учесть, что места в выводе распределены следующим образом,

• Колонка: Значения в positionMatrix.
• Строка: столбец positionMatrix.
• Кусочек: ряд в positionMatrix.

Если мы называем выходную матрицу map, это как применить accumarray:

[slices,rows] = ndgrid(1:size(positionMatrix,1),1:size(positionMatrix,2)); 
map = accumarray([rows(:) positionMatrix(:) slices(:)],ones(1,numel(rows))) 
map(:,:,1) = 
    1  0  0 
    0  1  0 
    0  0  1 
map(:,:,2) = 
    1  0  0 
    0  0  1 
    0  1  0 
map(:,:,3) = 
    0  1  0 
    1  0  0 
    0  0  1 

При необходимости, вы можете положить сторону три ломтика бок с map = reshape(map,size(map,1),[],1);.

Answer 4

Она также может быть сделано с ndgrid:

positionMatrixTr = positionMatrix.'; 
[M N] = size(positionMatrixTr); 
L = max(positionMatrixTr(:)); 
[jj kk] = ndgrid(1:M,1:N); 
array = zeros(M,L,N); 
array(sub2ind([M L N],jj(:),positionMatrixTr(:),kk(:))) = 1; 

Как и другие ответы, это дает результат в 3D-массиве.