function [out]=mgen(g,state,N)
%輸入g:m序列生成多項(xiàng)式（10進(jìn)制輸入）
%state:寄存器初始狀態(tài)（10進(jìn)制輸入）
%     N:輸入序列長(zhǎng)度
%test g=11;state=3;N=15;
gen=dec2bin(g)-48;
M=length(gen);
curState=dec2bin(state,M-1)-48;

for k=1:N
    out(k)=curState(M-1);
    a=rem(sum(gen(2:end).*curState),2);
    curState=[a curState(1:M-2)];
end