MATLAB提取灰度直方图特征

clc 
clear  all
I=imread （‘E:yangben/0/1.bmp‘）; 
J=I; 
%计算灰度图象的直方图数据 
L=256; %灰度级 
Ps = zeros（L1）; %统计直方图结果数据 
nk=zeros（L1）; 
[rowcol]=size（I）; 
n=row*col; %总像素个数 
for i = 1:row 
for j = 1:col 
num = double（I（ij））+1; %获取像素点灰度级 
nk（num） = nk（num）+1; %统计nk 
end 
end 
%计算直方图概率估计 
mu = 0;
for i=1:L 
Ps（i）=nk（i）/n;
mu = mu + i * Ps（i）;
end 

delta = 0;
for i=1:L
    delta