一段关于matlab的小程序 指点一下

[lfchun]

想对三维图像进行立方体覆盖 假设A是该图像的矩阵 用尺寸为a 的正方体覆盖图像,求用到的立方体总个数。
[r,l]=size(A);
a=2;%a是可以变的
x1=floor(r/a)+1;
y1=floor(l/a)+1;
for i=1:x1
for j=1:y1
x=(i-1)*a:i*a;
y=(j-1)*a:j*a;
h(i,j)=A(x,y);
end
end
我的目的是想把 矩阵A 按2*2的网格分割 然后把每一个小网格内的矩阵值赋给h(i,j)进而求h(i,j)的最大最小值. 可是运行老是 提示 我最后一句话错误
有没能人指点一下 该怎么修改 感激不尽

[thinkli]

假设A为4X4矩阵，按你这个程序，r=4， l=4， a=2, ==> x1=3, y1=3 ==> i最大值为3，j最大值为3， ==> x最大值为6，y最大值为6，于是访问A(x,y)时超范围了，就会error "index out of bounds"。

所以求x1,y1那里不该加1。