submmse_bmud.m

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

function submmse_bmud()
%
%2005.12.5

load goldseq;
goldseq = 1/sqrt(31)*goldseq;
s1=goldseq(1,:)';
load b;
r1=goldseq(1,:)'*b(1,:);
r2=goldseq(2,:)'*b(2,:);
r3=goldseq(3,:)'*b(3,:);
r4=goldseq(4,:)'*b(4,:);
r5=goldseq(5,:)'*b(5,:);
r6=goldseq(6,:)'*b(6,:);


for SNR=0:16
    r=r1+r2+r3+r4+r5+r6;
    r=awgn(r,SNR);
    [Us,Ds,Un,Dn] = subspace(r);
    k=size(Ds,1);
    sigma=sum(diag(Dn))/size(Dn,1);
    afa=s1'*Us*(Ds^-1)*Us'*s1;
    w1=(Us*(Ds^-1)*Us'*s1)/afa;
    b1=sign(w1'*r);
    [n,brate(SNR+1)]=symerr(b1,b(1,:));
end
SNR=0:16;
semilogy(SNR,brate,'-*');
title('基于子空间盲MMSE多用户检测');
xlabel('SNR(dB)');
ylabel('BER(dB)');