annotate tests/pointset_test/lshP.m @ 436:e43f8a7aca93 api-inversion

really zap all unnecessary sequenceLength from query.cpp
author mas01cr
date Wed, 24 Dec 2008 10:56:01 +0000
parents 8fcc7c590a0e
children
rev   line source
mas01mc@287 1 function P2 = lshP(w,c,k)
mas01mc@287 2 if nargin<3, k=1;end
mas01mc@287 3 if nargin<2, c=1;end
mas01mc@287 4 if nargin<1, w=4;end
mas01mc@287 5
mas01mc@287 6 P2 = 1 - 2*normcdf(-w./c) - 2./(sqrt(2*pi)*(w./c)) .* ( 1-exp(-w.^2./(2*c.^2)) );
mas01mc@287 7 if(k~=1)
mas01mc@287 8 P2 = P2.^k;
mas01mc@287 9 endif
mas01mc@287 10 endfunction