function mseq= m_seq(prim_poly);
a= de2bi( oct2dec(prim_poly));% 8進(jìn)制本原多項(xiàng)式直接轉(zhuǎn)換成2進(jìn)制。
a= a(end-1:-1:1);% 2進(jìn)制本原多項(xiàng)式位數(shù)順序顛倒。
n= length(a);
N= 2^n-1;
register= ones(1,n); % n級(jí)移位寄存器賦初值全“1”。
mseq= zeros(1,N);
mseq(1)= register(n);
for i= 2:N;
    newregister(1)=mod(sum(a.*register),2);
    for j=2:n
         newregister(j)=register(j-1);
    end
        register=newregister;
        mseq(i)=register(n);
end