Find sum of max values in matrix subject to constraints

Question

Seamus 2015-7-15

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/229927-find-sum-of-max-values-in-matrix-subject-to-constraints

评论： bio lim 2015-7-17

Hi I would greatly appreciate help on the following question. I have two matrices, a and b: a = [0 0 3 ; 0 7 2 ; 5 0 1 ] ; , with a representing yields; and b = [1 2 3 ; 1 2 100 ; 1 100 100 ] ; representing tenors.

I would like to find the highest combined value from a, subject to the sum of the coinciding values in b not exceeding three.

So the answer here would be 12, because the sum of 7+5 = 12, and the coinciding sum of b values would be 2+1 = 3.

Like i say any help at all would be greatly appreciated.

2 个评论
显示无隐藏无

Grzegorz Lippe 2015-7-15

编辑：Grzegorz Lippe 2015-7-15

在 MATLAB Online 中打开

Hello Seamus,

a nice riddle :)

A = [0 0 3 ; 0 7 2 ; 5 0 1 ] ;
B = [1 2 3 ; 1 2 100 ; 1 100 100 ] ;
siz = size(A) ;
b = B(:) ;
a = A(:) ;
[b1, b2] = meshgrid(b, b) ;
[a1, a2] = meshgrid(a, a) ;
% Sum up all possiblities, but only once
AA = triu(a1) + triu(a2) ;
BB = triu(b1) + triu(b2) ;
SIZ = size(AA) ;
% The sum of B must be less 4
validIdsofBB = triu(BB < 4) ;
% Find the maximum of all combined values in A, but only valid in B
[~, myID] = max(AA(:) .* validIdsofBB(:)) ;
%%The answer to the question:
AA(myID)
%%The Indizes to the original vectors:
[I,J] = ind2sub(SIZ,myID) ;
a(I)
a(J)
%%The Indizes to the original matrix
[i1,i2] = ind2sub(siz,I) ;
[j1,j2] = ind2sub(siz,J) ;
A(i1, i2)
A(j1, j2)

Seamus 2015-7-15

Gregor, that's absolute genius. Really, I thought i would have to loop and do this iteratively looking for something that met the right conditions. Thank you very much!

请先登录，再进行评论。

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

Answer 1

bio lim 2015-7-15

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/229927-find-sum-of-max-values-in-matrix-subject-to-constraints#answer_186108

在 MATLAB Online 中打开

 a = [0 0 3 ; 0 7 2 ; 5 0 1 ] ;
 b = [1 2 3 ; 1 2 100 ; 1 100 100 ] ; 
 [first_max, ind] = max(a(:));
 [m1, n1] = ind2sub(size(a),ind);
 c = a;
 c(m1,n1) = 0;
 [second_max, ind] = max(c(:));
 [m2,n2] = ind2sub(size(c),ind);
 sum1 = a(m1,n1) + a(m2,n2)
 if b(m1,n1) < 3 & b(m2,n2) < 3
 sum2 = b(m1,n1) + b(m2,n2)
 end

2 个评论
显示无隐藏无

Seamus 2015-7-15

Hi coffee murun, thanks a million for this. I had one follow up question on this. Is there an else statement that you can write which would repeat the initial part of the code in order to satisfy the below 3 constraint?

bio lim 2015-7-17