dctmax.m

混沌生成的很多形式

function [Ix,Iy]=dctmax(x,long)
x=abs(x);
[m,n]=size(x);
Mcol=[];
Im1=[];
In1=[];
for i=1:n,
   Mcol=[Mcol;x(:,i)];
end
Mcol=sortrows(Mcol);
Mcut=Mcol(m*n-long+1);
[Im In]=find(x>=Mcut);
if size(Im)>long,
   [Im1 In1]=find(x==Mcut);
   lh=size(Im)-long;   
   Ix=[Im;Im1(1:lh)];
   Iy=[In;In1(1:lh)];
end
Ix=Im;
Iy=In;