[求助]有哪位高手帮我看看这个程序
我搞了很久都没搞懂,程序如下面,特别是后半段,哪位高手帮我详细的解一下,谢谢了.clcclear
%读入原图象
trueImage=imread('autumn.tif');
alfa=.1;
LENGTH=2500;
subplot(2,3,1);
imshow(trueImage);
title('原始图象');
%对原图象进行DCT变换
dctF1=rgb2gray(trueImage);
subplot(2,3,2);
imshow(dctF1);
title('灰度图象');
I=dct2(dctF1);
subplot(2,3,3)
imshow(log(abs(I)),[ ]);
title('DCT cofficient matrix');
[m,n]=size(I);
%产生水印序列并对其排序
z=1:LENGTH;
randn('state',10);
watermark1=randn(LENGTH,1);
subplot(2,3,4)
[Y0,I0]=sort(watermark1);
plot(I0,Y0,'k')
title('watermark seqence');
%找出水印嵌入位置(幅值较大的n个频域成分)
A=dctF1(:);
[Y1,I1]=sort(A);
x=m*n;
k=LENGTH;
M=zeros(x,1);
%修改幅值较大的n个频域成分的幅值,嵌入水印(因为两个问题不同,所以有两个注释符)
for i=1:x
if k>=1
M(x)=Y1(x).*(1+alfa.*Y0(k));
k=k-1;
else
M(x)=Y1(x);
end
x=x-1;
end
N=zeros(x,1);
x=m*n;
for i=1:x
N(I1(i))=M(i);
end
a=1;
for j=1:n
for i=1:m
dctF2(i,j)=N(a);
a=a+1;
end
end
%DCT反变换,得到嵌入水印的图象
idctF1=idct2(dctF2);
subplot(2,3,5);
imshow(idctF1,[ ]);
title('嵌入水印后的图象');
end