pndemo.m

3种pn序列的产生代码

%n=5（对应5级移位寄存器）的m序列发生器
%File:PNdemo.m
function [pndata]=PNdemo(pntaps); %registers为反馈线状态c1-cn
pntaps=[0 1 0 0 1];
N=length(pntaps);%寄存器个数
pninitial=[0 1 0 0 0 ]; %寄存器的初始状态
pndata=zeros(1,2^N-1); %寄存器输出向量，长度为2^N-1，初始值设定为0
pnregister=pninitial; %设定寄存器的状态为初始状态
n=0;
kk=0;
while kk==0
      n=n+1;
      pndata(1,n)=pnregister(1,1);
      feedback=rem((pnregister*pntaps'),2); %计算反馈值，对应反馈系数为c0=1
      pnregister=[feedback,pnregister(1,1:N-1)];
      if pnregister==pninitial; %如果寄存器状态回到了初始状态，循环终止
      kk=1;
    end
end
pndata=2*pndata-1; %PN序列取值从{0，1}变化为{+1，=1}
Rm(1)=sum(pndata.*pndata) %求解PN序列发生器输出的自相关函数
for k=1:n-1
     Rm(k+1)=sum(pndata.*circshift(pndata,[0,k]))
end
subplot(2,1,1);
x=[0:n-1];
stem(x,pndata);   
title('m发生器的输出')                     
axis([0 30 -1.5 1.5]);
subplot(2,1,2);
plot(x,Rm);
title('m序列的自相关函数')
axis([0 30 -1.5 30]);
%程序结束