mseq4.m

根据M序列的产生原理

function  mseq4(len)
% k--模2加法门的输入（来自第k个寄存器的输出）
% num--寄存器的数目
% len--需要输出序列的长度
% 测试知有5个寄存器时，k=3时得到此时的m的序列;
% 测试知有6个寄存器时，k=5时得到此时的m的序列,测试发现要产生m序列，k的取值没有规律性
a=ones(1,4);  %寄存器初始化(全为1)
m=zeros(1,len);  %输出的m序列初始化
for i=1:len
    m(i)=a(4);       
    if(a(3)+a(4))==2|(a(3)+a(4))==0
           temp= 0;
      else
            temp=1;
     end
    for j=4:-1:1
        if j-1==0
            a(j)=temp;        
        else
        a(j)=a(j-1);
        end
    end
end
 m
 for  i=1:len
     if m(i)==1
         m(i)=-1;
     else m(i)=1;
     end
 end
         
 t=0:len-1;
 axis([0,len-1,-1.5,1.5]),hold on
 stairs(t,m),grid on