%定義小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;