function B=find_max(A,n)
%該函數查找任意維數矩陣中絕對值最大的 n 個元素
%矩陣A只能是n*n的方陣,n是元素絕對值最大的個數
%該函數能處理具有相同元素的矩陣A

%maked by Li_Ming_Quan on 2007-9-22 
tic                            %測定運行時間開始  
B=zeros(size(A));              %生成要存放的矩陣B
L=abs(A);                      %對矩陣A的每個元素取絕對值
M=reshape(L,1,numel(A));       %將矩陣L變為行向量
for i=1:numel(M)-1             %將相同的元素清零,只留下一個值
for j=i+1:numel(M)
if M(i)==M(j)
M(j)=0;
end
end
end
M=sort(M,'descend');           %處理后按降序排列
R=M(n);                        %取出絕對值最大n個元素的最小值
V=M(find(M>=R));               %確定絕對值最大的n個元素
for i=1:n;
[S,T]=find(V(i)==L);           %取出這n個元素的位置
for j=1:size(S)
B(S(j),T(j))=A(S(j),T(j));     %取出絕對值最大n個元素賦給矩陣
                               %B,其余為零
end
end
toc                            %測定運行時間結束