mseq6.m

根据M序列的产生原理

%function  mseq
a=ones(1,6);  %寄存器初始化(全为1)
mseq=zeros(1,252);
for i=1:252
   mseq(i)=a(6);       
    if(a(5)+a(6))==2|(a(5)+a(6))==0
           temp= 0;
      else
            temp=1;
     end
    for j=6:-1:1
        if j-1==0
            a(j)=temp;        
        else
        a(j)=a(j-1);
        end
    end
end
 for  i=1:252
     if mseq(i)==1
         mseq(i)=-1;
     else mseq(i)=1;
     end
 end
 time=1:252;
 simin=[time' mseq'];
 
 u=mseq(64:252);  %输入
 z=simout(64:252);%输出
 
 Ruz(1)=u*z/189;
 for i=2:189   %产生U(m-i) 
     temp=u(1); 
     u(1)=u(189);    
     for k=189:-1:3
         u(k)=u(k-1);
     end
      u(2)=temp;       
   Ruz(i)=u*z/189;
 end
 Ruz;
 c= Ruz(189)

g=zeros(1,189);
 for k=1:189
      g(k)=(Ruz(k)+c)*63/64;
 end
g;
t=1:189;
plot(t,g)