%%本函數用來生成m序列        不知道哪里錯了
%輸入格式為  seq=mrand（【1 0 0 1 0】）
%參數說明 connection為連接多項式，1表示相連  0表示斷開
%seq為輸出的結果m序列
function sep=mrand2(connection)
m=length(connection)
L=2^m-1
registers=[zeros(1,m-1) 1]
seq(1)=registers(m);
for i=2:L
    new_reg_cont(1)=connection(1)*seq(i-1);
    for j=2:m
        new_reg_cont(j)=bitxor(registers(j-1),connection(j)*seq(i-1));
    end
    registers=new_reg_cont;
    seq(i)=registers(m);
end