m_sequence1.asv

信号检测作业程序

%1+x+x^2+x^3+x^5
m1=zeros(1,31);
a0=zeros(1,31);
a1=zeros(1,31);
a2=zeros(1,31);
a3=zeros(1,31);
a4=zeros(1,31);
a0(1,1)=0;
a1(1,1)=0;
a2(1,1)=0;
a3(1,1)=0;
a4(1,1)=1;
m1(1,1)=a0(1,1);
for i=2:31
    a0(1,i)=a1(1,i-1);
    a1(1,i)=a2(1,i-1);
    a2(1,i)= a3(1,i-1);
    a3(1,i)=a4(1,i-1);
    a4(1,i)=mod((a3(1,i-1)+a2(1,i-1)+a1(1,i-1)+a0(1,i-1)),2);
    m1(1,i)=a0(1,i);
end
for i=1:31
    if (m1(1,i)==0)
        m1(1,i)=-1;
    end
end