view general/numerical/pdev.m @ 61:eff6bddf82e3 tip

Finally implemented perceptual brightness thing.

function N=pdev(p,A,D)
% pdev - p-deviation, generalisation of stddev for exponents other than 2
%
% pdev :: nonneg, [[D1]], I:natural -> [[D2]].
% pdev :: nonneg, [[N D]] -> [[ 1 D]].
%
% where D2 = set(D1,I,1), ie the Ith dimension is
% collapsed. I defaults to 1. The p-deviation is defined as 
%    pdev(p,x,i) = mean(abs(x).^p,i).^(1/p)
%
% See also pnorm.

if nargin<2, D=1; end;

if p==2
	N=sqrt(mean(abs(A).^2,D));
else
	N=real(mean(abs(A).^p,D).^(1/p));
end