g=imread (00.jpg);
I=rgb2gray(g);
J=I;
%計算灰度圖象的直方圖數據
L=256; %灰度級 
Ps = zeros(L,1); %出現頻率%統計直方圖結果數據
nk=zeros(L,1);%出現次數
[row,col]=size(I);
n=row*col; %總像素個數
for i = 1:row
for j = 1:col
num = double(I(i,j))+1; %獲取像素點灰度級
nk(num) = nk(num)+1; %統計nk 
end
end
%計算直方圖概率估計
for i=1:L
Ps(i)=nk(i)/n; 
end
figure;
subplot(3,1,1);imshow(J),title('原圖');
subplot(3,1,2),plot(nk),title('直方圖（nk）');
subplot(3,1,3),plot(Ps),title('直方圖（Ps）');