linearfeedbackshiftregister.m

来自「CDMA的Matlab例程」· M 代码 · 共 21 行

%File: LinearFeedbackShiftRegister.m
function [y,outstate]=...
    LinearFeedbackShiftRegister(R,generator,instate,N)
if max(generator)>R
    error(['The degree of the generetor polynomial,',...
        int2str(max(generator)),', cannot exceed R,',int2str(R),'.']);
end
if length(instate)>R
    error(['The length of the input state vector,',...
        int2str(length(instate)),', cannot exceed R,',int2str(R),'.']);
end
a=sort(generator);P=length(generator);M=length(instate)+1;
for k=1:N
    fee=instate((generator(2)));
    for q=3:P
        fee=bitxor(fee,instate(generator(q)));
    end
    instate=[fee,instate];y(k)=instate(1);instate(M)=[];
end
outstate=instate;
%End of function file.