帮帮忙：matlab基于dct数字水印程序运行有问题

Question

以下是我在网上搜的一个关于数字水印matlab实现dct域的程序，但运行时报错，不知道出错原因在哪里，程序源码如下：
clear all;
clc;
start_time=cputime;
%%%%%%%%%%%% 读取水印图像 %%%%%%%%
I=imread('mark.bmp');
I=rgb2gray(I);
I=double(I)/255;
I=ceil(I);
%%%%%%%%%%显示水印图像%%%%%%%%%%%%%
figure(1);
subplot(2,3,1);
imshow(I),title('水印图像')
dimI=size(I);
rm=dimI(1);cm=dimI(2);
%%%%%%%%%%%%%%%5 以下生成水印信息 %%
mark=I;
alpha=50,
k1=randn(1,8);
k2=randn(1,8);
a0=imread('lena.bmp');
psnr_cover=double(a0);
subplot(2,3,2),imshow(a0,[]),title('载体图像');
[r,c]=size(a0);
cda0=blkproc(a0,[8,8],'dct2');
%%%%%%%%%%%%%%%%%%%%% 嵌入 %%%%%%%%%%
cda1=cda0; % cda1 = 256_256
for i=1:rm % i=1:32
for j=1:cm % j=1:32
x=(i-1)*8;y=(j-1)*8;
if mark(i,j)==1
k=k1;
else
k=k2;
end
cda1(x+1,y+8)=cda0(x+1,y+8)+alpha*k(1);
cda1(x+2,y+7)=cda0(x+2,y+7)+alpha*k(2);
cda1(x+3,y+6)=cda0(x+3,y+6)+alpha*k(3);
cda1(x+4,y+5)=cda0(x+4,y+5)+alpha*k(4);
cda1(x+5,y+4)=cda0(x+5,y+4)+alpha*k(5);
cda1(x+6,y+3)=cda0(x+6,y+3)+alpha*k(6);
cda1(x+7,y+2)=cda0(x+7,y+2)+alpha*k(7);
cda1(x+8,y+1)=cda0(x+8,y+1)+alpha*k(8);

end
end
%%%%% 嵌入水印后图像 %%%%%%%%%%%%%%
。。。。。。。。
报错为On line 40 ==> cda1(x+1,y+8)=cda0(x+1,y+8)+alpha*k(1);
写不下那么多，程序运行结果只显示了水印图像和源图像，并未显示嵌入水印后图像和提取的水印，希望高手们指教下，小妹很着急
??? Index exceeds matrix dimensions.

Error in ==> D:\matlab\work\last\watermark.m
On line 40 ==> cda1(x+1,y+8)=cda0(x+1,y+8)+alpha*k(1);
你分析的挺对的，是矩阵越界的问题，但不知道怎么改，载体图像是按照8X8分块的。。。

Answer 1

拜托,把matlab的错误提示全部粘上来
你看下cda1(x+1,y+8)和cda0(x+1,y+8)的下标越界没有,这个可能性很大
---------------------------------------------------------------------
我不知道你用的一些函数到底是怎么写的,所以只能给点建议:
首先,你看下cda0是几X几的矩阵.
for i=1:rm % i=1:32
for j=1:cm % j=1:32
x=(i-1)*8;y=(j-1)*8;这3句话是错误的根源,如果cda0是8X8的矩阵,而你在32X32的区间索引,那么必定出错.你需要做的是,仔细计算X,Y的值,让他们正好落在cda0的范围内句可以了.

Answer 2

有木有提取部分的程序啊，我弄了半天不成功诶