m伪随机序列Matlab源代码

%5阶m序列

% 在MATLAB命令窗口输入以下：

% fbconnection=[0 1 0 0 1];

% mseq=m_sequence(fbconnection);

% mseq

function[mseq]=m_sequence(fbconnection)

n=length(fbconnection);

N=2^n-1;

register=[zeros(1,n-1) 1];  %移位寄存器的初始状态

mseq(1)=register(n);        %m序列的第一个输出码元

for i=2:N

    newregister(1)=mod(sum(fbconnection.*register),2);

    for j=2:n,

        newregister(j)=register(j-1);

    end;

    register=newregister;

    mseq(i)=register(n);

end