function  mseq(k,num,len)
% k--模2加法門的輸入（來(lái)自第k個(gè)寄存器的輸出）
% num--寄存器的數(shù)目
% len--需要輸出序列的長(zhǎng)度
% 測(cè)試知有5個(gè)寄存器時(shí)，k=3時(shí)得到此時(shí)的m的序列;
% 測(cè)試知有6個(gè)寄存器時(shí)，k=5時(shí)得到此時(shí)的m的序列,測(cè)試發(fā)現(xiàn)要產(chǎn)生m序列，k的取值沒(méi)有規(guī)律性
a=ones(1,num);  %寄存器初始化(全為1)
m=zeros(1,len);  %輸出的m序列初始化
for i=1:len
    m(i)=a(num);       
    if(a(k)+a(num))==2|(a(k)+a(num))==0
           temp= 0;
      else
            temp=1;
     end
    for j=num:-1:1
        if j-1==0
            a(j)=temp;        
        else
        a(j)=a(j-1);
        end
    end
end
 m
 t=0:len-1;
 axis([0,len-1,0,2]),hold on
 stairs(t,m)