% coeff: 本原多项式的系数
% reg: 寄存器的初始值
% mseq: 生成的m序列
function [mseq] = m_sequence(coeff, reg)
n = length(coeff);
N = 2^n-1;
register = reg; %定义移位寄存器的初始状态
newregister = zeros(1, n);
mseq = zeros(1, N);
for i = 1:N
mseq(i) = register(n);
newregister(1)= mod(sum(coeff.*register),2);
newregister(2:n) = register(1:n-1);
register = newregister;
end