A=randi(500,10);%数据矩阵，这里用随机的，换你自己的数据[a b]=meshgrid(1:numel(A));B=find(sum(A(a)==A(b))>1);[i j]=ind2sub(size(A),B');disp([i j]);
