% designed by superdont % blog.youkuaiyun.com/superdont clc; clear all; o=imread('E:\工作\科研\奇异值\image\lena.bmp'); od=double(o); w=imread('E:\工作\科研\奇异值\image\w.bmp'); %对分好的块进行svd变换 [u1 d1 v1]=svd(od); d12=d1+10*w; [u2 d2 v2]=svd(d12); ow=u1*d2*v1'; ow=uint8(ow); figure,imshow(ow); imwrite(ow,'E:\工作\科研\奇异值\image\lenawatermarked.bmp','bmp'); owr=imread('E:\工作\科研\奇异值\image\lenawatermarked.bmp'); owr=double(owr); [tu1 td1 tv1]=svd(owr); td12=u2*td1*v2'; tw=td12-d1 for i=1:512 for j=1:512 if tw(i,j)>5 tw(i,j)=1; else tw(i,j)=0; end end end % tw=uint8(tw) figure,imshow(tw,[]);