n = 128;   %哈達(dá)碼矩陣的階數(shù)
n1 = log2(n);
W = zeros(n,n);
W(1,1) = 1;
for t = 1 : n1  %產(chǎn)生哈達(dá)碼矩陣
    for i = 1 : n
        for j = 1 : n
            if (i <= 2^(t-1) && j <= 2^(t-1))
                W(i ,j) = W(i ,j);
            elseif (i <= 2^(t-1) && j > 2^(t-1))
                W(i ,j) = W (i , j-2^(t-1));
            elseif (i > 2^(t-1) && j <= 2^(t-1))
                W(i ,j) = W (i-2^(t-1) , j);
            else
                W(i ,j) = -W (i-2^(t-1) , j-2^(t-1));
            end     
        end
    end   
end
x=randint(1,n);       %發(fā)送數(shù)據(jù)
t=10;   %分配給此用戶的為哈達(dá)碼矩陣的第t行
s=W(t,:).*x;     %產(chǎn)生發(fā)送的數(shù)據(jù)
snr=-10;          %信噪比
enn = 10^(snr/10);      
sigman = sqrt(var(s)/(2*enn));        %噪聲方差
x_noise = zeros(1,n);
x_noise = sigman*(randn(1,n));          %產(chǎn)生噪聲
r = x_noise + s;          %接收信號(hào)
for i=1:n
   v(i)=sum(r.*W(i,:));            %相關(guān)接收機(jī)處理
end
num=find(v==max(v));           %找到最大值所對(duì)應(yīng)的沃爾什碼