尝试这个:
x = unique(perms([0 0 2 4]),'rows'); %# we'll convert 2 to 10 and 4 to 100 later
xi = randperm(size(x,1)); %# permute rows randomly
n = 6;
xi = xi(1:n); %# k random rows
x = x(xi,:);
x = [ x ones(n,1) ]'; %# add 1s and transpose
%# convert to strings representing binary numbers, then to cell array
out = reshape(cellstr(dec2bin(x)), size(x,1), n);
%# split the cell array by columns
out = mat2cell(out, size(out,1), ones(1,n));
如果您实际上需要数字(0、10、100),而不是字符串,则代码会短一些:
x = unique(perms([0 0 10 100]),'rows');
xi = randperm(size(x,1));
n = 6;
xi = xi(1:n);
x = x(xi,:);
x = [x ones(n,1)]';
out = mat2cell(x,size(x,1),ones(1,n));
更新:
根据您的评论,这是新代码:
n = 6;
P = repmat({zeros(5,3)},1,n); %# output array preallocation
x = unique(perms([0 0 1 2]),'rows');
xi = randperm(size(x,1));
xi = xi(1:n);
x = x(xi,:);
c = [x ones(n,1)+2]'; %# column index
r = repmat((1:5)',1,n); %# row index
for k=1:n
P{k}(sub2ind( [5 3], r(c(:,k)>0,k), c(c(:,k)>0,k) )) = 1;
end
如果您需要逻辑矩阵数组,只需在第二行替换zeros
为:false
P = repmat({false(5,3)},1,n);