function [mix] = dwthide(yuan,shui )
%-----------------------读入图像-------------------------------------%
% 开发平台为win8+matlabR2012
% matlab源代码:(注意:代码在matlab2012上测试运行通过,如果是matlab2007的话请将图像转换为灰度图像,利用函数rgb2gray)
% 通过对话框读入原始图像、水印图像,图像格式为.bmp
% 对原始图像进行处理:double()处理,进行小波分解 haar小波基
global ori Cwr Swr a;
ori = double(yuan);
mark = double(shui);
[Cwr,Swr] = wavedec2(mark,1,'haar');
[Cr,Sr] = wavedec2(ori,2,'haar');
a = 0.06;
for k = 0:1:size(Cr,2)/size(Cwr,2)-1
Cr(1+size(Cr,2)/4+k*size(Cwr,2)/4 : size(Cr,2)/4+(k+1)*size(Cwr,2)/4) = Cr(1+size(Cr,2)/4+k*size(Cwr,2)/4 : size(Cr,2)/4+(k+1)*size(Cwr,2)/4)+a*Cwr(1+size(Cwr,2)/4:size(Cwr,2)/2);
Cr(1+size(Cr,2)/2+k*size(Cwr,2)/4 : size(Cr,2)/2+(k+1)*size(Cwr,2)/4) = Cr(1+size(Cr,2)/2+k*size(Cwr,2)/4 : size(Cr,2)/2+(k+1)*size(Cwr,2)/4)+a*Cwr(1+size(Cwr,2)/2:3*size(Cwr,2)/4);
Cr(1+3*size(Cr,2)/4+k*size(Cwr,2)/4 : 3*size(Cr,2)/4+(k+1)*size(Cwr,2)/4) = Cr(1+3*size(Cr,2)/4+k*size(Cwr,2)/4 : 3*size(Cr,2)/4+(k+1)*size(Cwr,2)/4)+a*Cwr(1+3*size(Cwr,2)/4:size(Cwr,2));
end;
Cr(1:size(Cwr,2)/4) = Cr(1:size(Cwr,2)/4) + a*Cwr(1:size(Cwr,2)/4);
% 利用小波基重构图像
ori = waverec2(Cr,Sr,'haar');
mix = uint8(round(ori));
end