学步园

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