代码

终于搞定啦，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(
  '
  DCT cofficient matrix
  '
  );
[m,n]
  =
  size(dctF1);

  %
  
  %
  产生水印序列并对其排序
  %
  c
  =
  rand(
  '
  copyright
  '
  ,
  10
  ); 
  %
  by dont
  %
  watermark1
  =
  radon(LENGTH,
  1
  );
watermark1
  =
  rand(LENGTH,
  1
  ); 
  %
  by dont
subplot(
  2
  ,
  2
  ,
  3
  );
  %
  imshow(watermark1,[]);
imshow(watermark1,[]);
title(
  '
  watermark seqence
  '
  )
[Y0,I0]
  =
  sort(watermark1);
  %
  
  %
  找出水印嵌入位置（幅值较大的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
  ,
  2
  ,
  4
  );
imshow(idctF1,[ ]);
title(
  '
  嵌入水印后的图象
  '
  );