dm_bmud.m

来自「这是盲信号的代码 都已经通过编译了 做这方面的同仁可以参考一下 我觉得蛮惯用的」· M 代码 · 共 45 行

function brate=dm_bmud()
%
%2005.11.30
%author xiaobingzhang

%加载数据
load goldseq;
load b;
%gold码归一化
goldseq = 1/sqrt(31)*goldseq;
s1=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;
%直接盲检测
temp=zeros(31);
for i=1:10000
    temp = temp+r(:,i)*r(:,i)';
end
Cr = temp/10000;
m1=Cr^-1*s1;
z=m1'*r;
%for i=2:size(b,2);
%y(i)=sign(z(i)*z(i-1));
%end
y=sign(z);
[bnum(SNR+1),brate(SNR+1)] = symerr(y,b(1,:));
end
SNR=0:16;
semilogy(SNR,brate,'-*');
title('直接方法盲多用户检测');
xlabel('SNR(dB)');
ylabel('BER(dB)');