对于任意大小的矩阵x
,如何在给定矩阵的每一行中找到最后一个非零元素的索引?
例如,对于矩阵
x = [ 0 9 7 0 0 0; 5 0 0 6 0 3; 0 0 0 0 0 0; 8 0 4 2 1 0 ]
应该获得向量[ 3 6 0 5 ]
。
答案 0 :(得分:10)
这是一个较短的版本,结合了find和accumarray
x = [ 0 9 7 0 0 0; 5 0 0 6 0 3; 0 0 0 0 0 0; 8 0 4 2 1 0 ];
%# get the row and column indices for x
[rowIdx,colIdx] = find(x);
%# with accumarray take the maximum column index for every row
v = accumarray(rowIdx,colIdx,[],@max)'
v =
3 6 0 5
答案 1 :(得分:4)
这是一个版本:
x = [ 0 9 7 0 0 0; 5 0 0 6 0 3; 0 0 0 0 0 0; 8 0 4 2 1 0 ];
c = arrayfun(@(k) find(x(k,:)~=0,1,'last'), 1:size(x,1), 'UniformOutput',false);
c( cellfun(@isempty,c) ) = {0};
v = cell2mat(c);
v =
3 6 0 5
修改强>: 考虑这个替代解决方案:
[m,v] = max( cumsum(x'~=0) );
v(m==0) = 0;
v =
3 6 0 5
答案 2 :(得分:2)
bsxfun
的单行解决方案:
result = max(bsxfun(@times, x~=0, 1:size(x,2)).');
或者使用max
的两个输出:
[val, result] = max(fliplr(x~=0).',[],1); %'
result = (size(A,2)+1-result).*val;
答案 3 :(得分:1)
我的回答有点扭曲,但也应该有效
x = [ 0 9 7 0 0 0; 5 0 0 6 0 3; 0 0 0 0 0 0; 8 0 4 2 1 0 ];
[~,pos] = max([fliplr(x~=0),ones(size(x,1))],[],2);
v = size(x,2)-pos' +1;