【求助】请教FOR循环程序优化

[purplesapple]

二维直方图是基于象素之二维联合分布密度定义得到的.设(i,j),(k,l)两个任意象素点上的灰度值分别为g(i,j)和g(k,l),则图像灰度值的联合分布密度表示为
P(a,b)=Pk{g(i,j)=a,g(k,l)=b}.其中a,b均为0到L-1之间的灰度整量电平.

我的实现方法:
[m1,n1]=size(ca2);
T=m1*n1;
for ii=1:256
for jj=1:256
I0(ii,jj)=0;
for i1=1:m1
for j1=1:n1
for k=1:m1
for l=1:n1
if (I2(i1,j1)==(ii-1))&(I2(k,l)==(jj-1))
I0(ii,jj)=I0(ii,jj)+1;

end
end
end
end
end
end
end;
其中m1=130,n1=151
此程序运行速度很慢,请教高手有没有什么优化的方法!

谢谢各位的帮助