Cat映射的周期性

function cattest2 % Cat 加密图象 clear all f=imread('cat.jpg'); % 读曲图象 A=double(f)/255; [m,n,l]=size(A); Cat=[321 40;8 1];N=124; m=m-rem(m,N); n=n-rem(n,N); A=A(1:m,1:n,:); subplot(2,3,1) imagesc(A,[0 1]); mk=m/N;nk=n/N; for k=1:5 % 将图象分割成124*124的数据块 for mki=0:mk-1 for nki=0:nk-1 % 处理N*N的数据块 for i=N*mki:N*mki+123 for j=N*nki:N*nki+123 c=[i j]'; c=Catmap(c,Cat,N); B(mki*N+c(1)+1,nki*N+c(2)+1,:)=A(i+1,j+1,:); end end end end subplot(2,3,k+1) imagesc(B,[0 1]); A=B; end

相关程序 Cat.m