original0?=?imread（‘/Users/yanning/Pictures/test.jpg‘）;
original?=?rgb2gray（original0）;

figure;

%output?the?original?picture
subplot（231）;
imshow（original）;
title（‘original‘）;


histogram0?=?zeros（1256）;%原圖灰度矩陣
histogram1?=?zeros（1256）;%均衡化之后的灰度矩陣
histogram2?=?zeros（1256）;%線性變換后的灰度矩陣

[rowcol]?=?size（original）;
for?i=1:row
????for?j=1:col
????????value?=?double（original（ij））+1;
????????histogram0（value）?=?histogram0（value）+1;
????end
end
subplot（234）;
plot（histogram0）;

%直方圖均衡化
cdf（1）?=?histogram0（1）;
for?i=2:256
????cdf（i）?=?cdf（i-1）+histogram0（i）;
end
for?i=1:row
????for?j=1:col
????????value?=?original（ij）;
????????original（ij）?=?round（cdf（value）*256/（row*col））;
????????if（original（ij）>255）
????????????original（ij）=255;
????????end
????????if（original（ij）<0）
????????????