直方图-----方差

clc
clear all
img=imread('3.jpg');
p=imhist(img);
L=length(p);
sum=0;
for i=1:L
    for j=1:L
      sum=sum+((p(i)-p(j))*(p(i)-p(j)));
    end
end
result=sum/(2*256*256)