m_pn1.m

来自「建立了信号从扩频、加干扰信号到解扩、求误码率的仿真模型系统。 经过小波镜像处理」· M 代码 · 共 30 行

function M_seq=M_pn(f)   
%产生M序列
% f=[1 1 0 0 1 1];f的格式为f＝[Cr Cr-1 ....C1],C0=1是必须的。
r=length(f);
L=2.^r;
M_seq(1)=1;
Reg=ones(1,r);
R(1,:)=Reg;
for i=1:(L-1)
 s=0;
 x=0;
   for p=1:(r-1)
     x=x|Reg(p);
   end;
   x=~x;
   for k=1:r;
     if f(k)==1
     s=xor(s,Reg(k));
     end;
   end;
   s=xor(s,x);
   for j=(r-1):-1:1
     Reg(j+1)=Reg(j);
   end;
M_seq(i+1)=Reg(r);
Reg(1)=s;
R(i+1,:)=Reg;
end;
% M_seq=2*M_seq-1;
% plot(xcorr(M_seq,M_seq));grid;