decon_mud.m

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

function brate=decon_mud()

%加载数据
load goldseq;
load b;
%gold码归一化
goldseq = 1/sqrt(31)*goldseq;
%六信号解相关检测
S=goldseq(1:6,:)';
R=S'*S;
iR=inv(R);
e1=zeros(6,1);
e1(1)=1;
d1=S*iR*e1;
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;

y1=d1'*r;

y1=sign(y1);
[bnum(SNR+1),brate(SNR+1)] = symerr(y1,b(1,:));
end
SNR=0:16;
semilogy(SNR,brate,'-*');
title('解相关检测');
xlabel('SNR(dB)');
ylabel('BER(dB)');