clear?all;
A?=?imread（‘girl.bmp‘）;
imshow（A）;
[m?n]?=?size（A）;

B?=?zeros（1256）;%用來繪制直方圖

for?i?=?1:m
????for?j?=1:n
????????B（A（ij）+1）?=?B（A（ij）+1）?+?1;
????end
end

for?k?=?1:256
????B（k）?=?B（k）/m/n;
end

figurebar（1:256B‘g‘）;?????????????????%繪制原直方圖
title（‘原圖像直方圖‘）
xlabel（‘灰度值‘）
ylabel（‘出現概率‘）
axis（[025600.025]）


C?=?zeros（1256）;
C?=?cumsum（B）;

D?=?zeros（1256）;
for?k?=?1:256
????D（k）?