mm.m

来自「硕士CDMA作业」· M 代码 · 共 43 行

%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,
        a1(i,j)=b1(1,j);
    end
end
for k=1:63,
    if (a1(k,7)==0),
        c1(1,k)=1;
    else
        c1(1,k)=-1; 
    end
end
for k=1:63,
    c1(1,k+63)=c1(1,k);
    c1(1,k+126)=c1(1,k);
end
    
for t=0:100,
   s1=0 
   for n=1:63,
        s1=s1+c1(1,n)*c1(1,n+t);
    end
    r1(1,t+101)=s1/63
end
for k=1:100,
    r1(1,k)=r1(1,202-k)
end
    plot((-100:100),r1)
    xlabel('1+x^5+x^6的自相关函数')