m_sequence.m

信号检测作业程序

m0=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;
m0(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((a0(1,i-1)+a3(1,i-1)),2);
    m0(1,i)=a0(1,i);
end
for i=1:31
    if (m0(1,i)==0)
        m0(1,i)=-1;
    end
end