代码如下：

I = imread('');
dark_I=double(I)/255; 
avg=mean2(dark_I);  %求图像均值
[m,n]=size(dark_I);
s=0;
for x=1:mfor y=1:ns=s+(dark_I(x,y)-avg)^2; %求得所有像素与均值的平方和。end
end%求图像的方差
a1=var(dark_I(:)); %第一种方法：利用函数var求得。
a2=s/(m*n-1); %第二种方法：利用方差公式求得
a3=(std2(dark_I))^2; %第三种方法：利用std2求得标准差，再平方即为方差。

三种方法计算结果：

三种计算方式的结果是一样的。