clear all;

Ximage=double(imread('C:\Documents and Settings\xdjsj\My Documents\zyh\program\Lena.jpg'));
Yimage=double(imread('C:\Documents and Settings\xdjsj\My Documents\zyh\program\test2.jpg'));
Ih=size(Ximage,1);%圖像的長度
Iw=size(Ximage,2);%圖像的寬度
Image=double(zeros(size(Ximage)));
Image=Ximage;
s=3;
renewing=0;
hiding=0;
                                      %                 隱            藏
key1=4;
key2=0.6;

for(t=1:s)
figure(t)
subplot(3,2,1)
imshow(Image,[]);
title('待隱藏圖像');
Image=uint8(Image);
figure(t)
subplot(3,2,2)
imhist(Image);
title('待隱藏圖像的直方圖');
Image=double(Image);
figure(t)
subplot(3,2,3)
imshow(Yimage,[]);
title('模板圖像');
Yimage=uint8(Yimage);
figure(t)
subplot(3,2,4)
imhist(Yimage);
title('模板圖像的直方圖');
Yimage=double(Yimage);

for(i=1:Ih)
        for(j=1:Iw)
        
               key2=(key1*key2*(1-key2)); 

            if( key2< 0.001)
                key2=key2+0.1;
                Image(i,j)=mod(key2*(Image(i,j)-Yimage(i,j))/2+Yimage(i,j),256);
            else
                Image(i,j)=mod(key2*(Image(i,j)-Yimage(i,j))/2+Yimage(i,j),256);
            end
             
        end
end
figure(t)
subplot(3,2,5)
imshow(Image,[]);
title('結(jié)果圖像');
Image=uint8(Image);
figure(t)
subplot(3,2,6)
imhist(Image);
title('結(jié)果圖像的直方圖');
Image=double(Image);
end

for(i=1:Ih)
        for(j=1:Iw)
            if(Image(i,j)~=Yimage(i,j))
                hiding=hiding+(Image(i,j)-Yimage(i,j))^2;
            end
        end
end
    
   hiding=hiding/(Ih*Iw);
   hiding=hiding^(1/2);
   
                                         %                  恢            復(fù)
key=0.01;
for(t=1:s)
figure(t+s)
subplot(3,2,1)
imshow(Image,[]);
title('結(jié)果圖像');
Image=uint8(Image);
figure(t+s)
subplot(3,2,2)
imhist(Image);
title('結(jié)果圖像的直方圖');
Image=double(Image);
figure(t+s)
subplot(3,2,3)
imshow(Yimage,[]);
title('模板圖像');
Yimage=uint8(Yimage);
figure(t+s)
subplot(3,2,4)
imhist(Yimage);
title('模板圖像的直方圖');
Yimage=double(Yimage);


 
for(i=1:Ih)
       for(j=1:Iw)
         key=1/key2;
               Image(i,j)=mod(key*(Image(i,j)-Yimage(i,j))/2+Yimage(i,j),256);
                
        end
  end
figure(t+s)
subplot(3,2,5)
imshow(Image,[]);
title('恢復(fù)圖像');
Image=uint8(Image);
figure(t+s)
subplot(3,2,6)
imhist(Image);
title('恢復(fù)圖像的直方圖');
Image=double(Image);
end


for(i=1:Ih)
        for(j=1:Iw)
            if(Image(i,j)~=Ximage(i,j))
                renewing=renewing+(Image(i,j)-Ximage(i,j))^2;
            end
        end
    end
    
   renewing=renewing/(Ih*Iw);
   renewing=renewing^(1/2);