Mercurial > hg > audiodb
annotate tests/pointset_test/lshP.m @ 557:d3c96cbb91e3
INDEXING: fixed reverse-lookup roundoff error when data gets very large. Computation of offset was lacking an essential cast to off_t.
| author | mas01mc |
|---|---|
| date | Fri, 20 Feb 2009 19:08:45 +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 |