#histcmpchi2.m#

来自「KPMstats is a directory of miscellaneous」· M# 代码 · 共 12 行

function d = histCmpChi2(h1, h2)
% Compare two histograms using chi-squared
% function d = histCmpChi2(h1, h2)
%
% d(i,j) = chi^2(h1(i,:), h2(j,:)) = sum_b  (h1(b

[N B] = size(h1);
d = zeros(N,N);
for i=1:N
  hh1 = repmat(h1(i,:), N, 1);
  numer = (hh1 .- h2).^2;