I=imread（‘cumtlogo.bmp‘）;%讀入原圖
II=im2double（I）;??%轉化為[01）double型??????????%II為原圖像
[mn]=size（II（::1））;??%原圖像大小
af=0.1;??%嵌入強度
[USV]=svd（II（::1））;??%進行奇異值分解
M=imread（‘logo.tif‘）;??%讀入水印圖像
W=im2double（M）;??%轉化為[01）double型
[m1n1]=size（W（::1））;
WW=zeros（mn）;
for?i=1:m1
????for?j=1:n1
????????????WW（ij）=W（ij）;
????end
end
S1=S+af*WW;%加入水印后的對角陣
[U1SSV1]=svd（S1）;?%再進行奇異值分解
CWI=U*SS*V‘;??%嵌入水印后圖像
subplot（221）;?imshow（II（::1））;?title（‘原圖像‘）;??%顯示原圖像
su