m2.m

来自「信号检测作业程序」· M 代码 · 共 29 行

%1+x^2+x^5
m2=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)=1;
a4(1,1)=0;
m2(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)+a2(1,i-1)),2);
    m2(1,i)=a0(1,i);
end
for i=1:31
    if (m2(1,i)==0)
        m2(1,i)=-1;
    end
end
save mm2 m2