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