mm.asv

硕士CDMA作业

%1+x^5+x^6%
a1=zeros(63,7)
b1=zeros(7,1)
c1=zeros(1,189)
r1=zeros(1,201)
b1=[0 0 0 0 0 1 0]
a(1,6)=1
for i=1:63,
    b1(1,7)=b1(1,6);
    b1(1,6)=b1(1,5);
    b1(1,5)=b1(1,4);
    b1(1,4)=b1(1,3);
    b1(1,3)=b1(1,2);
    b1(1,2)=b1(1,1);
    b1(1,1)=xor(b1(1,6),b1(1,7));
    for j=1:7,
        a(i,j)=b1(1,j);
    end
end
for k=1:63,
    if (a(k,7)==0),
        c(1,k)=1;
    else
        c(1,k)=-1; 
    end
end
for k=1:63,
    c(1,k+63)=c(1,k);
    c(1,k+126)=c(1,k);
end
    
for t=0:100,
   s=0 
   for n=1:63,
        s=s+c(1,n)*c(1,n+t);
    end
    r(1,t+101)=s/63
end
for k=1:100,
    r(1,k)=r(1,202-k)
end
    plot((-100:100),r)