m_sequence.m

在matlab下产生m序周期,初态,特征式可以在主程序中自己设计.

function [m_out]=m_sequence(fbconnection,m,chu)
len=length(fbconnection);
n=2^len-1;                %序列长度
register=chu;             %初态
fan=0;
for i=1:n*m               % i循环一次输出一个m_out
    m_out(i)=register(len);
    for k=1:len                  % k循环累加反馈
        fan=fan+fbconnection(k)*register(k);
    end
    fan=mod(fan,2);            % 反馈系数
    newregister(1)=fan;
    for j=2:len               % j循环移位
        newregister(j)=register(j-1);
    end
    register=newregister;
    fan=0;
end