msjgenerator.m

来自「卷积码程序很适用的」· M 代码 · 共 26 行

function [m]=msjgenerator(n,c)


%初始化
x=[1,zeros(1,n-1)];
q=2^n-1;               
m=[0,zeros(1,q-1)];

%产生m序列    
for i=2:q  
 %求模2和
 s=0;                 
 for t=1:n            
  s=c(t+1)*x(t)+s;
 end;
   k=mod(s,2);
   
 %移位  
 for r=n:-1:2          
  x(r)=x(r-1);
 end; 
 
 x(1)=k; 
 m(i)=x(n);            
end;