m_test.m

自己编的基于最抖下降法

connections=[1 1 1 1 0 1];

m=length(connections)-1;

L=2^m-1; %移位寄存器序列的长度
registers=[ones(1,m-1) 1]; %寄存器的初始内容
seq(1)=registers(m); %seq表示得到的m序列
for ii=2:L
    s=0;
    new_reg_cont(1)=0;
    for jj=2:m
        new_reg_cont(jj)=registers(jj-1); %移位
    end;
    for jj=1:m
      %  new_reg_cont(1)=mod(new_reg_cont(1)+connections(jj)*registers(jj),2);
      s=connections(jj+1)*registers(jj)+s;
    end
    new_reg_cont(1)=mod(s,2);
    registers=new_reg_cont; % 当前的寄存器内容
    seq(ii)=registers(m); %m序列的下一个元素
end; 
bb=[seq seq seq seq seq seq seq seq seq];
for loop=0:93
    temp=0;
    for loop2=1:31
        temp=((-1)^(seq(loop2)))*((-1)^(bb(93+loop2+loop)))+temp;
        %temp=temp+((-1)^(seq(xxx)))*((-1)^(seq(xxx)))
    end
    rp(loop+1)=temp/31;
end
plot(rp);