con_dec.m

这是盲信号的代码 都已经通过编译了 做这方面的同仁可以参考一下 我觉得蛮惯用的

function brate=con_dec()
load goldseq;
load b;
%gold码归一化
goldseq = 1/sqrt(31)*goldseq;
%产生六个幅度为1的用户
for SNR=0:16
r1=goldseq(1,:)'*b(1,:);
r1=awgn(r1,SNR);
r2=goldseq(2,:)'*b(2,:);
r2=awgn(r2,SNR);
r3=goldseq(3,:)'*b(3,:);
r3=awgn(r3,SNR);
r4=goldseq(4,:)'*b(4,:);
r4=awgn(r4,SNR);
r5=goldseq(5,:)'*b(5,:);
r5=awgn(r5,SNR);
r6=goldseq(6,:)'*b(6,:);
r6=awgn(r6,SNR); 
r=r1+r2+r3+r4+r5+r6;
%直接解扩
s1=goldseq(1,:)*r;
s1=sign(s1);
[bnum(SNR+1),brate(SNR+1)] = symerr(s1,b(1,:));
end
SNR=0:16;
semilogy(SNR,brate,'-*');
title('传统检测');
xlabel('SNR(dB)');
ylabel('BER(dB)');