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?