%%% m序列產生方案
%%模擬移位寄存器
function b=mxulie(a,c,n)

b=[];
t=2.^n-1;
for i=1:t 
    b(i)=a(n);
 d=mod(sum(c.*a),2)  ;  
     
  for j=n:-1:2  
      a(j)=a(j-1);
    end
       a(1)=d;
    
end
b;