1、统计各灰度级的像素数nk(1),nk(2),nk(3),...,nk(256)MATLAB中输入如下语句：%求nknk=zeros(1,256);for i=0:255 num=find(I==(i+1));%I为原图像的矩阵 nk(i+1)=length(num); end
2、求概率prMATLAB输入：%求pr（rk）=nk/MN pr=zeros(1,256); for i=1:256 pr(i)=nk(i)/(width*height); end
3、对概率从前至后依次累加s(1) s(2) ...s(256)MATLAB格式：s=zeros(1,25鲂番黟谊6); for i=1:L256 for j=1:i s(i)=s(i)+pr(j); end end
4、将分布函数值量化到[0,255]中对应的灰度级MATLAB语言：sk=zeros(1,256); for i=1:256 sk(i)=round((256-1)*sk(i)); end
5、求灰度级出现的概率MATLAB语言：for i=0:255 pr(i+1)=sum(s(find(sk==i))); end pr %显示pr值
6、将原图像各灰度级的点的灰度改为新的灰度MATLAB语言：J=I; for i=1:256 J(find(J==(i-1)))=sk(i); end
