7-20.m

来自「经典数字信号处理滤波器的源代码 重点是利用巴特沃斯模拟滤波器转而设计其它数字滤」· M 代码 · 共 22 行

%例程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