Add an extra loop to create a 3D matrix:
N = 5; % No. of matrices
P = 5; % The number of rows
T = 10; % The number of columns
% Duration of projects
D = [4; 5; 3; 2; 4];
excludedcount = D-1;
X = zeros(P,T);
FinalArray = nan(N,P,T) ; % pre-allocation
for MatrixCount = 1:5
for row = 1:P
rv = [1, zeros(1, T - excludedcount(row) - 1)];
X(row, 1 : (T - excludedcount(row))) = rv(randperm(numel(rv)));
end
FinalArray(MatrixCount,:,:) = X ; % store it
end
