代码

终于搞定啦，dct读取图像的时候如果是值，会发现dct始终产生一个[1,10]的数组。
终于发现了问题的所在。
论文的思路差不多可以进行模拟啦，
今天发现上次和导师讨论的思路好像有点问题，空出的点并不能都嵌入呀，因为如果嵌入会无法分辨嵌入点和原来的点，所以可以插入的实际点要少！
下面是今天调试的代码

clear
%
%读入原图象
trueImage=imread('C:oat.png');
alfa=.1;
LENGTH=2500;
subplot(2,2,1);
imshow(trueImage);
title('原始图象');
%
%对原图象进行DCT变换
%dctF1=dct2('C:9.jpg');
dctF1=dct2(trueImage);
subplot(2,2,2);
imshow(log(abs(dctF1)),[]);
title('DCTcofficientmatrix');
[m,n]=size(dctF1);

%
%产生水印序列并对其排序
%c=rand('copyright',10);%bydont
%watermark1=radon(LENGTH,1);
watermark1=rand(LENGTH,1);%bydont
subplot(2,2,3);
%imshow(watermark1,[]);
imshow(watermark1,[]);
title('watermarkseqence')
[Y0,I0]=sort(watermark1);
%
%找出水印嵌入位置（幅值较大的n个频域成分）
A=dctF1(:);
[Y1,I1]=sort(A);
x=m*n;
k=LENGTH;
M=zeros(x,1);
%
%修改幅值较大的n个频域成分的幅值，嵌入水印（因为两个问题不同，所以有两个注释符）
fori=1:x
ifk>=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;
fori=1:x
N(I1(i))=M(i);
end
a=1;
forj=1:n
fori=1:m
dctF2(i,j)=N(a);
a=a+1;
end
end
%
%DCT反变换，得到嵌入水印的图象
idctF1=idct2(dctF2);
subplot(2,2,4);
imshow(idctF1,[]);
title('嵌入水印后的图象');