zuoye1.m

硕士CDMA作业

clear all
%function zuoye1
r=[0,0,0,0,0,1];
c=zeros(63,6);
c1=zeros(63,6);
c(1,:)=r;
c1(1,:)=r;
for i=2:63
    c(i,2:6)=c(i-1,1:5);
    c(i,1)=mod((c(i-1,1)+c(i-1,6)),2);
    c1(i,2:6)=c1(i-1,1:5);
    c1(i,1)=mod((c(i-1,1)+c(i-1,2)+c(i-1,5)+c(i-1,6)),2);
end
for i=1:63
    for j=1:6
        c(i,j)=(-1)^(c(i,j));%把0,1变成1.-1
        c1(i,j)=(-1)^(c1(i,j));
    end
end

l1=50;
wu199=zeros(l1,1);
wu299=zeros(l1,1);
for t=1:l1
    l=10;
    for t1=1:l
d=200 ;
b1=randsn(d);
b2=randsn(d);
s1=zeros(1,1);
s2=zeros(1,1);
for m=1:d
        r1=c(:,1)*b1(1,m);
        r2=c1(:,1)*b2(1,m);
        s1=[s1;r1];
        s2=[s2;r2];
end
s1=s1(2:63*d+1,1);
s2=s2(2:63*d+1,1);
Noi1=normrnd(0,t,63*d,1);
Noi2=normrnd(0,t,63*d,1);
s1=2*s1+Noi1;
s2=s2+Noi2;
s=s1+s2;
for i=1:d
    for j=1:63
        y1((i-1)*63+j,1)=s((i-1)*63+j,1)*c(j,1);
        y2((i-1)*63+j,1)=s((i-1)*63+j,1)*c1(j,1);
    end
end
y=[y1';y2']';
R=[1/63,0;0,1/63];
y9=y*R;%求得的最后的目标信号
%9被换成了*号
x19=zeros(d,1);
x29=zeros(d,1);
for i=1:d
    for j=1:63
        x19(i,1)=x19(i,1)+y9((i-1)*63+j,1);
        x29(i,1)=x29(i,1)+y9((i-1)*63+j,2);
    end
    x19(i,1)=sign(x19(i,1));
    x29(i,1)=sign(x29(i,1));
end
wu1=b1'-x19;
wu2=b2'-x29;
wu19=0;
wu29=0;
for i=1:d
    wu19=wu19+wu1(i,1)^2;
    wu29=wu29+wu2(i,1)^2;
end
wu1990(t,1)=wu19/(4*d);
wu2990(t,1)=wu29/(4*d);
wu199(t,1)=wu199(t,1)+wu1990(t,1);
wu299(t,1)=wu299(t,1)+wu2990(t,1);
wu19=0;
wu29=0;
end
wu199(t,1)=wu199(t,1)/l;
wu299(t,1)=wu299(t,1)/l;
end
%wu299=wu299*100;
figure(1)
%plot(wu199);
L=zeros(l1,1);
for i=1:l1
    L(i,1)=1/i;
end
semilogy(wu199,L,'-r*');
grid on;
xlabel('NSR Value');
ylabel('BER');
hold on
semilogy(wu299,L,'-g*');
hold off
title('单径信道的解相关解测其误码性能')
legend('one user','the other user')