clear all;

Image=double(imread('e:\資料\圖象數(shù)據(jù)庫\monkey.jpg'));
Yimage=double(imread('e:\資料\圖象數(shù)據(jù)庫\lena.jpg'));
Ih=size(Image,1);%圖像的長度
Iw=size(Image,2);%圖像的寬度


Himage=double(zeros(size(Image)));
for(t=1:8)
figure(t)
subplot(3,2,1)
imshow(Image,[]);%顯示原始圖像
title('原始圖像1');
Image=uint8(Image);
figure(t)
subplot(3,2,2)
imhist(Image);
title('原始圖像1的直方圖');
Image=double(Image);
figure(t)
subplot(3,2,3)
imshow(Yimage,[]);%顯示原始圖像
title('原始圖像2');
Yimage=uint8(Yimage);
figure(t)
subplot(3,2,4)
imhist(Yimage);
title('原始圖像2的直方圖');
Yimage=double(Yimage);


% Image=double(zeros(size(Ximage)));
key1=3.9;
key2=0.9;
% key=0;

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



figure(t)
subplot(3,2,5)
imshow(Image,[]);
title('加密圖像');
Image=uint8(Image);
figure(t)
subplot(3,2,6)
imhist(Image);
title('加密圖像的直方圖');
Himage=Image;
Image=Yimage;
Yimage=Himage;
end