annotate toolboxes/FullBNT-1.0.7/KPMstats/histCmpChi2.m @ 0:cc4b1211e677 tip

initial commit to HG from Changeset: 646 (e263d8a21543) added further path and more save "camirversion.m"
author Daniel Wolff
date Fri, 19 Aug 2016 13:07:06 +0200
parents
children
rev   line source
Daniel@0 1 function d = histCmpChi2(h1, h2)
Daniel@0 2 % Compare two histograms using chi-squared
Daniel@0 3 % function d = histCmpChi2(h1, h2)
Daniel@0 4 %
Daniel@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))
Daniel@0 6
Daniel@0 7 [N B] = size(h1);
Daniel@0 8 d = zeros(N,N);
Daniel@0 9 for i=1:N
Daniel@0 10 h1i = repmat(h1(i,:), N, 1);
Daniel@0 11 numer = (h1i - h2).^2;
Daniel@0 12 denom = h1i + h2 + eps; % if denom=0, then numer=0
Daniel@0 13 d(i,:) = sum(numer ./ denom, 2);
Daniel@0 14 end
Daniel@0 15