Hi all,
I want to find the indexes of the maximum and the second maximum in 3D array for a certain dimension. My strategy is to use the MAX function for finding the indexes of the first maximum. Use these indexes to make the value of the 3D array –Inf and call again the MAX function in order to find the indexes of the second maximum of the 3D array
My question is how to make optimal use of the returned indexes from the MAX function (without for loops).
cube1 = 10*rand([5 6 7]); %create 3D array
cube2 = cube1; %make copy for com parison
[dummy i1st] = max(cube1,[],3); %find first max across 3th dimension
cube1(i1st)=-inf; % this wont work
% for loop does work but is not elegant, require recoding when dimension in changed
for r=1:size(i1st,1) for c=1:size(i1st,2) cube1(r,c,i1st(r,c)) = -inf; endend[dummy i2nd] = max(cube1,[],3); %find second max across 3th dimension
Does anyone has some ideas how to get rid of the for loop?
Patrick
Best Answer