%本原多項(xiàng)式產(chǎn)生127位m序列;
%本原多項(xiàng)式可選用1，211；3，217；5，235；7，367；9，277；11，325；13，203；19，313；21，345
%產(chǎn)生9個(gè)127位的m序列
clc;
clear;
%在每一連接后均存在一個(gè)1才能符合本原多項(xiàng)式的8進(jìn)制表示；
connection(1,:) = [1 0 0 0 1 0 0 ];                     %211
connection(2,:) = [1 0 0 1 1 1 0 ];                     %235
connection(3,:) = [1 0 0 0 1 1 1 ];                     %217
connection(4,:) = [1 1 1 1 0 1 1 ];                     %367
connection(5,:) = [1 0 1 1 1 1 1 ];                     %277
connection(6,:) = [1 1 0 1 0 1 0 ];                     %325
connection(7,:) = [1 0 0 0 0 0 1 ];                     %203
connection(8,:) = [1 1 0 0 1 0 1 ];                     %313
connection(9,:) = [1 1 1 0 0 1 0 ];                     %345
%對抽頭系數(shù)進(jìn)行翻轉(zhuǎn)
for i = 1:9
    for j = 1:7
        link(i,j) = connection(i,8-j);
    end;
end;

%求出9個(gè)127位的m序列
for m = 1:9
    x(1,1:6)=zeros(1,6);
    x(1,7) = 1;
	for i=1:127
        x(i+1,1) = mod(sum(and(x(i,:),link(m,:))),2);
        for j=1:6
            x(i+1,j+1)=x(i,j);
        end;
	end;
    pn(m,:) = x(:,7)';
end;