寻找灰度图象分割阈值的matlab程序，结果不对，不知道问题在哪？

P1(t)=h(z)的累加 from z=0 到 z=t （1）
P2(t)=h(z)的累加 from z=t+1 到z=L-1

U1(t)=h(z)*z的累加 from z=0 到z=t，累加后的结果再除以P1(t) （2）
U2(t)=h(z)*z的累加from z=t+1 到z=L-1，累加后的结果再除以P2(t)

K1(t)=(z-u1(t))* (z-u1(t))*(z-u1(t) 的累加 from z=0 到z=t，累加后的结果再除以P1(t)
K2(t)= (z-u2(t))* (z-u2(t))*(z-u2(t) 的累加from z=t+1 到z=L-1，累加后的结果再除以P2(t) （3）

K(t)=|k1(t)|+|k2(t)| （4）

这是一个寻找最佳阈值进行图像分割的方法，根据公式4求得以各个灰度值做为阈值计算得出的偏态值，存入到tab（）数组中，其中最小值tab（i）对应的i-1即为最佳阈值。可是运行的结果出来后，tab数组中256个值都是0。请问程序哪儿出错了呢？应该怎样改呢？

a=imread('YuChu.bmp');
[M,N]=size(a);
h=imhist(a); %h中存放直方图数据
for i=1:256
tab(i)=0;
end

p=0;u=0;
p1=0;u1=0;
for i=0:255
p=p+h(i+1);
u=u+h(i+1)*i;
end
for j=0:255
for i=0:j
p1=p1+h(i+1);
u1=u1+h(i+1)*i;
end
p2=p-p1;
u2=u-u1;
if (p1) %计算一阶统计矩，公式2
u1=u1/p1;
else
u1=0;
end
if (p2)
u2=u2/p2;
else

什么是灰度图象 “渗透阈值”的定义？ ps里的阈值 灰度图象和伪彩色图象（如256色图象）的颜色表有什么区别？ 800X600的8位灰度图象的大小的计算方法是什么? 生态阈值的定义是什么 阈强度(阈值)的概念? 寻找一个合适的文件分割器 图像的灰度是什么意思？？ 什么是显示器的灰度？