mgen.m

北邮《通信原理基于Matlab的计算机仿真》书的源码 m文件 好不容易弄来的

%m序列发生器
function [out,pState] = mgena(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
pState = bin2dec( char(curState+48) );