我编的16*16图像分块的程序，请大家帮忙看一下

[wryclY]

f=imread('pout.tif');
[m,n,l]=size(f);
bk=16;
k1=m/bk;
k2=n/bk;
m=k1*bk;
n=k2*bk;
g=zeros(bk,bk);
if l==3
f=rgb2ycbcr(f);
end
g(i,j)=f(1:bk,1:bk,1);
figure;
imshow(g(i,j));

我的循环没有编完，写出来总是不对，请大家帮帮忙帮我编完好吗？谢谢！！