daniele@160: function mu = cumcoherence(obj)
daniele@160: obj = normalize(obj);
daniele@160: [M N] = size(obj.phi);
daniele@160: mu = zeros(M,1);
daniele@160: for m=1:M
daniele@160:     c = zeros(N);
daniele@160:     for i=1:N
daniele@160:         c(:,i) = abs(obj.phi'*obj.phi(:,i));
daniele@160:         c(i,i) = 0;
daniele@160:     end
daniele@160:     c = sort(c,'descend');
daniele@160:     c = c(1:m,:);
daniele@160:     if m==1
daniele@160:         mu(m) = max(c);
daniele@160:     else
daniele@160:         mu(m) = max(sum(c));
daniele@160:     end
daniele@160: end
daniele@160: end