wu1.m

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

function [p] = wu1(r,k)

p = [];
Q = zeros(k);
row_sum = sum(r');

for i = 1:k
   
    for j = 1:k

        if j == i
	    Q(i,j) = (row_sum(i) - r(i,j))/(k-1);
        else
            Q(i,j) = r(i,j)/(k-1);
        end

    end
   
end

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

p = A \ b;

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

p = (p / sum(p))';

end