m_sgenerator.m

自己编的基于最抖下降法

function [m]=msgenerator(n,c)
%本程序用于产生m序列

%参数c为序列生成多项式的系数向量
%参数n为序列生成多项式的阶数

%初始化
x=[ones(1,n)];    %赋初值 x第一位为1，其余n-1位为0
q=2^n-1;               %m序列周期
m=[1,zeros(1,q-1)];    %m序列初值 m第一位为1，其余q-1位为0

%产生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;               %将中间变量的值赋给x(1)
 m(i)=x(n);            %移位寄存器最后一位取出生成m序列
end;
for i=1:q
    if m(i)==0
        m(i)=-1;
    end;
end;