%m序列生成函數msequence.m
%輸入：本原多項式表中的8進制數
%輸出：m序列
function mseq=msequence(prim);
data=de2bi(oct2dec(prim));%8進制本原多項式直接轉換成2進制。
data=data(end-1:-1:1);%2進制本原多項式位數順序顛倒并去掉最后一位（其恒為1）。
n=length(data);%該序列的級數
N=2^n-1;%該序列的長度
register=ones(1,n); %n級移位寄存器賦初值全“1”
mseq=zeros(1,N);
mseq(1)=register(n);%讀取寄存器的最后一位值，即寄存器的輸出
for i=2:N
  newregister=mod(sum(data.*register),2);%輸出反饋
  for j=2:n
    newregister(j)=register(j-1);%逐級右移
  end
  register=newregister;%寄存器賦新值
  mseq(i)=register(n);
end