function A = BC_wspec_reference(P)
	% determina i parametri del codice
	[K,N_minus_K] = size(P);
	N = K + N_minus_K;
	
	% alloca la matrice delle molteplicità
	A = zeros(1+K,1+N);
	
	% calcola le parole di codice ed i relativi pesi
	for ik = 0:2^K-1
		b = dec2bin(ik,K)-'0';
		c = [b,mod(b*P,2)];
		A(1+sum(b),1+sum(c)) = A(1+sum(b),1+sum(c)) + 1;
	end
	