%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   作者：王宇                  %
%   內容：置亂函數 18/06/07     %  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%matrix:要處理的圖像矩陣
%k,N:k控制參數，N矩陣大小，k,N控制周期
%count跌代次數
%chaosC變換后矩陣

function [chaosC,row,col]=chaos(matrix,k,N,count)

%[m,n]=size(matrix);
m=64;
n=64;
chaosC=zeros(N);
randmtx=randCL(10,10,0.5,3.6)
row=zeros(N);
col=zeros(N);
%row=ones(N);
%col=ones(N);

 for i=1:m
    for j=1:n
        row(i,j)=mod(i+j,N);
        col(i,j)=mod(k*i+(k+1)*j,N);
    end;
 end

 for l=1:count-1
      for i=1:m          
        for j=1:n
            p=row(i,j);
            row(i,j)=mod(row(i,j)+col(i,j),N);      
            col(i,j)=mod(k*p+(k+1)*col(i,j),N);
        end        
      end      
 end
for i=1:m
    for j=1:n
        chaosC(row(i,j),col(i,j))=matrix(i,j);
    end;
end;