wu2.m

The goal of SPID is to provide the user with tools capable to simulate, preprocess, process and clas

function [p] = wu2(r,k)

p =[];
Q = zeros(k);
sum_r = [];

for i = 1:k
    sum_r(i) = (norm(r(:,i),2))^2;
end

for i = 1:k
    
     for j = 1:k
        if j==i
            Q(i,j) = sum_r(i) - (r(i,j)^2);
        else
            Q(i,j) = -r(i,j)*r(j,i);
        end
     end
   
end

A = [Q ones(k,1); (ones(1,k)) 0 ];
b = [zeros(k,1);1;];

p = A \ b;

% otherwise sometimes zeros with a neg sign! 
p = abs(p);

p = p / sum(p(1:k));

p = (p(1:k))';

end