idamnjanovic@60: function [dist,ratio] = dictdist(approxD,D,epsilon)
idamnjanovic@60: %DICTDIST Distance between dictionaries.
idamnjanovic@60: %  [DIST,RATIO] = DICTDIST(APPROXD,D) computes the distance between the
idamnjanovic@60: %  approximate dictionary APPROXD and the true dictionary D, where APPROXD
idamnjanovic@60: %  is NxK and D is NxM.
idamnjanovic@60: %
idamnjanovic@60: %  The distance between the dictionary APPROXD and a single atom A of D is
idamnjanovic@60: %  defined as:
idamnjanovic@60: %
idamnjanovic@60: %      DIST(APPROXD,A) = min  { 1-abs(APPROXD(:,i)' * A) }
idamnjanovic@60: %                         i
idamnjanovic@60: %
idamnjanovic@60: %  The distance between the dictionaries APPROXD and D is defined as:
idamnjanovic@60: %
idamnjanovic@60: %      DIST(APPROXD,D) = sum { dist(APPROXD, D(:,k)) } / M
idamnjanovic@60: %                         k
idamnjanovic@60: %
idamnjanovic@60: %  Note that 0 <= DIST(APPROXD,D) <= 1, where 0 implies that all atoms in D
idamnjanovic@60: %  appear in APPROXD, and 1 implies that the atoms of D are orthogonal to
idamnjanovic@60: %  APPROXD.
idamnjanovic@60: %
idamnjanovic@60: %  The similarity ratio between APPROXD and D is defined as:
idamnjanovic@60: %
idamnjanovic@60: %      RATIO(APPROXD,D) = #(atoms in D that appear in APPROXD) / M
idamnjanovic@60: %
idamnjanovic@60: %  where two atoms are considered identical when DIST(A1,A2) < EPSILON with
idamnjanovic@60: %  EPSILON=0.01 by default. Note that 0 <= RATIO(APPROXD,D) <= 1, where 0
idamnjanovic@60: %  means APPROXD and D have no identical atoms, and 1 means that all atoms
idamnjanovic@60: %  of D appear in APPROXD.
idamnjanovic@60: %
idamnjanovic@60: %  [DIST,RATIO] = DICTDIST(DICT1,DICT2,EPSILON) specifies a different value
idamnjanovic@60: %  for EPSILON.
idamnjanovic@60: 
idamnjanovic@60: %  Ron Rubinstein
idamnjanovic@60: %  Computer Science Department
idamnjanovic@60: %  Technion, Haifa 32000 Israel
idamnjanovic@60: %  ronrubin@cs
idamnjanovic@60: %
idamnjanovic@60: %  October 2007
idamnjanovic@60: 
idamnjanovic@60: 
idamnjanovic@60: if (nargin < 3), epsilon = 0.01; end
idamnjanovic@60: 
idamnjanovic@60: [n,m] = size(D);
idamnjanovic@60: 
idamnjanovic@60: approxD = normcols(approxD*spdiag(sign(approxD(1,:))));
idamnjanovic@60: D = normcols(D*spdiag(sign(D(1,:))));
idamnjanovic@60: 
idamnjanovic@60: identical_atoms = 0;
idamnjanovic@60: dist = 0;
idamnjanovic@60: 
idamnjanovic@60: for i = 1:m
idamnjanovic@60:   atom = D(:,i);
idamnjanovic@60:   distances = 1-abs(atom'*approxD);
idamnjanovic@60:   mindist = min(distances);
idamnjanovic@60:   dist = dist + mindist;
idamnjanovic@60:   identical_atoms = identical_atoms + (mindist < epsilon);
idamnjanovic@60: end
idamnjanovic@60: 
idamnjanovic@60: dist = dist / m;
idamnjanovic@60: ratio = identical_atoms / m;