function [seq]=m_sequence(connections)
%MSRG（模件抽頭型）結構
connections =[1 0 0 1];
m=length(connections);%移位寄存器的級數
L=2^m-1;%m序列長度
registers=[zeros(1,m-1) 1];%寄存器初始化
seq(1)=registers(m);%m序列的第一位取移位寄存器移位輸出的值
for i=2:L,
    new_reg_cont(1)=connections(1)*seq(i-1);%新寄存器的第一位等于連接值乘寄存器最后一位
    for j=2:m,
        new_reg_cont(j)=rem(registers(j-1)+connections(j)*seq(i-1),2);%其他位等于前邊的寄存器值加上連接值乘寄存器最后一位
    end
    registers=new_reg_cont;
    seq(i)=registers(m);%經過一次循環寄存器輸出一位得到m序列的其他位
end