dist.m

I developed an algorithm for using local ICA in denoising multidimensional data. It uses delay embed

function err = dist(A,B,pot);
% err = dist(A,B)
%
% calculates L^pot distance between A and B  
%
% 1.6.03 initial version
% 1.9.04 rewrite
% Peter Gruber

error(nargchk(2,3,nargin))     % check no. of input args

if (nargin<3) pot=2; end

A=reshape(A,1,prod(size(A)));
B=reshape(B,1,prod(size(B)));

err=sum(sum(sqrt((A-B).*conj(A-B)).^pot).^(1/pot))/size(A,2);

return;