Mercurial > hg > camir-aes2014
annotate toolboxes/FullBNT-1.0.7/KPMstats/histCmpChi2.m @ 0:e9a9cd732c1e tip
first hg version after svn
| author | wolffd |
|---|---|
| date | Tue, 10 Feb 2015 15:05:51 +0000 |
| parents | |
| children |
| rev | line source |
|---|---|
| wolffd@0 | 1 function d = histCmpChi2(h1, h2) |
| wolffd@0 | 2 % Compare two histograms using chi-squared |
| wolffd@0 | 3 % function d = histCmpChi2(h1, h2) |
| wolffd@0 | 4 % |
| wolffd@0 | 5 % d(i,j) = chi^2(h1(i,:), h2(j,:)) = sum_b (h1(i,b)-h2(j,b)^2 / (h1(i,b) + h2(j,b)) |
| wolffd@0 | 6 |
| wolffd@0 | 7 [N B] = size(h1); |
| wolffd@0 | 8 d = zeros(N,N); |
| wolffd@0 | 9 for i=1:N |
| wolffd@0 | 10 h1i = repmat(h1(i,:), N, 1); |
| wolffd@0 | 11 numer = (h1i - h2).^2; |
| wolffd@0 | 12 denom = h1i + h2 + eps; % if denom=0, then numer=0 |
| wolffd@0 | 13 d(i,:) = sum(numer ./ denom, 2); |
| wolffd@0 | 14 end |
| wolffd@0 | 15 |