zuoye3.asv

硕士CDMA作业

clear all
%function zuoye3
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

l=25;
for t=1:l
d=100;
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);
Noi11=normrnd(0,t,63*d,1);
Noi12=normrnd(0,t,63*d,1);
Noi21=normrnd(0,t,63*d,1);
Noi22=normrnd(0,t,63*d,1);
w1=20;
w2=30;
w3=40;
w4=50;
s11=s1+Noi11;
s12=s1+Noi12;
s21=s2+Noi21;
s22=s2+Noi22;
zn=10;
for i=1:zn
    s110(:,i)=s11.*exp(-j*pi*sin(w1*pi/180));
    s120(:,i)=s12.*exp(-j*pi*sin(w2*pi/180));
    s130(:,i)=s21.*exp(-j*pi*sin(w3*pi/180));
    s140(:,i)=s22.*exp(-j*pi*sin(w4*pi/180));
end
s=zeros(63*d+3,zn);
%延迟分别取1、2、3
s(1,:)=s110(1,:);
s(2,:)=s110(2,:)+s120(1,:);
s(3,:)=s110(3,:)+s120(2,:)+s130(1,:);
s(4,:)=s110(4,:)+s120(3,:)+s130(2,:)+s140(1,:);
for i=5:63*d
    s(i,:)=s110(i,:)+s120(i-1,:)+s130(i-2,:)+s140(i-3,:);
end
s(63*d+1,:)=s120(63*d,:)+s130(63*d-1,:)+s140(63*d-2,:);
s(63*d+2,:)=s130(63*d,:)+s140(63*d-1,:);
s(63*d+3,:)=s140(63*d,:);
%取第一、二路第一个信号
for i=1:d
    for j=1:63
        for k=1:zn
        y1((i-1)*63+j,k)=s((i-1)*63+j,k)*c(j,1);
        %y2((i-1)*63+j,k)=s((i-1)*63+j+2,k)*c1(j,1);
        end
    end
end
Rxx=s'*s;
Ryy=y1'*y1;
[V,D]=eig(Rxx,Ryy);
Ruu=(63/62)*Rxx-(1/62)*Ryy;
w=D(:,1)\Ruu;
Z0=y1*w';
x19=zeros(d,1);
for i=1:d
    for j=1:63
        x19(i,1)=x19(i,1)+Z0((i-1)*63+j,1);
    end
    x19(i,1)=sign(x19(i,1));
end
wu=b1'-x19;
wu19=0;
for i=1:d
    wu19=wu19+wu(i,1)^2;
end
wu199(t,1)=wu19/(4*d)