Creating two random matrices

Question

loukil sana 2016-11-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/311289-creating-two-random-matrices

编辑： James Tursa 2016-11-8

Capture d'écran 2016-11-08 22.23.37.png

Hi. I have 2 matrices. M1(4,6) and M2(4,2). The generation of the two matrices is random. Attached you can find an explanation of the constraints. How can I formulate that with an ILP? Thanks

4 个评论
显示 2更早的评论隐藏 2更早的评论

James Tursa 2016-11-8

Does this mean the individual entries in M1 and M2 need to be integers?

loukil sana 2016-11-8

yes. https://www.mathworks.com/videos/mixed-integer-linear-programming-in-matlab-91541.html

请先登录，再进行评论。

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

Answer 1

James Tursa 2016-11-8

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/311289-creating-two-random-matrices#answer_242529

编辑：James Tursa 2016-11-8

在 MATLAB Online 中打开

E.g.,

% random matrix
M = rand(4,8);
% normalize per the 4-element sums
M(:,1:4) = bsxfun(@rdivide,M(:,1:4),sum(M(:,1:4),2));
M(:,5:8) = bsxfun(@rdivide,M(:,5:8),sum(M(:,5:8),2));
% scale per the desired sum
M(1:2,1:4) = M(1:2,1:4) * 0.50 * 2000;
M(1:2,5:8) = M(1:2,5:8) * 0.50 * 2000;
M(3:4,1:4) = M(3:4,1:4) * 0.70 * 2000;
M(3:4,5:8) = M(3:4,5:8) * 0.30 * 2000;
% divy up the results
M1 = M(:,[1:3 5:7]);
M2 = M(:,[4 8]);

EDIT: For required integer entries, you can use this ad-hoc modification (not quite exactly uniform since the fractional parts always go into M2, but maybe close enough for your purposes)

% random matrix
M = rand(4,8);
% normalize per the 4-element sums
M(:,1:4) = bsxfun(@rdivide,M(:,1:4),sum(M(:,1:4),2));
M(:,5:8) = bsxfun(@rdivide,M(:,5:8),sum(M(:,5:8),2));
% scale per the desired sum
M(1:2,1:4) = M(1:2,1:4) * 0.50 * 2000;
M(1:2,5:8) = M(1:2,5:8) * 0.50 * 2000;
M(3:4,1:4) = M(3:4,1:4) * 0.70 * 2000;
M(3:4,5:8) = M(3:4,5:8) * 0.30 * 2000;
% Adjust values to integers
F = floor(M);
D = M - F;
F(:,4) = F(:,4) + round(sum(D(:,1:4),2));
F(:,8) = F(:,8) + round(sum(D(:,5:8),2));
% divy up the results
M1 = F(:,[1:3 5:7]);
M2 = F(:,[4 8]);