灰度直方图：

灰度直方图描述了一幅图像的绘图统计信息，主要应用于图像分割和图像灰度变换等处理过程中。从数学上说它是一个关于灰度的函数，如令x表示灰度值(一般0≤x≤255)，则f(x)表示当x为特定灰度时，一幅图像上灰度值为x的像素的数量，要注意的是这里的函数f(x)是一个离散的函数。从图形上来说，灰度直方图就是一个二维图，横坐标表示灰度值（灰度级别），纵坐标表示具有各个灰度值或者灰度级别的像素在图像中出现的次数或者概率；

imhist

目的：获取图像的灰度直方图。

用法：imhist( img )；直接显示图像img的灰度直方图；imhist（img，n）；显示一个统计n个灰度级信息的直方图，默认n为256；[counts, x] = imhist( img ) ；获取直方图信息，x为灰度级向量（一维），记录着图像中灰度的所有值，而counts也是一个向量（一维），里面记录着x中对应灰度值出现的个数。

注：imhist只能统计灰度图像的直方图，而对于RGB图像则需要分别统计每个通道的直方图。故使用imhist时，要将rgb图转换为灰度图。

clc;
clear all;
i=imread(\'lena.jpg\');
if(numel(size(i))==3) %判断图像是否非灰度图像；
　　i=rgb2gray(i);
end
imhist(i); %用自带函数imhist实现；
title(\'imhist\');
[m,n]=size(i);
h=zeros(1,256); %生成一个1*256的0矩阵；
for a=1:m
　　for b=1:n
　　　　h(i(a,b)+1)=h(i(a,b)+1)+1; %统计每个像素的数值，由于i(a,b)代表像素数值有0数值，而MATLAB中矩阵下标从1开始，所以用i(a,b)+1代替下标；
　　end
end
figure;
bar(h); %用bar画图函数实现；
axis([0 255 0 3000]); %控制bar的坐标范围，axis([Xmin Xmax Ymin Ymax]);
title(\'histogram\');
xlabel(\'intensity\'); 
ylabel(\'frequency\');

实验结果如下：