matlab加停顿,求高手解释下这段程序，实在是搞不懂了，要死的节奏了。。。。...-CSDN博客

该楼层疑似违规已被系统折叠隐藏此楼查看此楼

r=0.06

for k=0:3

whr(k+1,:)=Cwr(1+size(Cwr,2)/4+k*size(Cwr,2)/16:...

size(Cwr,2)/4+(k+1)*size(Cwr,2)/16)-...

Cr(1+size(Cr,2)/4+k*size(Cr,2)/16:...

size(Cr,2)/4+(k+1)*size(Cr,2)/16);

wvr(k+1,:)=Cwr(1+size(Cwr,2)/2+k*size(Cwr,2)/16:...

size(Cwr,2)/2 +(k+1)*size(Cwr,2)/16)-...

Cr(1+size(Cr,2)/2+k*size(Cr,2)/16:...

size(Cr,2)/2+(k+1)*size(Cr,2)/16);

wdr(k+1,:)=Cwr(1+3*size(Cwr,2)/4+k*size(Cwr,2)/16:...

3*size(Cwr,2)/4+(k+1)*size(Cwr,2)/16)-...

Cr(1+3*size(Cr,2)/4+k*size(Cr,2)/16:...

3*size(Cr,2)/4+(k+1)*size(Cr,2)/16);

end

whr=(whr(1,:)+whr(2,:)+whr(3,:)+whr(4,:))/(4*r);

wvr=(wvr(1,:)+wvr(2,:)+wvr(3,:)+wvr(4,:))/(4*r);

wdr=(wdr(1,:)+wdr(2,:)+wdr(3,:)+wdr(4,:))/(4*r);

war=(Cwr(1:size(Cwr,2)/16)-Cr(1:size(Cr,2)/16))/4;

g=0.03;

for k=0:3

whg(k+1,:)=Cwg(1+size(Cwg,2)/4+k*size(Cwg,2)/16:...

size(Cwg,2)/4+(k+1)*size(Cwg,2)/16)-Cg(1+...

size(Cg,2)/4+k*size(Cg,2)/16:...

size(Cg,2)/4+(k+1)*size(Cg,2)/16);

wvg(k+1,:)=Cwg(1+size(Cwg,2)/2+k*size(Cwg,2)/16:...

size(Cwg,2)/2+(k+1)*size(Cwg,2)/16)-Cg(1+...

size(Cg,2)/2+k*size(Cg,2)/16:...

size(Cg,2)/2+(k+1)*size(Cg,2)/16);

wdg(k+1,:)=Cwg(1+3*size(Cwg,2)/4+k*size(Cwg,2)/16:...

3*size(Cwg,2)/4+(k+1)*size(Cwg,2)/16)-...

Cg(1+3*size(Cg,2)/4+k*size(Cg,2)/16:...

3*size(Cg,2)/4+(k+1)*size(Cg,2)/16);

end

whg=(whg(1,:)+whg(2,:)+whg(3,:)+whg(4,:))/(4*g);

wvg=(wvg(1,:)+wvg(2,:)+wvg(3,:)+wvg(4,:))/(4*g);

wdg=(wdg(1,:)+wdg(2,:)+wdg(3,:)+wdg(4))/(4*g);

wag=(Cwg(1:size(Cwg,2)/16)-Cg(1:size(Cg,2)/16))/g;

b=0.12;

for k=0:3

whg(k+1,:)=Cwg(1+size(Cwg,2)/4+k*size(Cwg,2)/16:...

size(Cwg,2)/4+(k+1)*size(Cwg,2)/16)-Cg(1+...

size(Cg,2)/4+k*size(Cg,2)/16:...

size(Cg,2)/4+(k+1)*size(Cg,2)/16);

wvg(k+1,:)=Cwg(1+size(Cwg,2)/2+k*size(Cwg,2)/16:...

size(Cwg,2)/2+(k+1)*size(Cwg,2)/16)-Cg(1+...

size(Cg,2)/2+k*size(Cg,2)/16:...

size(Cg,2)/2+(k+1)*size(Cg,2)/16);

wdg(k+1,:)=Cwg(1+3*size(Cwg,2)/4+k*size(Cwg,2)/16:...

3*size(Cwg,2)/4+(k+1)*size(Cwg,2)/16)-...

Cg(1+3*size(Cg,2)/4+k*size(Cg,2)/16:...

3*size(Cg,2)/4+(k+1)*size(Cg,2)/16);

end

whg=(whg(1,:)+whg(2,:)+whg(3,:)+whg(4,:))/(4*g);

wvg=(wvg(1,:)+wvg(2,:)+wvg(3,:)+wvg(4,:))/(4*g);

wdg=(wdg(1,:)+wdg(2,:)+wdg(3,:)+wdg(4,:))/(4*g);

wag=(Cwg(1:size(Cwg,2)/16)-Cg(1:size(Cg,2)/16))/g;

b=0.12;

for k=0:3

whb(k+1,:)=Cwb(1+size(Cwb,2)/4+k*size(Cwb,2)/16:...

size(Cwb,2)/4+(k+1)*size(Cwb,2)/16)-...

Cb(1+size(Cb,2)/4+k*size(Cb,2)/16:...

size(Cb,2)/4+(k+1)*size(Cb,2)/16);

wvb(k+1,:)=Cwb(1+size(Cwb,2)/2+k*size(Cwb,2)/16:...

size(Cwb,2)/2+(k+1)*size(Cwb,2)/16)-...

Cb(1+size(Cb,2)/2+k*size(Cb,2)/16:...

size(Cb,2)/2+(k+1)*size(Cb,2)/16);

wdb(k+1,:)=Cwb(1+3*size(Cwb,2)/4+k*size(Cwb,2)/16:...

3*size(Cwb,2)/4+(k+1)*size(Cwb,2)/16)-...

Cb(1+3*size(Cb,2)/4+k*size(Cb,2)/16:...

3*size(Cb,2)/4+(k+1)*size(Cb,2)/16);

end

whb=(whb(1,:)+whb(2,:)+whb(3,:)+whb(4,:))/(4*b);

wvb=(wvb(1,:)+wvb(2,:)+wvb(3,:)+wvb(4,:))/(4*b);

wdb=(wdb(1,:)+wdb(2,:)+wdb(3,:)+wdb(4,:))/(4*b);

wab=(Cwb(1:size(Cwb,2)/16)-Cb(1:size(Cb,2)/16))/b;

%--------------------------------------------------

%5重构水印图像

cwr=[war,whr,wvr,wdr];

swr(:,1)=[sqrt(size(war,2)),sqrt(size(war,2)),2*sqrt(size(war,2))];

swr(:,2)=[sqrt(size(war,2)),sqrt(size(war,2)),2*sqrt(size(war,2))];

wr=waverec2(cwr,swr,'haar');

cwg=[wag,whg,wvg,wdg];

swg(:,1)=[sqrt(size(wag,2)),sqrt(size(wag,2)),2*sqrt(size(wag,2))];

swg(:,2)=[sqrt(size(wag,2)),sqrt(size(wag,2)),2*sqrt(size(wag,2))];

wg=waverec2(cwg,swg,'haar');

cwb=[wab,whb,wvb,wdb];

swb(:,1)=[sqrt(size(wab,2)),sqrt(size(wab,2)),2*sqrt(size(wab,2))];

swb(:,2)=[sqrt(size(wab,2)),sqrt(size(wab,2)),2*sqrt(size(wab,2))];

wb=waverec2(cwb,swb,'haar');

%--------------------------------------------------------------------

%三色叠加

temp=size(wr);

pic=zeros(temp(1),temp(2),3);

for i=1:temp(1);

for j=1:temp(2)

pic(i,j,1)=wr(i,j);

pic(i,j,2)=wg(i,j);

pic(i,j,3)=wb(i,j);

end

output=uint8(round(pic));

%转化为uint8

watermark_image_uint8=uint8(output);

imwrite(watermark_image_uint8,'watermark,bmp','bmp');

figure(2);

imshow(watermark_image_uint8);

title('水印');