How do I create a matrix with the diagonals set to 10 and non-diagonals a random integer (0 to 2) where each row sums to 10 (without including diagonal)?

Question

Josh 2016-2-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/266330-how-do-i-create-a-matrix-with-the-diagonals-set-to-10-and-non-diagonals-a-random-integer-0-to-2-wh

评论： Star Strider 2016-2-2

I need to create a 15x15 matrix where the diagonals are all 10 and the numbers off the diagonal are a random integer between 0 and 2. I also need that non-diagonal numbers to sum to 10 for each row. See the below matrix, first row for an example.

M =

   1     0     0     2     0     2     1     1     0     2     0     0     1     0
  10     0     0     0     0     0     0     0     0     0     0     0     0     0
   0    10     0     0     0     0     0     0     0     0     0     0     0     0
   0     0    10     0     0     0     0     0     0     0     0     0     0     0
   0     0     0    10     0     0     0     0     0     0     0     0     0     0
   0     0     0     0    10     0     0     0     0     0     0     0     0     0
   0     0     0     0     0    10     0     0     0     0     0     0     0     0
   0     0     0     0     0     0    10     0     0     0     0     0     0     0
   0     0     0     0     0     0     0    10     0     0     0     0     0     0
   0     0     0     0     0     0     0     0    10     0     0     0     0     0
   0     0     0     0     0     0     0     0     0    10     0     0     0     0
   0     0     0     0     0     0     0     0     0     0    10     0     0     0
   0     0     0     0     0     0     0     0     0     0     0    10     0     0
   0     0     0     0     0     0     0     0     0     0     0     0    10     0
   0     0     0     0     0     0     0     0     0     0     0     0     0    10

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Kirby Fears 2016-2-2

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/266330-how-do-i-create-a-matrix-with-the-diagonals-set-to-10-and-non-diagonals-a-random-integer-0-to-2-wh#answer_208283

编辑：Kirby Fears 2016-2-2

在 MATLAB Online 中打开

Here's a solution where you can change the values of matrix size, diagonal value, and non diagonal sum. Values between 0 and 2 are drawn to fill the rows. If they sum to the desired value, the row is filled. Otherwise they are redrawn, which ensures randomness.

% settings
matSize = 15;
diagVal = 10;
nonDiagSum = 10;
% Initialize matrix with diagonal filled in
targetMatrix = zeros(matSize);
targetMatrix(1:(matSize+1):end) = diagVal;
% Fill in matrix with random row that sums to nonDiagSum
i = 1;
while i<=matSize,
    randValues = randi(3,1,matSize-1)-1;
    if (sum(randValues)==nonDiagSum),
        idx = true(1,matSize);
        idx(i) = false;
        targetMatrix(i,idx) = randValues;
        i = i+1;
     end
end

This executes quickly for your example, but it could take an exceedingly long time to execute if matSize is very different from nonDiagSum.

Hope this helps.

2 个评论
显示无隐藏无

Josh 2016-2-2

This works great! Thank you!!

Star Strider 2016-2-2

@Kirby — Write your code up as a function and contribute it to the File Exchange. I did with a similar Answer that I thought was a one-off, and it's already gotten a bunch of downloads. I was quite pleasantly surprised.

请先登录，再进行评论。

Answer 2

John D'Errico 2016-2-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/266330-how-do-i-create-a-matrix-with-the-diagonals-set-to-10-and-non-diagonals-a-random-integer-0-to-2-wh#answer_208282

在 MATLAB Online 中打开

Simple.

First, generate all possible partitions of the number 10, as a sum of the integer [0,1,2]. I'll use my partitions function to see how to do it.

partitions(10,[0 1 2],[],10)
ans =
  10     0
   8     1
   6     2
   4     3
   2     4
   0     5

Thus, we can write 10 as

ones
ones and a 2
ones, and a pair of twos
ones and 3 twos
ones and 4 twos
twos

Your matrix has 9 columns in it, if we exclude the diagonal, so we can exclude the case where we have 10 ones.

Once you have generated all possible rows, just select randomly from that set to fill each row.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

the cyclist 2016-2-2

编辑：the cyclist 2016-2-2

在 MATLAB Online 中打开

John, that syntax of partitions() gives the error

Error using partitions (line 143)
All members of candidate_set must be > 0

for me, using R2015b on a Mac. I have the latest update from here.

请先登录，再进行评论。

Answer 3

the cyclist 2016-2-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/266330-how-do-i-create-a-matrix-with-the-diagonals-set-to-10-and-non-diagonals-a-random-integer-0-to-2-wh#answer_208285

在 MATLAB Online 中打开

matSize = 10;
% Partitions of matSize; Use's John D'Errico's utility from FEX
p = partitions(matSize);
% Exclude those that have values 3-10 in them.
p = p(sum(p(:,3:end),2)==0,:);
% Exclude the partition that is ten 1's, because there are only 9 spots to fill
p(1,:) = [];
% Remove the spurious columns
p(:,all(p==0)) = [];
% Append the zero counts
p = [matSize-1 - sum(p,2), p];
numberRows = size(p,1);
% For each partition, construct all the permutations
allPermutations = [];
for nr = 1:numberRows
    basis = [zeros(1,p(nr,1)) ones(1,p(nr,2)) 2*ones(1,p(nr,3))];
    thesePermutations = unique(perms(basis),'rows');
    allPermutations = [allPermutations; thesePermutations];
end
numberAllPermutations = size(allPermutations,1);
M = matSize * eye(matSize);
for nr = 1:matSize,
    rowRand = allPermutations(randi(numberAllPermutations),:);
    M(nr,1:nr-1) = rowRand(1:nr-1);
    M(nr,nr+1:end) = rowRand(nr:end);
end

2 个评论
显示无隐藏无

the cyclist 2016-2-2

编辑：the cyclist 2016-2-2

This is vastly slower than Kirby's solutions, probably because I am lazily generating many more partitions than required. Combining John's own more efficient use of partitions() in his answer would be much better. Unfortunately, I could not get his syntax to work.

the cyclist 2016-2-2

Oops. Nope. I just profiled my code. The slowness is NOT the partitioning, even though I am using it dumbly. The slow part is generating the permutations.

请先登录，再进行评论。