ms.m

此函数用来生成最大长度线性移位寄存器序列(m序列)

function [seq]=ms(fx,gx,lx)
%***************************************************
% 此函数用来生成最大长度线性移位寄存器序列(m序列)
% fx为特征多项式向量，gx为移位寄存器初始值向量
% lx为所需序列长度
%***************************************************
fx=[1 0 0 1 1];
gx=[1 0 0 0];
n=length(gx);      %n为移位寄存器级数
N=2^n-1;     %N为生成的m序列的周期
lx=30;
fb=0;
for i=1:lx
    seq(i)=gx(n);
    for j=1:n
        fb=fb +fx(j)*gx(j);
        fb=(mod(fb,2));
    end
    for t=n:-1:2      %寄存器移位
        gx(t)=gx(t-1);
    end
    gx(1)=fb;
    fb=0;
end