view arrows/stats/stats_cov.m @ 61:eff6bddf82e3 tip

Finally implemented perceptual brightness thing.

function C=stats_cov(stats)
	m1=mstep(stats);
	C=m1.cov;
end

function M=mstep(S)
	M.mu=S.amu/S.ag;
	M.cov=msym(S.aco/S.ag - M.mu*M.mu'); 
end

function A=msym(X), A=(X+X')/2; end