www.pudn.com > projective-and-affine.rar > var2.m

% V = VAR2(MTX,MEAN)
%
% Sample variance of a matrix.
%  Passing MEAN (optional) makes the calculation faster.

function res = var2(mtx, mn)

if (exist('mn') ~= 1)
  mn =  mean(mean(mtx));
end

if (isreal(mtx))
  res = sum(sum(abs(mtx-mn).^2)) / (prod(size(mtx)) - 1);
else
  res = sum(sum(real(mtx-mn).^2)) + i*sum(sum(imag(mtx-mn).^2));
  res = res  / (prod(size(mtx)) - 1);
end