混沌置乱演示

function zhiluanAll 
%对原始图像进行置乱
clear; 
clc; 
%os原始图像的路径 
[filename pathname]=uigetfile('*.bmp;*.tiff;*.tif', '读入图像'); 
os=[pathname filename]; 
o=imread(os); 
% figure,imshow(o);
%计算原始图像的大小 
[m n]=size(o); 
%l=zeros(m,n); 
l=linspace(0,0,m*n); 
l(1)=0.98; 
for i=2:m*n 
    l(i)=1-2*l(i-1)*l(i-1); 
end 
% %对生成的混沌序列进行排序 
% [lsort lindex]=sort(l); 
% t=linspace(0,0,m*n); 
% t(1)=0.98; 
% for i=1:m*n 
%     t(i+1)=1-2*t(i)*t(i); 
% end 
% for i=1:m*n 
%     if (l(i)>=0) 
%         l(i)=1; 
%     else 
%         l(i)=0; 
%     end 
% end 
% l=~l;
% subplot(2,2,2),imshow(l),title('混沌图像'); 
[lsort,lindex]=sort(l); 
%ozl =original zhi luan 表示置乱后的原始图像 
ozl=zeros(m,n); 
for i=1:m*n 
    ozl(i)=o(lindex(i)); 
end 
% ozl=reshape(ozl,m,n);
%%oo =original fu yuan 
ofy=zeros(m,n); 
for i=1:m*n 
    ofy(lindex(i))=ozl(i); 
end 
h= waitbar(0,'程序处理中，请耐心等待。。。'); 
for i=1:100, % computation here % 
waitbar(i/100) 
end 
close(h) ; 
%%%%%%%%%%%%%%%%%%错误的初始密钥
l=linspace(0,0,m*n); 
l(1)=0.9800001; 
for i=2:m*n 
    l(i)=1-2*l(i-1)*l(i-1); 
end 
[lsort,lindex]=sort(l); 
ofyError=zeros(m,n); 
for i=1:m*n 
    ofyError(lindex(i))=ozl(i); 
end
%%%%%%%%%%%%%%%%%%错误的初始密钥
subplot(2,2,1),imshow(o);%,title('原始图像'); 
title('(a)原始图像','position',[33.5 80.28]); 
subplot(2,2,2),imshow(ozl,[]);%,title('加密图像'); 
title('(b)加密图像','position',[33.5 80.28]); 
subplot(2,2,3),imshow(ofy,[]);%,title('复原图像'); 
title('(c) 正确密钥解密','position',[33.5 80.28]); 
subplot(2,2,4),imshow(ofyError,[]);%,title('复原图像'); 
title('(d) 错误密钥解密','position',[33.5 80.28]);