m_seq127.m

来自「产生周期为1023=(2^10-1)的Gold序列 delay_numb:G1固」· M 代码 · 共 25 行

function [m_seq] = m_seq127(c_poly)
%产生周期为127=(2^7-1)的m序列
%指定的本原多项式:c_poly
%默认本原多项式:f(x)=313=[1 0 1 0 0 1 1]
%===============================================
%设置本原多项式
if exist('c_poly')==0,
    c_poly = [1 0 1 0 0 1 1];    %注意顺序
end;
%序列生成   
x0 = zeros(1,7);  x0(7) = 1;
m_seq = zeros(1,127);
for i = 1:127,
   m_seq(i) = x0(7);
   a = x0*c_poly';
   b = 1;
   if(mod(a,2)==0)
      b = 0;
   end;
   x0(2:7) = x0(1:6);    x0(1) = b;
   if(m_seq(i)==0)
      m_seq(i) = -1;
   end;
end;