小m序列的源代码小m序列的源代码小m序列的源代码

%n=5（对应5级移位寄存器）的m序列发生器 
%File:PNdemo.m 
function [pndata]=PNdemo（pntapspninitial）; %registers为反馈线状态c1-cn 
% pntaps=[0 1 0 0 1]; 
% pninitial=[0 1 0 0 0 ]; %寄存器的初始状态 
N=length（pntaps）;%寄存器个数 
pndata=zeros（12^N-1）; %寄存器输出向量，长度为2^N-1，初始值设定为0 
pnregister=pninitial; %设定寄存器的状态为初始状态 
n=0; 
kk=0; 
while kk==0 
      n=n+1; 
      pndata（1n）=pnregister（11）; 
      feedback=rem（（pnregister*pntaps‘）2）; %计算反馈值，对应反馈系数为c0=1 
      pnregister=[feedbackpnregis