clc?
clear??all
I=imread?（‘E:yangben/0/1.bmp‘）;?
J=I;?
%計(jì)算灰度圖象的直方圖數(shù)據(jù)?
L=256;?%灰度級(jí)?
Ps?=?zeros（L1）;?%統(tǒng)計(jì)直方圖結(jié)果數(shù)據(jù)?
nk=zeros（L1）;?
[rowcol]=size（I）;?
n=row*col;?%總像素個(gè)數(shù)?
for?i?=?1:row?
for?j?=?1:col?
num?=?double（I（ij））+1;?%獲取像素點(diǎn)灰度級(jí)?
nk（num）?=?nk（num）+1;?%統(tǒng)計(jì)nk?
end?
end?
%計(jì)算直方圖概率估計(jì)?
mu?=?0;
for?i=1:L?
Ps（i）=nk（i）/n;
mu?=?mu?+?i?*?Ps（i）;
end?

delta?=?0;
for?i=1:L
????delta?