mgen.m

来自「通信中一些基本实验的程序AM、FM、ASK、FSK、PSK等等」· M 代码 · 共 14 行

function [out]=mgen(g,state,N)
%输入g:m序列生成多项式（10进制输入）
%state:寄存器初始状态（10进制输入）
%     N:输入序列长度
%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