%本原多項式產生127位m序列;
%本原多項式可選用1，211；3，217；5，235；7，367；9，277；11，325；13，203；19，313；21，345
%產生9個127位的m序列
clc;
clear;
%在每一連接后均存在一個1才能符合本原多項式的8進制表示；
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
%對抽頭系數進行翻轉
for i = 1:9
    for j = 1:7
        link(i,j) = connection(i,8-j);
    end;
end;

%求出9個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;