I=imread('rs.jpg'); 
[m,n]=size(I);%獲得該圖的行列數(shù)。 
m1=zeros(1,256);%為直方圖申請(qǐng)矩陣。 
m2=zeros(1,256); 
for i=1:m 
 for j=1:n 
 m1(I(i,j)+1)=m1(I(i,j)+1)+1;%統(tǒng)計(jì)各級(jí)灰度的像素?cái)?shù)。 
end 
end 
m2(1)=0; 
for k=2:256 
 m3(k)=m1(k)/(m*n);%計(jì)算出灰度的概率。 
 m2(k)=m2(k-1)+m3(k);%累加直方圖。 
end 
M=zeros(m,n); 
for i=1:m 
 for j=1:n 
 M(i,j)=(255*m2(I(i,j))); %將處理后的灰度裝入新的矩陣。 
 end 
end 
P=uint8(M);%轉(zhuǎn)成unit8型。 
subplot(3,2,1);imshow(I);title('原始圖像');subplot(3,2,2);imhist(I);title('原始直方圖');%輸出圖像和直方圖，加以對(duì)比。