main_mud.m

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

function main_mud()
%2005.11.29
%加载数据
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;
%MMSE检测器
S=goldseq(1:6,:)';
R=S'*S;
A1=1;
A2=2;
A3=5;
A4=10;
A5=15;
A6=20;
A=[A1,A2,A3,A4,A5,A6];
iR=inv(R+sum(A)^2*0.1*diag(A)^-2);
e1=zeros(6,1);
e1(1)=1;
m1=S*iR*e1;
%s1=goldseq(1,:)';
%产生六个幅度为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=A1*r1+A2*r2+A3*r3+A4*r4+A5*r5+A6*r6;

%y1=s1'*r;
%y1=d1'*r;
y1=m1'*r;
y1=sign(y1);
[bnum(SNR+1),brate(SNR+1)] = symerr(y1,b(1,:));
end




%加载数据
load goldseq;
load b;
%gold码归一化
goldseq = 1/sqrt(31)*goldseq;
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,:);
r=r1+r2+r3+r4+r5+r6;
SNR=10;
r=awgn(r,SNR);
%子空间
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,:);
r=r1+r2+r3+r4+r5+r6;

k=size(Ds,1);
k=size(Ds,1);
sigma=sum(diag(Dn))/size(Dn,1);
afa=s1'*Us*(Ds-sigma*eye(k))^-1*Us'*s1;
w1=(Us*(Ds-sigma*eye(k))^-1*Us'*s1)/afa;
b1=sign(w1'*r);
[n,brate]=symerr(b1,b(1,:));