sig_dec.m

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

function brate=sig_dec()

load b
load goldseq
goldseq = 1/sqrt(31)*goldseq;


for SNR=0:16
r1=goldseq(1,:)'*b(1,:);
r=awgn(r1,SNR);
y=goldseq(1,:)*r;
y=sign(y);
[bnum(SNR+1),brate(SNR+1)]=symerr(y,b(1,:));
end
SNR=0:16;
semilogy(SNR,brate,'-o');
title('单用户检测');
xlabel('SNR(dB)');
ylabel('BER(dB)');