lmsequence.m

来自「实验发现」· M 代码 · 共 18 行

%定义小m序列，a为生成多项式（常数1不用输入），b为寄存器初始值
function f=lmsequence(a,b)
   r=length(a);
   n=2^(r)-1;
   for i=1:n
       f(i)=b(r);
       x=a.*b;
       x1=mod(sum(x),2);    %往第n个移位寄存器输入的值
       for m=r:-1:2
           b(m)=b(m-1);
       end
       b(1)=x1;
   end
   f;