Mercurial > hg > audiodb

view tests/pointset_test/lshP.m @ 770:c54bc2ffbf92 tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
update tags
author convert-repo
date Fri, 16 Dec 2011 11:34:01 +0000
parents 8fcc7c590a0e
children
line wrap: on
line source
function P2 = lshP(w,c,k)
if nargin<3, k=1;end
if nargin<2, c=1;end
if nargin<1, w=4;end

P2 = 1 - 2*normcdf(-w./c) - 2./(sqrt(2*pi)*(w./c)) .* ( 1-exp(-w.^2./(2*c.^2)) );
if(k~=1)
  P2 = P2.^k;
endif
endfunction