logistic.m

数字图象加密程序（编写的数字图象加密程序）

clear all;

Image=double(imread('e:\work\lena.jpg'));

Ih=size(Image,1);%图像的长度
Iw=size(Image,2);%图像的宽度
Biao=(Ih*Iw)/256;
Junfang=double(zeros(1,1));

% for(t=1:3)
figure(1)
subplot(2,1,1)
imshow(Image,[]);%显示原始图像
title('原始图像');
Image=uint8(Image);
figure(1)
subplot(2,1,2)
imhist(Image);
title('原始图像的直方图');


True=double(zeros(1,256));


Image=double(Image);

key1=4;
key2=0.4;
Kimage=double(zeros(size(Image)));
key=0;
% Himage=zeros(size(Image));


% Kimage(1)=key2;
for(i=1:Ih)
        for(j=1:Iw)
            key2=key1*(key2*(1-key2));
            key=mod((key2*256),256);
            Kimage(i,j)=round(key);            
        end
end

for(i=1:Ih)
        for(j=1:Iw)
             Image(i,j)=mod((Kimage(i,j)+Image(i,j)),256);
         end
end
     
% Image=mod((Kimage+Image),256);
figure(2)
subplot(2,1,1)
imshow(Image,[]);
title('加密图像');
Image=uint8(Image);
figure(2)
subplot(2,1,2)
imhist(Image);
title('加密图像的直方图');
% for(i=1:Ih)
%         for(j=1:Iw)
%             for(k=1:256)
%                 if(Image(i,j)==(k-1))
%                     True(k)=True(k)+1;
%                 end
% %             break
%             
%             end
%         end
%     end
%   for(k=1:256) 
%       True(k)=(True(k)-Biao)^2;
%       Junfang=Junfang+True(k);
%   end
%   Junfang=Junfang/(Ih*Iw);
% end
% for(k=1:3)
%     Junfang(k)=Junfang(k)/(Ih*Iw);
% end