m_xulie.asv

硕士CDMA作业

clear all
%function m_xulie()
c=63;
r=[0,0,0,0,1,1];
x=zeros(c,6);
x1=zeros(c,6);
z=zeros(128,6);
z1=zeros(128,6);
m=31;
y=zeros(1,m);
y1=zeros(1,2*m-1);
x(1,:)=r;
x1(1,:)=r;
for i=2:64
    x(i,2:6)=x(i-1,1:5);
    %x(i,1)=mod((x(i-1,1)+x(i-1,6)),2);
    x(i,1)=xor(x(i-1,1),x(i-1,6));
    x1(i,2:6)=x1(i-1,1:5);
    x1(i,1)=mod((x(i-1,1)+x(i-1,2)+x(i-1,5)+x(i-1,6)),2);
end
for i=1:64
    for j=1:6
        x(i,j)=(-1)^(x(i,j));%把0,1变成1.-1
    end
end
z=[x;x];
z1=[x1;x1];
for p=1:m
    for q=1:64
        y(1,p)=y(1,p)+z(q,1)*z(q+p-1,1);%求相关
    end
end
%把相关的对称部分求出来
a=1;
for i=2:31
    y1(1,a)=y(1,33-i);
    a=a+1;
end
for i=1:31
    y1(1,a)=y(1,i);
    a=a+1;
end
plot((-30:30),y1);