cl;
img=imread(\'mask.jpg\');
imshow(img);
[x,y]=size(img);
img_man=zeros(x,y);
img_com=zeros(x,y);

%% 直方图均衡化算法
Max=max(max(img));
Min=min(min(img));
Hist=zeros(1,256);
for i=1:x
    for j=1:y
        Hist(img(i,j)+1)=Hist(img(i,j)+1)+1;
    end
end
figure,plot(Hist);

p=zeros(1,256);
for i=1:256
    p(i)=Hist(i)/(x*y);
end
figure,plot(p);

c=zeros(1,256);
for i=1:256
    c(i)=sum(p(1:i));
end
figure,plot(c);

for i=1:x
    for j=1:y
        img_man(i,j)=c(img(i,j)+1)*(Max-Min)+Min;
    end
end

figure,imshow(uint8(img_man))

Hist2=zeros(1,256);
for i=1:x
    for j=1:y
        Hist2(img_man(i,j)+1)=Hist2(img_man(i,j)+1)+1;
    end
end
figure,plot(Hist2);

%% matlab直方图均衡化函数
img_com=histeq(img);     
figure,imshow(img_com)

Hist3=zeros(1,256);
for i=1:x
    for j=1:y
        Hist3(img_com(i,j)+1)=Hist3(img_com(i,j)+1)+1;
    end
end
figure,plot(Hist3);