7-20.m

matlab数字信号与应用源代码7

%例程7－20  m序列的线性移位寄存器

function [seq]=m_seq(connections)
%connections为移位寄存器连接关系
%connections=[0 1 0 0 1]
m=length(connections);
L=2^m-1;                                  %移位寄存器长度
registers=[zeros(1,m-1) 1];                    %寄存器初始值
seq(1)=registers(m);                          %序列第一位
for i=2:L
    temp=0;
    for j=1:m
        temp=temp+connections(j)*registers(j);
    end
    new_reg_cont(1)=mod(temp,2);
    for k=2:m
        new_reg_cont(k)=registers(k-1);
    end
    registers=new_reg_cont;                  %当前寄存器内容
    seq(i)=registers(m);                      %序列下一位
end